#from bf_sd.r pe<-3.29 # t-ratio for parameter estimate x<-seq(-10,10,0.01) zero<-(x==0) absx<-abs(x) cent<-(x-pe) like<-1/(exp(cent*cent/2)) bfsign<-0 bfpo<-0 bfno<-0 bf10<-0 shape<-6 # shape parameter for Gamma distribution rate<-0 for (i in 1:150){ rate[i]<-.05*i f<-dgamma(absx,shape,rate[i]) negf<-f*(x<0) posf<-f*(x>0) pos<-posf*like*.01 neg<-negf*like*.01 null<-zero*like bfsign[i]<-sum(pos)/sum(neg) bfpo[i]<-sum(pos)/sum(null) bfno[i]<-sum(neg)/sum(null) bf10[i]<-(sum(pos)+sum(neg))/(2*sum(null)) } all<-cbind(rate,bf10,bfpo,bfno,bfsign)