Re: [R] Lda and Qda

?lda explains the object produced: please do study it. Hint: you asked for leave-one-out cross-validation, and what is the output from cross-validation of a classifer? The predicted class for each observation. How many observations do you have? You are using software from a contributed

Re: [R] warning on gamma option in par(args) or calling par(= new)?

You are calling par() *before* opening and (in on.exit) *after* closing the pdf() device. par() applies on a per-device basis, and if no device is open the default device will be opened. As for why you get the warnings about 'gamma' and 'new', see ?par and read their entries. Since you

Re: [R] How to catch data from the different dataframes and lm problem?

1) I have a dataframe. Based on “formulation” and “subject”, a dataframe is split into 4 dataframes. The example is as follows. Moreover, I want to calculate “test” value for these 4 dataframes. My question is that the “test” values not correct and I do not know where the problem is.

Re: [R] Conditionally incrementing a loop counter: Take 2

Mike Jones wrote: Hi, I am trying a for loop from 1 to 10 by 1. However, if a condition does not get met, I want to throw away that iteration. So if my loop is for (i in 1:10) and i is say, 4 and the condition is not met then I don't want i to go up to 5. Is there a way to do that? I can't

[R] Return Value of TCl/Tk window in R

Hello, I have the TCl/Tk command tkmessageBox(titel=,message=x,icon=question,type=okcancel) in my R script. Now I want to perform some operation in relation to the user's choice, something like if (okpressed) xxx else yyy What values does this command give and how are they used? Thank you,

Re: [R] Return Value of TCl/Tk window in R

Is it so hard to find out? Your tcl documentation will tell you what tk_messageBox returns, and as it is tcl string, you need to call tclvalue() on the value of tkmessageBox() to get an R character vector. On Fri, 28 Dec 2007, Richard Müller wrote: Hello, I have the TCl/Tk command

Re: [R] Return Value of TCl/Tk window in R

Richard Müller wrote: I have the TCl/Tk command tkmessageBox(titel=,message=x,icon=question,type=okcancel) in my R script. Now I want to perform some operation in relation to the user's choice, something like if (okpressed) xxx else yyy What values does this command give and how are they

Re: [R] Return Value of TCl/Tk window in R

This webpage has been very helpful for me: http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/ -Samu Richard Müller kirjoitti: Hello, I have the TCl/Tk command tkmessageBox(titel=,message=x,icon=question,type=okcancel) in my R script. Now I want to perform some operation in relation to

[R] index question

I was hoping for some advice regarding indexing, From a dataframe there are 27 variables of interest, with the prefix of pre. [7] Decision MHCDate pre01 pre0 pre012pre013 [13] pre02 pre02111 pre02114 pre0211 pre0212 pre029 [19] pre03apre0311 pre0312 pre03

Re: [R] index question

From a dataframe there are 27 variables of interest, with the prefix of pre. [7] Decision MHCDate pre01 pre0 pre012 pre013 [13] pre02 pre02111 pre02114 pre0211 pre0212 pre029 [19] pre03apre0311 pre0312 pre03 pre04 pre05 [25] pre06 pre07 pre08

[R] encoding question again

Hi, I'm running the actual R version in JGR (version 1.5-8 ). Sys.getlocale(category = LC_ALL) yields [1] LC_COLLATE=German_Germany.1252;LC_CTYPE=German_Germany.1252;LC_MONETARY=German_Germany.1252;LC_NUMERIC=C;LC_TIME=German_Germany.1252 I want to write some HTML-Code enhanced by

[R] Accesing the value

Hi R, x=A A=5 I need to get the value of A using x only. How do I do this? Thanks in advance, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}} __ R-help@r-project.org mailing list

Re: [R] Accesing the value

get(x) This is a FAQ: http://cran.cnr.berkeley.edu/doc/FAQ/R-FAQ.html#How-can-I-turn-a-string-into-a-variable_003f --sundar Shubha Vishwanath Karanth said the following on 12/28/2007 8:38 AM: Hi R, x=A A=5 I need to get the value of A using x only. How do I do this?

Re: [R] SAS to R - if you don't have a SAS license

