偽陽性率

• 今、検査をしている検査センターがある
• 偽陽性率がpという
• 検査センターでは、検査を続けている
• このセンターで、N回の検査をしたとき、そのうち、１件でも偽陽性である確率は

である

• このセンターでどんなに検査の精度をよくして(pを小さくして)も、Nが多くなれば偽陽性は避けられない
• では、個々の検査を受け取る、被検者にとっての偽陽性率はというと、pのまま

K<-10
as<-0:(-K)
lN<-0:K
# 偽陽性率は10^as[i]
# 検査件数は10^lN[j]

res<-matrix(0,length(as),length(lN))
for(i in 1:length(as)){
	p<-10^(as[i])
	N<-10^lN
	res[i,]<-1-(1-p)^N
}
matplot(t(res),type="l")