全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
8309 7
2012-12-20
想用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
求高手详细讲解怎么拟合
二维码

扫码加我 拉你入群

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

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

全部回复
2012-12-20 15:43:01
求助啊
二维码

扫码加我 拉你入群

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

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

2012-12-23 22:22:19
我有几点不是很清楚:
1) y=a+b*x1+c*log10(x2+d)+e*x3+f*x4,中的c*log10(x2+d)究竟是什么意思。LOG10是以10为底?还是其他什么意思?
2)"c为随机效应"是什么意思?
二维码

扫码加我 拉你入群

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

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

2012-12-24 09:35:58
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  表示单一爆破事件内的残差
     希望我能把问题表述明白了,我对这块理解很浅,请您给予指点。要是我哪里没有表示清楚,希望您能尽快联系我。若可以我当面向您请教

二维码

扫码加我 拉你入群

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

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

2012-12-24 09:45:57
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
...        …        …        …        ...





二维码

扫码加我 拉你入群

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

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

2012-12-24 11:24:23
哪位大神指点下
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群