新版 実用SAS生物統計ハンドブック正誤表頁

2021/8/31
誤正vi http://ace-jp.org/book/favor.html

https://ace-jp.org/book/favor.html

ix http://ace-jp.org/book/favor.html

https://ace-jp.org/book/favor.html

xi Ⅰ．対応のない２群の平均の差の検定としてt 検定を行う。（Student Ⅰ．対応のない２群の平均の差の検定としてt 検定を行う。（Studentの検定、Aspin-Welch の検定）のt 検定、Aspin-Welch の検定）

20 Ⅰ．対応のない２群の平均の差の検定としてt 検定を行う。（Student Ⅰ．対応のない２群の平均の差の検定としてt 検定を行う。（Studentの検定、Aspin-Welch の検定）
34 d010203.what<-list(num=0, toyo_0=0, toyo_9=0)のt 検定、Aspin-Welch の検定）d010203.what<-list(num=0, toyo_0=0, toyo_9=0)

tmp1<-scan(file="./d010203.csv",what=d010203.what,sep="")tmp1<-scan(file="./d010203.csv",what=d010203.what,sep="")d010203<-data.frame(tmp1)

d010203<-data.frame(tmp1)

d010203$diff<-d010203$toyo_9 - d010203$toyo_0

d010203$diff<-d010203$toyo_9 - d010203$toyo_0

t.test(d010203$toyo_0, d010203$toyo_9, paired=T)t.test(d010203$toyo_0, d010203$toyo_9, paired=T)#浮動小数点誤差の影響をなくすため、一度10 倍して整数にする。

wilcox.test(d010203$toyo_0, d010203$toyo_9, paired=T, correct =wilcox.test(d010203$toyo_0*10, d010203$toyo_9*10, paired = T,FALSE)

correct = FALSE)

wilcox.test(d010203$toyo_0, d010203$toyo_9, paired=T)
34 Wilcoxon signed rank test

wilcox.test(d010203$toyo_0*10, d010203$toyo_9*10, paired = T )Wilcoxon signed rank test

data: d010203$toyo_0 and d010203$toyo_9

data: d010203$toyo_0 * 10 and d010203$toyo_9 * 10

V = 66.5, p-value = 0.6273

V = 75.5, p-value = 0.9596

alternative hypothesis: true location shift is not equal to 0

alternative hypothesis: true location shift is not equal to 078 Wilcoxon signed rank test with continuity correction

Wilcoxon signed rank test with continuity correction

data: d010203$toyo_0 and d010203$toyo_9

data: d010203$toyo_0 * 10 and d010203$toyo_9 * 10

V = 66.5, p-value = 0.6446

V = 75.5, p-value = 0.9798

alternative hypothesis: true location shift is not equal to 0

alternative hypothesis: true location shift is not equal to 0259 【参考文献】【参考文献】

・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大

学出版会．

学出版会．

・宮原英夫，丹後俊郎 編．（1995）．「医学統計学ハンドブッ・宮原英夫，丹後俊郎 編．（1995）．「医学統計学ハンドブッ

ク」，朝倉書店．

ク」，朝倉書店．
・大橋靖雄，浜田知久馬，魚住龍史 著．（2016）．「生存時間解析応用編」，東京大学出版会．

265 【参考文献】【参考文献】

・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大

学出版会．

学出版会．
・大橋靖雄，浜田知久馬，魚住龍史 著．（2016）．「生存時間解析応用編」，東京大学出版会．

270 【参考文献】【参考文献】

・浜田知久馬．（2013）．「SAS 生存時間解析プロシジャの最新の・浜田知久馬．（2013）．「SAS 生存時間解析プロシジャの最新の機能拡張」，SUGI-J 論文集'13，pp.3-72．

機能拡張」，SUGI-J 論文集'13，pp.3-72．

・Klein J.P., Moeschberger M.L. (1997). Survival Analysis:

・Klein J.P., Moeschberger M.L. (1997). Survival Analysis:

Techniques for Censored and Truncated Data, Springer-Verlag.

Techniques for Censored and Truncated Data, Springer-Verlag.

・Westfall P.H., Tobias R.D., Rom D., Wolfinger R.D., Hochberg Y.

