Re: [R] Accessing individual variables from summary of lm

Hi Ben, > If you wouldn't mind, how do I access the individual components inside > coefficients matrix? What you want to know about is ?attributes ## attributes(model) model$coefficients model$coefficients[1] model$coefficients[2:4] model$coefficients[c(1,5)] HTH, Mark. ascentnet wrote: > >

Re: [R] Accessing individual variables from summary of lm

Hi Ben, Sorry (still a little out-of-tune), perhaps what you really need to know about is ?"[" HTH, Mark. Mark Difford wrote: > > Hi Ben, > >> If you wouldn't mind, how do I access the individual components inside >> coefficients matrix? > > What you want to know about is ?attributes > > #

Re: [R] adding the mean and standard deviation to boxplots

Dear users This is a message I was directing to Harold Baize but because I pressed the wrong button the message got lost g!!! So I’m doing it all over again: Lets suppose I have three batches of data: a <- rnorm(50,2500,300) b <- rnorm(50,3500,250) c <- rnorm(50,4000,200) # Now I want to

Re: [R] Accessing individual variables from summary of lm

Thanks for your help! If you wouldn't mind, how do I access the individual components inside coefficients matrix? thanks, Ben. Paul Smith wrote: > > On Tue, Jul 22, 2008 at 5:44 AM, ascentnet <[EMAIL PROTECTED]> wrote: >> I am trying to access the R2, intercept, and slope out of the summary >>

Re: [R] F test

--- On Tue, 7/22/08, Peter Dalgaard <[EMAIL PROTECTED]> wrote: From: Peter Dalgaard <[EMAIL PROTECTED]> Subject: Re: [R] F test To: [EMAIL PROTECTED] Cc: r-help@r-project.org Date: Tuesday, July 22, 2008, 11:51 PM Vicki Meng wrote: > Dear R users, > > I need to do a F test on the hypothesis t

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Dear all, Since you guys are frank, let me be frank as well. I did not ask anyone to impose on me their point of view on bootstrap. It's my impression that this is what you guys are trying to do - that's sad. Some of your emails in this discussion are worth less than junk mail - particularly the o

[R] rpart confidence intervals?

hi, I am using rpart with a continuous response variable. Is there a way to extract 95% confidence intervals for the predicted mean at each node? cheers. -- View this message in context: http://www.nabble.com/rpart-confidence-intervals--tp18601587p18601587.html Sent from the R help mailing lis

Re: [R] Plotting Multiple lines on one plot

You can set the range for each axis on your plot using the ylim and xlim parameters (use help to see how to implement them). You can find the range or minimum and maximum for each variable (or a number of variables at once) using the min and max commands. If you're ranges vary too much though i

[R] average replicate probe values

Hi, Could somebody tell me how I can average expression values of replicate probe sets in an data frame? Thanks Pal Kaposi-Novak MD PhD PIRT Fellow University of Pittsburgh Department of Pathology BST S408, 200 Lothrop Str Pittsburgh, PA , 15261 Tel: (412) 383-7748 [EMAIL PROTECTED] __

[R] rpart node number

