向各位大神求解：

模型采用的公式在图片里（不知道怎么直接显示在界面）

现有数据两组数据如下：

我用的代码是

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的方差估计值将每一组的衰减规律找出来。