・Westfall P.H., Tobias R.D., Rom D., Wolfinger R.D., Hochberg Y.

(1999). Multiple Comparisons and Multiple Tests Using the SAS(1999). Multiple Comparisons and Multiple Tests Using the SAS

System, SAS Institute Inc.

System, SAS Institute Inc.

・Hsu J.C. (1992). The Factor Analytic Approach to Simultaneous・Hsu J.C. (1992). The Factor Analytic Approach to Simultaneous

Inference in the General Linear Model, Journal of Computational

Inference in the General Linear Model, Journal of Computational

and Graphical Statistics, 1, pp.151-168.

and Graphical Statistics, 1, pp.151-168.

・Kramer C.Y. (1956). Extension of Multiple Range Tests to Group・Kramer C.Y. (1956). Extension of Multiple Range Tests to Group

Means with Unequal Numbers of Replications, Biometrics, 12,

Means with Unequal Numbers of Replications, Biometrics, 12,pp.307-310.

pp.307-310.

・Edwards D., Berry J.J. (1987). The Efficiency of Simulation-Based ・Edwards D., Berry J.J. (1987). The Efficiency of Simulation-BasedMultiple Comparisons, Biometrics, 43, pp.913–928.

Multiple Comparisons, Biometrics, 43, pp.913–928.

・Westfall P.H., Young S.S. (1993). Resampling-Based Multiple

・Westfall P.H., Young S.S. (1993). Resampling-Based Multiple

Testing: Examples and Methods for p-Value Adjustment, John

Testing: Examples and Methods for p-Value Adjustment, John

Wiley & Sons, Inc.

Wiley & Sons, Inc.

・Klein J.P., Moeschberger M.L. (2013). Survival Analysis:

・Klein J.P., Moeschberger M.L. (2013). Survival Analysis:

Techniques forCensored and Truncated Data, Springer-Verlag.

Techniques forCensored and Truncated Data, Springer-Verlag.

・永田靖，吉田道弘 著．（1997）．「統計的多重比較法の基礎」，

・永田靖，吉田道弘 著．（1997）．「統計的多重比較法の基礎」，

サイエンティスト社．

サイエンティスト社．
・大橋靖雄，浜田知久馬，魚住龍史 著．（2016）．「生存時間解析応用編」，東京大学出版会．

278 【参考文献】【参考文献】

・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大・大橋靖雄，浜田知久馬 著．（1995）．「生存時間解析」，東京大学出版会．

学出版会．

・宮原英夫，丹後俊郎 編．（1995）．「医学統計学ハンドブッ・宮原英夫，丹後俊郎 編．（1995）．「医学統計学ハンドブッ

ク」，朝倉書店．

ク」，朝倉書店．
・大橋靖雄，浜田知久馬，魚住龍史 著．（2016）．「生存時間解析応用編」，東京大学出版会．

283 【参考文献】【参考文献】

・浜田知久馬．（2013）．「SAS 生存時間解析プロシジャの最新の・浜田知久馬．（2013）．「SAS 生存時間解析プロシジャの最新の機能拡張」，SUGI-J'13 論文集，pp.3-72．

機能拡張」，SUGI-J'13 論文集，pp.3-72．

・Kalbfleisch J.D., Prentice R.L. (1980). The Statistical Analysis of・Kalbfleisch J.D., Prentice R.L. (1980). The Statistical Analysis of

Failure Time Data, John Wiley & Sons, Inc.

Failure Time Data, John Wiley & Sons, Inc.
・大橋靖雄，浜田知久馬，魚住龍史 著．（2016）．「生存時間解析応用編」，東京大学出版会．

358 do I=0 to 100;

do I=0 to 100;

C=(A+B)/2;

C=(A+B)/2;

if ((C-A)<=1.0E-7 or (B-C)<=1.0E-7) then continue;

if ((C-A)<=1.0E-7 or (B-C)<=1.0E-7) then continue;

if (I<IMAX) then do;

if (I<IMAX) then do;

C=(A*FB-B*FA)/(FB-FA);

C=(A*FB-B*FA)/(FB-FA);

