R/s_1.R

In fast: Implementation of the Fourier Amplitude Sensitivity Test (FAST)

`s` <-
function(m, factor=1,cukier=TRUE){
if(cukier){
	#min_runs_cukier75 <- c(NA,NA,38,78,142,182,334,486,630,806,974,1158,1374,1814,2038,2446,2734,3310,3838,4174,4702,5542,6174,6854,7110,8182,8934,9590,11358,11526,13014,14206,15046,16702,18374,19334,20422,21550,22678,14934,25782,27478,29486,31486,33950,36550,37854,39814,41518,43606)/2
	min_runs <- min_runs_cukier75
}else{
	#min_runs_mcrae82 <- c( 0,15,27,47,79,99,175,251,323,411,495,587,695,915,1027)
	min_runs <- min_runs_mcrae82
}
r=round(min_runs[m]*factor)
s=pi/r*(2*(1:r)-r-1)/2
return(s)
}