[R] Rounding error in seq(...)

2009-10-25 Thread David Airey
>> can you explain that a little more detailed? >> Perhaps I miss the background knowledge - but it seems just absurd to me. >> >> 0.1+0.1+0.1 is 0.3 - there is no rounding involved, is there? >> >> why is >> x <- 0.1 + 0.1 +0.1 >> not equal to >> y <- 0.3 > > Remember that this is in BINARY a

[R] variance components models with zero estimates

2008-06-07 Thread David Airey
When a variance components mixed model is run in Stata, if some of the variance components are zero, the model may not converge, for rational reasons. However, when the same model is run in SAS, the models with variance components that estimate to zero nonetheless converge. If I'm intereste

[R] latex{Hmisc} cannot find xdvi on MAC OS X

2008-01-24 Thread David Airey
. I also have an issue on the Mac getting latex(describe(iris)) to work. I get the .tex file but this error: > latex(describe(iris)) This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6) %&-line parsing enabled. entering extended mode (/var/folders/AJ/AJjttUSSFdSi1HiInJZyYTI/-Tmp-//Rtm

Re: [R] Two-way non-parametric ANOVA?

2008-01-23 Thread David Airey
. Just goes to show that Frank Harrell beats Google for statistics any day of the week! I have got to try that PO model!! On Jan 23, 2008, at 4:45 PM, Frank E Harrell Jr wrote: > David Airey wrote: >> While I defer to Frank as far as expertise goes, not having access >> t

Re: [R] Two-way non-parametric ANOVA?

2008-01-23 Thread David Airey
l Jr wrote: > David Airey wrote: >>> We need a two-way non-parametric ANOVA in order to analysis >>> properly some ecological data, do you know any reference in R? or >>> how to do it? Thank you very much All the best diana >>> >> A coup

Re: [R] Two-way non-parametric ANOVA?

2008-01-23 Thread David Airey
:) I was curious what you would say! Frank was kind enough to send this off list, but it's worth posting to the thread. On Jan 23, 2008, at 10:52 AM, Frank E Harrell Jr wrote: > David Airey wrote: >>> We need a two-way non-parametric ANOVA in order to analysis >>>

[R] Two-way non-parametric ANOVA?

2008-01-23 Thread David Airey
> We need a two-way non-parametric ANOVA in order to analysis properly > some ecological data, do you know any reference in R? or how to do > it? Thank you very much All the best diana > A couple more references here are below. I liked the Scheirer reference. /* Scheirer CJ, Ray WS, Hare N (1

Re: [R] extension to nlme self start SSmicmen?

2008-01-22 Thread David Airey
. Just to close the query, here are the data and my R code to use nlme to test for a strain difference in receptor binding using nlme and the self start function SSmicmen. Cheers, Dave "rownames","y","x","strain","pool" "1",1,54.0423393249512,0.169009417295456,"dba","5_10" "2",2,96.19451904

Re: [R] extension to nlme self start SSmicmen?

2008-01-22 Thread David Airey
> Dear list, > > Has anyone created a version of SSmicmen that allows testing for > group differences? The basic Michaelis-Menten equation is: > > (Bmax * X) / (Kd + X). > > The nlme package allows modeling of random effects for Bmax and Kd > as needed, but I curious how I can build in group di

[R] extension to nlme self start SSmicmen?

2008-01-22 Thread David Airey
Dear list, Has anyone created a version of SSmicmen that allows testing for group differences? The basic Michaelis-Menten equation is: (Bmax * X) / (Kd + X). The nlme package allows modeling of random effects for Bmax and Kd as needed, but I curious how I can build in group differences? I ha

[R] Input appreciated: R teaching idea + a way to improve R-wiki

2007-10-22 Thread David Airey
. I would love more quality online documentation around the same level as at UCLA (http://www.ats.ucla.edu/stat/; the Stata pages are fantastic), but I think I would like the first draft responsibility to fall to the well qualified instructor (Hi Matt!), and those with at least a Masters i

[R] saturation binding in nlme

2007-10-20 Thread David Airey
To estimate saturation binding parameters Bmax and Kd in a receptor saturation binding experiment, I use the following nonlinear equation and the nls() function: bmax*X*dummy + ns*X + background = total binding kd+X where X is concentration, and dummy is an indicator to allow sh