Wensui Liu wrote: while I move data between SAS and R all the time, personally I don't think your recommendation is very practical. Instead, I feel SAS transport file is much better than csv. Plus, the sas dataset created on unix can't be opened by sas viewer on windows. It is even undoable

Re: [R] SAS to R - if you don't have a SAS license

On Fri, 2007-12-28 at 11:21 -0600, Frank E Harrell Jr wrote: Wensui Liu wrote: while I move data between SAS and R all the time, personally I don't think your recommendation is very practical. Instead, I feel SAS transport file is much better than csv. Plus, the sas dataset created on

[R] Pause loop

Hi list How can I pause a loop (or any other function) for a defined time? E.g. for(i in 1:5) {print(1:i); 'function to pause for 10 seconds'} Thanks for help. -- Armin Goralczyk, M.D. -- Universitätsmedizin Göttingen Abteilung Allgemein- und Viszeralchirurgie Rudolf-Koch-Str. 40 39099

Re: [R] SAS to R - if you don't have a SAS license

Marc Schwartz wrote: On Fri, 2007-12-28 at 11:21 -0600, Frank E Harrell Jr wrote: Wensui Liu wrote: while I move data between SAS and R all the time, personally I don't think your recommendation is very practical. Instead, I feel SAS transport file is much better than csv. Plus, the sas

Re: [R] Pause loop

On Fri, 2007-12-28 at 19:07 +0100, Armin Goralczyk wrote: Hi list How can I pause a loop (or any other function) for a defined time? E.g. for(i in 1:5) {print(1:i); 'function to pause for 10 seconds'} Thanks for help. At least on Linux, see: ?Sys.sleep Note the caveats regarding

[R] unit attribute to list elements

Hi, I've started my own (first) package, part of which consists in listing common physical constants (Planck's constant, the speed of light in vacuum, etc). I'm wondering what would be a good way of dealing with pairs of value/unit. constants - list( cel = 2.99792458e8 , #m/s Z0 =

Re: [R] logistic mixed effects models with lmer

On Dec 28, 2007 9:35 AM, Sharon Goldwater [EMAIL PROTECTED] wrote: I have a question about some strange results I get when using lmer to build a logistic mixed effects model. I have a data set of about 30k points, and I'm trying to do backwards selection to reduce the number of fixed effects

[R] SAS to R - if you have SAS 8.2+

Hi there, the attached R function uses the SAS Integrated Object Model (IOM) and it can deal with SAS dates and long variable names. All you need to provide is the folder where the SAS data file is and the data file name without the extension. The function requires the rcom package. This is

[R] AIC and BIC in mixed effects model

