Articles of regresión

Especificando la fórmula en R con glm sin statement explícita de cada covariable

Me gustaría forzar variables específicas en regresiones glm sin especificar completamente cada una. Mi conjunto de datos real tiene ~ 200 variables. No he podido encontrar muestras de esto en mi búsqueda en línea hasta el momento. Por ejemplo (con solo 3 variables): n=200 set.seed(39) samp = data.frame(W1 = runif(n, min = 0, max = […]

Rápida regresión lineal simple por pares entre variables en un dataframe

He visto la regresión lineal simple por pares o general emparejada muchas veces en Desbordamiento de stack. Aquí hay un conjunto de datos de juguetes para este tipo de problema. set.seed(0) X <- matrix(runif(100), 100, 5, dimnames = list(1:100, LETTERS[1:5])) b <- c(1, 0.7, 1.3, 2.9, -2) dat <- X * b[col(X)] + matrix(rnorm(100 * […]

lme4 :: lmer informa que “la matriz del modelo de efectos fijos tiene una clasificación deficiente”, ¿necesito una solución y cómo hacerlo?

Estoy tratando de ejecutar un modelo de efectos mixtos que prediga F2_difference con el rest de las columnas como predictores, pero recibo un mensaje de error que dice: la matriz del modelo de efectos fijos tiene un rango deficiente, por lo que cae 7 columnas / coeficientes. Desde este enlace, el modelo de efectos fijos […]

predict.lm () con un nivel de factor desconocido en los datos de prueba

Estoy ajustando un modelo para factorizar datos y predecir. Si newdata en predict.lm() contiene un único nivel de factor desconocido para el modelo, todo predict.lm() falla y devuelve un error. ¿Hay una buena manera de que predict.lm() devuelva una predicción para los niveles de factores que conoce el modelo y NA para los niveles de […]

Extraer valores del coeficiente de regresión

Tengo un modelo de regresión para algunos datos de series de tiempo que investigan la utilización de medicamentos. El objective es ajustar una spline a una serie temporal y calcular un 95% de IC, etc. El modelo es el siguiente: id <- ts(1:length(drug$Date)) a1 <- ts(drug$Rate) a2 <- lag(a1-1) tg <- ts.union(a1,id,a2) mg <-lm (a1~a2+bs(id,df=df1),data=tg) […]

¿Cómo forzar a R a usar un nivel de factor específico como referencia en una regresión?

¿Cómo puedo decirle a R que use un cierto nivel como referencia si uso variables explicativas binarias en una regresión? Solo está usando algún nivel por defecto. lm(x ~ y + as.factor(b)) con b {0, 1, 2, 3, 4} . Digamos que quiero usar 3 en lugar del cero que usa R.

¿Cómo calcula predict.lm () el intervalo de confianza y el intervalo de predicción?

Ejecuté una regresión: CopierDataRegression <- lm(V1~V2, data=CopierData1) y mi tarea era obtener una Intervalo de confianza del 90% para la respuesta media dada V2=6 y Intervalo de predicción del 90% cuando V2=6 . Use el siguiente código: X6 <- data.frame(V2=6) predict(CopierDataRegression, X6, se.fit=TRUE, interval="confidence", level=0.90) predict(CopierDataRegression, X6, se.fit=TRUE, interval="prediction", level=0.90) y obtuve (87.3, 91.9) y […]

Montaje de un modelo lineal con múltiples LHS

Soy nuevo en R y quiero mejorar el siguiente script con una función *apply (He leído sobre apply , pero no pude lograr usarlo). Quiero usar la función lm en múltiples variables independientes (que son columnas en un dataframe). solía for (i in (1:3) { assign(paste0(‘lm.’,names(data[i])), lm(formula=formula(i),data=data)) } Formula(i) se define como formula=function(x) { as.formula […]

Regresión lineal y agrupar por en R

Quiero hacer una regresión lineal en R usando la función lm() . Mis datos son series temporales anuales con un campo por año (22 años) y otro por estado (50 estados). Quiero ajustar una regresión para cada estado para que al final tenga un vector de respuestas de lm. Puedo imaginar hacer un ciclo para […]

¿Cómo se puede depurar el error “los contrastes solo se pueden aplicar a factores con 2 o más niveles”?

Aquí están todas las variables con las que estoy trabajando: str(ad.train) $ Date : Factor w/ 427 levels “2012-03-24″,”2012-03-29”,..: 4 7 12 14 19 21 24 29 31 34 … $ Team : Factor w/ 18 levels “Adelaide”,”Brisbane Lions”,..: 1 1 1 1 1 1 1 1 1 1 … $ Season : int 2012 […]