search for: gee

I was trying to learn cross building an R pcakage for windows on my linux machine. I picked a relatively small library gee to test. I downloaded the cross-tools and put them in my path, jyan at ludwig:/a3/jyan/src/R-1.2.1/src/gnuwin32$ echo $PATH /home/jyan/cross-tools/bin:/home/jyan/cross-tools/i386-mingw32msvc/bin: /usr/lib:/usr/local/bin:/usr/bin:/bin:/usr/bin/X11:/usr/games Following the instructions, I was able...

...or the purposes of my analysis. 2) The subjects were not sampled randomly. Some of the subjects had a higher probability of selection, and I want to incorporate inverse-probability weights into my analysis to account for this. (The inverse-probability weights are already calculated). I know that GEE is an appropriate technique to deal with Issue #1, and I've toyed with the gee pack in R. R> library("gee") http://cran.r-project.org/web/packages/gee/gee.pdf My question is -- how can I incorporate the sampling weights into the GEE code? I don't see a spot for it based on...

Dear all, I need to fit a gee model with an auto-regressive correlation structure and I faced some problems. I attach a simple example: ####################################################### library(gee) library(geepack) # I SIMULATE DATA FROM POISSON DISTRIBUTION, 10 OBS FOR EACH OF 50 GROUPS set.seed(1) y <- rpois(500,...

Dear List: I found an error when I called the 'gee' function. I cannot solve and explain it. There are no errors when I used the 'geeglm' function. Both functions fit the gee model. The project supervisor recommends me to use the 'gee' function. But I cannot explain to him why this error happens. Would you help me solve t...

Hi, I'm attempting a comparative analysis using the function compare.gee, package (ape), which uses phylogeny as a correlation matrix in gee (package gee), in R version 2.2.0 on a Windows XP PC. I'm looking to model the relationship between a continuous explanatory variable and a binary response variable for 600 species, taking into account the phylogeny of those s...

Hi, I am fitting a GEE model using gee R package, but I am not sure how to get OR and its CI? Could anyone give me some hints? Here are some output: > gee.obj <- gee(Affection~Sibsex+Probandsex,id = FAMID,family = binomial,corstr = "independence",data =seldata) Beginning Cgee S-function, @(#) geeform...

Hi, I am trying to install gee on our server but I get the error below. I do not have root on this machine so no control on how R was installed itself. It looks like it cannot find blas libs, the only ones i can find on the machine are: /usr/lib64/libblas.so.3 -> libblas.so.3.0.3 /usr/lib64/libblas.so.3.0 -> libblas.so.3....

I'm running some GEE models but when I request the summary(pcb.gee) all I get are rows and rows of intercorelations and they fill up the screen buffer so I can not even scroll back to see what else might be in the summary. How do I get the summary function to NOT print the intercorrelations? Thanks, -- Dean Sonne...

Hello, I am interested in running generalized estimating equation models in R. Currently there are two main packages for doing so in R, geepack and gee. I understand that even though one can obtain similar to almost identical results using either of the two, that there are differences between the packages. The paper that introduces the geepack package ( https://www.jstatsoft.org/article/view/v015i02) states that the three features tha...

I'm not sure if the gee package is still actively maintained, but I for one find it extremely useful. However, I've come across a few infelicities that I'm hoping could be resolved for future versions. Hope it's okay to list them all in one post! They are: (1) AR(1) models don't fit when clustsize = 1 for...

windows Vista R 2.10.1 Is it possible to get p values from gee? Summary(geemodel) does not appear to produce p values.: > fit4<- gee(y~time, id=Subject, data=data.frame(data)) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) time 1.1215614 0.8504413 > summary(fit4) G...

Does anyone know how could I get the significant codes from mixed model output fitted with a GEE package? The output I got is the following: GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA gee S-function, version 4.13 modified 98/01/27 (1998) Model: Link: Logit Variance to Mean Relation: Binomial Correlation Structure: Exchangeable Call: gee(formula = b...

...and take into account the fact that I have repeated measurements for every subject. If I do this with repeated measures ANOVA (in which the patient is considered a random effect) I got the following results: library(nlme) data<-read.table("http://146.9.88.18/uploads/dataGEE.txt",header=TRUE) res<-lme(Y~T,random=~1|P,data=data) summary(res) So the p-value for significance of the treatment effect is 0.069. I would like to use also as a variant analysis a Generalized Estimation Equation Model, like library(gee)...

I'm using the function gee from the library(gee) gee(Y~X,id=clust.id,corstr="exchangeable",b=tmc$coef,family=binomial(link=logit),silent=T) Every time it runs, it dutifully prints out Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 user's initial regression estimate [,1] [1,] -...

Dear all, I don't understand the following behaviour: Running compar.gee (in library ape ) with and without the option 'data', it give me different results Example: .... Start R .... > load("eiber.RData") > ls() [1] "gee.na" "mydata" "mytree" > library(ape) > # runnig with the option data= mydata > com...

Hi all, I have fitted a gee model with the gee package and included restricted cubic spline functions. Here is the model: chol.g <- gee(SKIN ~ rcs(CHOLT, 3), id=ID, data=chol, family=binomial(link="logit"), corstr="exchangeable") To extract the log odds I use: predict.glm(chol.g, type = "link&...

Hello, I would like to try to fit a GEE with user-specified link function. I read through a couple of thread on the list, for example http://tolstoy.newcastle.edu.au/R/help/04/12/9768.html#start and http://tolstoy.newcastle.edu.au/R/help/06/04/25298.html. I noticed that they are all 6 or more years old and the answer is very clear for...

...ffect of this variable is normatively linear). I fit three models, but I am observing differences between the three. The first model is a quasibinomial without any subject effects using glm. The second is a random-effects model using lmer. The third model is a generalized estimating equation using gee from the gee package in which I cluster for the subject using an unstructured correlation matrix. The results of the first and the third model almost coincide, but the second, using lmer, shows an insginficant coefficient where I would expect a significant one. The other 2 models show the coefficie...

Hi, The standard error from logistic regression is slightly different from the naive SE from GEE under independence working correlation structure. Shouldn't they be identical? Anyone has insight about this? Thanks, Qiong a<-rbinom(1000,1) b<-rbinom(1000,2,0.1) c<-rbinom(1000,10,0.5) summary(gee(a~b, id=c,family="binomial",corstr="independence"))$coef summary...

I am trying to learn the gee function in R. So I try to generate some data and use this function. I have the following lines: ######################################## Gee # Generating lny=10+2*Si-Si^2+eta # eta ~ N(0,1) # Si ~ U(0,11) eta <- vector(mode="numeric",100) eta <- rnorm(100) Si <- vector(mod...