Re: [R] Use of survreg.distributions

2008-04-25 Thread Abdus Sattar
xactsurv, 1])) : NaNs produced 2: In log(y) : NaNs produced Note, the data file y.data does not contain any missing data! Do you know why it is giving me such an error message please? Thank you again for your helpful comment/suggestion.   Best Regards,   Abdus Sattar [EMAIL PROTECTED] ---

Re: [R] Use of survreg.distributions

2008-04-24 Thread Abdus Sattar
n these codes then I got the following error message, Error in survreg(Surv(y, y >= -5, type = "left") ~ x + : Invalid distribution object Does anybody can help me in identifying the error(s) in these code please?

Re: [R] gradient and Hessian from lmer2

2007-10-29 Thread Abdus Sattar
Dear R-user: Could you please tell me how to extract gradient and Hessian from the following example model of lmer2 package? (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) I would really appreciate your kind help. Best, Sattar

[R] estfun & df

2007-10-08 Thread Abdus Sattar
Hello EVERYONE, I need an URGENT help from you please! How can I see the "estfun" (empirical estimating function) and "df" (degree of freedom) from the following mixed-model please? (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Many thanks in advance for your kind help.

[R] Extracting df (degree of freedom) & estfun (estimating function) from model built in lmer or lmer2

2007-10-05 Thread Abdus Sattar
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 th

Re: [R] Robust or Sandwich estimates in lmer2

2007-09-19 Thread Abdus Sattar
Message From: Thomas Lumley <[EMAIL PROTECTED]> To: "Doran, Harold" <[EMAIL PROTECTED]> Cc: Sundar Dorai-Raj <[EMAIL PROTECTED]>; Abdus Sattar <[EMAIL PROTECTED]>; [EMAIL PROTECTED]; [EMAIL PROTECTED] Sent: Wednesday, September 19, 2007 1:13:20 PM Subject: Re: [

[R] Robust or Sandwich estimates in lmer2

2007-09-19 Thread Abdus Sattar
Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question i