=== 用量反応解析 ===                                               [DANS V7.4]

データ名:用量反応解析例数不揃いの時 (doser_05.dns 2023-05-13)

用量項目:用量
反応項目:反応

    用量          x=用量      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       0               0         5       36.66     2.93138     1.31095
     0.1             0.1         6     38.1167     2.57714     1.05211
       1               1         8     39.3375      2.2058    0.779867
      10              10        10        41.5     2.34284     0.74087
     100             100         8        43.6     2.33972    0.827216
    1000            1000         6       44.25     2.48495     1.01448
----------------------------------------------------------------------
    全体         160.665        43     40.8372     3.46829     0.52891

反応ロジット
    用量          x=用量      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       0               0         5    -6.90675           0           0
     0.1             0.1         6    -5.34182     2.57178     1.04993
       1               1         8    -4.55373     3.31423     1.17176
      10              10        10   -0.850982      5.1153      1.6176
     100             100         8     3.19654     4.07255     1.43986
    1000            1000         6     3.78973     3.69132     1.50698
----------------------------------------------------------------------
    全体         160.665        43    -1.47008     5.31738    0.810893

・用量反応直線:y=40.1023+0.00457442x
  寄与率:r2乗=0.365002    回帰係数の標準誤差=0.00109425

・反応(y)=40.8372の時の用量:x=160.665
  95%信頼区間:下限=-1089.09  上限=1410.42

                  用量反応直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         104.021         1       104.021     17.4758    0.00017086***
ズレ         180.966         4       45.2414      7.6007   0.000141825***
--------------------------------------------------------------------------
用量         284.986         5       56.9973     9.57572   6.32614e-06***
残差         220.234        37       5.95227
--------------------------------------------------------------------------
全体          505.22        42

・用量反応ロジット直線:ly=ly0+c・x=-2.61187+0.00710659x
  logit(y)=ly=ly0+c・x=ln[(y0-ym)/{M-(y0-ym)}]+c・x=ln[(y-ym)/{M-(y-ym)}]
  反応:y=ym+M/{1+exp(-ly)}  範囲外データ近似数=26
  最小値:ym=40.419  最大値:ym+M=40.419+2.95265=43.3717
  パラメータ初期値:y0=43.3717  ly0=0  c=-2.87034
  パラメータ推定値:y0=40.6209  ly0=-2.61187  c=0.00710659
  y=ym+M/2=41.8954の時の用量:D50=x=-ly0/c=367.527
  重寄与率:R2乗=0.381633  重相関係数:R=0.617764

・反応:y=41.8954(ly=0)の時の用量:D=x=367.527
  95%信頼区間:下限=-822.981  上限=1684.5

              用量反応ロジット直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         251.055         1       251.055     17.5369   0.000167179***
ズレ          406.79         4       101.698     7.10386   0.000238577***
--------------------------------------------------------------------------
用量         657.846         5       131.569     9.19047   9.43699e-06***
残差         529.685        37       14.3158
--------------------------------------------------------------------------
全体         1187.53        42

・用量反応曲線:y=ym+M/[1+[{M-(y0-ym)}/(y0-ym)]exp(-c・x)]
                =ym+M/{1+exp(-ly0-c・x)}=ym+M/{1+exp(-ly)}
  logit(y)=ly=ly0+c・x=ln[(y0-ym)/{M-(y0-ym)}]+c・x=ln[(y-ym)/{M-(y-ym)}]
  計算方法:滑降Simplex+Newton法  反復回数=23
  最小値(固定):ym=40.419  最大値(固定):ym+M=40.419+2.95265=43.3717
  パラメーター初期値:y0=40.6209  ly0=-2.61187  c=0.00710659
  パラメーター収束値:y0=40.6209  ly0=-2.61187  c=0.00710659
  ※エラー発生！


