Displaying 2 results from an estimated 2 matches for "rsm1".
Did you mean:
rsm
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...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+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...
2012 Apr 02
2
graficar superficie de respuesta con los valores de las variables originales ...
estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y
hago un grafico 2D de la superficie de respuesta con contour(), pero como
el argumento de esta función es la superficie de respuesta ajustada y ésta
se ha ajustado codificando las variables, de acuerdo a la recomendación de
la documentación del paquete rsm, obtengo un gráfico con los ejes en
unidades codificadas