Я новичок в R. У меня есть набор данных с 245 переменными и 32 регионами или географическими c областями. Мне нужно получить коэффициенты модели ols для каждого из 32 регионов. Переменная, содержащая информацию для 32 состояний, - это ubica_geo_enh. Набор данных включает более 25 635 наблюдений за один год.
Я думал, что этот код будет работать:
y<- c(esc_inf) #dependent variable
x1<- c(esc_pad) #independent variable
x2<- c(esc_mad) #independent variable
edo<- c(ubica_geo_enh) #This variable goes fron 1:32
predictorlist<- list("x1","x2")
for (i in edo){
model[i] <- lm(y[i] ~ x1[i] + x2[i])
}
summary(model_[i])
Он работает, но возвращает результат для всего набора данных, а не каждый отдельный регион (1:32) Буду очень признателен, если кто-то может мне помочь. Заранее спасибо