R程序解决非线性混合模型

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

想用R软件拟合一个非线性混合效应模型
y=a+b*x1+c*log10(x2+d)+e*x3+f*x4
其中a,b,c,d,e,f,是拟合系数，拟合过程中，c为随机效应，其余参数为固定效应
数据结构是这样的（由于数据很长，只是写了几个作为说明）
y       x1        x2     x3    x4
0.1    4           0      0       1
02.    5         10      0       0
2.1    8          3      1        0
5       6       100     1        0
2      6.8      50      0        1
6       7        30      0        0
3      7.3      80       1       0


数据文件为1.txt   拟合的初始参数为（-1，0.6，-2，10，-0.2，0.3）——分别对应a,b,c,d,e,f
求高手详细讲解怎么拟合

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：混合模型 非线性 R程序 混合效应模型 随机效应 程序 混合

求助啊

我有几点不是很清楚：
1） y=a+b*x1+c*log10(x2+d)+e*x3+f*x4，中的c*log10(x2+d)究竟是什么意思。LOG10是以10为底？还是其他什么意思？
2）"c为随机效应"是什么意思？

TaskShare 发表于 2012-12-23 22:22
我有几点不是很清楚：
1） y=a+b*x1+c*log10(x2+d)+e*x3+f*x4，中的c*log10(x2+d)究竟是什么意思。LOG10是 ...

1）是的，c*log10(x2+d)  其中，c，d是要拟合的系数，log10表示以10为底，
2）y=a+b*x1+c*log10(x2+d)+e*x3+f*x4
     这里我说c为随机效应是我表述有错误，我对混合效应模型的理解有问题
    数据是这样的，例如，在A点发生了一个爆破，其大小用为x1表示；该爆破在空间各处都有一个影响，y表示在x2点的爆破影响大小，x3，x4的有三种组合，（0，0）代表该爆破效果较差，（0，1）代表一般，（1，0）代表完美
     现在我想区分残差项,将其表示为2项，事件间的和事件内的---数学表示如下：      yij={a+b*xi+c*log10(xij+d)+e*xi1+f*xi2}+ui+εij   ui表示不同爆破事件间的残差，εij  表示单一爆破事件内的残差
     希望我能把问题表述明白了，我对这块理解很浅，请您给予指点。要是我哪里没有表示清楚，希望您能尽快联系我。若可以我当面向您请教

TaskShare 发表于 2012-12-23 22:22
我有几点不是很清楚：
1） y=a+b*x1+c*log10(x2+d)+e*x3+f*x4，中的c*log10(x2+d)究竟是什么意思。LOG10是 ...

还有一点没说到，不同的爆破有不同的x1，对某一个爆破（给定一个x1）会有多个x2（该爆破在空间多点都有影响）,每一个爆破有对应的x3,x4； （0，0）代表该爆破效果较差，（0，1）代表一般，（1，0）代表完美，     数据形式如下 ：
y        x1        x2        x3        x4
3        4          0.5        0        0
4        4          1           0        0
5        4          2           0        0
6        4          3           0        0
…        …        …          …        …
5        5         0           1        0
6        5        1            1        0
7        5        1.5         1        0
8        5        2.5         1        0
…        …        …        …        …
4        6        2            0        1
5.2     6        3            0        1
5.3     6        3.2         0        1
6        6        3.5         0        1
7        6        8            0        1
...        …        …        …        ...

哪位大神指点下

请问楼主知道怎么做了么？我也遇到同样的问题？

请问楼主拟合非线性混合模型的R包要怎么才能找到？