Displaying 1 result from an estimated 1 matches for "fitm2".
Did you mean:
fit2
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...runif(n,0,1); # For random variate generation
t <- (-1*log(r)) / (lambda*v*exp(age*beta[1]+tc*beta[2]+hdlc*beta[3]+sbp*beta[4]+diab*beta[5]+smok*beta[6]))
fitM1=coxph(Surv(t,c)~age+tc+hdlc+sbp+diab+smok); # Model 1, no frailty
cf[i,1,] <- round(fitM1$coef,6) # Model 1 coefficients
fitM2=coxph(Surv(t,c)~age+tc+hdlc+sbp+diab+smok+
frailty(id,dist="gamma",sparse=T,method="em")) # Model 2, frailty
cf[i,2,] <- round(fitM2$coef,6) # Model 2, coefficients
# store estimated Variance of random effect
varre=fitM2$history$frailty$history[length(fitM2$history...