Dear R Users: I am trying to compare several structures of the within-patient covariance such as unstructured, Autoregressive, and spatial by using the MIXED effects model. Can AIC, BIC be negative ? If yes, then in what situations they may be negative. Thanks a lot. [[alternative HTML

Re: [R] Return Value of TCl/Tk window in R

res - tkmessageBox(title = test,message = Continue?, +icon =question, type = okcancel) if (tclvalue(res) == ok) 1 else 2 Happy new year! Philippe Grosjean Richard Müller wrote: Hello, I have the TCl/Tk command tkmessageBox(titel=,message=x,icon=question,type=okcancel) in my R

[R] gplot function - error

Hi all, i installed the ggplot2 package via install.packages(), but the gplot function was not recognized in R console command. Is there any paths to configure? The error message reports that the function was not found. Thanks Ricardo

Re: [R] Conditionally incrementing a loop counter: Take 2

Dear Mike, You could use a repeat loop and manage the index yourself: i - 0 repeat{ x - runif(1) if (x .1){ i - i + 1 cat(x = , x, \n) } if (i == 10) break } But if your example problem reflects your actual application, why not just

Re: [R] groupedData function not found

Hi Andrea, did you try require(nlme) ? Adding a package is not necessarily the same as loading it. Cheers Andrew On Thu, Dec 27, 2007 at 10:42:27PM +, andrea previtali wrote: Hello, I'm trying to use the groupedData function and R is giving me the message: Error: can not find

Re: [R] unit attribute to list elements

Is this what you want? Lines - cel 3.0 m/s + Z0 367 ohm + eps0 8.9e-12 F/m + Constants.DF - read.table(textConnection(Lines), as.is = TRUE) Constants - as.list(Constants.DF[[2]]) names(Constants) - Constants.DF[[1]] for(i in seq_along(Constants)) comment(Constants[[i]]) - Constants.DF[i, 3]

Re: [R] unit attribute to list elements

On 28 Dec 2007, at 20:06, Gabor Grothendieck wrote: Is this what you want? Perfect! Thanks a lot! Lines - cel 3.0 m/s + Z0 367 ohm + eps0 8.9e-12 F/m + Constants.DF - read.table(textConnection(Lines), as.is = TRUE) Constants - as.list(Constants.DF[[2]]) names(Constants) -

Re: [R] gplot function - error

You changed the intial letter: the function is qplot and not gplot. domenico Ricardo Perrone wrote: Hi all, i installed the ggplot2 package via install.packages(), but the gplot function was not recognized in R console command. Is there any paths to configure? The error message reports

[R] Forcing virtual digits

Hi there, I'm using XP with R 2.6.1 I've got the following question: Is there a way to force virtual digits? I mean: In the external csv-file, I read in some numbers, e.g. 33 which is written into a variable. I want to add this variable into a plot with an additional digit: 33.0.

[R] two plots on the same page

I'd like to know why I cannot get a plot and the QQnorm in the same sheet. The commands are simple but: library(nlme) glmod1 - gls(upfmla,correlation=corAR1(),method=ML) summary(glmod1) par(mfrow = c(2,1)) plot(glmod1, main=GLS Residuals vs. GLS Fitted) qqnorm(glmod1) No matter what (I

[R] print ggplot2 summary on graphic

Hi all: The code below works well to print annotations on base graphics but I can't make it work using the ggplot2 package. I tried to use the capture.output function but only prints the summary on the R console and no on the graphic. Any ideas on how to print output on ggplot2 graphs? For a

[R] FYI: Package installation problem with windows

Hello, maybe this info is helpful for someone else, at least it is a reminder for me if the problem should reoccur: The last hours R drove me nuts because it wouldn't install new packages, or more irritating still: It would chose randomly which packages to install when given a list (unable to

Re: [R] AIC and BIC in mixed effects model

Chuanjun Zhang zhangchu at umkc.edu writes: Dear R Users: I am trying to compare several structures of the within-patient covariance such as unstructured, Autoregressive, and spatial by using the MIXED effects model. Can AIC, BIC be negative ? If yes, then in what situations they may be

[R] How can I change the language back to english?

I have reset my windows and re-installed R (I explicitly asked for english) and although the windows XP version is in english I am having R display in french. like: R est un logiciel libre livré sans AUCUNE GARANTIE. Vous pouvez le redistribuer sous certaines conditions. Tapez 'license()' ou

Re: [R] two plots on the same page

-Original Message- From: [EMAIL PROTECTED] on behalf of Maura E Monville Sent: Fri 12/28/2007 2:18 PM To: [EMAIL PROTECTED] Subject: [R] two plots on the same page I'd like to know why I cannot get a plot and the QQnorm in the same sheet. The commands are simple but:

Re: [R] Forcing virtual digits

?sprintf x - 33 plot(0, main=sprintf(Value = %.1f, x)) On Dec 28, 2007 5:13 PM, Thomas Schwander [EMAIL PROTECTED] wrote: Hi there, I'm using XP with R 2.6.1 I've got the following question: Is there a way to force virtual digits? I mean: In the external csv-file, I read in some

Re: [R] FYI: Package installation problem with windows

On Fri, 28 Dec 2007, Owe Jessen wrote: Hello, maybe this info is helpful for someone else, at least it is a reminder for me if the problem should reoccur: The last hours R drove me nuts because it wouldn't install new packages, or more irritating still: It would chose randomly which

Re: [R] How can I change the language back to english?

This is rw-FAQ Q3.2, Q3.4. But the underlying problem is that Windows is telling R it is in French, and perhaps you need to sort that out. On Fri, 28 Dec 2007, francogrex wrote: I have reset my windows and re-installed R (I explicitly asked for english) and although the windows XP version