similar to: problems with function geese in geepack

Good evening, In my research I am studying the marginal models, where the main goal is on the structure of the association. My practical example has cluster with up to 600 observations and with this database, the function geese() return me the following message: This application has requested the runtime to terminate it is on unusual way. Please contact the application’s support team

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091126/7134fc17/attachment-0001.pl>

Hi there, I just installed the package under windows, and it works very well. However, as I tried to install it on a cluster which uses linux OS, "Linux lionxc.aset.psu.edu 2.6.18-194.el5 #1 SMP Tue Mar 16 21:52:39 EDT 2010 x86_64 x86_64 x86_64 GNU/Linux" I always failed. The message during the installation process is as follows: > install.packages("geepack") Warning

Hi all I am analyzing a data set containing information about the behaviour of marine molluscs on a vertical wall. Since I have replicate observations on the same individuals I was thinking to use the geepack library. The data are organised in a dataframe with the following variables Date = date of sampling, Size = dimensions (mm) Activity duration of activity (min) Water = duration of

Hi, I downloaded both gee and geepack, and I am trying to understand the differences between the two libraries. I used the same data and estimated the same model, with a correlation structure autoregressive of order 1. Surprisingly for me, I found very different results. Coefficients are slightly different in value but sometimes opposite in sign. Moreover, the estimate of rho (correlation

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 <-

Dear all, We have got data (response and predictor variables) for each country of the world; I started by fitting standard GLM and tested for spatial correlation using variogram models (geoR) fitted to the residuals of the GLM. Spatial autocorrelation is significant. Therefore, I think about using general estimation equations (geeglm or geese in geepack) allowing for residual spatial

We are currently analyzing data on children clustered in day care-centers (DCC). We have tried to use geepack and gee libraries to estimate an overall incidence rate for absences (=number of absences/risk time) by specifying geese(number.absences ~ offset(log(risktime)), id=day.care.id, family=poisson("log"), data=dcc, corstr="exch",

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 A quick question. Standard errors reported by gee/yags differs from the ones in geeglm (geepack). require(gee) require(geepack) require(yags) mm <- gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable") mm2 <- geeglm(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable", std.err = "san.se") mm3 <- yags(breaks ~

I am working with geese from geepack, does anyone know how to get the predicted values from this program. I tried the predict function and got the error Error in predict(results.e, newdat = prog) : no applicable method for "predict" Thank you for any help you can send, Suzette ================================= Suzette Blanchard, Ph.D. UCSD-PPRU

Hi, I am trying to fit a ordinal GEE model using ordgee {geepack}. In order to check the validity of the function, I specified the correlation structure as independence (i.e. constr = "independence") and compared the result with that using polr {MASS}. Because a GEE model with an independent working correlation structure is equivalent to an ordinary GLM model, we would expect the same

Dear R users, I fitted a GEE model using the function 'geese' (or 'geeglm') with user defined correlation matrix. I want to get the var-cov matrix of the regression coefficients. But the output provides only limited information. I would be very much thankful if you could kindly let me know how to get it..since I am struggling lot getting this. Thanks -- View this message in

  I am trying to used "ordgee" from "geepack" for an ordinal dataset. When I write the code it returns  "Warning message:In binomial(link) : use of binomial(link=link) is deprecated" , but the program runs. Even when I run your example for "ohio" and "respdis", it returns the same error. Please guide me [[alternative HTML version deleted]]

what is exactly difference between "gee", "geese", and "ordgee" [[alternative HTML version deleted]]

Hello, I want to fit a GEE with a user-defined link function. For the user-defined link-function I still read http://finzi.psych.upenn.edu/R/Rhelp01/archive/6555.html and http://finzi.psych.upenn.edu/R/Rhelp02a/archive/25727.html. Only for testing purposes I added a new link function (corlogit) in make.link (as well as in binomial) with exactly the same code as logit before using my code.

event.nab.2 is 0/1 and I dichotomized va to get va.2 to see if I could get geeglm to work. glm has no problem with the data but geeglm chokes. Each subject (patient.id) has at most 2 observations and more than 3/4 of the subjects have 2 observations. I have even worse problems trying to use glmmPQL from MASS and worse still trying to use lmer from lme4. But I figured a marginal model would work.

This is a summary and extension of the thread "GLMM (lme4) vs. glmmPQL output" http://maths.newcastle.edu.au/~rking/R/help/04/01/0180.html In the new revision (#Version: 0.4-7) of lme4 the standard errors are close to those of the 4 other methods. Thanks to Douglas Bates, Saikat DebRoy for the revision, and to G?ran Brostr?m who run a simulation. In response to my first posting, Prof.

Hello, I'm trying to fit a special kind of proportional odds model from: Whitehead et al. (2001). Meta-analysis of ordinal outcome using individual patient data. Statistics in medicine 20: 2243-2260. (model 2) The data are as follows: library(nlme) library(geepack) library(Design) library(MASS) options(contrasts=c("contr.SAS","contr.poly")) counts <-

Hello, I have the following r-codes for solving a quasilikelihood estimating equation: >library(geepack) >fit<-geese(y~x1+x2+x3,jack=TRUE,id=id,scale.fix=TRUE,data=dat,mean.link = "logit", corstr="independence") Now my question is how can I calculate the confidence interval of the parameters of the above model "fit"? [[alternative HTML version deleted]]