[R] Autogenerate tags in tag=value pairs for

2007-03-06 Thread Jon Olav Vik
Dear list, Is there a way to programmatically specify tag names for the ... (ellipsis) part of the argument list to a function? In other words, a way to do this: x - data.frame(A=1:5) if the name A was not hardcoded but given by a variable, and without resorting to: x - data.frame(1:5)

[R] How to test for significance of random effects?

2006-05-06 Thread Jon Olav Vik
36.50451 .0001 ## Various log likelihoods: logLik(lm1,REML=FALSE) logLik(lm1,REML=TRUE) logLik(lme1,REML=FALSE) logLik(lme1,REML=TRUE) Any help is highly appreciated. Best regards, Jon Olav Vik __ R-help@stat.math.ethz.ch mailing list https

[R] Robust measures of goodness of fit?

2006-03-30 Thread Jon Olav Vik
as possibly contaminated by some other phenomenon than the one I'm interested in. What I think I want is a robust estimator of the R-squared which would be obtained in the absence of such contamination. Sincerely, Jon Olav Vik __ R-help@stat.math.ethz.ch

[R] Why is AIC from lmer 2 less than that from lme?

2006-03-28 Thread Jon Olav Vik
on Windows XP SP2 (details below), lmer() from Matrix 0.995-8, and lme() from nlme 3.1-68.1 Thank you in advance for your help. Sincerely, Jon Olav Vik library(nlme) summary(lme(travel~1,data=Rail,random=~1|Rail)) Linear mixed-effects model fit by REML Data: Rail AIC BIC