# draw 9,999 random draws from the prior given sample configuration of Stone et al. data from 18 species pairs of parisitoid wasps codistributed across Iberian peninsula and Hungary # see 2012 Stone, G.N., K. Lohse, J. A. Nicholls, P. Fuentes-Utrilla, F. Sinclair, K. Schönrogge, G. Csóka, G. Melika, J-L. Nieves-Aldrey, J. Pujade-Villar, M. Tavakoli, R. R. Askew and M. J. Hickerson. Reconstructing Community Assembly in Time and Space Reveals Enemy Escape in a Western Palearctic Insect Community. Current Biology; 22: 532-537 # batch file with priors and sample configuration is HU_IB_batch_pars.txt #cd $HOME/msbayes/trunk/src #./msbayes.pl -c HU_IB_batch_pars.txt -r 9999 -o prior_IB_HUPARSnewPsi0_10K #R script Prior <- scan(("prior_IB_HUPARSnewPsi0_10K"),nlines=9999) PRIOR<-matrix(Prior,ncol=383,byrow=T) pdf("Psi_vs_Omega.pdf") plot(PRIOR[,2],PRIOR[,5],xlim=c(0,18),ylim=c(0,0.7),lty=2,lwd=0.5) dev.off()