在使用library( systemfit )包建立非线性联立方程组时
程序如下：
library( systemfit )
rd<-read.csv("2.csv")
ABOVE.formula<-ABOVE~1/(1+r1*DBH^r2*H^r3)* TOTAL
UNDER.formula<-UNDER~1/(1+(1/r1*DBH^-r2* H^-r3))*TOTAL
labels <- list( "ABOVE.formula","UNDER.formula")
inst <- ~ DBH+H+ABOVE+ UNDER+TOTAL
start.values<-c(r1=0.1,r2=0.3, r3=-0.5)
model<-list(ABOVE.formula,UNDER.formula)
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
结果的残差图存在方差异性
想通过权重取出异方差，假如ABOVE.formula公式的权重是：1/DBH^1.93;UNDER.formula的为：1/DBH^1.49
那要怎么计算呢？可以把程序发给我吗？非常感谢数据2在附件中