source("MCanalysisNPI-latter-plot-PIRAT.R")
./clean-up.sh 
q()
means*1000
means*1000
q()
