[R-es] r

[Carlos Ortega]

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