我想要以FID编号来分组进行回归，也就是FID编号=1的数据做一次回归，FID编号=2的数据做一次回归，以此类推。并且希望能把每个回归分析的变量系数，每个系数的P值，以及R2放到一个文件里。
form<-LN.yield.~year+Tave+Tave2+I+I2+Rad
for (i in 1:2408) {
lrm<-lm(form, data=subset(data_all,data_all$FID编号==i))  
}
这样写报错了，求大神指导