[R-es] para que incluya las tres
jbetancourt
jbetancourt en iscmc.cmw.sld.cu
Mar Jun 30 15:43:07 CEST 2015
Estimados
Gracias por la ayuda, le agregue para calcular IT que era de nuestro interes también y funciona bien
Mis saludos
library(R0)
vacunar <- function(pop.size, incid, durEst, ...){
resEst <- est.R0.AR(pop.size=pop.size, incid=incid,...)
ro <- resEst$R
vac <- 1 - (1/ro)
IT <- ro/durEst
return(data.frame(ro=ro,vac=vac,IT=IT))
}
vacunar(pop.size = 500, incid = c(1,5,12,13,15,19,22,34,41,53,70), durEst = 7)
Ya, pero puedes conseguir las salida de los dos resultados de la función de Jorge fácilmente:
#------------
vacunar <- function(pop.size, incid, durEst, ...){
resEst <- est.R0.AR(pop.size=pop.size, incid=incid,...)
ro <- resEst$R
vac <- 1 - (1/ro)
return(data.frame(ro=ro,vac=vac))
}
vacunar(pop.size = 500, incid = c(1,5,12,13,15,19,22,34,41,53,70), durEst = 7)
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es