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) #------------ > 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) ro vac 1 1.480649 0.3246206 Saludos, Carlos Ortega www.qualityexcellence.es El 29 de junio de 2015, 13:04, josebetancourt.cmw < josebetancourt....@infomed.sld.cu> escribió: > EStimados , > > la opción que propone carlos si da todas las soluciones, la segunda solo > da una salida > > saludos y gracias > > josé > > > > *De:* Carlos Ortega [mailto:c...@qualityexcellence.es] > *Enviado el:* sábado, 27 de junio de 2015 12:30 > *Para:* jbetancourt > *CC:* Lista R > *Asunto:* Re: [R-es] pr > > > > Hola, > > Puedes hacerlo de esta forma: > > #----------------- > library(R0) > > durEst <-7 #7 dias estimacion duracion > resEst <- est.R0.AR <http://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 <http://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 > > > > > > > > > > la otra propuesta solo que da salida a un elemento > > 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 <http://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.- > > > > > > > > 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 <http://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 > > -- 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