Displaying 3 results from an estimated 3 matches for "ddimer".
Did you mean:
dimer
2007 Sep 19
3
Robust or Sandwich estimates in lmer2
..., I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of the standard error? If anybody could offer me a suggestion I would greatly appreciate it. Thank you.
Model-1:
> p.mle<-lmer2(ddimer~race+steroid+psi+sofa+apache + (apache|subject), method="ML", data=final, robust=TRUE, cluster="id", weights=final$w)
> beta=fixef(p.mle)
> Vcov=vcov(p.mle, useScale=FALSE)
> se=sqrt(diag(Vcov))
> beta
(Intercept) race steroid psi sofa...
2007 Feb 08
0
How to get p-values, seperate vectors of regression coefficients and their s.e. from the "yags" output?
...tion please? Also, how can I get the standard errors of these reg. coefficients as a seperate vector? Notice, as it is highlighted below, beta1=coef(wee) giving me "NULL" and also summary(wee) is not giving me nothing!
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...
2007 Oct 05
0
Extracting df (degree of freedom) & estfun (estimating function) from model built in lmer or lmer2
Hello R-users:
Could you please tell me how can I extract the "df (degree of freedom)" and "estfun (estimating functions)" for the following lmer (or lmer2) model?
wtd.mixed<-lmer(ddimer~race+steroid+psi+sofa+apache + (1|subject), method="ML", data=final, cluster="id", weights=w)
I tried the following codes:
- for the degree of freedom (erorr message in the red):
>so <- summary(wtd.mixed)
> so$df[1:6]
NULL
Warning message:
$ operator not defined...