He entrenado un modelo de regresión lineal, utilizando un conjunto de variables / características. Y el modelo tiene un buen rendimiento. Sin embargo, me he dado cuenta de que no hay una variable con una buena correlación con la variable predicha. ¿Como es
Tengo datos de pacientes tratados con 2 tipos diferentes de tratamientos durante la cirugía. Necesito analizar su efecto sobre la frecuencia cardíaca. La medición del ritmo cardíaco se toma cada 15 minutos. Dado que la duración de la cirugía puede ser diferente para cada paciente, cada paciente...
Se puede realizar una regresión logit en R usando dicho código: > library(MASS) > data(menarche) > glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, + family=binomial(logit), data=menarche) > coefficients(glm.out) (Intercept) Age -21.226395 1.631968 Parece que el algoritmo de...
Realicé una evaluación por computadora de diferentes métodos para ajustar un tipo particular de modelo utilizado en las ciencias paleo. Tenía un conjunto de entrenamiento de gran tamaño y, por lo tanto, al azar (muestreo aleatorio estratificado) dejé de lado un conjunto de prueba. Ajusté diferentes...
Considere el modelo lineal simple: yy=X′ββ+ϵyy=X′ββ+ϵ\pmb{y}=X'\pmb{\beta}+\epsilon donde ϵi∼i.i.d.N(0,σ2)ϵi∼i.i.d.N(0,σ2)\epsilon_i\sim\mathrm{i.i.d.}\;\mathcal{N}(0,\sigma^2) y X∈Rn×pX∈Rn×pX\in\mathbb{R}^{n\times p} ,p≥2p≥2p\geq2 yXXX contiene una columna de constantes. Mi pregunta es, dado...
¿Existe tal paquete que proporcione la estimación del modelo de efectos mixtos binomiales negativos inflados a cero en R? Con eso quiero decir: Inflación cero donde puede especificar el modelo binomial para inflación cero, como en la función zeroinfl en el paquete pscl: zeroinfl (y ~ X | Z, dist...
Digamos que está en la biblioteca de su departamento de estadística y que se encuentra con un libro con la siguiente imagen en la portada. Probablemente pensarás que este es un libro sobre cosas de regresión lineal. ¿Cuál sería la imagen que te haría pensar en modelos lineales...
Estamos trabajando con algunas regresiones logísticas y nos hemos dado cuenta de que la probabilidad estimada promedio siempre es igual a la proporción de unos en la muestra; es decir, el promedio de los valores ajustados es igual al promedio de la muestra. ¿Alguien puede explicarme la razón o...
Actualización : Perdón por otra actualización, pero he encontrado algunas posibles soluciones con polinomios fraccionales y el paquete de riesgos de la competencia con el que necesito ayuda. El problema No puedo encontrar una manera fácil de hacer un análisis de coeficiente dependiente del...
La regresión beta (es decir, GLM con distribución beta y generalmente la función de enlace logit) a menudo se recomienda para tratar la respuesta, también conocida como variable dependiente que toma valores entre 0 y 1, como fracciones, razones o probabilidades: regresión para un resultado...
Pregunta simple: ¿Cómo especificar una distribución lognormal en el argumento de la familia GLM en R? No pude encontrar cómo se puede lograr esto. ¿Por qué lognormal (o exponencial) no es una opción en el argumento familiar? En algún lugar de los R-Archives leí que uno simplemente tiene que usar...
He estado aprendiendo sobre el modelo de riesgo proporcional de Cox. Tengo mucha experiencia ajuste modelos de regresión logística, y los modelos de manera de construir la intuición que he estado comparando ajuste usando coxphDel R "supervivencia" con modelos de regresión logística aptos utilizando...
¡Estoy realmente sorprendido por el hecho de que Poisson GLM acepta números no enteros! Mira: Datos (contenido de data.txt): 1 2001 0.25 1 1 2002 0.5 1 1 2003 1 1 2 2001 0.25 1 2 2002 0.5 1 2 2003 1 1 R script: t <- read.table("data.txt") names(t) <- c('site', 'year', 'count',...
Después de instalar un modelo cox, es posible hacer predicciones y recuperar el riesgo relativo de nuevos datos. Lo que no entiendo es cómo se calcula el riesgo relativo para un individuo y a qué se relaciona (es decir, el promedio de la población). ¿Alguna recomendación sobre recursos para ayudar...
Tengo datos de densidad de peces que estoy tratando de comparar entre varias técnicas de recolección diferentes, los datos tienen muchos ceros y el histograma parece vaugley apropiado para una distribución de Poisson, excepto que, como densidades, no son datos enteros. Soy relativamente nuevo en...
Tengo seis variables dependientes (datos de conteo) y varias variables independientes, veo que en un MMR el script es así: my.model <- lm(cbind(DV1,DV2,DV3,DV4,DV5,DV6) ~ IV1 + IV2 + ... + IVn) Pero, dado que mis datos son recuentos, quiero usar un modelo lineal generalizado y probé...
Estoy siguiendo AW van der Vaart, estadísticas asintóticas (1998). Habla de experimentos estadísticos, alegando que son diferentes de un modelo estadístico, pero no define ninguno. Mi pregunta: ¿Qué es (1) un experimento estadístico, (2) un modelo estadístico y (3) cuál es el ingrediente clave que...
Tengo un modelo lineal clásico, con 5 posibles regresores. No están correlacionados entre sí y tienen una correlación bastante baja con la respuesta. Llegué a un modelo donde 3 de los regresores tienen coeficientes significativos para su estadística t (p <0.05). Agregar una o las dos variables...
Tengo problemas para entender el resultado de mi lmer()modelo. Es un modelo simple de una variable de resultado (Apoyo) con interceptaciones estatales / efectos aleatorios estatales variables: mlm1 <- lmer(Support ~ (1 | State)) Los resultados de summary(mlm1)son: Linear mixed model fit by...
He estado usando el MCMCglmmpaquete recientemente. Estoy confundido por lo que se refiere en la documentación como estructura R y estructura G. Estos parecen relacionarse con los efectos aleatorios, en particular especificando los parámetros para la distribución previa sobre ellos, pero la...