Displaying 6 results from an estimated 6 matches for "alphainit".
Did you mean:
alphainfo
2007 Feb 08
0
How to get p-values, seperate vectors of regression coefficients and their s.e. from the "yags" output?
...thing!
The following is the output from the yags analysis:
> yf=formula(Ddimer~newrace+steroid+treatment+SOFA+PSI)
> wee=yags(yf, id=Subject, data=final, cor.met=as.double(rep(0:6, 872)), family=gaussian, corstruct="exchangeable", control=yags.control(), weights=w, betainit=NULL, alphainit=.1, subset=NULL)
> wee
YAGS (yet another GEE solver) $Date: 2004/10/22 18:49:23 $
Call:
yags(formula = yf, id = Subject, cor.met = as.double(rep(0:6,
872)), family = gaussian, corstruct = "exchangeable", control = yags.control(),
weights = w, betainit = NULL, alphainit = 0.1,...
2008 Oct 29
2
call works with gee and yags, but not geepack
...quot;\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 following error:
Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor, corp, :
nrow(zsca) and length(y) no...
2004 Mar 19
0
yags, GEEs, and GLMMs
...tion (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 the parameter estimates, std errors, t- and p-values etc.?
And if someone is willing to answer questions about the glmmPQL function, plea...
2004 Mar 19
0
yags, GEEs and GLMMs
...tion (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 the parameter estimates, std errors, t- and p-values etc of the model?
And if someone is willing to answer questions about the glmmPQL f...
2009 Dec 08
0
Difference in S.E. gee/yags and geeglm(/geese)
...sion, id=wool, data=warpbreaks,
corstr="exchangeable")
mm2 <- geeglm(breaks ~ tension, id=wool, data=warpbreaks,
corstr="exchangeable", std.err = "san.se")
mm3 <- yags(breaks ~ tension, id=wool, data=warpbreaks,
corstr="exchangeable", alphainit=0.)
# S.E.
# gee:
sqrt(diag(mm$robust.variance))
#(Intercept) tensionM tensionH
# 5.77 7.46 3.73
# geeglm:
sqrt(diag(mm2$geese$vbeta))
# [1] 8.17 10.56 5.28
# yags:
sqrt(diag(slot(mm3, "robust.parmvar")))
# [1] 5.77 7.46 3.73
Any explanation of this behavior i...
2009 Jan 22
0
detecting the source of memory consumption (example provided)
...s <- 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(abs(z.gee), lower.tail = FALSE);
pvalues[Var_num] <- signif(pval,3);
}