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

データ名:非線形最小2乗法-ロジスティック関数テストデータ (nlra_04.dns 2023-05-15)

用量項目:相対経過日数
反応項目:累積感染者数

    用量          x=用量      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       0               0         1          20           0           0
       1               1         1          61           0           0
       2               2         1          64           0           0
       3               3         1          64           0           0
       4               4         1          70           0           0
       5               5         1         135           0           0
       6               6         1         175           0           0
       7               7         1         219           0           0
       8               8         1         218           0           0
       9               9         1         285           0           0
      10              10         1         355           0           0
      11              11         1         454           0           0
      12              12         1         542           0           0
      13              13         1         542           0           0
      14              14         1         634           0           0
      15              15         1         634           0           0
      16              16         1         634           0           0
      17              17         1         691           0           0
      18              18         1         691           0           0
      19              19         1         691           0           0
      20              20         1         705           0           0
      21              21         1         705           0           0
      22              22         1         705           0           0
      23              23         1         705           0           0
      24              24         1         706           0           0
      25              25         1         706           0           0
----------------------------------------------------------------------
    全体            12.5        26     438.885     266.654     52.2952

反応ロジット
    用量          x=用量      例数      平均値    標準偏差    標準誤差
----------------------------------------------------------------------
       0               0         1    -6.90675           0           0
       1               1         1    -3.08424           0           0
       2               2         1    -2.98376           0           0
       3               3         1    -2.98376           0           0
       4               4         1    -2.80656           0           0
       5               5         1    -1.71618           0           0
       6               6         1    -1.31851           0           0
       7               7         1   -0.966188           0           0
       8               8         1   -0.973551           0           0
       9               9         1     -0.5215           0           0
      10              10         1  -0.0991504           0           0
      11              11         1    0.491772           0           0
      12              12         1     1.09686           0           0
      13              13         1     1.09686           0           0
      14              14         1     2.03761           0           0
      15              15         1     2.03761           0           0
      16              16         1     2.03761           0           0
      17              17         1     3.42059           0           0
      18              18         1     3.42059           0           0
      19              19         1     3.42059           0           0
      20              20         1     4.48855           0           0
      21              21         1     4.48855           0           0
      22              22         1     4.48855           0           0
      23              23         1     4.48855           0           0
      24              24         1     4.63168           0           0
      25              25         1     4.63168           0           0
----------------------------------------------------------------------
    全体            12.5        26    0.842981     3.13029      0.6139

・用量反応直線:y=22.1966+33.335x
  寄与率:r2乗=0.914241    回帰係数の標準誤差=2.08403

・反応(y)=438.885の時の用量:x=12.5
  95%信頼区間:下限=7.42915  上限=17.5708

                  用量反応直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰     1.62517e+06         1   1.62517e+06     255.855   2.65811e-14***
ズレ          152446        24       6351.92
--------------------------------------------------------------------------
用量     1.77761e+06        25

・用量反応ロジット直線:ly=ly0+c・x=-4.13374+0.398138x
  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=31.1808  最大値:ym+M=31.1808+681.391=712.572
  パラメータ初期値:y0=41.3754  ly0=-4.18721  c=0.402415
  パラメータ推定値:y0=41.9265  ly0=-4.13374  c=0.398138
  y=ym+M/2=371.876の時の用量:D50=x=-ly0/c=10.3827
  重寄与率:R2乗=0.946356  重相関係数:R=0.972808

・反応:y=371.876(ly=0)の時の用量:D=x=10.3827
  95%信頼区間:下限=6.42655  上限=14.2958

              用量反応ロジット直線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰         231.826         1       231.826      423.39   9.38919e-17***
ズレ         13.1411        24      0.547548
--------------------------------------------------------------------------
用量         244.967        25

・用量反応曲線: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法  反復回数=25
  最小値(固定):ym=31.1808  最大値(固定):ym+M=31.1808+681.391=712.572
  パラメーター初期値:y0=41.9265  ly0=-4.13374  c=0.398138
  パラメーター収束値:y0=43.1608  ly0=-4.02316  c=0.406093
  y=ym+M/2=371.876の時の用量:D50=x=-ly0/c=9.90701
  残差平方和(Q)=8536.14  赤池の情報量基準(AIC)=239.354
  重寄与率:R2乗=0.995494  重相関係数:R=0.997744

・反応:y=371.876の時の用量:D=x=9.90701
  95%信頼区間:下限=6.06919  上限=13.7448

                  用量反応曲線の分散分析表(ANOVA table)
要因          平方和    自由度    平均平方和         F値   有意確率p値
--------------------------------------------------------------------------
回帰     1.88582e+06         1   1.88582e+06     5302.12             0***
ズレ         8536.14        24       355.673
--------------------------------------------------------------------------
用量     1.89436e+06        25


