联立方程组求解错误

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

library( systemfit )
rd<-read.csv("123.csv")
STEM.formula<-STEM~ a1*DBH^b1/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)* (a0*DBH^b0)
SKIN.formula<-SKIN~ a2*DBH^b2/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)* (a0*DBH^b0)
LEAF.formula<-LEAF~ a3*DBH^b3/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)* (a0*DBH^b0)
BRANCH.formula<-BRANCH~ a4*DBH^b4/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)* (a0*DBH^b0)
UNDER.formula<-UNDER~ a5*DBH^b5/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)* (a0*DBH^b0)
TOTAL.formula<-UNDER~ (a0*DBH^b0)^2/( a1*DBH^b1+ a2*DBH^b2+
a3*DBH^b3+ a4*DBH^b4+ a5*DBH^b5)
labels <- list( " STEM.formula ", " SKIN.formula ", " LEAF.formula ", " BRANCH.formula ",
" UNDER.formula ", " TOTAL.formula ")
inst <- ~ DBH+STEM+SKIN+LEAF+BRANCH+UNDERL+TOTAL
start.values <- c(a0=0.3,b0=2.3,a1=0.1,b1=2.3,a2=0.02,b2=2.0,a3=0.06,b3=1.4,
a4=0.003,b4=2.8,a5=0.02,b5=2.4)
model<-list(STEM.formula,SKIN.formula,LEAF.formula,BRANCH.formula,UNDER.formula, TOTAL.formula)
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
时出错警告信息：
1: In nlm(knls, startvals, typsize = abs(startvals), iterlim = maxiter, :
NA/Inf被换成最大的正值
2: In nlm(knls, startvals, typsize = abs(startvals), iterlim = maxiter, :
NA/Inf被换成最大的正值
错误于kronecker(X, Y) :
在为'kronecker'函数选择方法时评估'X'参数出了错: 错误于diag(diag(qr.solve(S, tol = solvtol))) :
在为'diag'函数选择方法时评估'x'参数出了错: 错误于diag(qr.solve(S, tol = solvtol)) :
在为'diag'函数选择方法时评估'x'参数出了错: 错误于qr.solve(S, tol = solvtol) : solve里的'a'是奇异矩阵：
请问是什么原因？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：联立方程组 联立方程 方程组 Library System library 方程组

数来指点一下呀

请高手指点一下

xingzhaoh 发表于 2012-4-18 08:21
请高手指点一下

123.csv 请传上来

epoh 发表于 2012-4-18 15:55
123.csv 请传上来

有道理。

epoh 发表于 2012-4-18 15:55
123.csv 请传上来

文件已上传，
还要请您用这个方程代码计算一下，这个对小弟来说比较重要：
谢谢
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)* (a0*DBH^b0)
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)* (a0*DBH^b0)
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)* (a0*DBH^b0)
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)* (a0*DBH^b0)
UNDER.formula<-UNDER~ 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)* (a0*DBH^b0)
TOTAL.formula<-TOTAL~(a0*DBH^b0)^2/( 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 ",
" UNDER.formula ",  " TOTAL.formula ")
inst <- ~ DBH +STEM+SKIN+LEAF++H+CW+ BRANCH+UNDERL+TOTAL
start.values <- c(a0=0.19,b0=2.5,a1=0.06,b1=2.5,a2=0.03,b2=2,a3=0.4,b3=1.5,c3=-0.9,d3=1,
a4=0.2,b4=2.8,c4=-1.8,d4=1.2,a5=0.01,b5=2.6)
model<-list(STEM.formula, SKIN.formula, LEAF.formula,BRANCH.formula,
UNDER.formula ,TOTAL.formula)            
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )

xingzhaoh 发表于 2012-4-19 08:08
文件已上传，
还要请您用这个方程代码计算一下，这个对小弟来说比较重要：
谢谢

多给几组start.values
这个差太远了,几乎没有iteration

给的越接近才有机会求解

只是初始值的问题吗

xingzhaoh 发表于 2012-4-19 13:31
只是初始值的问题吗

初步显现出来的是初始值问题
有出现其他问题再说
理论上这是你的模型
还有其他问题,你最清楚

请问初始值可以影响其相关系数，还是只要能求解出来其相关系数都一样