Displaying 2 results from an estimated 2 matches for "treatp".
Did you mean:
treat
2010 Nov 16
1
glmer, Error: Downdated X'X is not positive definite,49
Dear list,
I am new to this list and I am new to the world of R. Additionally I am not
very firm in statistics either but have to deal. So here is my problem:
I have a dataset (which I attach at the end of the post) with a binomial
response variable (alive or not) and three fixed factors
(trapping,treat,sex). I do have repeated measures and would like to include
one (enclosure) random factor.
I
2005 Dec 14
3
Memory shortage running Repeated Measures (nlme)
...location of 254Mb: see help(memory.size)".
here is part of the log:
***********************************************************
> aphids<-read.table("aphid.txt",header=T)
> attach(aphids)
> names(aphids)
[1] "site" "time" "shade" "treatp" "aphid"
> library(nlme)
> aphids<-groupedData(aphid~time|site/shade/treatp,data=aphids,outer=~shade*treatp)
> model<-lme(aphid~shade*treatp,random=~time|site)
Error in logLik.lmeStructInt(lmeSt, lmePars) :
Calloc could not allocate (89359209 of 8) memory
In ad...