search for: gee

Displaying 20 results from an estimated 509 matches for "gee".

Did you mean: ge
2001 Feb 23
0
cross building
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...
2012 Jun 12
2
GEE with Inverse Probability Weights
...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...
2009 Feb 09
1
gee with auto-regressive correlation structure (AR-M)
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,...
2008 Sep 07
1
an error to call 'gee' function in R
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...
2005 Nov 18
0
Likely cause of error (code=1) in compar.gee/gee
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...
2010 Sep 08
0
How to get OR and CI from GEE R package
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...
2011 Jan 26
1
Compilation errors when installing gee
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....
2005 Sep 28
1
gee models summary
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...
2008 Dec 08
0
gee niggles
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...
2010 Sep 10
2
gee p values
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...
2011 Aug 15
1
Get significant codes from a model output fit with GEE package
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...
2006 Apr 10
5
p values for a GEE model
...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)...
2009 Oct 13
2
gee: suppress printout
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,] -...
2004 Feb 08
1
APE: compar.gee( )
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...
2008 Dec 01
1
gee + rcs
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&...
2010 Oct 12
1
GEE with user-specified link function
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...
2008 Jul 07
1
GLM, LMER, GEE interpretation
...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...
2008 Sep 09
2
naive variance in GEE
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...
2004 Aug 18
1
Gee
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...
2003 Oct 24
1
gee and geepack: different results?
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 some...