Antecedentes Tengo datos de un estudio de campo en el que hay cuatro niveles de tratamiento y seis réplicas en cada uno de los dos bloques. (4x6x2 = 48 observaciones) Los bloques están separados aproximadamente 1 milla, y dentro de los bloques, hay una cuadrícula de 42 parcelas de 2m x 4m y una...
Tengo una pregunta sobre si usar o no un desplazamiento. Suponga un modelo muy fácil, donde desea describir el número (general) de goles en el hockey. Entonces tienes goles, número de juegos jugados y un "delantero" ficticio variable que es igual a 1 si el jugador es delantero y 0 en caso...
Tengo una pregunta sobre la validez del uso de RMSE (Root Mean Squared Error) para comparar diferentes modelos logísticos. La respuesta es bien 0o 1y las predicciones son las probabilidades entre 0- 1? ¿La forma aplicada a continuación también es válida con las respuestas binarias? # Using...
Tengo el siguiente modelo lineal: Para abordar la heteroscedasticidad residual, he tratado de aplicar una transformación logarítmica en la variable dependiente como Iniciar sesión( Y+ 1 )Iniciar sesión(Y+1)\log(Y + 1) pero todavía veo el mismo efecto de abanico en los residuos. Los valores de...
Estoy revisando la regresión lineal. El libro de texto de Greene dice: Ahora, por supuesto, habrá otros supuestos sobre el modelo de regresión lineal, como E(ϵ|X)=0E(ϵ|X)=0E(\epsilon|X)=0 . Esta suposición combinada con la suposición de linealidad (que en efecto define ϵϵ\epsilon ), le da...
El modelo de regresión de Poisson inflado a cero se define para una muestra por y además supone que los parámetros \ mathbf {\ lambda} = (\ lambda_1, \ dots, \ lambda_n) y \ textbf {p} = (p_1, \ dots, p_n) satisfacen(y1,…,yn)(y1,…,yn)(y_1,\ldots,y_n)λ = ( λ 1 , ... , λ n )Yi={0kwith...
Tengo un GLMM de la forma: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Cuando lo uso drop1(model, test="Chi"), obtengo resultados diferentes a los que uso Anova(model, type="III")del paquete del automóvil o summary(model). Estos dos...
Supongamos que tengo 10 estudiantes, cada uno de los cuales intenta resolver 20 problemas de matemáticas. Los problemas se califican como correctos o incorrectos (en datos largos) y el rendimiento de cada alumno se puede resumir mediante una medida de precisión (en subjdatos). Los modelos 1, 2 y 4...
Estoy enseñando una clase sobre integración de funciones de varias variables y cálculo de vectores este semestre. La clase está compuesta por la mayoría de los estudiantes de economía y de ingeniería, con un puñado de matemática y física. Enseñé esta clase el semestre pasado, y descubrí que muchos...
Una forma de encontrar la precisión del modelo de regresión logística usando 'glm' es encontrar el gráfico AUC. ¿Cómo verificar lo mismo para el modelo de regresión encontrado con la variable de respuesta continua (familia = 'gaussiano')? ¿Qué métodos se utilizan para verificar qué tan bien se...
¿Alguien sabe cómo calcular (o extraer) el apalancamiento y las distancias de Cook para un merobjeto de clase (obtenido a través del lme4paquete)? Me gustaría trazar estos para un análisis de
Estoy tratando de ajustar un modelo de tiempo discreto en R, pero no estoy seguro de cómo hacerlo. He leído que puede organizar la variable dependiente en diferentes filas, una para cada observación de tiempo, y usar la glmfunción con un enlace logit o cloglog. En este sentido, tengo tres...
En el siguiente código realizo una regresión logística en datos agrupados usando glm y "a mano" usando mle2. ¿Por qué la función logLik en R me da una probabilidad de registro logLik (fit.glm) = - 2.336 que es diferente del logLik (fit.ml) = - 5.514 que obtengo a mano? library(bbmle) #successes...
Simplemente estoy tratando de recalcular con dnorm () la probabilidad de registro proporcionada por la función logLik de un modelo lm (en R). Funciona (casi perfectamente) para un gran número de datos (por ejemplo, n = 1000): > n <- 1000 > x <- 1:n > set.seed(1) > y <- 10 +...
Vi esta notación para los mínimos cuadrados ordinarios aquí . minw∥Xw−y∥22minw‖Xw−y‖22 \min_w \left\| Xw - y \right\|^2_2 Nunca he visto las barras dobles y los 2 en la parte inferior. ¿Que significan esos símbolos? ¿Tienen terminología específica para
Tengo una muestra de 1,449 puntos de datos que no están correlacionados (r-cuadrado 0,006). Al analizar los datos, descubrí que al dividir los valores de las variables independientes en grupos positivos y negativos, parece haber una diferencia significativa en el promedio de la variable...
Tanto bayesglm()(en el paquete arm R) como varias funciones en el paquete MCMCpack están destinadas a hacer una estimación bayesiana de modelos lineales generalizados, pero no estoy seguro de que realmente estén calculando lo mismo. Las funciones MCMCpack utilizan la cadena de Markov Monte Carlo...
Supongamos que tengo un NNN -vector YYY de variables dependientes, y un NNN -vector XXX de variable independiente. Cuando YYY se traza contra 1X1X\frac{1}{X} , veo que hay una relación lineal (tendencia al alza) entre los dos. Ahora, esto también significa que hay una tendencia a la baja lineal...
¿Cuáles son las diferencias en el uso de modelos lineales generalizados, como la determinación automática de relevancia (ARD) y la regresión de Ridge, en comparación con modelos de series de tiempo como Box-Jenkins (ARIMA) o suavizado exponencial para el pronóstico? ¿Existen reglas generales sobre...
Estoy ejecutando glms en R (modelos lineales generalizados). Pensé que conocía los valores, hasta que vi que llamar a un resumen para un glm no le da un valor de referencia superior del modelo en su conjunto, al menos no en el lugar donde lo hacen los modelos lineales. Me pregunto si esto se da...