qoiqpwqr 发表于 2012-4-24 08:02 
怎么有两个f1(x)和g1(x)
能不能帮我看看这个问题的原因
library( systemfit )
rd<-read.
delim("11.txt")
STEM.formula<-STEM~ a1*DBH^b1/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3*H^c3*CW^d3 + a4*DBH^b4*H^c4*CW^d4)*( a5*DBH^b5)
SKIN.formula<-SKIN~ a2*DBH^b2/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3*H^c3*CW^d3 + a4*DBH^b4*H^c4*CW^d4)*( a5*DBH^b5)
LEAF.formula<-LEAF~ a3*DBH^b3*H^c3*CW^d3/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3*H^c3*CW^d3 + a4*DBH^b4*H^c4*CW^d4)*( a5*DBH^b5)
BRANCH.formula<-BRANCH~ a4*DBH^b4*H^c4*CW^d4/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3*H^c3*CW^d3 + a4*DBH^b4*H^c4*CW^d4)*( a5*DBH^b5)
ABOVE.formula<- ABOVE~ a5*DBH^b5/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3*H^c3*CW^d3 + a4*DBH^b4*H^c4*CW^d4)*( a5*DBH^b5)
labels<-list("STEM.formula","SKIN.formula",“LEAF.formula","BRANCH.formula","ABOVE.formula")
start.values <- c(a1=0.06,b1=2.5,a2=0.03,b2=2.0,a3=0.4,b3=1.5,c3=-0.9,d3=1.0,
a4=0.2,b4=2.8,c4=-1.8,d4=1.2,a5=0.1,b5=2.4)
inst <- ~ DBH+STEM+SKIN+LEAF+BRANCH+ABOVE
model<-list(STEM.formula, SKIN.formula, LEAF.formula,BRANCH.formula, ABOVE.formula)
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
出现错误于eqnlabels[
] : 类别为'closure'的对象不可以取子集
请您帮我分析一下,这个一直没做出来呀
这个比较重要,请您帮我分析一下原因