...
 
Commits (3)
......@@ -15,7 +15,7 @@ set xlabel '# Exposures'
set key top right
set bars 0
set term post size 5in, 9in eps enhanced color
set term post size 5in, 5in eps enhanced color
#set term png size 700,400
set out "Ecount_extinction.eps"
......
......@@ -15,7 +15,7 @@ set xlabel 'Time to first Exposure (days)'
set key bottom right
set bars 0
set term post size 5in, 9in eps enhanced color
set term post size 5in, 5in eps enhanced color
#set term png size 700,400
set out "Efirst_extinction.eps"
......
......@@ -15,12 +15,15 @@ set xlabel 'Hit probability'
set key bottom right
set bars 0
set term post size 5in, 9in eps enhanced color
set term post size 5in, 5in eps enhanced color
#set term png size 700,400
set out "Zhit_extinction.eps"
plot "out/Zhit_extinction.txt" using 3:1:ytic(2) w p t "", \
#using (0):($1-0.3):3:($1+0.3):ytic(2) w boxxyerrorbars t "", \
set out "Zhit_outbreak.eps"
plot "out/Zhit_outbreak.txt" using 3:1:ytic(2) w p t "", \
......
No preview for this file type
......@@ -3,10 +3,10 @@ histo <- hist(d$V1, breaks=40, plot=FALSE)
hdata = data.frame(histo$mids, histo$density)
write.table(hdata, file="out/finished_extinction_t_histo.txt", col.names = FALSE, row.names = FALSE)
d=read.table("out/finished_hit_stop_t.txt")
histo <- hist(d$V1, breaks=40, plot=FALSE)
hdata = data.frame(histo$mids, histo$density)
write.table(hdata, file="out/finished_hit_stop_t_histo.txt", col.names = FALSE, row.names = FALSE)
#d=read.table("out/finished_hit_stop_t.txt")
#histo <- hist(d$V1, breaks=40, plot=FALSE)
#hdata = data.frame(histo$mids, histo$density)
#write.table(hdata, file="out/finished_hit_stop_t_histo.txt", col.names = FALSE, row.names = FALSE)
......
local NRUNS = assert(tonumber(arg[1]), 'missing count parameter')
local E_PLOT_COUNT = 70
local E_PLOT_COUNT = 40
local mathutil = require('lib.mathutil')
......
......@@ -15,7 +15,7 @@ set term post size 5in, 2in eps enhanced color
set out "states.eps"
filename = 'out/states1.log'
filename = 'out/states8.log'
stats filename nooutput
......