[R] mixed models

[Jeff Newmiller]

Please keep the mailing list in the loop by using reply-all.

I don't think there is a requirement that the number of levels is equal, but there may be problems if you don't have the minimum number of records corresponding to each combination of levels specified in your model. 

You can change the csv extension to txt and attach for the mailing list. Or, better yet, you can use the dput function to embed the data directly in your sample code. 

Also,  please learn to post plain text email to avoid corruption of R code by the HTML formatting. 
-- 
Sent from my phone. Please excuse my brevity.

On May 25, 2016 2:26:54 PM PDT, James Henson <jfhenson1 at gmail.com> wrote:
>Good afternoon Jeff,
>The sample sizes for levels of the factor "Irrigation" are not equal.
>If
>'nlme' requires equal sample sizes this may be the problem. The same
>data
>frame runs in 'lme4' without a problem.
>
>Best regards,
>James
>
>
>On Wed, May 25, 2016 at 3:41 PM, James Henson <jfhenson1 at gmail.com>
>wrote:
>
>> Good afternoon Jeff,
>>
>> When working with this data frame, I just open the .csv file in R
>Studio.
>> But, we should not send .csv file to R_help.  What should I send?
>>
>> Best regards,
>> James
>>
>> On Wed, May 25, 2016 at 2:52 PM, Jeff Newmiller
><jdnewmil at dcn.davis.ca.us>
>> wrote:
>>
>>> You forgot to show the commands to us that you used to read the data
>in
>>> with (your example is not "reproducible"). This step can make all
>the
>>> difference in the world as to whether your analysis commands will
>work or
>>> not.
>>> --
>>> Sent from my phone. Please excuse my brevity.
>>>
>>> On May 25, 2016 11:59:06 AM PDT, James Henson <jfhenson1 at gmail.com>
>>> wrote:
>>>
>>>> Greetings R community,
>>>>
>>>> My aim is to analyze a mixed-effects model with temporal
>pseudo-replication
>>>> (repeated measures on the same experimental unit) using ‘nlme’. 
>However,
>>>> my code returns the error message “Error in na.fail.default’, even
>though
>>>> the data frame does not contain missing values. My code is below,
>and the
>>>> data file is attached as ‘Eboni2.txt.
>>>>
>>>> library("nlme")
>>>>
>>>> str(Eboni2)
>>>>
>>>> head(Eboni2)
>>>>
>>>> model1 <- lme(preDawn ~ Irrigation, random=~season_order|treeNo,
>>>> data=Eboni2)
>>>>
>>>> I am genuinely confused.  Hope someone can help.
>>>>
>>>> Best regards,
>>>>
>>>> James F. Henson
>>>>
>>>> ------------------------------
>>>>
>>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>>> 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.
>>>>
>>>>
>>

	[[alternative HTML version deleted]]