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

データ名:用量反応解析 (doser_04.dns 2011-5-24)

層別内容:用量=1-4

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

    用量     x=log(用量)      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       1               0        10       791.5     39.7555     12.5718
       2         0.30103        10       861.9     20.1243     6.36387
       4         0.60206        10         926     25.6775     8.11993
----------------------------------------------------------------------
    全体         0.30103        30       859.8     62.7861     11.4631

反応ロジット
    用量     x=log(用量)      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       1               0        10     1.40471    0.248084   0.0784511
       2         0.30103        10     1.92496    0.179394   0.0567294
       4         0.60206        10      3.0272     1.38291    0.437313
----------------------------------------------------------------------
    全体         0.30103        30     2.11896     1.04687    0.191132

・対数用量反応直線:y=792.55+223.4x
  寄与率:r2乗=0.999269    回帰係数の標準誤差=22.0553

・反応(y)=859.8の時の用量:D=2(x=log(D)=0.30103)
  95%信頼区間:下限=1.04219(x=0.0179461)  上限=3.83808(x=0.584114)

                  用量反応直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         90451.2         1       90451.2     102.598   1.07754e-10***
ズレ           66.15         1         66.15   0.0750334      0.786229
--------------------------------------------------------------------------
用量         90517.4         2       45258.7     51.3366   6.30864e-10***
残差         23803.4        27       881.607
--------------------------------------------------------------------------
全体          114321        29

・対数用量反応ロジット直線:ly=ly0+c・x=1.30771+2.69491x
  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)}  範囲外データ近似数=1
  最小値:ym=0  最大値:ym+M=0+989=989
  パラメータ初期値:y0=353.5  ly0=-0.586529  c=2.81227
  パラメータ推定値:y0=778.471  ly0=1.30771  c=2.69491
  y=ym+M/2=494.5の時の用量:D50=0.327152(x=log(D50)=-ly0/c=-0.485251)
  重寄与率:R2乗=0.958875  重相関係数:R=0.979222

・反応:y=494.5(ly=0)の時の用量:D=0.327152(x=log(D)=-0.485251)
  95%信頼区間:下限=0.028197(x=-1.5498)  上限=1.41718(x=0.151426)

              用量反応ロジット直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         13.1625         1       13.1625     19.6831   0.000138312***
ズレ        0.564527         1      0.564527    0.844191      0.366337
--------------------------------------------------------------------------
用量          13.727         2       6.86351     10.2637    0.00048381***
残差         18.0554        27      0.668719
--------------------------------------------------------------------------
全体         31.7824        29

・対数用量反応曲線: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法  反復回数=29
  最小値(固定):ym=0  最大値(固定):ym+M=0+989=989
  パラメーター初期値:y0=778.471  ly0=1.30771  c=2.69491
  パラメーター収束値:y0=789.091  ly0=1.37302  c=2.02211
  y=ym+M/2=494.5の時の用量:D50=0.20941(x=log(D50)=-ly0/c=-0.679002)
  残差平方和(Q)=24752.7  赤池の情報量基準(AIC)=307.501
  重寄与率:R2乗=0.999768  重相関係数:R=0.999884

・反応:y=494.5の時の用量:D=0.20941(x=log(D)=-0.679002)
  95%信頼区間:下限=0.0181397(x=-1.74137)  上限=2.4175(x=0.383366)

                  用量反応曲線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰     4.09289e+06         1   4.09289e+06     4642.53             0***
ズレ         949.297         1       949.297     1.07678      0.308623
--------------------------------------------------------------------------
用量     4.09384e+06         2   2.04692e+06      2321.8             0***
残差         23803.4        27       881.607
--------------------------------------------------------------------------
全体     4.11764e+06        29


