=== プロビット法による平行線検定 ===                               [DANS V7.4]

データ名:プロビット直線の平行線検定用テストデータ (parapl_01.tab 2006-3-12)

用量項目(縦):用量
反応項目(横):反応

・群 1:薬剤 (0:標準薬 1:被験薬)=0

用量＼反応    無(  %  )    有(  %  )  合計(  %  )
-------------------------------------------------
250           11( 73.3)     4( 26.7)    15(100.0)
300            8( 53.3)     7( 46.7)    15(100.0)
360            4( 26.7)    11( 73.3)    15(100.0)
-------------------------------------------------
合計          23( 51.1)    22( 48.9)    45(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
     250     2.398   0.266667(    4/15   )    0.266667     4.37707     4.39624
     300     2.477   0.466667(    7/15   )    0.466667     4.91635     4.96972
     360     2.556   0.733333(   11/15   )    0.733333     5.62293     5.54321
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=-14.5007+7.85996x
  反応率=50%(y=5)の時の用量:D50=302.708(log(D50)=2.48102)
共通対数用量-プロビット直線:y=-12.9712+7.24267x
  反応率=50%(y=5)の時の用量:D50=302.902(log(D50)=2.4813)

・群 2:薬剤 (0:標準薬 1:被験薬)=1

用量＼反応    無(  %  )    有(  %  )  合計(  %  )
-------------------------------------------------
300           12( 80.0)     3( 20.0)    15(100.0)
360           10( 66.7)     5( 33.3)    15(100.0)
440            6( 40.0)     9( 60.0)    15(100.0)
-------------------------------------------------
合計          28( 62.2)    17( 37.8)    45(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
     300     2.477        0.2(    3/15   )         0.2     4.15838       4.061
     360     2.556   0.333333(    5/15   )    0.333333     4.56927     4.63448
     440     2.643        0.6(    9/15   )         0.6     5.25335     5.26568
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=-12.3658+6.65277x
  反応率=50%(y=5)の時の用量:D50=407.672(log(D50)=2.61031)
共通対数用量-プロビット直線:y=-13.88+7.24267x
  反応率=50%(y=5)の時の用量:D50=404.362(log(D50)=2.60677)
※反復計算回数=3

                       平行線検定表(分散分析表)
要因                      平方和    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------------
群 1の直線性              6.5551         1        6.5551     0.0104584*
群 1の異質性(ズレ)     0.0417899         1     0.0417899      0.838021
群 2の直線性             4.91433         1       4.91433     0.0266348*
群 2の異質性(ズレ)     0.0718264         1     0.0718264      0.788696
--------------------------------------------------------------------------
共通直線性               11.3904         1       11.3904   0.000738261***
異質性(ズレ)合計        0.113616         2      0.113616      0.944775
非平行性               0.0790686         1     0.0790686      0.778564
--------------------------------------------------------------------------
全体                      11.583         4

○効力比Rとその95%信頼区間
・群2/群1:R=0.749086 (log(R)=-0.125468)
  信頼区間:下限=0.573083(log(R)=-0.241783)  上限=0.899065(log(R)=-0.0462091)

・修正切片差のBonferroni型多重比較(Bonferroni type multiple comparison)
群         - 群              χ2乗値  自由度   有意確率p値
-------------------------------------------------------------
2          - 1              7.25299       1    0.00707831**
-------------------------------------------------------------

・修正切片差のBonferroni型95%同時信頼区間(simultaneous confidence interval)
群      - 群        修正切片差    絶対精度        下限        上限
------------------------------------------------------------------
2       - 1          -0.908724    0.661335    -1.57006   -0.247389
------------------------------------------------------------------


=== プロビット法による平行線検定 ===                               [DANS V7.4]

データ名:生物統計ハンドブック1-8-IV (parapl_02.tab 2016-11-30)

用量項目(縦):群(用量)
反応項目(横):効果

・群 1:軽度

群(用量)＼効果  無効(  %  )  有効(  %  )  合計(  %  )
-----------------------------------------------------
0                 27( 67.5)    13( 32.5)    40(100.0)
10                25( 65.8)    13( 34.2)    38(100.0)
20                19( 70.4)     8( 29.6)    27(100.0)
30                 8( 42.1)    11( 57.9)    19(100.0)
-----------------------------------------------------
合計              79( 63.7)    45( 36.3)   124(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
       0         -      0.325(   13/40   )           0
      10         1   0.342105(   13/38   )   0.0253411     3.04584    -1.85172
      20     1.301   0.296296(    8/27   )           0          ↓     2.27201
      30     1.477   0.578947(   11/19   )    0.376218     4.68457     4.68424
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=-18.32+15.5739x
  反応率=50%(y=5)の時の用量:D50=31.4322(log(D50)=1.49737)
共通対数用量-プロビット直線:y=-15.5505+13.6987x
  反応率=50%(y=5)の時の用量:D50=31.6353(log(D50)=1.50017)

・群 2:中等度

群(用量)＼効果  無効(  %  )  有効(  %  )  合計(  %  )
-----------------------------------------------------
0                 24( 70.6)    10( 29.4)    34(100.0)
10                20( 64.5)    11( 35.5)    31(100.0)
20                27( 75.0)     9( 25.0)    36(100.0)
30                12( 46.2)    14( 53.8)    26(100.0)
-----------------------------------------------------
合計              83( 65.4)    44( 34.6)   127(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
       0         -   0.294118(   10/34   )           0
      10         1   0.354839(   11/31   )   0.0860215     3.63433    -1.93192
      20     1.301       0.25(    9/36   )           0          ↓     2.19181
      30     1.477   0.538462(   14/26   )    0.346154     4.60427     4.60404
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=-18.3336+15.5288x
  反応率=50%(y=5)の時の用量:D50=31.813(log(D50)=1.5026)
共通対数用量-プロビット直線:y=-15.6306+13.6987x
  反応率=50%(y=5)の時の用量:D50=32.0646(log(D50)=1.50603)

・群 3:重度

群(用量)＼効果  無効(  %  )  有効(  %  )  合計(  %  )
-----------------------------------------------------
0                 12( 75.0)     4( 25.0)    16(100.0)
10                15( 71.4)     6( 28.6)    21(100.0)
20                20( 74.1)     7( 25.9)    27(100.0)
30                25( 55.6)    20( 44.4)    45(100.0)
-----------------------------------------------------
合計              72( 66.1)    37( 33.9)   109(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
       0         -       0.25(    4/16   )           0
      10         1   0.285714(    6/21   )    0.047619     3.33161    -2.18127
      20     1.301   0.259259(    7/27   )   0.0123457      2.7538     1.94246
      30     1.477   0.444444(   20/45   )    0.259259     4.35437     4.35468
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=9.4246-3.43251x
  反応率=50%(y=5)の時の用量:D50=19.4549(log(D50)=1.28903)
共通対数用量-プロビット直線:y=-15.88+13.6987x
  反応率=50%(y=5)の時の用量:D50=33.4371(log(D50)=1.52423)
※反復計算回数=12

                       平行線検定表(分散分析表)
要因                      平方和    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------------
群 1の直線性            0.039171         1      0.039171       0.84311
群 1の異質性(ズレ)     0.0506823         1     0.0506823       0.82188
群 2の直線性           0.0384996         1     0.0384996      0.844443
群 2の異質性(ズレ)      0.550538         1      0.550538      0.458098
群 3の直線性         0.000409222         1   0.000409222       0.98386
群 3の異質性(ズレ)      0.142857         1      0.142857      0.705457
--------------------------------------------------------------------------
共通直線性             0.0667839         1     0.0667839      0.796078
異質性(ズレ)合計        0.744077         3      0.744077      0.862791
非平行性               0.0112959         2     0.0112959      0.994368
--------------------------------------------------------------------------
全体                    0.822157         6

○効力比Rとその95%Bonferroni型同時信頼区間
・群2/群1:R=0.986611 (log(R)=-0.00585419)
  信頼区間:下限=0.757605(log(R)=-0.120557)  上限=1.28484(log(R)=0.108849)
・群3/群1:R=0.946113 (log(R)=-0.0240571)
  信頼区間:下限=0.543893(log(R)=-0.264486)  上限=1.64578(log(R)=0.216372)
・群3/群2:R=0.958952 (log(R)=-0.018203)
  信頼区間:下限=0.622424(log(R)=-0.205914)  上限=1.47743(log(R)=0.169508)

・修正切片差のBonferroni型多重比較(Bonferroni type multiple comparison)
群         - 群              χ2乗値  自由度   有意確率p値
-------------------------------------------------------------
2          - 1            0.0191298       1             1
3          - 1             0.376702       1             1
3          - 2             0.265872       1             1
-------------------------------------------------------------

・修正切片差のBonferroni型95%同時信頼区間(simultaneous confidence interval)
群      - 群        修正切片差    絶対精度        下限        上限
------------------------------------------------------------------
2       - 1          -0.080195     1.38807    -1.46827     1.30788
3       - 1          -0.329552     1.28542    -1.61498    0.955871
3       - 2          -0.249357     1.15773    -1.40709     0.90837
------------------------------------------------------------------


=== プロビット法による平行線検定 ===                               [DANS V7.4]

データ名:tckプロビット法による共分散分析デモデータ (parapl_03.tab 2017-11-27)

用量項目(縦):群(用量)
反応項目(横):効果

・群 1:男

群(用量)＼効果  無効(  %  )  有効(  %  )  合計(  %  )
-----------------------------------------------------
0                 20( 69.0)     9( 31.0)    29(100.0)
10                18( 60.0)    12( 40.0)    30(100.0)
20                16( 53.3)    14( 46.7)    30(100.0)
30                11( 37.9)    18( 62.1)    29(100.0)
-----------------------------------------------------
合計              65( 55.1)    53( 44.9)   118(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
       0         -   0.310345(    9/29   )           0
      10         1        0.4(   12/30   )        0.13     3.87361     3.87308
      20     1.301   0.466667(   14/30   )    0.226667     4.25013     4.42549
      30     1.477    0.62069(   18/29   )        0.45     4.87434     4.74863
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=1.50343+2.23389x
  反応率=50%(y=5)の時の用量:D50=36.7486(log(D50)=1.56524)
共通対数用量-プロビット直線:y=2.03803+1.83505x
  反応率=50%(y=5)の時の用量:D50=41.125(log(D50)=1.61411)

・群 2:女

群(用量)＼効果  無効(  %  )  有効(  %  )  合計(  %  )
-----------------------------------------------------
0                 22( 81.5)     5( 18.5)    27(100.0)
10                20( 71.4)     8( 28.6)    28(100.0)
20                18( 64.3)    10( 35.7)    28(100.0)
30                15( 55.6)    12( 44.4)    27(100.0)
-----------------------------------------------------
合計              75( 68.2)    35( 31.8)   110(100.0)
                                                           観測      期待プロ
実数用量  対数用量     反応率(反応数/総数)  補正反応率  プロビット   ビット(y)
------------------------------------------------------------------------------
       0         -   0.185185(    5/27   )           0
      10         1   0.285714(    8/28   )    0.123377     3.84173      3.6784
      20     1.301   0.357143(   10/28   )    0.211039     4.19718      4.2308
      30     1.477   0.444444(   12/27   )    0.318182     4.52721     4.55394
------------------------------------------------------------------------------
群別対数用量-プロビット直線:y=2.35872+1.45208x
  反応率=50%(y=5)の時の用量:D50=65.9113(log(D50)=1.81896)
共通対数用量-プロビット直線:y=1.84335+1.83505x
  反応率=50%(y=5)の時の用量:D50=52.5047(log(D50)=1.7202)
※反復計算回数=5

                       平行線検定表(分散分析表)
要因                      平方和    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------------
群 1の直線性             2.38004         1       2.38004      0.122894
群 1の異質性(ズレ)      0.241973         1      0.241973      0.622785
群 2の直線性              1.0473         1        1.0473       0.30613
群 2の異質性(ズレ)     0.0296309         1     0.0296309       0.86333
--------------------------------------------------------------------------
共通直線性               3.27863         1       3.27863     0.0701876+
異質性(ズレ)合計        0.271604         2      0.271604      0.873016
非平行性                0.148715         1      0.148715      0.699767
--------------------------------------------------------------------------
全体                     3.69895         4

○効力比Rとその95%信頼区間
・群2/群1:R=0.783264 (log(R)=-0.106092)
  信頼区間:下限=0.335184(log(R)=-0.474716)  上限=1.83034(log(R)=0.262533)

・修正切片差のBonferroni型多重比較(Bonferroni type multiple comparison)
群         - 群              χ2乗値  自由度   有意確率p値
-------------------------------------------------------------
2          - 1             0.348636       1      0.554886
-------------------------------------------------------------

・修正切片差のBonferroni型95%同時信頼区間(simultaneous confidence interval)
群      - 群        修正切片差    絶対精度        下限        上限
------------------------------------------------------------------
2       - 1          -0.194684    0.646237    -0.84092    0.451553
------------------------------------------------------------------


