R help es - Jun 2015 - pr

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

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 <jbetancourt en
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 en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

?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 <cof en 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 <jbetancourt en
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 en 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 en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
	[[alternative HTML version deleted]]