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]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
