Lặp qua các biến số trong hồi quy bằng cách sử dụng R

Question

Tôi đang cố gắng chạy 96 lần hồi quy và lưu kết quả dưới dạng 96 đối tượng khác nhau. Để làm phức tạp mọi thứ, tôi muốn chỉ số trên một trong các biến số trong mô hình cũng thay đổi 96 lần. Tôi đã gần như giải quyết vấn đề nhưng tôi đã không may trúng một bức tường. Mã cho đến nay là,

for(i in 1:96){ 

    assign(paste("z.out", i,sep=""), lm(rMonExp_EGM~ TE_i + Month2+Month3+Month4+Month5+Month6+Month7+Month8+Month9+ 
    Month10+Month11+Month12+Yrs_minus_2004 + 
    as.factor(LGA),data=Pokies)) 

} 

này hoạt động ở phía bên tạo đối tượng (ví dụ tôi có z.out1 - z.out96) nhưng tôi dường như không thể có được subscript trên covariate thay đổi là tốt.

Tôi có 96 biến được gọi là TE_1, TE_2 ... TE_96 trong tập dữ liệu. Như vậy, chỉ số trên TE_, "i" cần thay đổi để tương ứng với từng đối tượng tôi tạo. Đó là, z.out1 nên giữ kết quả từ mô hình này:

z.out1 <- lm(rMonExp_EGM~ TE_1 + Month2+Month3+Month4+Month5+Month6+Month7+Month8+Month9+ 
    Month10+Month11+Month12+Yrs_minus_2004 + as.factor(LGA),data=Pokies) 

Và z.out96 nên là:

z.out96 <- lm(rMonExp_EGM~ TE_96+ Month2+Month3+Month4+Month5+Month6+Month7+Month8+Month9+ 
    Month10+Month11+Month12+Yrs_minus_2004 + as.factor(LGA),data=Pokies) 

Hy vọng rằng điều này làm cho ý nghĩa. Tôi rất biết ơn mọi lời khuyên/lời khuyên.

Answer 1

tôi sẽ đưa các kết quả trong một danh sách và tránh những for loop và assign báo cáo

Bạn có thể sử dụng một sự kết hợp của reformulate và update để tạo ra công thức của bạn

orig_formula <- MonExp_EGM~ Month2+Month3+Month4+Month5+Month6+Month7+Month8+Month9+ 
Month10+Month11+Month12+Yrs_minus_2004 + as.factor(LGA) 


te_variables <- paste0('TE_', 1:96) 
# Or if you don't have a current version of R 
# te_variables <- paste('TE', 1:96, sep = '_') 

new_formula <- lapply(te_variables, function(x,orig = orig_formula) { 
    new <- reformulate(c(x,'.')) 
    update(orig, new)}) 
## it works!  
new_formula[[1]] 
## MonExp_EGM ~ TE_1 + Month2 + Month3 + Month4 + Month5 + Month6 + 
## Month7 + Month8 + Month9 + Month10 + Month11 + Month12 + 
## Yrs_minus_2004 + as.factor(LGA) 
new_formula[[2]] 
## MonExp_EGM ~ TE_2 + Month2 + Month3 + Month4 + Month5 + Month6 + 
## Month7 + Month8 + Month9 + Month10 + Month11 + Month12 + 
## Yrs_minus_2004 + as.factor(LGA) 


models <- lapply(new_formula, lm, data = pokies) 

Có bây giờ sẽ là 96 yếu tố trong danh sách models

Bạn có thể đặt tên cho chúng để phản ánh các tên được hoạch định ban đầu của bạn

names(models) <- paste0('z.out', 1:96) 
# or if you don't have a current version of R 
# names(models) <-paste('z.out', 1:96 ,sep = '') 

và sau đó truy cập vào một mô hình duy nhất bằng cách

models$z.out5 

vv

hoặc tạo bản tóm tắt của tất cả các mô hình

summaries <- lapply(models, summary) 

vv ....

# just the coefficients 
coefficients <- lapply(models, coef) 

# the table with coefficient estimates and standard.errors 

coef_tables <- apply(summaries, '[[', 'coefficients')