Tengo que estimar un modelo lineal simple con una variable dependiente que es un índice (ir) y una variable independiente que son precipitaciones de una estación meteorológica. En total son 18 estaciones y tengo que obtener los r cuadrado de cada una de esas regresiones. La pregunta entonces es, ¿cómo puedo estimar todas esas regresiones dejando fija la variable dependiente y cambiando sólo la independiente (las estaciones) ? Lo que se me ocurre es un for que recorra un array que contenga los nombres de las estaciones.
df <- data.frame(data)
vec <- data.frame(colnames(df))
for (i in vec) {
random_model = lm(data = df, formula = ir ~ i )
summary(random_model)
}
pero me arroja un error :
> for (i in vec) {
+
+ random_model = lm(data = df, formula = ir ~ i )
+ summary(random_model)
+ }
Error in model.frame.default(formula = ir ~ i, data = df, drop.unused.levels = TRUE) :
variable lengths differ (found for 'i')