if (FALSE) # Time consuming code part
data(Schizo)
fit.ica <- ICA.BinCont.BS(Dataset = Schizo, Surr = BPRS, True = PANSS_Bin, nb = 10,
Theta.S_0=c(-10,-5,5,10,10,10,10,10), Theta.S_1=c(-10,-5,5,10,10,10,10,10),
Treat=Treat, M=50, Seed=1)
fit.spf <- SPF.BinCont(fit.ica, a=-5, b=5)
summary(fit.spf)
plot(fit.spf)
Run the code above in your browser using DataLab