[R-es] r

josebetancourt.cmw josebetancourt.cmw en infomed.sld.cu
Lun Jun 29 13:04:18 CEST 2015


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



Más información sobre la lista de distribución R-help-es