if (C<=A) then C=A+0.6*2.22E-16*abs(A);if (C<=A) then C=A+0.6*2.22E-16*abs(A);else if (C>-B) then C=B-0.6*2.22E-16*abs(B);else if (C>-B) then C=B-0.6*2.22E-16*abs(B);

end;

end;

FC=(&P-(C-&DELTA))/((C-&DELTA)*(1-C+&DELTA))+(&P-FC=(&P-(C-&DELTA))/((C-&DELTA)*(1-C+&DELTA))+(&P-C)/(C*(1-C));

C)/(C*(1-C));

if FC=0 then goto ED;

if FC=0 then goto ED;

if (FC>0) and (FA>0) then do;

if (FC>0) and (FA>0) then do;

A=C; FA=FC;

A=C; FA=FC;

end;

end;

else do;else do;

B=C; FB=FC;

B=C; FB=FC;

end;end;

end;

end;

360 proc power;

proc power;

twosamplefreq test = FM

twosamplefreq test = FM

sides

sides

=L=U

groupproportions = (0.75 0.75)

groupproportions = (0.75 0.75)

nullpdiff= 0.10

nullpdiff= -0.10alpha= 0.05alpha= 0.05power= 0.80power= 0.80npergroup

npergroup= .;360 分布
手法

= .;

run;

run;
Asymptotic normal

分布

Asymptotic normal

Normal approximation

手法

Normal approximation

裾の数L

帰無仮説の比率差
有意水準

裾の数
0.1

0.05

U

帰無仮説の比率差
有意水準-0.10.05

グループ 1 の比率

0.75

グループ 1 の比率

グループ 2 の比率

0.75

グループ 2 の比率

0.75

名目の検出力0.8

名目の検出力0.8

362 for (ix in 0:100) {0.75for (ix in 0:100) {C <- (A+B)/2C <- (A+B)/2

if((C-A)<=1.0E-7 | (B-C)<=1.0E-7) {next}

if((C-A)<=1.0E-7 | (B-C)<=1.0E-7) {next}

if(ix<imax) {if(ix<imax) {C <- (A*fb-B*fa)/(fb-fa)

C <- (A*fb-B*fa)/(fb-fa)

if(C<=A) {if(C<=A) {C <- A+0.6*2.22E-16*abs(A)C <- A+0.6*2.22E-16*abs(A)} else if (C > -B) {} else if (C > -B) {C <- B-0.6*2.22E-16*abs(B)C <- B-0.6*2.22E-16*abs(B)}}}}fc <- (P-(C-DELTA))/((C-DELTA)*(1-C+DELTA))+(P-C)/(C*(1-fc <- (P-(C-DELTA))/((C-DELTA)*(1-C+DELTA))+(P-C)/(C*(1-C))C))

if(fc==0) {if(fc==0) {endsw <- 1

endsw <- 1
breakbreak}}if((fc>0)&(fa>0)) {if((fc>0)&(fa>0)) {A <- CA <- C

fa <- fc

fa <- fc} else {} else {B <- C

B <- Cfb <- fc

fb <- fc}}}364}【参考文献】【参考文献】
・緒方宏泰，鮫島政義 編．（1992）．「医薬品のバイオアベイラビ
リティと生物学的同等性試験第 2 版」，薬業時報社．
・Fleiss J.L., Levin B., Paik M.C. (2003). Statistaical Methods forRates and Proportions, Third Edition, John Wiley & Sons, Inc.
・奥村晴彦．（1991）．「C 言語による最新アルゴリズム事典」，技術評論社．
・丹後俊郎．（1994）．「医学への統計学」，朝倉書店．

・緒方宏泰，鮫島政義 編．（1992）．「医薬品のバイオアベイラビ
リティと生物学的同等性試験第 2 版」，薬業時報社．
・Fleiss J.L., Levin B., Paik M.C. (2003). Statistical Methods forRates and Proportions, Third Edition, John Wiley & Sons, Inc.
・奥村晴彦．（1991）．「C 言語による最新アルゴリズム事典」，技術評論社．
・丹後俊郎．（1994）．「医学への統計学」，朝倉書店．
・魚住龍史，飯塚政人，浜田知久馬．（2015）．「非劣性試験における割合の差の信頼区間と例数設計」，SUGI-J'15論文集．pp.422431．