Hello, I have the following design, counts were collected at different transects, different depths and different sites at different times. Time is continuous and assumed to be random, all the others are categorical fixed where transect is nested within depth which is nested within site. I would like an expert opinion about the following code where intercept is modeled as random (I am not sure if this is the right way). a<-factor(transect) b<-factor(depth) c<-factor(site) g<-lmer(count~(1|time)+(time|a:b:c), family="poisson") Thanks, Y [[alternative HTML version deleted]]
This should definitely be posted on the r-sig-mixed-models list, not here. -- Bert On Thu, Jul 19, 2012 at 6:16 AM, Yolande Tra <yolande.tra@gmail.com> wrote:> Hello, > > I have the following design, counts were collected at different transects, > different depths and different sites at different times. Time is continuous > and assumed to be random, all the others are categorical fixed where > transect is nested within depth which is nested within site. > > I would like an expert opinion about the following code where intercept is > modeled as random (I am not sure if this is the right way). > > a<-factor(transect) > b<-factor(depth) > c<-factor(site) > g<-lmer(count~(1|time)+(time|a:b:c), family="poisson") > > Thanks, > Y > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >-- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm [[alternative HTML version deleted]]
Maybe Matching Threads
- Error in as.xts
- loop searching the id corresponding to the given index (timestamp)
- lmer: Estimated variance-covariance is singular, false convergence
- readRDS, In as.double.xts(fishReport$count) : NAs introduced by coercion
- loop for inserting rows in a matrix