search for: peventm1

...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+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(...