Estimado Jorge
Perdóneme que lo moleste de nuevo, hay otra condición además de que
sum(y)=1 y es que
y[1] tiene que dar 0
en el ejemplo
y<- c (0.0000000000, 0.6321985783, 0.2325728597, 0.0855587737, 0.0314753138,
0.0115791209, 0.0042597205, 0.0015670636, 0.0005764905,
0.0002120790)
y[1]=0
sum(y)=1
esto se utiliza en el paquete EpiEstim donde se estiman las probabilidades
de esa manera
saludos
De: Dr. José A Betancourt Bethencourt [mailto:jbetancourt@iscmc.cmw.sld.cu]
Enviado el: miércoles, 18 de diciembre de 2013 09:15
Para: 'Jorge I Velez'
CC: 'Lista R'
Asunto: RE: [R-es] rgamma
Gracias y feliz fin de año a Ud. y a los demás miembros de la red
De: Jorge I Velez [mailto:jorgeivanvelez@gmail.com]
Enviado el: miércoles, 18 de diciembre de 2013 07:11
Para: Dr. José A Betancourt Bethencourt
CC: Lista R
Asunto: Re: [R-es] rgamma
Hola José,
Intenta lo siguiente:
x <- rgamma(25, 1, 2)
y <- x/sum(x)
y
sum(y)
par(mfrow = c(1, 2))
hist(x)
hist(y)
Saludos,
Jorge.-
2013/12/18 Dr. José A Betancourt Bethencourt <
<mailto:jbetancourt@iscmc.cmw.sld.cu> jbetancourt@iscmc.cmw.sld.cu>
Estimados
Quisiera tener la sintaxis generar una serie de números entre 0 y 1 con
distribución gamma de 25 puntos y que la suma de todos los puntos debe dar
1.
Por ejemplo:
c (0.000, 0.001, 0.012, 0.043, 0.078, 0.104, 0.117, 0.116, 0.108, 0.094,
0.078, 0.063, 0.049,
0.038, 0.028, 0.021, 0.015, 0.011, 0.008, 0.005, 0.004, 0.003,
0.002, 0.001, 0.001)
No hemos podido construirla a partir de la ayuda
rgamma(n, shape, rate = 1, scale = 1/rate)
Saludos
José
--
Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
Nacional de Salud. La persona que envia este correo asume el compromiso de
usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org <mailto:R-help-es@r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-help-es
De: Dr. José A Betancourt Bethencourt [mailto:jbetancourt@iscmc.cmw.sld.cu]
Enviado el: miércoles, 18 de diciembre de 2013 09:15
Para: 'Jorge I Velez'
CC: 'Lista R'
Asunto: RE: [R-es] rgamma
Gracias y feliz fin de año a Ud. y a los demás miembros de la red
De: Jorge I Velez [mailto:jorgeivanvelez@gmail.com]
Enviado el: miércoles, 18 de diciembre de 2013 07:11
Para: Dr. José A Betancourt Bethencourt
CC: Lista R
Asunto: Re: [R-es] rgamma
Hola José,
Intenta lo siguiente:
x <- rgamma(25, 1, 2)
y <- x/sum(x)
y
sum(y)
par(mfrow = c(1, 2))
hist(x)
hist(y)
Saludos,
Jorge.-
2013/12/18 Dr. José A Betancourt Bethencourt <jbetancourt@iscmc.cmw.sld.cu
<mailto:jbetancourt@iscmc.cmw.sld.cu> >
Estimados
Quisiera tener la sintaxis generar una serie de números entre 0 y 1 con
distribución gamma de 25 puntos y que la suma de todos los puntos debe dar
1.
Por ejemplo:
c (0.000, 0.001, 0.012, 0.043, 0.078, 0.104, 0.117, 0.116, 0.108, 0.094,
0.078, 0.063, 0.049,
0.038, 0.028, 0.021, 0.015, 0.011, 0.008, 0.005, 0.004, 0.003,
0.002, 0.001, 0.001)
No hemos podido construirla a partir de la ayuda
rgamma(n, shape, rate = 1, scale = 1/rate)
Saludos
José
--
Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
Nacional de Salud. La persona que envia este correo asume el compromiso de
usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org <mailto:R-help-es@r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece
Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de
Salud. La persona que envia este correo asume el compromiso de usar el servicio
a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/
[[alternative HTML version deleted]]