Rで厳密な連検定。 exact runs test

Rで厳密な連検定。 exact runs test
正規分布近似を使わない場合。
両側検定にするのが微妙にむずい。
https://ncss-wpengine.netdna-ssl.com/wp-content/themes/ncss/pdf/Procedures/NCSS/Analysis_of_Runs.pdf
の５ページを参考にする

install.packages("randtests")
library(randtests)
runs.test2<-function(x){
		
	r<-length(rle(c(unlist(x)))$length)
	if(r==1){
		return(NA)
	}
	
	n1<-sum(x)
	n2<-sum(1-x)
	n<-length(x)
	rmin<-2
	er<-2*n1*n2/n+1
	rmax<-2*min(n1,n2)+1
	
		rmin<-1
		rrs<-er-abs(r-er)
		rrf<-er+abs(r-er)
		
		ans1<-sum(druns(rmin:rrs,n1,n2))
		ans2<-sum(druns(rrf:rmax,n1,n2))
		ans<-ans1+ans2	
	
	ans
}