La aplicación directa de los métodos de bootstrap a la prueba de hipótesis es estimar el intervalo de confianza del estadístico de prueba al calcularlo repetidamente en las muestras de bootstrap (deje que el estadístico muestreado de bootstrap se llame ). Rechazamos si el parámetro hipotético (que...
Tengo una pregunta sobre la distribución posterior de Dirichlet. Dada una función de probabilidad multinomial, se sabe que el posterior es , donde es el número de veces que hemos visto observación.Dir(αi+Ni)Dir(αi+Ni)Dir({\alpha_i + N_i})NiNiN_iithithi^{th} ¿Qué sucede si comenzamos a disminuir s...
Quiero verificar si realmente entendí el análisis factorial [clásico, lineal] (FA), especialmente los supuestos que se hacen antes (y posiblemente después) de FA. Algunos de los datos deben estar inicialmente correlacionados y existe una posible relación lineal entre ellos. Después de hacer el...
Me gustaría usar lme4para ajustar una regresión de efectos mixtos y multcompcalcular las comparaciones por pares. Tengo un conjunto de datos complejo con múltiples predictores continuos y categóricos, pero mi pregunta se puede demostrar utilizando el ChickWeightconjunto de datos incorporado como...
La documentación de R para cualquiera no arroja mucha luz. Todo lo que puedo obtener de este enlace es que usar cualquiera de los dos debería estar bien. Lo que no entiendo es por qué no son iguales. Hecho: La función de regresión gradual en R, step()usa extractAIC(). Curiosamente, ejecutar un...
De Wikipedia Una Red Bayesiana Dinámica (DBN) es una Red Bayesiana que relaciona variables entre sí en pasos de tiempo adyacentes. Esto a menudo se llama un BN de dos tiempos porque dice que en cualquier momento T, el valor de una variable se puede calcular a partir de los regresores internos y...
Digamos que tenemos una pregunta simple de "sí / no" a la que queremos saber la respuesta. Y hay N personas "votando" por la respuesta correcta. Cada votante tiene un historial: una lista de 1 y 0, que muestra si tenían razón o no sobre este tipo de preguntas en el pasado. Si asumimos la historia...
Para analizar los recuentos de aves inflados a cero, me gustaría aplicar modelos de recuento inflados a cero usando el paquete R pscl . Sin embargo, al observar el ejemplo proporcionado en la documentación para una de las funciones principales ( ? Zeroinfl ), empiezo a dudar de cuál es la verdadera...
La implementación randomForest no permite el muestreo más allá del número de observaciones, incluso cuando se realiza un muestreo con reemplazo. ¿Por qué es esto? Funciona bien: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris,...
Encontré esta distribución en un juego de computadora y quería aprender más sobre su comportamiento. Viene de la decisión de si un determinado evento debe ocurrir después de un número determinado de acciones del jugador. Los detalles más allá de esto no son relevantes. Parece aplicable a otras...
Fondo biologico Con el tiempo, algunas especies de plantas tienden a duplicar sus genomas completos, obteniendo una copia adicional de cada gen. Debido a la inestabilidad de esta configuración, muchos de estos genes se eliminan y el genoma se reorganiza y se estabiliza, listo para duplicarse...
Un ejemplo de una buena medida de la separabilidad de clases en estudiantes discriminantes lineales es la relación discriminante lineal de Fisher. ¿Existen otras métricas útiles para determinar si los conjuntos de características proporcionan una buena separación de clases entre las variables de...
Me gustaría describir el "pico" y el "peso" de la cola de varias funciones de densidad de probabilidad sesgadas. Las características que quiero describir, ¿se llamarían "curtosis"? ¿Solo he visto la palabra "curtosis" usada para distribuciones
Probablemente estoy lidiando con un problema que probablemente se haya resuelto cientos de veces antes, pero no estoy seguro de dónde encontrar la respuesta. Cuando mediante regresión logística, dadas muchas características tratando de predecir un valor categórico binario y , estoy interesado en...
Estoy entrenando un proceso gaussiano con un núcleo ARD con muchos parámetros maximizando la capacidad marginal de los datos, en lugar de la validación cruzada. Sospecho que es demasiado adecuado. ¿Cómo puedo probar esta sospecha en un contexto
Con respecto a la ponderación de puntaje de propensión (IPTW) al realizar el modelado de riesgo proporcional de Cox de los datos de supervivencia del tiempo hasta el evento: Tengo datos de registro prospectivos en los que estamos interesados en observar el efecto del tratamiento de un...
En el documento seminal "Filtrado de partículas Rao-Blackwellised para redes bayesianas dinámicas" por A. Doucet et. Alabama. Se propone un filtro secuencial Monte Carlo (filtro de partículas), que utiliza una subestructura lineal en un proceso de Markov x k = ( x L k , x N k ) . Al marginar esta...
Estoy intentando modelar una variable de respuesta que está teóricamente limitada entre -225 y +225. La variable es la puntuación total que obtuvieron los sujetos al jugar un juego. Aunque en teoría es posible que los sujetos obtengan un puntaje de +225. A pesar de esto porque el puntaje no solo...
Estoy impartiendo un curso sobre estadísticas básicas y estamos haciendo la prueba t para dos muestras independientes con variaciones desiguales (prueba de Welch). En los ejemplos que he visto, los grados ajustados de libertad utilizados por la prueba de Welch son siempre menores o iguales que ....
Para los datos de recuento que he recopilado, uso la regresión de Poisson para construir modelos. Hago esto usando la glmfunción en R, donde uso family = "poisson". Para evaluar posibles modelos (tengo varios predictores) utilizo el AIC. Hasta aquí todo bien. Ahora quiero realizar una validación...