刚接触哑变量线性模型
数据在附件中：
code如下：
rd<-read.csv("11.csv")
attach(rd)
model1<-lm(H~(a1*f1+a2*f2+a3*f3+a4*f4+a5*f5+a6*f6+a7*f7+a8*f8+a9*f9+a10*f10+a11*f11+a12*f12+a13*f13+a14*f14+a15*f15+a16*f16+a17*f17+a18*f18+a19*f19+a20*f20+a21*f21+a22*f22+a23*f23+a24*f24
+ +a25*f25+a26*f26+a27*f27+a28*f28+a29*f29+a30*f30+a31*f31+a32*f32+a33*f33+a34*f34+a35*f35+a36*f36+a37*f37+a38*f38+a39*f39+a40*f40+a41*f41+a42*f42+a43*f43+a44*f44+a45*f45+a46*f46+a47*f47+a48*f48
+ +a49*f49+a50*f50+a51*f51+a52*f52+a53*f53+a54*f54+a55*f55+a56*f56+a57*f57+a58*f58+a59*f59+a60*f60+a61*f61+a62*f62+a63*f63+a64*f64+a65*f65+a66*f66+a67*f67+a68*f68+a69*f69+a70*f70+a71*f71+a72*f72
+ +a73*f73+a74*f74+a75*f75+a76*f76+a77*f77+a78*f78+a79*f79+a80*f80+a81*f81+a82*f82+a83*f83+a84*f84+a85*f85+a86*f86+a87*f87+a88*f88+a89*f89+a90*f90+a91*f91+a92*f92+a93*f93+a94*f94+a95*f95+a96*f96
+ +a97*f97+a98*f98+a99*f99+a100*f100+a101*f101+a102*f102+a103*f103+a104*f104+a105*f105+a106*f106+a107*f107+a108*f108+a109*f109+a110*f110+a111*f111+a112*f112+a113*f113+a114*f114+a115*f115+a116*f116
+ +a117*f117+a118*f118+a119*f119+a120*f120+a121*f121+a122*f122+a123*f123+a124*f124+a125*f125+a126*f126+a127*f127+a128*f128+a129*f129+a130*f130+a131*f131+a132*f132+a133*f133+a134*f134+a135*f135
+ +a136*f136+a137*f137+a138*f138+a139*f139)+(b1*f1+b2*f2+b3*f3+b4*f4+b5*f5+b6*f6+b7*f7+b8*f8+b9*f9+b10*f10+b11*f11+b12*f12+b13*f13+b14*f14+b15*f15+b16*f16+b17*f17+b18*f18+b19*f19+b20*f20+b21*f21
+ +b22*f22+b23*f23+b24*f24+b25*f25+b26*f26+b27*f27+b28*f28+b29*f29+b30*f30+b31*f31+b32*f32+b33*f33+b34*f34+b35*f35+b36*f36+b37*f37+b38*f38+b39*f39+b40*f40+b41*f41+b42*f42+b43*f43+b44*f44+b45*f45
+ +b46*f46+b47*f47+b48*f48+b49*f49+b50*f50+b51*f51+b52*f52+b53*f53+b54*f54+b55*f55+b56*f56+b57*f57+b58*f58+b59*f59+b60*f60+b61*f61+b62*f62+b63*f63+b64*f64+b65*f65+b66*f66+b67*f67+b68*f68+b69*f69
+ +b70*f70+b71*f71+b72*f72+b73*f73+b74*f74+b75*f75+b76*f76+b77*f77+b78*f78+b79*f79+b80*f80+b81*f81+b82*f82+b83*f83+b84*f84+b85*f85+b86*f86+b87*f87+b88*f88+b89*f89+b90*f90+b91*f91+b92*f92+b93*f93
+ +b94*f94+b95*f95+b96*f96+b97*f97+b98*f98+b99*f99+b100*f100+b101*f101+b102*f102+b103*f103+b104*f104+b105*f105+b106*f106+b107*f107+b108*f108+b109*f109+b110*f110+b111*f111+b112*f112+b113*f113
+ +b114*f114+b115*f115+b116*f116+b117*f117+b118*f118+b119*f119+b120*f120+b121*f121+b122*f122+b123*f123+b124*f124+b125*f125+b126*f126+b127*f127+b128*f128+b129*f129+b130*f130+b131*f131+b132*f132
+ +b133*f133+b134*f134+b135*f135+b136*f136+b137*f137+b138*f138+b139*f139)*A, data=rd)

错误于eval(expr, envir, enclos) : 找不到对象'a1'

请问怎么设置变量a1-a139和b1-b139，才能求出
a1-a139和b1-b139

不胜感激