similar to: gee

Hi there, Please help me this out. > m [,1] [,2] [1,] 1.1 1.0 [2,] 1.0 1.1 > chol2inv(m) [,1] [,2] [1,] 1.5094597 -0.7513148 [2,] -0.7513148 0.8264463 > CT -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Any hints on how to do a significance test using gee? Is the "robust z" the thing to use? There is apparently another implementation of gee called yags (http://www.biostat.harvard.edu/~carey), but I can't figure out how to install it. I think it does significance tests. -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Home page:

Dear all, It looks like the parameters var-cov matrix returned by gee() is not adjusted for the scale parameter: > fm1 <- gee(nbtrp ~ strate * saison + offset(log(surf)), family = poisson, data = Eff2001, + id = loc, tol = 1e-10, corstr = "exchangeable") [1] "Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27" [1] "running glm to get initial

A bug (misfeature) has been found in the S gee library (and thus in the R gee library). The problem, which is shared by nearly all gee implementations, involves the calculation of working correlations when some clusters have only one observation. For compatibility reasons nearly everyone uses the computing formula from the first SAS macro by Karim, rather than the formula from the original GEE

Hi all, 1. Can nonlinear models be used in gee? For example, I have a dataset which contains 2 variables x and y, I wrote data(ex) atttach(ex) a<-100 b<- -0.5 c<-4.5 d<-20 Then: a. y~gee(y~d+(a-d)/(1+(x/c)^b)) Error in terms.formula(formula, data = data) : invalid power in formula b. y~gee(y~d+(a-d)/(1+(x/c)) Error in model.frame(formula, rownames, variables, varnames, extras,

This is the first announcement of the Comprehensive R Archive Network (CRAN) CRAN is a collection of sites which carry identical material, consisting of the R&R R distribution(s), the contributed extensions, documentation for R, and binaries. The CRAN master site can be found at the URL ftp://ftp.ci.tuwien.ac.at/pub/R (Austria) and is currently being mirrored daily at

This is the first announcement of the Comprehensive R Archive Network (CRAN) CRAN is a collection of sites which carry identical material, consisting of the R&R R distribution(s), the contributed extensions, documentation for R, and binaries. The CRAN master site can be found at the URL ftp://ftp.ci.tuwien.ac.at/pub/R (Austria) and is currently being mirrored daily at

Attached is an updated version of the CRAN announcement, also to be sent out as soon as 0.50 is out. Again, please have a look if possible. -k ************************************************************************ This is the first announcement of the Comprehensive R Archive Network (CRAN) CRAN is a collection of sites which carry identical material, consisting of the

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

Dear Lister I need to perform Generalized estimating equations on some data. Has 'R' any function or routine to do it? (GEE is same as GLM but with correlated y's) Thanks in advance Mostafa _________________________________________________________________ Protect your PC - get McAfee.com VirusScan Online http://clinic.mcafee.com/clinic/ibuy/campaign.asp?cid=3963

`update.packages' fails when a package with incorrect Version field is found in R_LIBS: > update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain' length 49161 bytes opened URL .......... .......... .......... .......... ........ downloaded 48Kb Error in if (a[k] > b[k]) return(TRUE) else if (a[k] < b[k]) return(FALSE) :

I am using the 0.50 alpha version of R packaged (in 3 parts) by Kurt for Debian Linux. Package: r-base Status: install ok installed Priority: optional Section: local Maintainer: Kurt Hornik <Kurt.Hornik at ci.tuwien.ac.at> Version: 0.50b7-1 Depends: libc5, xlib6, libreadline2 Description: R, a language not entirely unlike the language S. Package: r-contrib Status: install ok

Hi, I'm trying to fit a marginal model via GEE but I'm getting strange results and few problems. If I set the working correlation as exchangeable I'm getting the same fitting when I set as independent. Comparing to SAS results it shouldn't happen. If I try to use another working correlation (like AR-M or stat_M_dep), R just exits without giving any error message. Another doubt

In lm and glm, the weights command should only be used when the variances of the observations are inversely proportional to the weights. Recently, however, a question came up regarding how one could estimate lm and glm models with weights where the weights refer to number of observations (counts). Because lm and glm do not handle this case, SE values will be wrong if one uses the weight command

Dear all: Y=a+bX1+cX2 In the Gee (Generalized Estimation Equation) arguments: The arument Corstr has sveral choices: "independence" "fixed" "stat_M_dep" "non_stat_M_dep" "exchangeable" "AR-M" "unstructured" What does each term mean? How do I choose among them? How do I know the correlation structure of

Full_Name: Emmanuel Charpentier Version: R 1.1.1 OS: Windows 98 FE / Linux (Debian 2.2) Submission from: (NULL) (164.2.255.244) This bug report is pro forma. I've already reported it informally on R-help and R-devel. The current (0.7) version of RODBC has a problem with missing values. It correctly gets missing "character" values (as blanks), but puts the last non-missing value of

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, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate

Dear R team, I had a strange problem that appeared suddenly, and which I'd like to to mention in case it is of interest (although I suspect that I may have unwittingly done something dumb: if so please accept my apology) > library(ctest) ## or any package other than "oz" that I tried Error in as.character(table) : Object ".Dyn.libs" not found > traceback() [1]

Binary distributions of R-0.49 for NEXTSTEP (Intel and M68k) are now available on CRAN: http://www.ci.tuwien.ac.at/R/bin/i386-nextstep/R.0.49.I.b.tar.gz http://www.ci.tuwien.ac.at/R/bin/m68k-nextstep/R.0.49.N.b.tar.gz Because NEXTSTEP doesn't support dynamic loading, I've compiled in a number of additional functions (mostly from the user contributed section of CRAN) which I find

Versions of KernSmooth VR5.3pl037 boot gee rpart sm (delete src/Makefile: I'll get the version on CRAN changed) tree for R-0.64 are now on CRAN. Please note that spatial (in the VR bundle) and sm really do need 0.64 as we use the platform-independent file handling, and sm uses chull. This will be the last version of the VR5.x software for the second edition, as a new more `R