Displaying 1 result from an estimated 1 matches for "peventm2".
Did you mean:
peventm1
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...;- 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+cf[i,1,4]*sbp+cf[i,1,5]*diab+cf[i,1,6]*smok)
peventM1 <- exp(-10*lambda)^(exp(rsM1))
rsM2<- (cf[i,1,1]*age+cf[i,1,2]*tc+cf[i,1,3]*hdlc+cf[i,1,4]*sbp+cf[i,1,5]*diab+cf[i,1,6]*smok)
peventM2 <- exp(-10*lambda)^(estv*exp(rsM1))
# Proportion of more accurate predictions
pred <- sum(abs(pevent-peventM2) < abs(pevent-peventM1))/n
print(paste("Theta = ",l," proportion pred. M2 more accurate than M1 = ",pred));
}
corresponding output:
[1] "Theta = 1...