楼主: xingzhaoh
1272 3

[问答] R做非线性分析时可以直接除以权函数吗 [推广有奖]

副教授

21%

还不是VIP/贵宾

-

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

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
用R作非线性分析时:
nls.sol<-nls(TOTAL~a*DBH^b,data=rd,start=list(a=0.01,b=2))

现在知道他们的权重了,可以直接除以权重吗?

nls.sol<-nls(TOTAL*1/DBH^1.86~a*DBH^b*1/DBH^1.86,data=rd,start=list(a=0.01,b=2))
二维码

扫码加我 拉你入群

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

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

关键词:非线性 权函数 nls SOL

沙发
waveland 发表于 2013-1-22 16:54:36 |只看作者 |坛友微信交流群
嗯,可以。。。

使用道具

藤椅
xingzhaoh 发表于 2013-1-22 18:50:10 |只看作者 |坛友微信交流群
waveland 发表于 2013-1-22 16:54
嗯,可以。。。
在使用library( systemfit )包建立非线性联立方程组时
程序如下:
library( systemfit )
rd<-read.csv("5.csv")
LEAF.formula<- LEAF~ 1/( 1+ r1*DBH^r2*H^r3*CW^r4)* CWB
BRANCH.formula<-BRANCH~1/( 1+(1/ r1*DBH^-r2*H^-r3*CW^-r4))* CWB
labels <- list("LEAF.formula", "BRANCH.formula")
inst <- ~ DBH+H+CW+LEAF+BRANCH
start.values<-c(r1=0.4,r2=0.9, r3=-0.2,r4=-1.0)
model<-list(LEAF.formula, BRANCH.formula)
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )
我用权重后精度怎么降低了
相关系数成负数了,请问怎么回事

library( systemfit )
rd<-read.csv("4444.csv")
LEAF.formula<- LEAF*1/DBH^0.23~ 1/( 1+ r1*DBH^r2*H^r3*CW^r4)* CWB*1/DBH^0.23
BRANCH.formula<-BRANCH*1/DBH^0.96~1/( 1+(1/ r1*DBH^-r2*H^-r3*CW^-r4))* CWB*1/DBH^0.96
labels <- list("LEAF.formula", "BRANCH.formula")
inst <- ~ DBH+H+CW+LEAF+BRANCH
start.values<-c(r1=0.4,r2=0.9, r3=-0.2,r4=-1.0)
model<-list(LEAF.formula, BRANCH.formula)
model.ols <- nlsystemfit( "OLS", model, start.values, data=rd, eqnlabels=labels )

使用道具

板凳
xingzhaoh 发表于 2013-1-22 19:17:12 |只看作者 |坛友微信交流群
xingzhaoh 发表于 2013-1-22 18:50
在使用library( systemfit )包建立非线性联立方程组时
程序如下:
library( systemfit )
https://bbs.pinggu.org/thread-2177875-1-1.html
指点一下这个问题可以吗?

使用道具

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

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-17 21:24