similar to: Variance component estimation in glmmPQL

Dear R users, I have been using the pdMat class "pdIdnot" (from the mgcv package)instead of "pdIdent" to avoid overflow in GLMM fits with the MASS package function glmmPQL, of the following form: fit1 <- glmmPQL(fixed=y0~-1+xx0, random=list(gp=pdIdent(~-1+zz0)), family=binomial) # vulnerable to overflow fit2 <- glmmPQL(fixed=y0~-1+xx0,

Hi all, I have an experiment where 5 raters assessed the quality of 24 web sites. (each rater rated each site once). I want to come up with a measure of reliability of the ratings for the web sites ie to what extent does each rater give the same (or similar) rating to each web site. My idea was to fit a random effects model using lme and from that, calculate the intraclass correlation as a

Hi all - R2.0.1, OS X Perhaps while there is some discussion of lme going on..... I am trying to execute a glmm using glmmPQL from the MASS libray, using the example data set from McCullagh and Nelder's (1989, p442) table 14.4 (it happens to be the glmm example for GENSTAT as well). The data are binary, representing mating success (1,0) for crosses between males and females from two

Hi, all, I have a question about random effects model. I am dealing with a three-factor experiment dataset. The response variable y is modeled against three factors: Samples, Operators, and Runs. The experimental design is as follow: 4 samples were randomly chosen from a large pool of test samples. Each of the 4 samples was analyzed by 4 operators, randomly selected from a group of

Dear All, I would like to ask you to help me with my memeory. I remember using some function that would list all the possible methods I could apply to an object. Say, if I had an object of class=lme, it would tell me that that I could do stuff like qqnorm(myobjct), or VarCorr(myobject). In general, a very complete list. I though this list of all possible methods would pop out by typing

To those who might be interested -- following is the solution to my previous post regarding reproducing output from SAS Proc Mixed for a two-factor crossed random effects ANOVA model. I am graciously endebted to the kind replys from two statisticians for this solution whose names I will refrain from mentioning for the sake of privacy. I hope this helps someone?! -- Phil Turk > hw7 <-

Douglas Bates wrote: > > The GLMM function in the lme4 package allows you to specify crossed > random effects within the random argument without the need for the > pdBlocked and pdIdent constructions. Simply ensure that your grouping > factors are defined in such a way that each distinct group has a > different level in the grouping factor (this is usually not a problem

Hi R helpers, I'm trying to fit a rather complex model to some simulated data using lme and am not getting the correct results. It seems there might be some identifiability issues that could possibly be dealt with by specifying starting parameters - but I can't see how to do this. I'm comparing results from R to those got when using GenStat... The raw data are available on the

hello! this is a question, how can i specify the random part in the GLMM-call (of the lme4 library) for compound matrices just in the the same way as they defined in the lme-Call (of the nlme library). For example i would just need random=list(my.Subject=pdBlocked(list(pdIdent(~... , ...),pdIdent(~... , ...)))) this specification , if i also attach library(nlme) , is not

hello! this is a question, how can i specify the random part in the GLMM-call (of the lme4 library) for compound matrices just in the the same way as they defined in the lme-Call (of the nlme library). For example i would just need random=list(my.Subject=pdBlocked(list(pdIdent(~... , ...),pdIdent(~... , ...)))) this specification , if i also attach library(nlme) , is not

Consider the Assay data where block, sample within block and dilut within block is random. This model can be fitted with (where I define Assay2 to get an ordinary data frame rather than a grouped data object): Assay2 <- as.data.frame(Assay) fm2<-lme(logDens~sample*dilut, data=Assay2, random=list(Block = pdBlocked(list(pdIdent(~1), pdIdent(~sample-1),pdIdent(~dilut-1))) )) Now, block

Dear all, I am trying to fit a model with crossed random effects using lme. In this experiment, I have been measuring oxygen consumption (mlmin) in bird nestlings, originating from three different treatments (treat), in a respirometer with 7 different channels (ch). I have also measured body mass (mass) for these birds. id nest treat year mlmin mass ch hack 1EP51711 17

Hello, I work with a mixed model with 4 predictor variables Time, Size, Charge, Density and Size, Charge, Density are factors, all with two levels. Hence I want to put their interactions with Time into the model. But, I have two data sets (Replication 1 and 2) and I want that Replication is random effect. Here is my code: knots <- default.knots(Time) z <- outer(Time, knots, "-")

Hi, I have a data set on germination and plant growth with the following variables: dataset=fm mass (response) sub (fixed effect) moist (fixed effect) pop (fixed effect) mum (random effect nested within population) iheight (covariate) plot (random effect- whole plot factor for split-plot design). I want to see if moist or sub interacts with mum for any of the pops, but I am getting an error

I have been using lme from nlme to do a 3-way anova with all the effects treated as random. I was wondering if someone could direct me to an example of how to do this using lmer from lme4. I have 3 main effects, tim, trt, ctr, and all the interaction effects tim*trt*ctr. The response variable is ge. Here is my lme code: dat <-

Try this data(Assay) as1 <- lme(logDens~sample*dilut, data=Assay, random=pdBlocked(list( pdIdent(~1), pdIdent(~sample-1), pdIdent(~dilut-1)))) update(as1,random=pdCompSymm(~sample-1)) update(as1,random=pdCompSymm(~sample-1)) update(as1,random=pdCompSymm(~sample-1)) update(as1,random=pdCompSymm(~sample-1)) I'm

Hello, I have used nlme to fit a model, the R syntax is like fmla0<-as.formula(paste("~",paste(colnames(ldata[,9:13]),collapse="+"),"-1")) > fmla1<-as.formula(paste("~",paste(colnames(ldata[,14:18]),collapse="+"),"-1")) >

Dear R-help group, How can I define a lme with 3 factors(a,b,c), where c is nested in b, and a is crossed with b/c? I think that: lme(response ~ ..., data = Data, random = pdBlocked(list(pdIdent(~ a - 1), pdIdent(~ b - 1)))) is one part of the answer and: lme(response~..., data=Data, random=~1|b/c) is the other part of the answer but how can I combine them?? Could anybody please help

Hello. I have some trouble with mixed effects in R, similar to problems that other people had with not nested models and lme, as I understand from the mailing list archive. Unfortunately, I could not understand the solutions that were proposed... I have a data set with response variable (y) and two explanatory variables x1 and x2 (x1 - fixed factor, x2 - random factor). Fixed factor x1 is

---------- Forwarded message ---------- From: Freddy Gamma <freddy.gamma@gmail.com> Date: 2011/1/21 Subject: TRADUCING lmer() syntax into lme() To: r-sig-mixed-models@r-project.org Dear Rsociety, I'd like to kingly ask to anyone is willing to answer me how to implement a NON NESTED random effects structure in lme() In particular I've tried the following translation from lmer to