search for: yag

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 'http://www.stats.ox.ac.uk/...

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

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

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

...k Crawley, when are you going to write a chapter on mixed models with binary response variables?......). I have used the glmmPQL function (about which I also have a whole list of questions....) as mentioned in MASS, but now I was told that GEE models might be better suited for my data. So I got the yags package and tried a model like this: > gee1<-yags(alive~sex*weight, family=binomial,id=id,corstr='exchangeable', alphainit=0.) But then the summary only gives me this: > summary(gee1) Length Class Mode 0 yagsResult list How do I get R to give me...

...k Crawley, when are you going to write a chapter on mixed models with binary response variables?......). I have used the glmmPQL function (about which I also have a whole list of questions....) as mentioned in MASS, but now I was told that GEE models might be better suited for my data. So I got the yags package and tried a model like this: > gee1<-yags(alive~sex*weight, family=binomial,id=id,corstr='exchangeable', alphainit=0.) But then the summary only gives me this: > summary(gee1) Length Class Mode 0 yagsResult list How do I get R to give m...

Dear R users, I search the R archives and noted that the same problem has been posted but without solution. I know there IS instructions by the author of yags, but I just couldn't figured out. I know gee and geepack can also perform generalized estimating equation, but the reason why I need yags is I want to perform a Small Sample Adjustments by Michael P. Fay and Barry I. Graubard in Biometrics, 2001, 57: 1198-1206. I'm working on some de...

...hlighting 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 <- yags(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat,alphainit=0.05) However, I am making a mistake with: library(geepack) fit.geese <- geese(score ~ chem + time, id=id, family=gaussian,corstr="exch",data=dat) I obtain the follow...

Full_Name: Markus Jantti Version: 1.5.1 OS: GNU/Linux Unstable (Sid) Submission from: (NULL) (193.166.59.44) I was trying to get the package yags 3.20 by Vince Carey <stvjc@channing.harvard.edu> to run, but - after successfull installation -- loading the library failed with the message > library(yags) Error in stop(paste("Duplicated slot names: ", paste(slots[duplicated(slo : unused argument(s) (sep ...) Error...

...0,80),50,replace=TRUE) id <- rep(1:25,each=2) var1 <- rnorm(50); var2 <- rnorm(50); var3 <- rnorm(50); myData <- data.frame(response,x1,age,id,var1,var2,var3) numVars <- ncol(myData)-4; pvalues <- rep(-1,numVars); names(pvalues) <- colnames(myData)[5:ncol(myData)]; library(yags) for (Var_num in 1:numVars) { fit.yags <- yags(myData$response ~ myData$age+myData$x1*myData[,(Var_num+4)], id=myData$id, family=gaussian,corstr="exchangeable",alphainit=0.05) z.gee <- fit.yags at coefficients[5]/sqrt(fit.yags at robust.parmvar[5,5]); pval <- 2 * pnorm(...

markus.jantti@iki.fi wrote: > > Full_Name: Markus Jantti > Version: 1.5.1 > OS: GNU/Linux Unstable (Sid) > Submission from: (NULL) (193.166.59.44) > > I was trying to get the package yags 3.20 by Vince Carey > <stvjc@channing.harvard.edu> > to run, but - after successfull installation -- loading the library failed with > the > message > > > library(yags) > Error in stop(paste("Duplicated slot names: ", paste(slots[duplicated(slo > : >...

...39;t require actually installing recent versions of all those old packages? I suppose I could just use R 2.10.1, but suppose at some point I want to use both an old package and a new package that was built under R 2.12.1 in the same program? That has happened by the way. I wanted to use deSolve and yags. Since I don't have an old version of yags, I had to install the current version on CRAN, and it won't work under 2.10.1. My second question is, if not, should the R developers reconsider their strategic decision to invalidate packages just because they were built under early versions of...

Hi, I'm using 'geepack' to run Generalized Estimating Equations. I'm aware that I can use anova to compare two models, but would it be possible to test QIC on R? It seems that there were similar questions a couple of years ago, but the question has not been answered yet. I'd appreciate if someone could show me the code! Thank you, Sachi [[alternative HTML version

...ested random effect, or so I believe I read while scanning backlogs of the R-Help list. glmmPQL() is in the ballpark of what I need, but once again seems to lack the "outer" variable specification that nlme() has, and which PROC GENMOD also has---and which I need. I read someplace of yags() that apparently uses GEE to estimate parameters of nonlinear models including GLIMs/mixed models, just the way PROC GENMOD (and many another program) does. But on trying to install it (either v4.0-1.zip or v4.0-2.tar.gz from Carey's site, or Ripley's Windows port) from a local, down...

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

Hello Everybody: I am searching for an R package for fitting Generalized Estimating Equations (GEE) with weights (i.e. Weighted Estimating Equations). From the R documentation I found "geese(geepack)" for fitting Generalized Estimating Equations. In this documentation, under the paragraph “weights” it has been written, “an optional vector of weights to be used in the fitting process.

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

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

Default arguments are evaluated in the function frame, not in the calling environment (nor in the same place as explicit arguments). > Which to me reads that a with statement as above is equivalent to > > > attach(data) ; aov.SS1(y=Obs) ; detach(data) > > Or is that just wishful thinking?? The latter. On Thu, 7 Oct 2004, RenE J.V. Bertin wrote: > Hello, > >

I have would like to use logistic regression to analyze the percentage of correct responses in a 2 alternative forced choice task. The question is whether one needs to take into account the fact expected probabilities for the percentage of correct responses ranges between 0.5 and 1 in this case and how to adjust the link function accordingly in R (see details below). Gabriel Subjects were asked