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 R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es