similar to: GLMMpql: reporting on main effects

Dear All, I have some code that works in S-Plus for writing saving a graphics file to disk :- graphsheet(type = "auto", format = "WMF", file = "G:\\north0l.wmf", pages = "auto", print.background = F, orientation="landscape", color.style="color") plot(x,y) dev.off() This works fine in S-Plus. I have tried playing

Dear all, I an running R on a Windows 2000 machine (1.5Gb RAM) and am trying to load the lme4 package, however, whenever I attempt to load the library "lme4" I get the following error message I have installed lme4 and Matrix locally from the zip file with no problems. R version 2.6.2 (2008-02-08) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0

Dear All, I am trying to put some Confidence intervals on some regressions from a linear model with no luck. I can extract the fitted values using 'predict', but am having difficulty in getting at the confidence intervals, or the standard errors. Any suggestions would be welcome Cheers Guy Using Version 2.1.0 (2005-04-18) on a PC vol.mod3 <-

Dear All, Using:- R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 and Jose Pinheiro, Douglas Bates, Saikat DebRoy and Deepayan Sarkar (2005). nlme: Linear and nonlinear mixed effects models. R package version 3.1-65. on a WINDOWS 2000 machine I am trying to run the scripts from the Mixed Models book and am running into some

Hello together, Is there a tool to test the statistical differences between parameter estimates of a nlsList fit? I fitted degradation data using the nlsList method and want to find out whether derived rate constants are significantly different from each other at the grouping factors soil and temperature. Regards Frank Scherr

Hi all. I have an issue that I cannot resolve. I am trying to read in lots of data that are stored in xml files. But after I read them in and copy the relevant data, then remove the document etc, it doesn't free up the memory. When I monitor it in windows task manager the memory usage just climbs with each iteration until R crashes. I can replicate the problem with the small example:

Hi, It feels like I should be able to do something like: R CMD INSTALL lib='/usr/lib64/R/library' repos='http://proxy.url/cran' package We have a bunch of servers (compute nodes in a Rocks cluster) in an isolated subnet, there is a basic pass-through proxy set up on the firewall (the head node) which just passes HTTP requests through to our nearest CRAN mirror. when using

I'm a bit puzzled by how to write out additive random effects in glmmPQL. In my situation, I have a factorial design on two (categorical) random factors, A and B. At each combination, I have a binary response, y, and two binary fixed covariates, C and D. If everything were fixed, I would use glm(y ~ A + B + C + D, family = binomial) My first thought was to use glmmPQL(y ~ A + B, random

Hello, I am using the R function glmmPQL to fit a logistic GLMM, with weights. I am finding that I get fairly different parameter estimates in glmmPQL from fitting the full dataset (with no "weight" statement) and an equivalent, shorter dataset with the weights statement. I am using the weights statement in the 'glmmPQL' function exactly as in the 'glm' function. I

Hi ! I try to compare some models obtained from glmmPQL. model1 <- glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 +I(freq8_4^2), random=~1|num, binomial); model2 <- glmmPQL(y~red*yellow+I(red^2)+I(yellow^2)+densite8+I(densite8^2)+freq8_4 , random=~1|num, binomial); anova(model1, model2) here is the answer : Erreur dans anova.lme(model1, model2) : Objects must

Dear all, for a cross-validation I have to use predict.glmmPQL() , where the formula of the corresponding glmmPQL call is not given explicitly, but constructed using as.formula. However, this does not work as expected: x1<-rnorm(100); x2<-rbinom(100,3,0.5); y<-rpois(100,2) mydata<-data.frame(x1,x2,y) library(MASS) # works as expected model1<-glmmPQL(y~x1, ~1 | factor(x2),

I have come across the previous communication on this list in September (copied below) because I had received the same error message. I understand from Brian Ripley's reply that anova should not be used with glmmPQL because it is not an adequate method, and that this is now shown with an error message. My question is, what method *should* be used? Using summary does not give me the result

Dear R-help, I recently tried implementing glmmPQL in 2.3.1, and I discovered a few differences as compared to 2.2.1. I am fitting a regression with fixed and random effects with Gamma error structure. First, 2.3.1 gives different estimates than 2.2.1, and 2.3.1, takes more iterations to converge. Second, when I try using the anova function it says, "'anova' is not available

Hi, In glmmPQL in the MASS library, the function uses repeated calls to the function lme(), using ML. Does anyone know how you can change this to REML? I know that in lme(), the default is actually set to REML and you can also specify this as 'method=REML' or 'method'ML' but this isn't applicable to glmmPQL(). I'd appreciate any help or advice! Thanks, Emma

hello all, i'm an R newbie struggling a bit with the glmmPQL function (from the nlme pack). i think i've managed to run the model successfully, but can't seem to plot the resulting function. plot(glmmPQL(...)) plots the residuals rather than the model... i know this should be basic, but i can't seem to figure it out. many thanks in advance. j -- View this message in context:

Hello, currently I am trying to fit a generalized linear mixed model using the glmmPQL function in the MASS package. I am working with the data provided by the book from Heck, Thomas and Tabata (2012) - https://www.routledge.com/Multilevel-Modeling-of-Categorical-Outcomes-Using-IBM-SPSS/Heck-Thomas-Tabata/p/book/9781848729568 I was wondering, which variance-covariance structure the glmmPQL

Hola! I am encountering the following problem, in a multilevel analysis, using glmmPQL from MASS. This occurs with bothj rw1062 and r-devel, respectively with nlme versions 3.1-38 and 3.1-39 (windows XP). > S817.mod1 <- glmmPQL( S817 ~ MIEMBROScat+S901+S902A+S923+URBRUR+REGION+ + S102+S103+S106A+S108+S110A+S109A+S202+S401+S557A+S557B+ + YHOGFcat,

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,

Dear all, I have two questions concerning model simplification in GlmmPQL, for for random and fixed effects: 1. Fixed effects: I don't know if I can simply specify anova(model) and trust the table that comes up with the p value for each variable in the fixed effects formula. I have read that the only way to test for fixed effects is to do approximate wald tests based on the standard errors

I fit the following model using glmmPQL from MASS: fit.glmmPQL <- glmmPQL(ifelse(class=="Disease",1,0)~age+x1+x2,random=~1|subject,family=binomial) summary(fit.glmmPQL) The response is paired (pairing denoted by subject), although some subjects only have one response. Also, there is a perfect positive correlation between the paired responses. x1 and x2 can and do differ within each