Estoy tratando de ajustar una distribución gamma inversa de tres parámetros a mis datos en R o Python. Me gustaría hacer esto usando la estimación de máxima verosimilitud (MLE). El pdf del gamma inverso de tres parámetros viene dado por: Donde Γ es la función gamma, ρ es la forma, α es la...
Antecedentes : Israel (y el Medio Oriente en general) está lleno de protestas. Tengo curiosidad, cuando se me da una imagen, de estimar cuántas personas hay en esa imagen (a menudo una imagen de una gran multitud). ¿Qué modelado puede ofrecer alguna solución a este problema? (y, por supuesto, se...
Como institución financiera, a menudo nos encontramos con el análisis de datos de series temporales. Muchas veces terminamos haciendo regresión usando variables de series de tiempo. A medida que esto sucede, a menudo encontramos residuos con una estructura de series de tiempo que viola la...
Tengo el plmpaquete y me gustaría ejecutar pruebas de raíz unitaria en algunas variables. Obtuve el siguiente error: > purtest(data$tot.emp) Error in data.frame(baldwin = c(59870, 61259, 60397, 58919, 57856, 57227, : arguments imply differing number of rows: 14, 19, 11, 12, 1, 20, 18, 10,...
¿Cómo completo el cuadrado desde el punto donde lo dejé y es correcto hasta ahora? Tengo un previo normal para ββ\beta de la forma p(β|σ2)∼N(0,σ2V)p(β|σ2)∼N(0,σ2V)p(\beta|\sigma^2)\sim \mathcal{N}(0,\sigma^2V),
Tengo enormes dificultades para evaluar con el paquete en R. En mi caso, los valores de , , son siempre números reales positivos. Aun así, la función hipergeométrica es increíblemente sensible a sus valores. No estoy buscando precisión extrema; Puedo usar Excel para obtener una estimación...
Cuando calculo modelos de regresión con R, utilizo regularmente la función de relevel para que mi modelo también me dé resultados para el otro nivel. Noté que a veces, pero no con frecuencia, esto cambió el modelo en el sentido de que los niveles de otros factores que eran significativos antes del...
Estoy analizando un determinado conjunto de datos y necesito entender cómo elegir el mejor modelo que se ajuste a mis datos. Estoy usando R. Un ejemplo de datos que tengo es el siguiente: corr <- c(0, 0, 10, 50, 70, 100, 100, 100, 90, 100, 100) Estos números corresponden al porcentaje de...
Dado que la estimación posterior de de una probabilidad normal y una gamma inversa antes de es:σ′2σ′2\sigma'^{2}σ2σ2\sigma^2 σ′2∼IG(α+n2,β+∑ni=1(yi−μ)22)σ′2∼IG(α+n2,β+∑i=1n(yi−μ)22)\sigma'^{2}\sim\textrm{IG}\left(\alpha + \frac{n}{2}, \beta +\frac{\sum_{i=1}^n{(y_i-\mu)^2}}{2}\right) que es...
Una variable de respuesta y es una función no lineal de varias variables predictoras X (en mis datos reales la respuesta está distribuida binomialmente, pero aquí estoy usando un valor normalmente distribuido para simplificar). Puedo modelar las relaciones entre los predictores y la respuesta...
Quiero que R muestre los datos que me da de la summary()función en una tabla para poder compartir esto fácilmente. Actualmente solo estoy haciendo summary()en la consola y luego tomando una captura de pantalla, pero prefiero que esto se genere como una buena tabla al igual que todos mis gráficos....
Necesito hacer algunas comparaciones de medias simples entre grupos (pruebas F ANOVA básicas) en datos con valores faltantes. Utilizo el paquete de ratones en R para la imputación múltiple, pero solo puedo agrupar resultados para los coeficientes del modelo lineal, o .R2R2R^2 ¿Alguien sabe cómo...
Estoy calculando algunas probabilidades condicionales y los intervalos de confianza del 95% asociados. Para muchos de mis casos, tengo recuentos directos de xéxitos fuera de los nensayos (de una tabla de contingencia), por lo que puedo usar un intervalo de confianza binomial, como se proporciona...
Supongamos que estoy construyendo un clasificador de regresión logística que predice si alguien está casado o soltero. (1 = casado, 0 = soltero) Quiero elegir un punto en la curva de recuperación de precisión que me brinde al menos un 75% de precisión, por lo que quiero elegir umbralest1t1t_1 y...
Tengo un problema de regresión múltiple, que intenté resolver usando una regresión múltiple simple: model1 <- lm(Y ~ X1 + X2 + X3 + X4 + X5, data=data) Esto parece estar explicando el 85% de la varianza (según R cuadrado) que parece bastante bueno. Sin embargo, lo que me preocupa es la trama...
Actualización : dado que ahora sé que mi problema se llama separación casi completa , actualicé la pregunta para reflejar esto (gracias a Aaron). Tengo un conjunto de datos de un experimento en el que 29 participantes humanos (factor code) trabajaron en un conjunto de ensayos y responsefue 1 o...
Tengo 2 factores Ay B(5 × 3) y una covariable Xen un diseño dentro de la asignatura. Así es como especifico mi modelo general: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Mi interpretación es que estoy mirando un gráfico y~x, donde la pendiente cambia debido a la covariable, y la...
Realmente no sé qué es posible, y me gustaría un puntero en la dirección correcta. Tengo medidas de tiempo y posición que pueden ser cualquier cosa, desde una persona caminando, un vehículo en una carretera, un estacionamiento o una impresora en una oficina. Necesito calcular los tiempos de viaje...
Tengo el modelo que necesito estimar, con y .Y=π0 0+π1X1+π2X2+π3X3+ ε ,Y=π0 0+π1X1+π2X2+π3X3+ε, Y = \pi_0 + \pi_1 X_1 + \pi_2 X_2 + \pi_3 X_3 + \varepsilon, ∑kπk= 1 para k ≥ 1∑kπk=1 para k≥1\sum_k \pi_k = 1 \text{ for }k \geq 1πk≥ 0 para k ≥ 1πk≥0 0 para k≥1\pi_k\ge0 \text{ for }k \geq 1 La...
Un rápido vistazo a la página de ayuda del generador de números aleatorios de espectáculos R que se puede elegir entre 7 generadores predefinidos ( Wichmann-Hill, Marsaglia-Multicarry, Super-Duper, Mersenne-Twister, Knuth-TAOCP-2002, Knuth-TAOCP, L'Ecuyer-CMRG). ?Random El valor predeterminado es...