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

データ名:医学統計ハンドブック表4.3 (doser_02.dns 2017-11-21)

用量項目:用量(mg/kg)
反応項目:薬剤投与後の相対肝重量

    用量     x=log(用量)      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       0  (計算から除外)         6     33.3667     5.47272     2.23423
    12.5         1.09691         6     40.8167     10.2384     4.17983
      50         1.69897         6     38.8667     2.35683    0.962173
     200         2.30103         6     44.1333     8.33323     3.40203
----------------------------------------------------------------------
    全体         1.69897        18     41.2722      7.6088     1.79341

反応ロジット
    用量     x=log(用量)      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
    12.5         1.09691         6   -0.651403      3.7641     1.53669
      50         1.69897         6    -1.50609     0.65461    0.267244
     200         2.30103         6   -0.693956     1.60737    0.656207
----------------------------------------------------------------------
    全体         1.69897        18   -0.950482     2.28405    0.538355

・対数用量反応直線:y=36.5925+2.75443x
  寄与率:r2乗=0.387875    回帰係数の標準誤差=3.7122

・反応(y)=41.2722の時の用量:D=50(x=log(D)=1.69897)
  95%信頼区間:下限=3.49721e-05(x=-4.45628)  上限=7.14855e+07(x=7.85422)

                  用量反応直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         33.0008         1       33.0008    0.550555      0.469551
ズレ         52.0803         1       52.0803    0.868859      0.366033
--------------------------------------------------------------------------
用量         85.0811         2       42.5406    0.709707      0.507577
残差         899.115        15        59.941
--------------------------------------------------------------------------
全体         984.196        17

・対数用量反応ロジット直線:ly=ly0+c・x=-0.890441-0.0353393x
  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=33.3667  最大値:ym+M=33.3667+27.9333=61.3
  パラメータ初期値:y0=17.25  ly0=0  c=0.500237
  パラメータ推定値:y0=41.4958  ly0=-0.890441  c=-0.0353393
  y=ym+M/2=47.3333の時の用量:D50=6.35465e-26(x=log(D50)=-ly0/c=-25.1969)
  重寄与率:R2乗=0.00195143  重相関係数:R=0.044175

・反応:y=47.3333(ly=-5.55112e-16)の時の用量:D=6.35465e-26(x=log(D)=-25.1969)
  95%信頼区間:下限=0(x=-1892.46)  上限=inf(x=1842.06)

              用量反応ロジット直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰      0.00543222         1    0.00543222  0.00094855      0.975836
ズレ         2.77828         1       2.77828    0.485131       0.49676
--------------------------------------------------------------------------
用量         2.78371         2       1.39185     0.24304      0.787269
残差         85.9029        15       5.72686
--------------------------------------------------------------------------
全体         88.6867        17

・対数用量反応曲線: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法  反復回数=16
  最小値(固定):ym=33.3667  最大値(固定):ym+M=33.3667+27.9333=61.3
  パラメーター初期値:y0=41.4958  ly0=-0.890441  c=-0.0353393
  パラメーター収束値:y0=34.0746  ly0=-3.64949  c=-6.24478
  y=ym+M/2=47.3333の時の用量:D50=0.260371(x=log(D50)=-ly0/c=-0.584407)
  残差平方和(Q)=2109.09  赤池の情報量基準(AIC)=141.772
  重寄与率:R2乗=0  重相関係数:R=0

・反応:y=47.3333の時の用量:D=0.260371(x=log(D)=-0.584407)
  95%信頼区間:下限=0.0375228(x=-1.4257)  上限=1.80672(x=0.256891)

                  用量反応曲線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰               0         1             0           0             1
ズレ         746.348         1       746.348     12.4514    0.00303948**
--------------------------------------------------------------------------
用量         746.348         2       373.174     6.22569     0.0107511*
残差         899.115        15        59.941
--------------------------------------------------------------------------
全体         1645.46        17


