[R] LINEAR MIXED EFFECT

[ANARPCG]

Milton's point is dead-on, and I would highly encourage you to give the
posting guide a look.

That said... you might try "na.action = na.omit" in your call to...
actually, we don't know what function you are using (see first point). 
Regardless, sounds like you have missing data and na.action is set to
na.fail (ie, "fail" if any missing data).

cheers, Dave


milton ruser wrote:
> 
> Dear Ana Golveia,
> 
> It is completelly impossible someone realise what kind or help you need or
> what is happening. I suggest you give a look on the posting guide, mainly
> that part about a minimum reproducible code with self explaining
> information, etc.
> 
> Cheers
> 
> milton
> 
> On Wed, Nov 11, 2009 at 7:22 AM, ANARPCG <a.gouveia07 at imperial.ac.uk>
> wrote:
> 
>>
>> CAN ANYONE PLEASE HELP ME WITH THIS
>> i HAVE TO DO A MIXED EFFECT LINEAR MODEL WITH MY DATA DUE TO THE FACT
>> THAT
>> I
>> have pseudoreplication!
>> Although after reading and trying it for several times can get around due
>> to
>> "Error in na.fail.default(list(date = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 
>> :
>> missing values in object"
>> I uploaded my data file
>> Thank you so much
>> Kind regards
>> AG
>> http://old.nabble.com/file/p26300394/rawoctobercalciumexperiment2.txt
>> rawoctobercalciumexperiment2.txt
>>
>> --
>> View this message in context:
>> http://old.nabble.com/LINEAR-MIXED-EFFECT-tp26300394p26300394.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> R-help at 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<http://www.r-project.org/posting-guide.html>
>> and provide commented, minimal, self-contained, reproducible code.
>>
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.
> 
> sorry I am new on this things 
> the point is after trying to run the model all that happens is:
>> exp<-read.table(file=file.choose(),header=T)
>> attach(exp)
>> names(exp)
> [1] "group"     "date"      "diet"      "weight"    "thickness" "length"   
> [7] "width"     "liplength" "lipwidth" 
>> exp$diet=factor(exp$diet,levels=c("zeropercent","tenpercent","twentypercent","thirtypercent","fortypercent","cuttleprecent"))
>> exp=na.omit(exp)
>> library(nlme)
>> random=~date
>> model<-lme(weight~date*diet)
> Error in na.fail.default(list(date = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L,  : 
>   missing values in object
> 
>  I have pseudoreplications due to the fact that the measurements of the
> replicates have different starting points so I was advised to do lme. I
> never used it before and I cant get arround with it! the help I wanted
> from you is to help me to understand how to do a "runable" model! 
> 
> 
> 
http://old.nabble.com/file/p26315302/rawoctobercalciumexperiment2.txt
rawoctobercalciumexperiment2.txt 
-- 
View this message in context: http://old.nabble.com/LINEAR-MIXED-EFFECT-tp26300394p26315302.html
Sent from the R help mailing list archive at Nabble.com.