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.cmw en 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:cof en 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 en 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]]