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