Re: [R] Converting list of data frame to data frame

2007-04-29 Thread Ajit Pawar
Douglas/R-help, Thanks for your reply. I did try the solution but the result is not what I expect and I also get the following warning message: --- Warning message: number of columns of result is not a multiple of vector length (arg 1) in: rbind(1, c(6, 9, 10,

Re: [R] digest readability/ thunderbird email client/ R-help

2007-04-29 Thread Michael Dewey
At 21:00 28/04/2007, Martin Maechler wrote: Hi Sam, I posted a request for information on how various mailers handled this so I could summarise the results and put it on a helpful webpage. Nobody was willing to sing the praises of his/her favourite mailer so the only one I can suggest is the

Re: [R] Converting list of data frame to data frame

2007-04-29 Thread jim holtman
since you didn't supply a reproducible example, here is a test that I ran with varying number of rows that seems to work fine. You might want to provide at leat an str of the structure that you are doing the rbind on. x - lapply(1:10, function(z){ + data.frame(a=runif(z), b=rnorm(z),

Re: [R] randomForest gives different results for formula call v. x, y methods. Why?

2007-04-29 Thread Gavin Simpson
On Sat, 2007-04-28 at 21:13 -0400, David L. Van Brunt, Ph.D. wrote: Just out of curiosity, I took the default iris example in the RF helpfile... but seeing the admonition against using the formula interface for large data sets, I wanted to play around a bit to see how the various options

Re: [R] weight

2007-04-29 Thread Douglas Bates
On 4/28/07, John Kane [EMAIL PROTECTED] wrote: IIRC you have a yes/no smoking variable scored 1/2 ? It is possibly being read in as a factor not as an integer. try class(df$smoking.variable) to see . Good point. In general I would recommend using str(df) to check on the class or

[R] Too slow to execute!

2007-04-29 Thread Usman Shehu
Greetings, I have the following simple function but what worries me is that it takes about 5 or more minutes to execute. My machine runs on windows with 1.8GHz and 256 Ram. Re=NULL for(i in 1:10){ + x=rnorm(20) + Re[i]=(x-2*10)/20 + Re + } I would appreciate any help on how to make it

Re: [R] Too slow to execute!

2007-04-29 Thread Gavin Simpson
On Sun, 2007-04-29 at 06:56 -0700, Usman Shehu wrote: Greetings, I have the following simple function but what worries me is that it takes about 5 or more minutes to execute. My machine runs on windows with 1.8GHz and 256 Ram. Re=NULL for(i in 1:10){ + x=rnorm(20) +

Re: [R] Too slow to execute!

2007-04-29 Thread Chuck Cleland
Usman Shehu wrote: Greetings, I have the following simple function but what worries me is that it takes about 5 or more minutes to execute. My machine runs on windows with 1.8GHz and 256 Ram. Re=NULL for(i in 1:10){ + x=rnorm(20) + Re[i]=(x-2*10)/20 + Re + } I would appreciate

Re: [R] randomForest gives different results for formula call v. x, y methods. Why?

2007-04-29 Thread David L. Van Brunt, Ph.D.
In the words of Simpson (2007), D'OH! I knew it had to be something simple! On 4/29/07, Gavin Simpson [EMAIL PROTECTED] wrote: On Sat, 2007-04-28 at 21:13 -0400, David L. Van Brunt, Ph.D. wrote: Just out of curiosity, I took the default iris example in the RF helpfile... but seeing the

Re: [R] Example of mcmcsamp() failing with lmer() output

2007-04-29 Thread Douglas Bates
On 4/27/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi, I would appreciate help with the following model 1= gunload - read.table(hh('datasets/gunload.dat'), header = T) gunload$method - factor(gunload$method, labels = c('new', 'old')) gunload$physique - factor(gunload$group, labels =

[R] probl with optimization

2007-04-29 Thread fc502
hi everyone! i've got a problem in optimizing the following function fun - function(theta){ theta - rbind(beta,lambda) fun - sum(exp(h(beta)%*%lambda)) fun } where h(beta) is h - function(beta,...){ g1 - matrix(0,b,M) g2 - matrix(0,b,M) h1 - matrix(0,b,1) h2 - matrix(0,b,1) for(f in

Re: [R] Too slow to execute!

2007-04-29 Thread rolf
Chuck Cleland wrote: snip Why use a loop here at all? It seems like this would do it pretty fast: Re - rnorm(n = 2e+06, mean = - 1) system.time(Re - rnorm(n = 2e+06, mean = - 1)) [1] 0.77 0.01 0.78 NA NA I'm not sure what you expect Re to be (a list with 10 vectors

[R] thanks

2007-04-29 Thread Sam . McClatchie
Thanks Martin and Michael The plain digest option works well with Thunderbird. Best fishes Sam Sam McClatchie, Fisheries oceanographer Southwest Fisheries Science Center, NOAA, 8604 La Jolla Shores Drive La Jolla, CA 92037-1508, U.S.A. email [EMAIL PROTECTED] Cellular: 858 752 8495

Re: [R] probl with optimization

2007-04-29 Thread Prof Brian Ripley
On Sun, 29 Apr 2007, [EMAIL PROTECTED] wrote: hi everyone! i've got a problem in optimizing the following function fun - function(theta){ theta - rbind(beta,lambda) fun - sum(exp(h(beta)%*%lambda)) fun } Well, what are beta and lambda? Trying to optimize a function that does not

Re: [R] Too slow to execute!

2007-04-29 Thread jim holtman
Yet another way of doing it: system.time({ + x - replicate(10, (rnorm(20) - 20) / 20) + }) [1] 3.44 0.06 3.71 NA NA str(x) num [1:20, 1:10] -1.041 -0.980 -0.987 -1.062 -1.012 ... On 4/29/07, Usman Shehu [EMAIL PROTECTED] wrote: Greetings, I have the following simple

Re: [R] probl with optimization

2007-04-29 Thread fc502
ok...beta and lambda are the parameters that i ve to find. beta is a scalar and lambda is a 2x1 vector i ll have a look to the debugging tools, thnx! f On Apr 29 2007, Prof Brian Ripley wrote: On Sun, 29 Apr 2007, [EMAIL PROTECTED] wrote: hi everyone! i've got a problem in optimizing the

Re: [R] Greek symbols in xtable rows

2007-04-29 Thread Roger Bivand
On Tue, 17 Apr 2007, James W. MacDonald wrote: Roger Bivand wrote: Andris Jankevics andza at osi.lv writes: Thank You for answer. I Tried code that You provided in reply to my question on both PC's with Linux and Widows OS. Yes, I have also been bitten by the upgrade to xtable, which

Re: [R] The confidence level of p-value of ks.boot

2007-04-29 Thread Jasjeet Singh Sekhon
Hi Gala, The default p-value is the bootstrap p-value for the ks-test. Bootstrapping is highly recommended because the bootstrapped Kolmogorov-Smirnov test, unlike the standard test, provides correct coverage even when there are point masses in the distributions being compared. The bootstrap

[R] impossible to open SPSS file

2007-04-29 Thread Frank Thomas
Hi, I try to import a SPSS file with foreign but don't succeed. I looked into the help pages but find no answer. The file has been produced with SPSS 15 but I got the same error with a file in SPSS portable format. mydata-read.spss(C:\Rdata\ESS1_ICT.SAV) Erreur dans read.spss(C:RdataESS1_ICT.SAV)

Re: [R] impossible to open SPSS file

2007-04-29 Thread David Barron
Have a look at FAQ 2.16 in the R for Windows FAQ. On 29/04/07, Frank Thomas [EMAIL PROTECTED] wrote: Hi, I try to import a SPSS file with foreign but don't succeed. I looked into the help pages but find no answer. The file has been produced with SPSS 15 but I got the same error with a file in

[R] Going from Coda Files to R2WinBUGS

2007-04-29 Thread Jarrett Byrnes
I'm currently using JAGS as my Bayesian program of choice due to working off of an older mac running OSX. I'd like to utilize some of the functions from R2WinBUGS, however. As such, I'm attempting to write something to go from coda output dumped by JAGS into the bugs object format - I've

[R] adjusting a sigmoidal shaped regression

2007-04-29 Thread Milton Cezar Ribeiro
Hi R-friends How can I adjust a sigmondal shaped regression in R? My data looks like x-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16) y-c(10,11,11,14,15,35,42,49,80,120,130,138,149,150,151,150) plot(y~x) modlin-glm(y~x) yest-predict(modlin) lines(yest~x) thanks a lot Miltinho Brazil

[R] windows completion in 2.5?

2007-04-29 Thread David Hugh-Jones
Can anyone tell me if the new R 2.5 has autocompletion for the default windows shell? I know they integrated rcompletion but wasn't sure if it was just Unix. Cheers David __ R-help@stat.math.ethz.ch mailing list

Re: [R] Converting list of data frame to data frame

2007-04-29 Thread John Kane
I just tried an artifical example and the approach seems to work okay # Joining data.frames stored in a list # From Douglas Bates aa - 1:4 bb - 2:5 cc - 3:6 dd - 4:7 ee - 5:8 ff - c(rep(NA,4)) gg - 6:9 lst - list(data.frame(aa,bb), data.frame(cc,dd), data.frame(ee,ff,gg)) tatiana - names(lst)

Re: [R] weight

2007-04-29 Thread John Kane
--- Douglas Bates [EMAIL PROTECTED] wrote: On 4/28/07, John Kane [EMAIL PROTECTED] wrote: IIRC you have a yes/no smoking variable scored 1/2 ? It is possibly being read in as a factor not as an integer. try class(df$smoking.variable) to see . Good point. In general I would

Re: [R] windows completion in 2.5?

2007-04-29 Thread Duncan Murdoch
On 29/04/2007 8:06 PM, David Hugh-Jones wrote: Can anyone tell me if the new R 2.5 has autocompletion for the default windows shell? I know they integrated rcompletion but wasn't sure if it was just Unix. Yes, it's in Windows too. Hit tab and it will either auto-complete the current

Re: [R] adjusting a sigmoidal shaped regression

2007-04-29 Thread Gabor Grothendieck
Try this: library(drc) plot(y ~ x) mod - multdrc(y ~ x) lines(x, fitted(mod)) On 4/29/07, Milton Cezar Ribeiro [EMAIL PROTECTED] wrote: Hi R-friends How can I adjust a sigmondal shaped regression in R? My data looks like x-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16)

[R] A problem with svIDE in Tinn-R?

2007-04-29 Thread Kenneth Cabrera
Hi R users: I want to know if any one of you had the problem with Tinn-R, when you call the library svIDE on the new R 2.5.0, (because in the old R 2.4.1 works with out any problem). I got this message: library(svIDE) Loading required package: tcltk Loading Tcl/Tk interface ... done Warning

Re: [R] A problem with svIDE in Tinn-R?

2007-04-29 Thread apjaworski
Similar thing here. Happens in 2.5.0 patched and 2.6.0 (both April-25-2007) and svIDE 0.9-5 The only difference is that the warnings are: 1: '\A' is an unrecognized escape in a character string 2: unrecognized escape removed from ;for Options\AutoIndent: 0=Off, 1=follow language scoping and

[R] Quadratcount() plotting in R spatstat

2007-04-29 Thread Adrian Baddeley
John [EMAIL PROTECTED] writes: I am trying to plot a quadratcount object over a ppp object in the spatstat package. For queries about a particular package, please email the package authors. I ultimately want to get something like this

Re: [R] A problem with svIDE in Tinn-R?

2007-04-29 Thread Prof Brian Ripley
It is a warning, and a warning about a long-standing error in svIDE. This is not to do with Tinn-R except perhap that the message gets mangled. I get from command-line R 1: '\A' is an unrecognized escape in a character string 2: unrecognized escape removed from ;for Options\AutoIndent: 0=Off,