请选择 进入手机版 | 继续访问电脑版
楼主: Nothing_LF
2595 4

[有偿编程] 加权 非线性混合效应模型 [推广有奖]

  • 0关注
  • 0粉丝

本科生

64%

还不是VIP/贵宾

-

威望
0
论坛币
6 个
通用积分
12.0005
学术水平
1 点
热心指数
0 点
信用等级
0 点
经验
539 点
帖子
15
精华
0
在线时间
198 小时
注册时间
2015-4-4
最后登录
2023-9-23

Nothing_LF 学生认证  发表于 2019-4-23 11:47:06 |显示全部楼层 |坛友微信交流群
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
如题,利用nlme()语句进行非线性混合效应模型的加权回归,权函数是固定的,如下:\[g(x)=D^1.2\]
具体程序如下:
  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) :\n  longer object length is not a multiple of shorter object length",  : argument 2 (type 'list') cannot be handled by 'cat'


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


希望大家帮忙解决,有偿也可。部分数据在附件
data-wls.xlsx (12.73 KB)




二维码

扫码加我 拉你入群

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

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

关键词:非线性混合效应模型 线性混合效应模型 混合效应模型 线性混合效应 非线性混合 加权 nlme weights nonlinear mixed-effects model

Nothing_LF 学生认证  发表于 2019-4-23 11:49:04 |显示全部楼层 |坛友微信交流群
权函数错了,是g(x)=D^1.2

使用道具

Nothing_LF 学生认证  发表于 2019-4-24 17:53:12 |显示全部楼层 |坛友微信交流群
有懂的请联系我企鹅:2448687875,有偿请教

使用道具

花花花ya 发表于 2021-10-5 16:17:03 |显示全部楼层 |坛友微信交流群
你好 请问您的问题解决了吗?我也遇到了同样的问题。

使用道具

Nothing_LF 学生认证  发表于 2021-10-24 11:22:00 |显示全部楼层 |坛友微信交流群
花花花ya 发表于 2021-10-5 16:17
你好 请问您的问题解决了吗?我也遇到了同样的问题。
改改权函数形式就行

使用道具

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

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

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

GMT+8, 2024-3-29 04:00