similar to: 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

Cordial saludo. Estoy trabajando con una densidad conjunta algo elaborada (suma y cociente de dos variables aleatorias positivas usando la cópula Frank) y requiero encontrar la distribución marginal de cada una de las variables. A manera de ejemplo (no es el caso específico, que tiene mucho más código), si partimos de la distribución exponencial bivariada #----- theta<-0.5

Hola comunidad, tengo el siguiente problema con la funcion stat_density2d (estimacion de densidad 2d.) del paquete ggplot2 stat_density2d (mapping = NULL, data = NULL, geom = "density2d", position = "identity", na.rm = FALSE, contour = TRUE, n = 100, ...) de esta función requiero ver el cálculo númerico de la estimación de la densidad, cuento con el siguiente codigo y

Estimados una consulta me encuentro graficando un histograma cuyos datos no provienen de una distribución clásica como la normal exponenial, poisson, etc, Lo que necesito es colocar una curva no paramétrica que permita evidenciar el ajuste de los datos a esa curva ya que son muchos (alrededor de 80000). Muchas gracias [[alternative HTML version deleted]]

Estimados amigos y expertos del R, Les escribo para hacerles una pregunta que parece un poco sencilla pero me ha costado mucho. Estoy tratando de graficar los datos correspondientes a tasas de consumo de algunos organismos cuando están expuestos o no expuestos a un químico sobre tres tipos de presa. Por definición, debía ajustar los datos a un glm con distribución gama. Las gráficas pueden ser

Hola a tod en s, tengo una duda que se relaciona con alternativas al uso de variables globales. En principio, si se quiere usar un generador de v.a con la librería Runuran sólo se permite definir las funciones de densidad (o el núcleo de las mismas) con funciones con un único argumento en (x). Sin embargo, necesito pasar a las funciones más argumentos que van cambiando en las iteraciones de

Marta, Puedes ver en el código source de ggplot stat_density si es que lo tienes descargado como sorce o puedes ir a GitHub https://github.com/hadley/ggplot2/blob/master/R/stat-density-2d.r y ver como calcula los datos. Por lo que veo usa la función de la librería MASS kde2d, para ver que hace esta escribe kde2d en la consola. Espero te sirva. Daniel El 18 de enero de 2014, 22:35, marta

Hola Alguien me podría decir como hacer una grafica del tipo persp() de la densidad de una distribucion normal bivariante estandarzada con correlacion 0.5?? gracias [[alternative HTML version deleted]]

Hola, Manuel, No entiendo tu pregunta (la repito aqui para que sea mas explicito): hay alguna forma de comparar la probabilidad de supervivencias (en este caso anual) entre grupos sin utilizar un chi-cuadrado y un valor de P. Entiendo que lo que hace survdiff es comparar las curvas de supervivencia, pero yo quiero comparar la probabilidad de supervivencia entre grupos al final del estudio. Con

hola todosnecesito hacer un gráfica de este tipo, alguna idea ??? gracias Patricio F. ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130603/f3fc33e3/attachment-0001.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: densidad.jpg Type: image/jpeg

Hello, I have this code to plot a certain normal distribution and represent the pnorm value for a certain "x": x<-300 xx <- seq(2.5,7.5, by=0.1) yy <- dnorm(xx,5.01,0.77) d<-signif(pnorm(log(x), 5.01,0.77),4) xpts <- round(exp(0:8)) par(bg = "antiquewhite") plot(xx,yy, type="l", col="blue", lwd=2, xaxt="n",

Hola amigos, Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en el ejemplo de la ayuda: DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) summary(fm1DNase1) Aquí se está modelando la densidad óptica de un ensayo relacionada de forma no lineal (logística) con (el logaritmo) de la concentración de una proteína.

Creo que el problema se podría simplificar y aún obtener un resultado aceptable si simplemente haces un ajuste lineal en un entorno del punto. Como ya tienes los valores de los puntos de cada curva como dices Carlos, creo que con esa información podrías identificar un entorno adecuado en el cual ajustar linealmente. Cuéntanos como te va, Saludos !! Eric. On 24-08-20 12:52, Carlos Ortega

Buenas tardes, tengo una variable bimodal (*var)*, de presencias y ausencias (1s y 0s) y otra variable, *prob*, con las probabilidades que le asigna un modelo (entre 0 y 1). Con: *ggplot(Preds, aes(x=prob, fill= var )) + geom_density(alpha=.3)* obtengo la distribución de las presencias y de las ausencias, por separado, en función del valor de probabilidad asignado. Las dos curvas se cruzan en un

Hola, Estoy haciendo un PCA con el paquete ade4. En el gráfico, las etiquetas de las especies suelen quedar sobrepuestas unas sobre las otras y no se pueden distinguir individualmente. ¿Hay alguna manera de generar un poco de espacio entre ellas para poder visualizarlas todas? gracias Yésica [[alternative HTML version deleted]]

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150128/9e23caec/attachment.html>

Estimada comunidad, nuevamente con algo que no se como hacer en R, pero aprendiendo .... El problema es el siguiente: 1. tengo un conjunto de 10 mil datos (n) cuyo valor va entre 0 y 10.000 a los que aplico una funcion sm.density() para obtener una especie de histograma pero con una curva que parece continua. 2. Una vez que obtengo esa curva necesito calcular el area bajo la curva en ciertos

Por favor ver comentario con las las letras mayúsculas data("CPS1985", package = "AER") cps <- CPS1985 cps [1:10,] dim( cps) names( cps) library("quantreg") #regresion usual cps_lm <- lm(log(wage) ~ experience + I(experience^2) + + education, data = cps) #regresion por quantiles, tau son los cuantiles

Hola, Sí, aquí tienes la solución: http://stackoverflow.com/questions/3297552/how-to-plot-a-set-of-densities-in-3d-using-r/3297815#3297815 Saludos, Carlos Ortega www.qualityexcellence.es El 3 de junio de 2013 18:14, Patricio Fuenmayor Viteri < patricio.fuenmayor@outlook.com> escribió: > hola todos > necesito hacer un gráfica de este tipo, alguna idea ??? > > gracias >

Matlab hace ese tipo de gráficos. Y no sé hasta qué punto son o no buenos para representar el fenómeno subyacente. Yo me decantaría por una alternativa al estilo de http://www.bankofengland.co.uk/publications/PublishingImages/inflationreport/cpimktmay13large.gif Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 3 de junio de 2013 18:14, Patricio Fuenmayor Viteri