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

Responder a