搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  data-wls.xlsx
资料下载链接地址: https://bbs.pinggu.org/a-2797015.html
附件大小:
12.73 KB   举报本内容
如题,利用nlme()语句进行非线性混合效应模型的加权回归,权函数是固定的,如下:[LaTex]g(x)=D^1.2[/LaTex]
具体程序如下:
  1. m1 <- nlme(model = WL~ exp(a)*D^b
  2. ,data = d1
  3. ,fixed = a+b~1
  4. ,random = a+b~1|Pn
  5. ,start = c(-2,1.5)
  6. ,control = nlmeControl(maxIter = 200,msMaxIter = 200))

  7. m3 <- update(m1,weights=varPower(form = ~D,fixed = 1.2))
复制代码


错误:
> warnings(m3)
警告信息:
1: In conLin$Xy * varWeights(object) :
longer object length is not a multiple of shorter object length Error in cat("1: In conLin$Xy * varWeights(object) :\nlonger object length is not a multiple of shorter object length",: argument 2 (type 'list') cannot be handled by 'cat'


如果用默认的varPower()选项,form = ~fitted(.),程序可以正常运行,但是这个fitted(.)在介绍中只说明是拟合值,并未说明是怎么来的,现在要输出也无法输出这个值


希望大家帮忙解决,有偿也可。部分数据在附件







    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2026-1-9 12:20