.... Un código más optimizado para la aproximación Monte Carlo de la
distribución de IH
N=10000 # tamaño simulación Monte Carlo
n=5 # numero de uniformes
cdf.IH <-function(x,n,N) {
z=replicate(N,sum(runif(n)))
apply(outer(x,z,">="),1,mean)
}
x=seq(0,5,.1)
y=cdf.IH(x,n=n,N=N)
plot(x,y,type="l")
Un saludo. Olivier
----- Mensaje original -----
De: "Genaro Llusco" <gellusco en gmail.com>
Para: r-help-es en r-project.org
Enviados: Martes, 21 de Abril 2015 23:07:14
Asunto: [R-es] distribucion de IRWIN HALL
estimados
estoy considerando programar la funcion de distribucion de Irwin hall.
lamentablemente no he tenido exito, pido que alguien me pueda colaborar con
aquello, les quedo agradecido de antemano.
--
atte. Lic. Genaro Llusco Silvestre
gellusco en gmail.com
Telf: 74028671
blog personal es:
http://www.cientificest.blogspot.com
[[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
_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
Alternativamente, uno puede implementar la fórmula que aparece en El día 22 de abril de 2015, 12:13, Olivier Nuñez <onunez en unex.es> escribió:> .... Un código más optimizado para la aproximación Monte Carlo de la distribución de IH > > N=10000 # tamaño simulación Monte Carlo > n=5 # numero de uniformes > > cdf.IH <-function(x,n,N) { > z=replicate(N,sum(runif(n))) > apply(outer(x,z,">="),1,mean) > } > x=seq(0,5,.1) > y=cdf.IH(x,n=n,N=N) > plot(x,y,type="l") > > Un saludo. Olivier > > > ----- Mensaje original ----- > De: "Genaro Llusco" <gellusco en gmail.com> > Para: r-help-es en r-project.org > Enviados: Martes, 21 de Abril 2015 23:07:14 > Asunto: [R-es] distribucion de IRWIN HALL > > estimados > > estoy considerando programar la funcion de distribucion de Irwin hall. > lamentablemente no he tenido exito, pido que alguien me pueda colaborar con > aquello, les quedo agradecido de antemano. > -- > atte. Lic. Genaro Llusco Silvestre > gellusco en gmail.com > Telf: 74028671 > blog personal es: > http://www.cientificest.blogspot.com > > [[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 > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es
http://en.wikipedia.org/wiki/Irwin%E2%80%93Hall_distribution (Perdón, por error se me ha cortado el mensaje). Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 22 de abril de 2015, 15:34, Carlos J. Gil Bellosta <cgb en datanalytics.com> escribió:> Alternativamente, uno puede implementar la fórmula que aparece en > > > El día 22 de abril de 2015, 12:13, Olivier Nuñez <onunez en unex.es> escribió: >> .... Un código más optimizado para la aproximación Monte Carlo de la distribución de IH >> >> N=10000 # tamaño simulación Monte Carlo >> n=5 # numero de uniformes >> >> cdf.IH <-function(x,n,N) { >> z=replicate(N,sum(runif(n))) >> apply(outer(x,z,">="),1,mean) >> } >> x=seq(0,5,.1) >> y=cdf.IH(x,n=n,N=N) >> plot(x,y,type="l") >> >> Un saludo. Olivier >> >> >> ----- Mensaje original ----- >> De: "Genaro Llusco" <gellusco en gmail.com> >> Para: r-help-es en r-project.org >> Enviados: Martes, 21 de Abril 2015 23:07:14 >> Asunto: [R-es] distribucion de IRWIN HALL >> >> estimados >> >> estoy considerando programar la funcion de distribucion de Irwin hall. >> lamentablemente no he tenido exito, pido que alguien me pueda colaborar con >> aquello, les quedo agradecido de antemano. >> -- >> atte. Lic. Genaro Llusco Silvestre >> gellusco en gmail.com >> Telf: 74028671 >> blog personal es: >> http://www.cientificest.blogspot.com >> >> [[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 >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es
Pues, si hay una formula cerrada sencilla, no dudarlo y olvidar mi aproximación Monte Carlo. Gracias Carlos. Un saludo. Olivier ----- Mensaje original ----- De: "Carlos J. Gil Bellosta" <cgb en datanalytics.com> Para: "Olivier Nuñez" <onunez en unex.es> CC: "Genaro Llusco" <gellusco en gmail.com>, "r-help-es" <r-help-es en r-project.org> Enviados: Miércoles, 22 de Abril 2015 15:34:26 Asunto: Re: [R-es] distribucion de IRWIN HALL Alternativamente, uno puede implementar la fórmula que aparece en El día 22 de abril de 2015, 12:13, Olivier Nuñez <onunez en unex.es> escribió:> .... Un código más optimizado para la aproximación Monte Carlo de la distribución de IH > > N=10000 # tamaño simulación Monte Carlo > n=5 # numero de uniformes > > cdf.IH <-function(x,n,N) { > z=replicate(N,sum(runif(n))) > apply(outer(x,z,">="),1,mean) > } > x=seq(0,5,.1) > y=cdf.IH(x,n=n,N=N) > plot(x,y,type="l") > > Un saludo. Olivier > > > ----- Mensaje original ----- > De: "Genaro Llusco" <gellusco en gmail.com> > Para: r-help-es en r-project.org > Enviados: Martes, 21 de Abril 2015 23:07:14 > Asunto: [R-es] distribucion de IRWIN HALL > > estimados > > estoy considerando programar la funcion de distribucion de Irwin hall. > lamentablemente no he tenido exito, pido que alguien me pueda colaborar con > aquello, les quedo agradecido de antemano. > -- > atte. Lic. Genaro Llusco Silvestre > gellusco en gmail.com > Telf: 74028671 > blog personal es: > http://www.cientificest.blogspot.com > > [[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 > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es