Quiero hacer lo siguiente: 1) Regresión OLS (sin término de penalización) para obtener coeficientes beta ; representa las variables utilizadas para la regresión. Hago esto porb∗jbj∗b_{j}^{*}jjj lm.model = lm(y~ 0 + x) betas = coefficients(lm.model) 2) Regresión de lazo con un término de...
Cuando se modelan proporciones continuas (p. Ej., Cobertura vegetal proporcional en los cuadrantes de la encuesta, o proporción de tiempo dedicado a una actividad), la regresión logística se considera inapropiada (p. Ej., Warton y Hui (2011). El arcoseno es estúpido: el análisis de proporciones en...
Estoy intentando el análisis factorial confirmatorio (CFA) usando lavaan. Me está costando interpretar el resultado producido por lavaan. Tengo un modelo simple: 4 factores, cada uno respaldado por elementos de los datos recopilados de la encuesta. Los factores están en línea con lo que miden los...
He construido un modelo glm en R y lo he probado usando un grupo de prueba y entrenamiento, así que estoy seguro de que funciona bien. Los resultados de R son: Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -2.781e+00 1.677e-02 -165.789 < 2e-16 *** Coeff_A 1.663e-05...
Estoy tratando de especificar una fórmula para un modelo lineal de efectos mixtos (con lme4) para mi diseño experimental, pero no estoy seguro de que lo esté haciendo bien. El diseño: básicamente estoy midiendo un parámetro de respuesta en plantas. Tengo 4 niveles de tratamiento y 2 niveles de...
Estoy usando el RandomForestpaquete R y estoy confundido sobre cómo interpretar los valores del eje Y en sus gráficos de dependencia parcial. Los documentos de ayuda indican que la gráfica es una "representación gráfica del efecto marginal de una variable en la probabilidad de la clase". Sin...
Estoy buscando un método de agrupamiento óptimo (discretización) de una variable continua con respecto a una variable binaria de respuesta (objetivo) dada y con un número máximo de intervalos como parámetro. ejemplo: tengo un conjunto de observaciones de personas con variables de "altura" (numeral...
Tengo un marco de datos que contiene valores en 4 columnas: Por ejemplo: ID, price, click count,rating Lo que me gustaría hacer es "dividir" este marco de datos en N grupos diferentes donde cada grupo tendrá el mismo número de filas con la misma distribución de precios, recuento de clics y...
Estoy leyendo el documento de teoría de Doug Bates sobre el paquete lme4 de R para comprender mejor la esencia de los modelos mixtos, y encontré un resultado intrigante que me gustaría entender mejor, sobre el uso de la máxima verosimilitud restringida (REML) para estimar la varianza . En la...
Esto es similar a Bootstrap: la estimación está fuera del intervalo de confianza Tengo algunos datos que representan recuentos de genotipos en una población. Quiero estimar la diversidad genética usando el índice de Shannon y también generar un intervalo de confianza usando bootstrapping. Sin...
Estoy leyendo el libro de Kevin Murphy: Aprendizaje automático: una perspectiva probabilística. En el primer capítulo, el autor explica la maldición de la dimensionalidad y hay una parte que no entiendo. Como ejemplo, el autor afirma: Considere que las entradas están distribuidas uniformemente a...
Me pregunto cómo puedo simular una muestra de n vidas de distribución de Weibull que incluya observaciones censuradas a la derecha de Tipo I. Por ejemplo, tengamos n = 3, forma = 3, escala = 1 y la tasa de censura = .15, y el tiempo de censura = .88. Sé cómo generar una muestra de Weibull pero no...
Estoy revisando la sección LAB §6.6 sobre Regresión de cresta / lazo en el libro 'Una introducción al aprendizaje estadístico con aplicaciones en R' de James, Witten, Hastie, Tibshirani (2013). Más específicamente, estoy tratando de aplicar el Ridgemodelo scikit-learn al conjunto de datos...
Identifiqué varios lugares en los libros de texto donde el GLM se describe con 5 distribuciones (a saber, gamma, gaussiano, binomial, gaussiano inverso y Poisson). Esto también se ejemplifica en la función familiar en R. Ocasionalmente me encuentro con referencias al GLM donde se incluyen...
Me gustaría simular un proceso de excursión browniana (un movimiento browniano que está condicionado siempre será positivo cuando 0<t<10<t<10 \lt t \lt 1 a 000 en t=1t=1t=1 ). Dado que un proceso de excursión browniana es un puente browniano que está condicionado a ser siempre positivo,...
Estoy tratando de modelar las intensidades medias de parásitos que afectan a un huésped en R usando un modelo binomial negativo. Sigo recibiendo 50 o más advertencias que dicen: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 ¿Cómo puedo lidiar con esto? Mi código se ve así:...
Cuando me enteré del previo de Jeffreys en mi clase de inferencia estadística de posgrado, mis profesores hicieron que pareciera interesante, principalmente por razones históricas y no porque alguien lo usaría alguna vez. Luego, cuando tomé el análisis de datos bayesianos, nunca se nos pidió que...
Estoy leyendo un artículo cuyo método se basa completamente en la prueba de razón de probabilidad. El autor dice que la prueba LR contra alternativas unilaterales es UMP. Él procede afirmando que "... incluso cuando [la prueba LR] no se puede demostrar que es uniformemente más potente, la prueba...
Estoy analizando una tabla de 2x2 de un pequeño conjunto de datos de 30 pacientes. Estamos intentando retrospectivamente encontrar algunas variables que den una pista sobre qué tratamiento elegir. Las variables (obs normal / extraño) y la decisión de tratamiento (A / B) son de especial interés y,...
Lo siguiente es un extracto de la Introducción a las estadísticas bayesianas de Bolstad . Para todos los expertos, esto puede ser trivial, pero no entiendo cómo el autor concluye que no tenemos que hacer ninguna integración para calcular la probabilidad posterior de algún valor de . Entiendo la...