全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3912 5
2012-04-23
请高手指点如何用R求解附件中的方程

求解出f1(x),f2(x),f1(x),f3(x),f4(x)

          g1(x),g1(x),g2(x),g3(x),g4(x)

就是f1(x),f2(x),f1(x),f3(x),f4(x)和 g1(x),g1(x),g2(x),g3(x),g4(x)之间的关系



请大侠指点


附件列表

求解方程.doc

大小:24.5 KB

 马上下载

请大侠指教

二维码

扫码加我 拉你入群

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

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

全部回复
2012-4-23 20:28:39
没看到里面有f5和g5
二维码

扫码加我 拉你入群

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

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

2012-4-24 08:01:32
qoiqpwqr 发表于 2012-4-23 20:28
没看到里面有f5和g5
是我多写了,求f1(x),f2(x),f1(x),f3(x),f4(x)

          g1(x),g2(x),g3(x),g4(x)之间的关系

麻烦你了



二维码

扫码加我 拉你入群

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

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

2012-4-24 08:02:50
xingzhaoh 发表于 2012-4-24 08:01
是我多写了,求f1(x),f2(x),f1(x),f3(x),f4(x)

          g1(x),g1(x),g2(x),g3(x),g4(x)之间的关系 ...
怎么有两个f1(x)和g1(x)
二维码

扫码加我 拉你入群

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

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

2012-4-24 08:38:42
xingzhaoh 发表于 2012-4-24 08:01
是我多写了,求f1(x),f2(x),f1(x),f3(x),f4(x)

          g1(x),g2(x),g3(x),g4(x)之间的关系
呵呵,这下改好了
二维码

扫码加我 拉你入群

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

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

2012-4-24 08:46:42
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'的对象不可以取子集

请您帮我分析一下,这个一直没做出来呀


这个比较重要,请您帮我分析一下原因

附件列表

11.txt

大小:5.75 KB

 马上下载

二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群