similar to: QIC for GEE

Hello, Is it possible to derive Pan's QIC (2001 Biometrics 57:120) from either a fitted gee() object in the gee package or from a geese() fit in the geepack package? If so, would anyone be kind enough to provide me with code to do so? I realize that QIC is part of the output from yags() but I would like to use one of the other functions. Thanks. Richard

Hi, I'm running GEE using geepack. I set corstr = "ar1" as below: > m.ar <- geeglm(L ~ O + A, + data = firstgrouptxt, id = id, + family = binomial, corstr = "ar1") > summary(m.ar) Call: geeglm(formula = L ~ O + A, family = binomial, data = firstgrouptxt, id = id, corstr = "ar1") Coefficients:

Hi, I'm analyzing my data using GEE, which looks like below: > interact <- geeglm(L ~ O + A + O:A, + data = data1, id = id, + family = binomial, corstr = "ar1") > summary(interact) Call: geeglm(formula = lateral ~ ontask + attachment + ontask:attachment, family = binomial, data = firstgroupnowalking, id = id, corstr = "ar1") Coefficients:

Dear R list members, I have a question regarding xyplot. I managed to make a xyplot of all the IDs by using the syntax below: xyplot(PA ~ CRPC + CRPT | ID, data = redinteract) Now, I'd like to make a graph of a specific ID number (e.g., only ID number 301). I thought I could use "subset", but it seems to be not working. Could anyone let me know how I can make a graph of a

Dear R list members, I'd like to make a graph of coefficients of the intercept, variable 1, and variable 2 (and possibly the interaction between variable 1 and variable 2). When I use the lmList function as attached below, it shows a nice coefficient graph. > PACRP.lis <- lmList(PAffect ~ CRPC + CRPT + CINT | ID, redinteract) > coef(PACRP.lis) > PACRPlis.coef <-

Hi, Probably this is a very simple question for all the programmers, but how do you change from 32-bit builds (default) to 64-bit builds? I've been trying to run Anova to compare two models, but I get the following error message: Error: cannot allocate vector of size 1.2 Gb R(3122,0xa0ab44e0) malloc: *** mmap(size=1337688064) failed (error code=12) *** error: can't allocate region ***

Hi, I'm actually trying to rank a set of candidate models with an information criterion (AICc, QIC, BIC). The problem I have is that I use mixed-effect cox regression only available with the package {coxme} (see the example below). #Model1 >spring.cox <- coxme (Surv(start, stop, Real_rand) ~ strata(Paired)+R4+R3+R2+(R3|Individual), spring) I've already found some explications in

Dear R users, Using R lme function, I found that both fixed and random effects of variable A on variable B are significant. Now, I'd like to analyze what variables are predicting differences in the slope. In other words, I'd like to know what variables (e.g., variable C) are predicting individual differences in the effects of A on B. I have many data points for A and B for each

Hi all, I've been running loglinear models for three-way tables: one of the variables having three levels, and the other two having two levels each. An example looks like below: > yes.no <- c("Yes","No") > switch <- c("On","Off") > att <- c("BB","AA","CC") > L <- gl(2,1,12,yes.no) > T <-

I'm using the package geepack to fit GEE models. Does anyone know of methods for add1 and drop1 for a 'geeglm' model object, or perhaps a method for extractAIC based on the QIC of Pan 2001? I see there has been some mention of this on R-help a few years ago (RSiteSearch("QIC")). The package does provide an anova method for its model objects, and update() seems to work:

I have included data at the bottom of this email. It can be read in by highlighting the data and then using this command: dat <- read.table("clipboard", header = TRUE,sep="\t") I can obtain solutions with both of these: library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <-

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

Hello R-users: I am using "yags" for fitting GEE which is giving me the same result as "Proc GENMOD". Now I have couple of questions related to yags output. (By the way, someone told me to run the geeglm for the same analysis and I did run but did not get the same result as of genmod and don't know how to correct the geeglm codes so that all three will be same!)

Hi , I use library(gee),library(geepack),library(yags) perform GEE data analysis , but all of them cannot do variable selection! Both step and stepAIC can do variable selection based on AIC criterion under linear regression and glm, but they cannot work when model is based on GEE. I want to ask whether any variable selection function or package under GEE model avaliable now? Thanks! Best,

--On Monday, May 11, 2015 02:26:17 PM -0700 John R Pierce <pierce at hogranch.com> wrote: > never met a unix that didn't come with Perl already installed, or as a > base option SunOS-4 :) Didn't have emacs, either, nor an ANSI-C compiler. And the OS came on QIC-150 tape (ie: 150 MB total capacity). Not that that defeats the argument ... Devin

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL

Hello, R exits every time I want to fit a GEE with AR-1 correlation structure. It just closes itself without any error. My Data has about 3300 observation (using the big data with about 20 000 observation it doesn't work at all). I'm a Windows 2000 User of R 2.0.1 but it was the same problem on a MacOSX. Is it my fault? Thank you for your help! Johanna Brandt

Hi all, My data are binary and within-subject correlated; there are three factors of interest and two of them are within-subject. So, I have considered modelling the data using a binomial regression with a GEE approach (which can be achieved using the function geeglm). One problem is that I got some empty cells, so I can't simply use logit as a link function. I was wondering if you know any

Here is an example given from ?yags library(methods) data(stackloss) Y1 <- yags(stack.loss~Air.Flow,id=1:21, data=stackloss) How can I access parts of the output. I tried: > str(Y1) Formal class 'yagsResult' [package "yags"] with 25 slots ..@ coefficients : num [1:2] -44.13 1.02 ..@ coefnames : chr(0) > Y1$coefnames Error in Y1$coefnames : $