楼主: xingzhaoh
3944 5

[程序分享] 如何用R求方程 [推广有奖]

副教授

20%

还不是VIP/贵宾

-

威望
0
论坛币
3407 个
通用积分
13.9638
学术水平
5 点
热心指数
6 点
信用等级
3 点
经验
583 点
帖子
421
精华
0
在线时间
703 小时
注册时间
2011-9-14
最后登录
2023-11-17

楼主
xingzhaoh 发表于 2012-4-23 20:19:31 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
请高手指点如何用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
下载链接: https://bbs.pinggu.org/a-1097794.html

24.5 KB

请大侠指教

沙发
qoiqpwqr 发表于 2012-4-23 20:28:39
没看到里面有f5和g5

藤椅
xingzhaoh 发表于 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)之间的关系

麻烦你了



板凳
qoiqpwqr 发表于 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)

报纸
xingzhaoh 发表于 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)之间的关系
呵呵,这下改好了

地板
xingzhaoh 发表于 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

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-7 13:10