Displaying 4 results from an estimated 4 matches for "yagsresult".
2008 Oct 29
1
how can I access parts of yags output
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 : chr(0)
> Y1$coefnames
Error in Y1$coefnames : $ operator not defined for this S4 class
Thanks,
Juliet
2004 Mar 19
0
yags, GEEs, and GLMMs
...ht 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, please let me know and I'll send you the most urgent questions, but I don't want to bother everyone with it.
Many tha...
2004 Mar 19
0
yags, GEEs and GLMMs
...ht 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 function, please let me know and I'll send you the most urgent questions, but I don't want to bother everyone with...
2007 Feb 08
0
How to get p-values, seperate vectors of regression coefficients and their s.e. from the "yags" output?
...?
Working correlation model: exchangeable
alpha est: 0.7344
NULL
Pan QIC(R): 7534.732
QLS: 56989.3
Rotnitzky-Jewell: 9.477, 143.987
yags/R: $Id: yags.R,v 1.5 2004/10/22 18:49:23 stvjc Exp $
> beta1=coef(wee)
> beta1
NULL
> summary(wee)
Length Class Mode
1 yagsResult S4
> summary(wee)
FYI, in the following geeglm analysis, I have gotten beta=coef(wgee) as a row vector(highlighted):
mf=formula(Ddimer~newrace+steroid+treatment+SOFA+PSI)
> wgee=geeglm(mf, id=Subject, data=na.omit(final), weights=w, family=gaussian("identity"), corstr=...