Displaying 2 results from an estimated 2 matches for "range_sp".
2013 Sep 16
1
Solapamiento de funciones de densidad
...o de dos
funciones de densidad (que creo que es lo que quieres hacer), es esta:
kernel.dist<-function(sp1,sp2) {
n.obs<-c(length(sp1),length(sp2))
sdevs<-c(sd(sp1,na.rm=TRUE),sd(sp2,na.rm=TRUE))
means<-c(mean(sp1,na.rm=TRUE),mean(sp2,na.rm=TRUE))
h<-1.06*sdevs*(n.obs^-0.2)
range_sp<-range(c(means[1]+5*sdevs[1],means[1]-5*sdevs[1],means[2]+5*sdevs[2
],means[2]-5*sdevs[2]))
length.grad<-seq(range_sp[1],range_sp[2],length.out=128)
pdf_x<-density(sp1,
bw=h[1],kernel="gaussian",from=range_sp[1],to=range_sp[2],n=128,na.rm=TRUE)$
x
pdf_sp1<-density(sp1,
b...
2013 Sep 15
1
Solapamiento de funciones de densidad
Hola a todos.
Estoy interesado en obtener el grado de solapamiento de dos parametros obtenidos a partir de una MCMC.
La gente que he consultado realiza inspecciones visuales (graficando ambas funciones de densidad de los valores obtenidos y viendo si hay mucho o poco solapamiento), pero me gustaría incluirlo en una serie de simulaciones con lo que no podría realizarlas.
Conoceia alguna manera