Hola a todos, Tomando el ejemplo de Carlos, podrías construir una funcion que realice lo que necesitas:
vacunar <- function(pop.size, incid, durEst, ...){ resEst <- est.R0.AR(pop.size=pop.size, incid=incid,...) ro <- resEst$R vac <- 1 - (1/ro) vac } vacunar(pop.size = 500, incid = c(1,5,12,13,15,19,22,34,41,53,70), durEst = 7) ## [1] 0.3246206 Saludos cordiales, Jorge.- 2015-06-28 2:30 GMT+10:00 Carlos Ortega <c...@qualityexcellence.es>: > Hola, > > Puedes hacerlo de esta forma: > > #----------------- > library(R0) > > durEst <-7 #7 dias estimacion duracion > resEst <- est.R0.AR(pop.size=500, incid=c(1,5,12,13,15,19,22,34,41,53,70)) > #el resultado es 1.48 > > ro <- resEst$R > > vac <-1-(1/ro) # (1-1/Ro) > vac > > IT <- ro/durEst # Ro/7 > IT > #----------------- > > "est.R0.AR" devuelve un objeto en el que uno de los elementos es "R" que > almacena la estimación. > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > > El 27 de junio de 2015, 16:01, jbetancourt <jbetanco...@iscmc.cmw.sld.cu> > escribió: > > > Estimados > > No tengo experiencia en programación y quisiera preparar el siguiente > > script para los estudiantes > > Una supuesta enfermedad, al parecer de transmisión aérea y que suele > durar > > 7 días, entra en una localidad de 500 habitantes, la incidencia diaria > es > > c(1,5,12,13,15,19,22,34,41,53,70) > > Al calcular Ro de la manera siguiente: > > > > library(R0) > > est.R0.AR(pop.size=500, incid=c(1,5,12,13,15,19,22,34,41,53,70)) #el > > resultado es 1.48 > > > > yo quisiera programar el rango de personas a vacunar (1-1/Ro) sin tener > > que escribir 1.48 sino que el script lo tome del resultado de la > > estimación realizada en est.R0.AR, lo mismo para calcular las personas > > a vacunar Ro/7 > > vac <-1-(1/1.48) # (1-1/Ro) > > vac > > IT<-1.48/7 # Ro/7 > > IT > > Mis saludos cordiales > > José > > > > [[alternative HTML version deleted]] > > > > > > _______________________________________________ > > R-help-es mailing list > > R-help-es@r-project.org > > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es