Displaying 1 result from an estimated 1 matches for "estv".
Did you mean:
est
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...heta = ",l, " variance in estimated ind. frailty = ", var(exp(fitM2$frail))))
print(paste("Theta=",l," org beta ",beta," M1: ",cf[i,1,], " M2: ", cf[i,2,]))
# Calculate the 'actual' 10-year risk and the risk estimated using M1 and M2
estv <- exp(rep(fitM2$frail,2)) # Individual frailty values from M2
# Original risk score
rs<- age*beta[1]+tc*beta[2]+hdlc*beta[3]+sbp*beta[4]+diab*beta[5]+smok*beta[6]
pevent <- exp(-10*lambda)^(exp(rs)) # Original 10-year risk
rsM1<- (cf[i,1,1]*age+cf[i,1,2]*tc+cf[i,1,3]*hdlc+c...