向各位大神求解: 模型采用的公式在图片里(不知道怎么直接显示在界面)
现有数据两组数据如下:
时间t
|
PCI值
|
交通量(N)
|
板厚H
|
0
|
100
|
130894
|
37
|
2
|
94
|
130894
|
37
|
3
|
90
|
130894
|
37
|
4
|
85
|
130894
|
37
|
6
|
80
|
130894
|
37
|
0
|
100
|
185823
|
40
|
3
|
97
|
185823
|
40
|
4
|
95
|
185823
|
40
|
5
|
90
|
185823
|
40
|
我用的代码是 data roads; inputroad year y H N;
if (road=1) then do; M=0; end; else do; M=1; end; datalines;
1 0.01 100 37 130894 1 2 94 37130894 1 3 90 37130894 1 4 85 37130894 1 6 80 37130894 2 0.01 100 40 185823 2 3 97 40185823 2 4 95 40185823 2 5 90 40185823 ; run; proc nlmixed data=roads; parms a0=1 b0=0.1 a1=-0.1 a2 =0.5 u0=-1.8 u1=0.1 u2=0.7 s2a3=0.0001 s2b3=0.01 s2e=4; bounds u1>0, b0>0.07,s2a3>0.00001,s2b3>0.008; A = a0+u0*M+(a1+u1*M)*log10(N)+(a2+u2*M)*log10(H)+a3; B = b0+b3; AI = 10**(A); BI = 10**(B);
ex = 10**(-(AI/year)**BI); den = 100*(1-ex); model y ~normal(den,s2e); random a3b3 ~ normal([0,0],[s2a3,0,s2b3]) subject=road; run;
我用的算法算出来的结果好像有问题,如何根据bi的方差估计值将每一组的衰减规律找出来。
|