similar to: glmmPQL questions

Hi, I have large data frame with many columns. A short example is given below: > dataH host ms01 ms31 ms33 ms34 1 cattle 4 20 9 6 2 sheep 4 3 4 5 3 cattle 4 3 4 5 4 cattle 4 3 4 5 5 sheep 4 3 5 5 6 goat 4 3 4 5 7 sheep 4 3 5 5 8 goat 4 3 4 5 9 goat 4 3 4 5 10 cattle

Dear group, Here is my data frame: futures <- structure(list(DESCRIPTION = c("CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "LIVE CATTLE Aug/10", "LIVE CATTLE Aug/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11

Dear group, I have these 2 following data frame: allcon <- structure(list(DESCRIPTION = structure(1:17, .Label = c("COFFEE C Jul/10", "COPPER May/10", "CORN Jul/10", "CORN May/10", "COTTON NO.2 Jul/10", "CRUDE OIL miNY May/10", "GOLD Jun/10", "HENRY HUB NATURAL GAS May/10", "ROBUSTA COFFEE (10)

Dear group, Here is the kind of data.frame I obtain every day with my function : futures <- structure(list(DESCRIPTION = c("CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "CORN Jul/10", "LIVE CATTLE Aug/10", "LIVE CATTLE Aug/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 Jul/10", "SUGAR

Hi all, I'm having trouble with the model generating 'dredge' function in the MuMIn 'Multi-model Inference' package. Here's the script: globalmodel<- glm(TB~lat+protocol+tested+ streams+goats+hay+cattle+deer, family="binomial") chat<- deviance(globalmodel)/59 #There we 59 residual degrees of freedom in this global model. models<- dredge(globalmodel,

Hi I discovered a similar problem early this week and a colleage on irc was told it was a copy and paste error. in man syslogd, the following can be found: Sucker rod def. -- 3/4, 7/8 or 1in. hardened steel rod, male threaded on each end. Primary use in the oil industry in Western North Dakota and other locations to pump 'suck' oil from oil wells. Sec- ondary

On 2/22/23 19:20, Andrey Drobyshev wrote: > Since commits b28cd1dc ("Remove requested_guestcaps / rcaps"), f0afc439 > ("Remove guestcaps_block_type Virtio_SCSI") support for installing > virtio-scsi driver is missing in virt-v2v. AFAIU plans and demands for > bringing this feature back have been out there for a while. E.g. I've > found a corresponding issue

Hi there I have created a multivariate logistic regression model looking at the presence/absence of disease on farms. I would like to plot the diagnostic plots recommended by Hosmer & Lemeshow to look particularly for any points of high influence. In order to do this I need to extract values for delta-beta. The function dfbeta gives a value for change in each coefficient but I am looking

http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the method in Mohsen's paper,but I found I really can get the same or close output as he did,so,any

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

could someone please give me the lowdown on samba and passwords? i keep trying to connect to any share name on either of my unix boxes, either from each other or from a win95 machine... every time i can't get anything due to passwords... does samba have it's own password file? from trial-and-error it seems that it doesn't read out your /etc/passwd file for user's home dirs...

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,