yoden_fun <- function(data,filename)

{

sum = 0

cut = sort(unique(data[,1]))

for(i in 1:length(sort(unique(data[,1]))))

{

dat = table(data[,2],data[,1]<cut[i])

if(ncol(dat)==1)

{

dat = cbind(dat,0)

}

tokuido =dat[1,2]/(dat[1,1]+dat[1,2])

kando = dat[2,1]/(dat[2,1]+dat[2,2])

if((kando+tokuido)>sum)

{

sum = kando + tokuido

dat2 = dat

cutpoint = cut[i]

}

}

ans <- fisher.test(dat2)

p_value <- ans[1]

tokuido =dat2[1,2]/(dat2[1,1]+dat2[1,2])

kando = dat2[2,1]/(dat2[2,1]+dat2[2,2])

N <- dat2[1]+dat2[2]+dat2[3]+dat2[4]

seitou = (dat2[1,2] + dat2[2,1])/N

result = cbind(colnames(data)[1],cutpoint,dat2[1,2],dat2[2,2],dat2[1,1],dat2[2,1],p_value,kando,tokuido,kando+tokuido,seitou,N)

colnames(result) = c(“名前”,“カットポイント”,“a”,“b”,“c”,“d”,“p値”,“感度”,“特異度”,“感度+特異度”,“正答率”,“N”)

write.table(result,filename,append=T,quote=F,sep=“,”,row.names=F,col.names=T,fileEncoding=“CP932”)

write.table(“”,filename,append=T,quote=F,sep=“,”,row.names=F,col.names=F,fileEncoding=“CP932”)

}

fn = “yoden.csv”

data = df[,c(2,4)]

yoden_fun(data,fn)