有一组不同文献中的自变量x和依变量y数据,分别用线性模型和非线性模型建立模型。线性模型为:Y=a+bX+c(X*X)+d(X*X*X);非线性模型为:Y=a*(1-exp(-exp(b)*(X^c)),其中的a、b、c和d为系数,想求试验数据的固定效应和随机效应。
其中线性模型Y=a+bX+c(X*X)+d(X*X*X)的分析固定效应和随机效应的SAS程序为:
PROC MIXED;
CLASS Study;
MODEL Y=X X*X X*X*X/Solution OUTP=Predictionset OUTPM=PredY;
RANDOM intercept X X*X X*X*X/TYPE=UN SUBJECT=Study;
RUN;
但非线性模型:Y=a*(1-exp(-exp(b)*(X^c)),不知用什么代码得出固定效应和随机效应,非常期待大家的帮助,谢谢。
数据大概如下:
| Study | X | Y |
1 | 80.15 | 839.00 |
2 | 30.54 | 558.10 |
2 | 51.28 | 859.20 |
2 | 69.88 | 790.50 |
2 | 89.88 | 775.70 |
3 | 68.63 | 1063.00 |
3 | 89.20 | 941.00 |
4 | 64.32 | 763.63 |
5 | 76.68 | 885.17 |
6 | 71.75 | 844.00 |
7 | 61.75 | 740.00 |
8 | 38.83 | 653.70 |
8 | 69.73 | 717.20 |
9 | 79.57 | 854.09 |
10 | 25.38 | 570.00 |
10 | 38.92 | 700.00 |
11 | 36.84 | 742.00 |
12 | 60.90 | 811.51 |
12 | 79.16 | 849.08 |
13 | 80.78 | 903.00 |
14 | 41.21 | 718.00 |
15 | 49.97 | 711.06 |
16 | 46.57 | 650.00 |
17 | 34.73 | 720.00 |
18 | 27.21 | 652.98 |
18 | 42.76 | 828.27 |
19 | 73.60 | 829.70 |
20 | 51.73 | 840.83 |
20 | 54.65 | 905.56 |