Articles of lm

Establecer uno o más de los coeficientes a un entero específico

Estoy usando un modelo de lm estándar y me gustaría establecer los coeficientes de una o más de mis variables en un entero específico. Por ejemplo, me gustaría que el coeficiente de mi clima y las variables de precio sean 647 y 15 respectivamente. Estoy usando la función lm con una fórmula estándar. Las cosas […]

Trazar la curva de regresión polinomial en R

Tengo una regresión polinómica simple que hago de la siguiente manera attach(mtcars) fit <- lm(mpg ~ hp + I(hp^2)) Ahora, trazo lo siguiente > plot(mpg~hp) > points(hp, fitted(fit), col=’red’, pch=20) Esto me da lo siguiente Quiero conectar estos puntos en una curva suave, usar líneas me da lo siguiente > lines(hp, fitted(fit), col=’red’, type=’b’) Que […]

¿Cómo extraigo solo el número de un número con nombre (sin el nombre)?

Estoy buscando solo el valor del coeficiente del modelo lineal B1 (newx), no el nombre. Solo quiero el valor 0.5. No quiero el nombre “newx”. newx <- c(0.5,1.5.2.5) newy <- c(2,3,4) out <- lm(newy ~ newx) out ve así: Call: lm(formula = newy ~ newx) Coefficients: (Intercept) newx 1.5 1.0 Llegué aquí. Pero ahora estoy […]

pasando un vector de variables a la fórmula lm ()

Estaba tratando de automatizar una parte de mi código para que la progtwigción sea menos tediosa. Básicamente, estaba tratando de hacer una selección gradual de variables usando fastbw() en el paquete rms . Me gustaría pasar la lista de variables seleccionadas por fastbw() en una fórmula como y ~ x1+x2+x3 , “x1” “x2” “x3” es […]

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 * […]

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) […]

¿Por qué el uso de update en una lm dentro de un data.table agrupado pierde sus datos de modelo?

Ok, esta es una extraña. Sospecho que esto es un error dentro de data.table , pero sería útil si alguien puede explicar por qué sucede esto, ¿qué está haciendo exactamente la update ? Estoy usando la list(list()) truco dentro de data.table para almacenar modelos ajustados. Cuando crea una secuencia de objetos lm cada uno para […]

¿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 […]