# set terminal canvas  rounded size 600,400 enhanced fsize 10 lw 1.6 fontscale 1 name "stats_2" jsdir "."
# set output 'stats.2.js'
set key fixed left top vertical Left noreverse enhanced autotitle nobox
set label 1 "min" at 2.07944, graph 0.1, 0 center norotate back nopoint offset character 0, -1, 0
set label 2 "max" at 2.30259, graph 0.9, 0 center norotate back nopoint offset character 0, 1, 0
set arrow 1 from 2.07944, graph 0.1, 0 to 2.07944, 2.20087, 0.00000 head back filled linewidth 1.000 dashtype solid
set arrow 2 from 2.30259, graph 0.9, 0 to 2.30259, 3.13972, 0.00000 head back filled linewidth 1.000 dashtype solid
set style increment default
set offsets 0, 0, 0.5, 0.5
set style data lines
set title "Use of stats command to find min/max/mean before plotting\nTwo data columns" 
set autoscale tfixmin
set autoscale tfixmax
set autoscale ufixmin
set autoscale ufixmax
set autoscale vfixmin
set autoscale vfixmax
set xrange [ * : * ] noreverse writeback
set autoscale xfixmin
set autoscale xfixmax
set x2range [ * : * ] noreverse writeback
set autoscale x2fixmin
set autoscale x2fixmax
set yrange [ * : * ] noreverse writeback
set autoscale yfixmin
set autoscale yfixmax
set y2range [ * : * ] noreverse writeback
set autoscale y2fixmin
set autoscale y2fixmax
set zrange [ * : * ] noreverse writeback
set autoscale zfixmin
set autoscale zfixmax
set cbrange [ * : * ] noreverse writeback
set autoscale cbfixmin
set autoscale cbfixmax
set rrange [ * : * ] noreverse writeback
set autoscale rfixmin
set autoscale rfixmax
f(x) = log(1+x)
A_records = 20
A_invalid = 0
A_headers = 0
A_blank = 1
A_blocks = 1
A_outofrange = 0
A_columns = 11
A_mean = 2.54077965
A_stddev = 0.222745079388137
A_ssd = 0.228531629485718
A_skewness = 0.978323019217487
A_kurtosis = 3.60561959806089
A_adev = 0.178523515
A_mean_err = 0.0498073139165462
A_stddev_err = 0.0352190894230769
A_skewness_err = 0.547722557505166
A_kurtosis_err = 1.09544511501033
A_sum = 50.815593
A_sumsq = 130.103532004915
A_min = 2.20087
A_max = 3.139718
A_median = 2.4609945
A_lo_quartile = 2.3865975
A_up_quartile = 2.6562155
A_index_min = 7
A_index_max = 9
B_records = 20
B_invalid = 0
B_headers = 0
B_blank = 1
B_blocks = 1
B_outofrange = 0
B_columns = 11
B_mean_x = 2.11678082303767
B_stddev_x = 0.792134394132666
B_ssd_x = 0.812712740322205
B_skewness_x = -1.11744391283858
B_kurtosis_x = 3.52503054370255
B_adev_x = 0.632964368155614
B_mean_err_x = 0.177126635259625
B_stddev_err_x = 0.125247444920837
B_skewness_err_x = 0.547722557505166
B_kurtosis_err_x = 1.09544511501033
B_sum_x = 42.3356164607535
B_sumsq_x = 102.16475902296
B_min_x = 0.0
B_max_x = 2.99573227355399
B_median_x = 2.35024018289621
B_lo_quartile_x = 1.70059869083108
B_up_quartile_x = 2.740319461671
B_index_min_x = 0
B_index_max_x = 19
B_mean_y = 2.54077965
B_stddev_y = 0.222745079388137
B_ssd_y = 0.228531629485718
B_skewness_y = 0.978323019217487
B_kurtosis_y = 3.60561959806089
B_adev_y = 0.178523515
B_mean_err_y = 0.0498073139165462
B_stddev_err_y = 0.0352190894230769
B_skewness_err_y = 0.547722557505166
B_kurtosis_err_y = 1.09544511501033
B_sum_y = 50.815593
B_sumsq_y = 130.103532004915
B_min_y = 2.20087
B_max_y = 3.139718
B_median_y = 2.4609945
B_lo_quartile_y = 2.3865975
B_up_quartile_y = 2.6562155
B_index_min_y = 7
B_index_max_y = 9
B_slope = -0.066943671163917
B_intercept = 2.68248472934352
B_slope_err = 0.0643729477667789
B_intercept_err = 0.145491947243556
B_correlation = -0.238067590736956
B_sumxy = 106.725360630742
B_pos_min_y = 2.07944154167984
B_pos_max_y = 2.30258509299405
## Last datafile plotted: "orbital_elements.dat"
plot 'orbital_elements.dat' index 1 using (f($0)):2 title "  Data" lw 2,      B_mean_y title "  Mean",      B_slope * x + B_intercept title "Linear fit"