Hello there. I have a tree with only 13 nodes. See my results below. How do I get node number that is up to 33? Also, what happen to the missing node number less than 13? > rpartModel <- rpart(y~., data=data.frame(y=y,x=x),method="anova", cp=0.005, + mi

Re: [R] Can't Load Text Files

Jia Ying Mei If I understand you correctly I think you need read.table(). Try: ?read.table and note the sep argument (sep='\t') HTH . Peter Alspach > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Jia Ying Mei > Sent: Wednesday, 23 July

Re: [R] Can't Load Text Files

setwd("c:\\temp") #you need to insert a ENTER after your last value > df1<-read.table("testinput.txt", head=T, sep="\t") > df1 Date France Germany 1 3/15/07 1 2 2 3/10/07 2 4 > df2<-read.table("testdate.txt", head=T, sep="\t") > df2 Date 1 3/15/07 2 3/14/07 3 3/13

[R] Can't Load Text Files

Hi, I've completely forgotten how to do R, and it seems I can't even load simple files correctly. I have two tab deliminated text files (attached, these were test files of the larger data files I'm trying to work with) that I am trying to load into R so that I can merge them. They were origi

Re: [R] rpart$where and predict.rpart

Brian, thanks for the prompt reply. I am quite new to R and have no clue on how to call internal function. Tried the following but failed: > .Internal("pred.rpart") Error in .Internal("pred.rpart") : invalid .Internal() argument > .Internal() Error in .Internal() : 0 arguments

Re: [R] scatter plot using ggplot

Thanks for this mail. However I am getting following error on that : Error in inherits(formula, "formula") : could not find function "aes" And regarding my 3rd query I actually wanted to get "bigger circles" as points, not intended to increase plot size. Regards, --- On Tue, 7/22/08, hadley w

Re: [R] an extremely simple question

Hi Edna, I am not sure the completely right answer, but I suspect that you use data(df) to invoke a data.frame from the sample datasets provided from the packages available on your R system. By other side, if you have a data.frame loaded on your environent, you can use "attach" do turn the columns

Re: [R] TextMate and R

Angelo, Assuming you are on a mac (TextMate), [EMAIL PROTECTED] is a better list for mac specific questions There is a bundle for R ( http://manual.macromates.com/en/bundles#getting_more_bundles ). Not sure if the R bundle these days comes with the default TextMate install. With the bun

[R] select significant variables

Hi, I want to do significant variables selection for continuous/discrete dependent variables. I am wondering what the best package or what the best method should be used. Thank you for your time. Best regards, Brandon _ [[elided

[R] an extremely simple question

Hi! What is the difference between using data(somedataset) vs. attach(somedataset), please? thanks, Edna Bell __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/po

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

On 23/07/2008, at 1:17 PM, Frank E Harrell Jr wrote: Michal Figurski wrote: Hmm... It sounds like ideology to me. I was asking for technical help. I know what I want to do, just don't know how to do it in R. I'll go back to SAS then. Thank you. -- Michal J. Figurski You don't understand

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Michal Figurski wrote: Hmm... It sounds like ideology to me. I was asking for technical help. I know what I want to do, just don't know how to do it in R. I'll go back to SAS then. Thank you. -- Michal J. Figurski You don't understand any of the theory and you are using techniques you don

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Michal Figurski wrote: Dear Marc and all, Thank you for all the due respect. I tried to explain as much explicitly as I could what I am trying to do in my first email. I did not invent this procedure, it was already published in the paper: T. Pawinski, M. Hale, M. Korecka, W.E. Fitzsimmons,

[R] deblur microscope image stacks

I'm new to working with microscope image stacks. I've been asked to use software to deblur the image volumes, eg: http://www.mediacy.com/index.aspx?page=AutoDeblurVisualize Is there a toolkit for this within R? I would like to try something that can take advantage of dual processors (64bit) and

[R] Plotting Multiple lines on one plot

Hi, I'm trying to plot multiple lines on one plot. I have a data frame raw, and i want to plot raw$date on the x-axis and raw$theta, raw$vega, raw$delta, and a few others on the y-axis, with a legend. However, I'm not sure what the scale of those data sets are, and I don't know which one will ha

Re: [R] Editor for Mac

Try Komodo Edit and SciviewsK. They work on Windows, Linux, and Mac. Shige On Wed, Jul 23, 2008 at 1:48 AM, Angelo Scozzarella < [EMAIL PROTECTED]> wrote: > Thanks for the answers, > but I'm looking for an editor like Tinn-R. > > Is there something like it for Mac? > > > Angelo Scozzarella > > _

Re: [R] F test

Vicki Meng wrote: Dear R users, I need to do a F test on the hypothesis that a 2 by 1 vector (X_1, X_2)' has the mean vector (M_1, M_2)'. Specifically, I would like to assume the X vector comes from a bivariate Normal distribution (M, Sigma). Then, given 1000 observations on X, I wanted to tes

Re: [R] Erro: cannot allocate vector of size 216.0 Mb

On 23/07/2008, at 10:54 AM, José Augusto Jr. wrote: My apologies, I used an improper language in my last email. My sincere apologies for anyone that became offended. I certainly wasn't offended! I was kidding. Sorry if I sounded as if I were criticizing you.

[R] F test

Dear R users,   I need to do a F test on the hypothesis that a 2 by 1 vector (X_1, X_2)' has the mean vector  (M_1, M_2)'. Specifically, I would like to assume the X vector comes from a bivariate Normal distribution (M, Sigma). Then, given 1000 observations on X, I wanted to test if the means of

Re: [R] Erro: cannot allocate vector of size 216.0 Mb

My apologies, I used an improper language in my last email. My sincere apologies for anyone that became offended. Sincerely, jamaj 2008/7/22 Rolf Turner <[EMAIL PROTECTED]>: > > On 23/07/2008, at 10:03 AM, José Augusto Jr. wrote: > > >> I'm writing this message to encourage people migrating to

Re: [R] .asp (not really an R question) but related

Its probably a Microsoft Active Server Page that is serving the pdf. The document itself is not likely an asp file. On Tue, Jul 22, 2008 at 3:24 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > I am trying to get the paper that explains the stl procedure it is a .asp > file. What can I use to rea

Re: [R] Erro: cannot allocate vector of size 216.0 Mb

On 23/07/2008, at 10:03 AM, José Augusto Jr. wrote: I'm writing this message to encourage people migrating to Ubuntu. W$ e M$ are shit. Such language in a young gentleman!!! I agree, but. I'm also amazed that your message got past the censorious censoring program `

Re: [R] Erro: cannot allocate vector of size 216.0 Mb

Dear all, I have some news. First, i sucessfully installed Ubuntu 8.04 LTS in my notebook with R too. The problem simply disappeared. Now, i can embed my 1.000.000 observations time-series even in 20 dimensions without problem, in the same notebook. I'm writing this message to encourage people

Re: [R] help with simulate AR(1) data

On 23/07/2008, at 9:50 AM, Peter Dalgaard wrote: Rolf Turner wrote: matrix(unlist(lapply(1:1000,function(x){arima.sim(list(ar=0.3), 100)+mean})),nrow=1000,byrow=TRUE) Whatever did replicate() ever do to you to deserve getting so rudely ignored? replicate(1000, arima.sim(list(ar=.3),100))

[R] TextMate and R

Hi, how can I use TextMate as editor for R? Thanks Angelo Scozzarella __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] Change font-face in title

See ?title where its mentioned that calls are ok and a formula is a call. Actually the first ~ is not needed since its already a formula. On Tue, Jul 22, 2008 at 1:12 PM, Mike Prager <[EMAIL PROTECTED]> wrote: > "Gabor Grothendieck" <[EMAIL PROTECTED]> wrote: >> >> plot(1, type="l", xlab="Wellenl

Re: [R] How to....

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Some idea of the the data looks like and what you want to do with it would be helpful. ?read.csv ?mean ?tapply On Tue, Jul 22, 2008 at 4:31 PM, Jia Ying

Re: [R] help with simulate AR(1) data

Rolf Turner wrote: matrix(unlist(lapply(1:1000,function(x){arima.sim(list(ar=0.3),100)+mean})),nrow=1000,byrow=TRUE) Whatever did replicate() ever do to you to deserve getting so rudely ignored? replicate(1000, arima.sim(list(ar=.3),100)) + mean (OK, so t(...) then) -- O__ Peter D

Re: [R] Decoding subscripts/superscripts from CSVs

Hi Nina, > Subscripts, superscripts, and italics are encoded in carats, and I was > wondering > if R can actually recognize those and print actual superscripts, etc. > Here's an example: I know that ladies are fond of diamonds (perhaps you mean carets?); though it isn't quite clear what you wan

Re: [R] How to filter a data frame?

Thank you all!! :-) rlearner309 wrote: > > I have a question about how to filter the data frame: > Suppose my data frame has variables like gender, age,... How to get a > subset of the data frame, with only female (or male) and/or age > 50...? > What is the typical syntax? I tried several

Re: [R] help with simulate AR(1) data

On 23/07/2008, at 7:52 AM, cathelf wrote: Hi, sorry for bothering your guys again. I want to simulate 100 AR(1) data with cor(x_t, x_t-1)=rho=0.3. The mean of the first 70 data (x_1 to x_70) is 0 and the mean of the last 30 data (x_71 to x_100) is 2. Can I do it in the following way? x <

Re: [R] Lattice: How to draw curves from given formulae

On 7/22/08, John Smith <[EMAIL PROTECTED]> wrote: > Thanks, Deepayan, > > I actually want to have a picture like the Figure 5.9 in "Mixed-Effects > Models in S and S-PLUS", which have a separate panel for each function and > has the function name been printed in strip. xyplot(linear(d) + logist

[R] How to....

Hi everyone, I am not new to R, but its been over a year since I've used it, so I don't remember how to do some things. I have a data set (in excel currently, but will be converted to text), that has date, country and price, with a price for every country for every date. I need to find the m

Re: [R] Cannot re-start R after bus error

thank you for the prompt help - found the file it was in a different directory R > Date: Tue, 22 Jul 2008 16:13:19 -0500 > From: [EMAIL PROTECTED] > Subject: RE: [R] Cannot re-start R after bus error > To: [EMAIL PROTECTED] > > i realize you said that

Re: [R] Cannot re-start R after bus error

On Tue, 22 Jul 2008 20:22:08 + "Ralph S." <[EMAIL PROTECTED]> wrote: > > Hi, > > I ran a program that was exhausting the (unix) server's available memory. I > got a bus error and then chose to save and quit workspace (option 4). > > That may not have been a smart idea - when I try to use R

Re: [R] rollmean and stl

On Tue, 22 Jul 2008, stephen sefick wrote: I need to investigate how rollmean and the trend returned from stl differ. I am trying to find out exactly what the trend part of stl is (I have just started coding in R and do not know fortran). I need to extract this because it will be used in furthe

[R] Cannot re-start R after bus error

Hi, I ran a program that was exhausting the (unix) server's available memory. I got a bus error and then chose to save and quit workspace (option 4). That may not have been a smart idea - when I try to use R now, I get the error message: Error in load(name, envir = .GlobalEnv) : ReadItem:

[R] Decoding subscripts/superscripts from CSVs

Hi, I have a CSV file with various biological reactions. Subscripts, superscripts, and italics are encoded in carats, and I was wondering if R can actually recognize those and print actual superscripts, etc. Here's an example: S-adenosyl-L-methionine + rRNA = S-adenosyl-L-homocysteine + rRNA c

Re: [R] Function Error

2008/7/22 Angelo Scozzarella <[EMAIL PROTECTED]>: > Hi, > > Why this function doesn't work? >return(x$breaks[cl]+ >(N/2-sum(x$counts[1:(cl-1)]))/(x$densit[cl]*N) >} else median(x) > } Have you read the error message? I get: Error: unexpected '}' in: ":(cl-1)]))/(x$densit[cl]*

Re: [R] Reading the data from specific columns

See argument colClasses in ?read.table. On Tue, Jul 22, 2008 at 5:28 PM, Vidhu Choudhary <[EMAIL PROTECTED]> wrote: > Hi All, > I have a huge tab separated file which a many columns. I want to read only > few specific columns in R object. There are no column names/headers in the > file. Can someon

[R] Reading the data from specific columns

Hi All, I have a huge tab separated file which a many columns. I want to read only few specific columns in R object. There are no column names/headers in the file. Can someone please suggest how can I do that. I checked the help of read.delim but wasnt able to find such option Thank you Reagards V

Re: [R] Function Error

At a first guess, because you are missing a right paren in thereturn() statement of the final else section. Otherwise, we'd need to know about what "doesn't work" means - what is x supposed to be? What do you want the function to do? What error messages are you getting? You haven't provided any o

[R] .asp (not really an R question) but related

I am trying to get the paper that explains the stl procedure it is a .asp file. What can I use to read this. thanks stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We

Re: [R] how to speed up this for loop?

> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Nordlund, > Dan (DSHS/RDA) > Sent: Tuesday, July 22, 2008 10:31 AM > To: r-help@r-project.org > Subject: Re: [R] how to speed up this for loop? > > > -Original Message- > > From: [EMAIL PROTEC

Re: [R] Function Error

Hi Angelo, please read the posting guide, before posting! Anyway, for general rule, if you want to see what's going on in a function at each step, it may be helpful to install library("debug") that has very useful functions as "mtrace" that my resolve your problem. Cheers Anna

Re: [R] Function Error

Please give commented, minimal, self-contained, reproducible code. In this case, not only the function definition, but a simple example showing *how* it does not work, and what you were expecting it to do! Angelo Scozzarella wrote: Hi, Why this function doesn't work? function (x) { if (

[R] help on R benchmarks

Hi, All, I am just wondering if I can find some R benchmarks that has mixed loop and function calls and both are computation expensive? I searched the R packages but could not find a statisfactory one since the examples provide small sample input data and the run time is short. A benchmark that ru

[R] help with simulate AR(1) data

Hi, sorry for bothering your guys again. I want to simulate 100 AR(1) data with cor(x_t, x_t-1)=rho=0.3. The mean of the first 70 data (x_1 to x_70) is 0 and the mean of the last 30 data (x_71 to x_100) is 2. Can I do it in the following way? x <- arima.sim(list=(ar=0.3), 100) mean <- c(rep(0, 70

[R] rollmean and stl

I need to investigate how rollmean and the trend returned from stl differ. I am trying to find out exactly what the trend part of stl is (I have just started coding in R and do not know fortran). I need to extract this because it will be used in further calculations, and it needs to be verified to

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Dear Marc and all, Thank you for all the due respect. I tried to explain as much explicitly as I could what I am trying to do in my first email. I did not invent this procedure, it was already published in the paper: T. Pawinski, M. Hale, M. Korecka, W.E. Fitzsimmons, L.M. Shaw. Limited Sam

Re: [R] Calculating Betweenness - Efficiency problem

Dear Gabor, Thank you very much for the insights. I have been using the igraph package for my computations. But I did not know about graph.data.frame(). Thanks again for that. So I did run my data using the steps you had provided. Weirdly, even though the .csv file has approximately 300,000 reco

Re: [R] Help with rearranging data

Yes, my apologies, I want to do this in R. > Date: Tue, 22 Jul 2008 11:33:01 -0700> From: [EMAIL PROTECTED]> To: [EMAIL > PROTECTED]> CC: r-help@r-project.org> Subject: Re: [R] Help with rearranging > data> > Hi,> > is your question how you do this in R? Note that this is a > mailing> list fo

Re: [R] Help with rearranging data

Hi, is your question how you do this in R? Note that this is a mailing list for R. There are probably other lists/groups better suited for Excel-specific problems. Cheers Henrik On Tue, Jul 22, 2008 at 10:59 AM, William Pepe <[EMAIL PROTECTED]> wrote: > > In Excel, I have data that looks like

[R] How to simulate heteroscedasticity (correlation)

Hi, I would like to generate two correlated variables. I found that funktion for doing that: a <- rmvnorm(n=1,mean=c(20,20),sigma=matrix(c(5,0.8*sqrt(50), 0.8*sqrt(50),10),2,2)) (using library(mvtnorm)) Now I also want to generate two correlated variables where the error variance vary

[R] Function Error

Hi, Why this function doesn't work? function (x) { if (is.factor(x)) { if (!is.ordered(x)) { warning("La mediana non si puo' calcolare!!!") return(NA) } me <- median(unclass(x)) if (me - floor(me) != 0) { warning("Mediana in

Re: [R] Lattice: How to draw curves from given formulae

Thanks, Deepayan, I actually want to have a picture like the Figure 5.9 in "Mixed-Effects Models in S and S-PLUS", which have a separate panel for each function and has the function name been printed in strip. Thanks On Tue, Jul 22, 2008 at 1:08 PM, Deepayan Sarkar <[EMAIL PROTECTED]> wrote:

Re: [R] Help with rearranging data

Here is a start that will get you a dataframe. You will have to decide how you want to convert the long to wide format since that did not clearly come through in the mail: > x <- readLines('/tempxx.txt') > # remove blank lines > x <- x[x != ""] > # create matrix > y <- matrix(x[-(1:6)], ncol=6, b

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

> install.packages('fortunes') > library(fortunes) > fortune(28) > -Original Message- > From: Marc Schwartz [mailto:[EMAIL PROTECTED] > Sent: Tuesday, July 22, 2008 1:29 PM > To: Michal Figurski > Cc: Doran, Harold; r-help@r-project.org; Frank E Harrell Jr; > Bert Gunter > Subject: Re:

[R] Help with rearranging data

In Excel, I have data that looks like this. Month Seg Prod Biller Cycle Sales 1/1/2008 A Table Phyllis 1 500 1/1/2008 A Table Phyllis 2 600 1/1/2008 A Table Phyllis 3 650 1/1/2008 A Table Doreen 1 345 1/1/2008 A Table Doreen 2 451 1/1/2008 A Table Doreen 6 550 1/1/2008 A Table Doreen

[R] Editor for Mac

Thanks for the answers, but I'm looking for an editor like Tinn-R. Is there something like it for Mac? Angelo Scozzarella __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-

Re: [R] how to speed up this for loop?

> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Rheannon > Sent: Tuesday, July 22, 2008 9:09 AM > To: r-help@r-project.org > Subject: Re: [R] how to speed up this for loop? > > > Here is the full code: > > R <- 5328 > #Temp is a temperature file w

Re: [R] data transformation

on 07/22/2008 11:24 AM Christian Hof wrote: Dear all, how can I, with R, transform a presence-only table (with the names of the species (1st column), the lat information of the sites (2nd column) and the lon information of the sites (3rd column)) into a presence-absence (0/1) matrix of specie

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Michal, With all due respect, you have openly acknowledged that you don't know enough about the subject at hand. If that is the case, on what basis are you in a position to challenge the collective wisdom of those professionals who have voluntarily offered *expert* level statistical advice t

Re: [R] randomForest Tutorial

The help pages for 'randomForest' are quite clear in stating how to use the 'randomForest()' function. Input is a data frame (not a csv-file) and a formula involving attribute names. There is a kind of tutorial for classification and clustering with Random Forests on Leo Breiman's web page

Re: [R] Does R have SQL interface in windows?

On Tue, 22 Jul 2008, rlearner309 wrote: Seems that the RmySQL package does not support windows... Not so (there is even a binary on CRANextras, but read the README), and RODBC does as well. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.s

[R] 2 Courses*** R/Splus Advanced Programming : August 14-15 in Seattle ! by XLSolutions Corp

XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our*** R/S Advanced Programming ***course in Seattle on August 14-15 and *** R/Splus Fundamentals and Programming Techniques*** course in New York City on July 28-29 Ask for group discount and reserve your seat Now - Earlybir

Re: [R] drawing segments through points with pch=1

David Epstein <[EMAIL PROTECTED]> wrote: > Please excuse me for asking such basic questions: > > Here is my code > > y=c(1.21,0.51,0.14,1.62,-0.8,0.72,-1.71,0.84,0.02,-0.12) > > ybar=mean(y) > > ll=length(y); > > ybarv=rep(ybar,ll) > > x=1:ll > > plot(x,ybarv,pch=1) > > segments(x[1],ybar,x[ll],y

Re: [R] Change font-face in title

"Gabor Grothendieck" <[EMAIL PROTECTED]> wrote: > > plot(1, type="l", xlab="Wellenlänge [nm]", col="darkslategray", > main = ~ Spektrum ~ italic("Deschampsia caespitosa")) That clearly works, but why? Does the tilde cause the rest of the text to be parsed as an expression? This looks ext

Re: [R] tklistbox and extracting selection to R

The <<- assignment can be very tricky and Tk windows have their own world. If you just copy and paste the below (with ntx and levels.tx defined, it would have been nice to have the example fully self running), then the OnOk function has been defined and linked to the button, but never run, so w

Re: [R] Lattice: How to draw curves from given formulae

On 7/21/08, John Smith <[EMAIL PROTECTED]> wrote: > Dear R Users: > > I have a list function as: > Flat: y = 0 > Linear: y = -(1.65/8)d > Logistic: y = 0.015 - 1.73/{1+exp[1.2(4-d)]} > Umbrella: y= -(1.65/3)d + (1.65/36)d^2 > Emax: y = -1.81d/(0.79+d) > Sigmoid Emax: y = -1.70d^5/(4^5+d^5)

Re: [R] S4 : setGeneric for classical methods

Prof Brian Ripley <[EMAIL PROTECTED]> a écrit : On Tue, 22 Jul 2008, [EMAIL PROTECTED] wrote: Hi the list (well, half of the list, only the one who are not on holidays...) I am trying to make an S4 package. When I run it on a console, everything seems ok. When I run R CMD check, I got an

Re: [R] normalised/transformed regressions

Many thanks, let me find that, Tolga "Greg Snow" <[EMAIL PROTECTED]> 22/07/2008 18:01 To "[EMAIL PROTECTED]" <[EMAIL PROTECTED]> cc "r-help@r-project.org" Subject RE: [R] normalised/transformed regressions There is a whole taskview on CRAN devoted to robust statistics (including varian

Re: [R] normalised/transformed regressions

There is a whole taskview on CRAN devoted to robust statistics (including variants of regression), I would suggest you look through that to see if any of those tools work for you. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111

Re: [R] S4 : setGeneric for classical methods

On Tue, 22 Jul 2008, [EMAIL PROTECTED] wrote: Hi the list (well, half of the list, only the one who are not on holidays...) I am trying to make an S4 package. When I run it on a console, everything seems ok. When I run R CMD check, I got an error --- 8< -- Error in setMeth

Re: [R] normalised/transformed regressions

Hi Greg, First of all, many thanks for your detailed response. I was trying to achieve robustness in the regressions without making any assumptions about any transformations that might be needed, as I have no prior apart from realised. This seemed like a neat way of doing it, but there may be

[R] tklistbox and extracting selection to R

Dear experts, I am trying to understand why is it that when I paste (into the R console) the following code to select an option from a list: require(tcltk) tt<-tktoplevel() tl<-tklistbox(tt,height=ntx,selectmode="single",background="white") tkgrid(tklabel(tt,text="Select the legend of

Re: [R] scatter plot using ggplot

On Tue, Jul 22, 2008 at 3:42 AM, Megh Dal <[EMAIL PROTECTED]> wrote: > I used ggplot to create a scatter plot : > > library(ggplot) > library(mnormt) > Sigma = matrix(c(1, 0.6, 0.6, 1), 2, 2) > x = rmnorm(20, c(0,0), Sigma) > xx = x[order(x[,1]),] > y = xx[,1] > z = xx[,2] > qplot(z, y, type="point

Re: [R] how to speed up this for loop?

Here is the full code: R <- 5328 #Temp is a temperature file with daily temp. values for each day of the year, for 5328 locations Temp <- rnorm(100) dim(Temp) <- c(365,5328) # transpose temp dataframe in order to perform rollmean (running mean) on proper data section Temp_T <- t(Temp) library(zoo

[R] data transformation

Dear all, how can I, with R, transform a presence-only table (with the names of the species (1st column), the lat information of the sites (2nd column) and the lon information of the sites (3rd column)) into a presence-absence (0/1) matrix of species occurrences across sites, as given in the

[R] S4 : setGeneric for classical methods

Hi the list (well, half of the list, only the one who are not on holidays...) I am trying to make an S4 package. When I run it on a console, everything seems ok. When I run R CMD check, I got an error --- 8< -- Error in setMethod("plot", "ClassX", function(x, y) { : no exis

Re: [R] Does R have SQL interface in windows?

Its hard to tell from your post why you believe it is not supported by Windows. Please read the posting guide prior to posting. Some possible problemsyou may have: - spelling is RMySQL, not RmySQL as in your post - maybe you simply have not found the binary? Set your mirror to the CRAN extras

Re: [R] normalised/transformed regressions

It is possible to write a function to do what you describe, but the real question is why would you want to do that? It looks like you are trying to force your data to fit a set of assumptions that are not needed. The normality assumption in regression models is that the residuals are normal, o

Re: [R] saving plot both as jpg and pdf

Just an advice: I recommend use the png device instead of jpeg for plotting purposes, it gives better results for graphs. Jpeg is better for photos. Off course IMHO :-) I wish a nice day Tomas stephen sefick wrote: you have to go in sequence- you are opening a pdf device then a jpeg device

Re: [R] saving plot both as jpg and pdf

Plotting functions only send the plot to the current device (in your case the 2nd one opened). You can see all the currently open devices with dev.list, you can see which is current with dev.cur, you can change which is current using dev.set (and possibly dev.next and dev.prev). Look at the he

Re: [R] How to filter a data frame?

Hi, rlearner309 wrote: I have a question about how to filter the data frame: Suppose my data frame has variables like gender, age,... How to get a subset of the data frame, with only female (or male) and/or age > 50...? What is the typical syntax? I tried several condition expressions, but non

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

The bootstrap **can** be used for bias correction. However, it may not be such a good thing to do. I quote from Efron and Tibshirani's AN INTRODUCTION TO THE BOOTSTRAP (p.138): "... bias estimation is usually interesting and worthwhile, but the exact use of a bias estimate is often problematic. Bi

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Hi Michal, This paper by John Fox may help you to precise what you are looking for and to perform your analyses http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-bootstrapping.pdf Nael On Tue, Jul 22, 2008 at 3:51 PM, Michal Figurski < [EMAIL PROTECTED]> wrote: > Dear all, > > I don't

Re: [R] saving plot both as jpg and pdf

Hi, Thank you both for the suggestions it seems i have to send the plot twice since i cannot have 2 devices active in the same time . I Knew about savePlot but this saves jpf only at 72 dpi and i wanted a better resolution for jpg files. I will look into dev.copy2pdf, maybe i can hav

Re: [R] saving plot both as jpg and pdf

you have to go in sequence- you are opening a pdf device then a jpeg device and ploting to the jpeg device. My suggestion would be to do pdf() plot(yourplot) dev.off() jpeg plot(yourplot) dev.off() this should work On Tue, Jul 22, 2008 at 9:50 AM, Monica Pisica <[EMAIL PROTECTED]> wrote: > > H

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Probably a good idea for you. The R help list is useful for both programming AND statistical advice for those who want it. > -Original Message- > From: Michal Figurski [mailto:[EMAIL PROTECTED] > Sent: Tuesday, July 22, 2008 10:44 AM > To: Doran, Harold; r-help@r-project.org > Subject:

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Hmm... It sounds like ideology to me. I was asking for technical help. I know what I want to do, just don't know how to do it in R. I'll go back to SAS then. Thank you. -- Michal J. Figurski Doran, Harold wrote: I think the answer has been given to you. If you want to continue to ignore tha

  1   2   >