Thanks! I made this correction
ONKELINX, Thierry wrote: > > > Please note that the case of characters is important in R. Hence "True" > differs from "TRUE". You need to use REML = TRUE instead of REML = True! > The only valid logical values are TRUE, FALSE, T and F. But I recommend > to use only TRUE and FALSE. Because you can overwrite T and F. You could > do something like T <- FALSE, which is ofcourse not a good idea. TRUE > and FALSE can not be overwritten. > > HTH, > > Thierry > > > ------------------------------------------------------------------------ > ---- > ir. Thierry Onkelinx > Instituut voor natuur- en bosonderzoek / Research Institute for Nature > and Forest > Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, > methodology and quality assurance > Gaverstraat 4 > 9500 Geraardsbergen > Belgium > tel. + 32 54/436 185 > thierry.onkel...@inbo.be > www.inbo.be > > To call in the statistician after the experiment is done may be no more > than asking him to perform a post-mortem examination: he may be able to > say what the experiment died of. > ~ Sir Ronald Aylmer Fisher > > The plural of anecdote is not data. > ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answer can be extracted from a given body of > data. > ~ John Tukey > > -----Oorspronkelijk bericht----- > Van: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > Namens boeinguy2 > Verzonden: donderdag 7 mei 2009 0:21 > Aan: r-help@r-project.org > Onderwerp: [R] Help with lme4 model specification > > > I am new to R and am trying to specify a model for mixed model analysis. > > When I run the following model I get an error: > > AAT<- lmer(Y ~ S + A + (1|S:A/H), data=AT, REML=True) > > The error looks like this: > > Error in Spar_loc:`:` : NA/NaN argument > In addition: Warning messages: > 1: In model.matrix.default(mt, mf, contrasts) : > variable 'Spar_loc' converted to a factor > 2: In Spar_loc:`:` : > numerical expression has 720 elements: only the first used > > I am having trouble specifying th random component. It should reflect > the > random term as H nested within the interaction of S & A. What am I > doing > wrong? > :confused: > -- > View this message in context: > http://www.nabble.com/Help-with-lme4-model-specification-tp23415908p2341 > 5908.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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, minimal, self-contained, reproducible code. > > Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver > weer > en binden het INBO onder geen enkel beding, zolang dit bericht niet > bevestigd is > door een geldig ondertekend document. The views expressed in this message > and any annex are purely those of the writer and may not be regarded as > stating > an official position of INBO, as long as the message is not confirmed by a > duly > signed document. > > ______________________________________________ > 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, minimal, self-contained, reproducible code. > > -- View this message in context: http://www.nabble.com/Help-with-lme4-model-specification-tp23415908p23428528.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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, minimal, self-contained, reproducible code.