[R-es] para que incluya las tres

[jbetancourt]

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]]