# # Exercise the "stats" command. # Skip this demo cleanly if the stats option was not configured. # if (!strstrt(GPVAL_COMPILE_OPTIONS,"+STATS")) { print "No support for stats command" } else { set title "Use of stats command to find min/max/mean before plotting\nOne data column" set style data line set offset 0,0,.5,.5 set autoscale fix set key left Left stats 'orbital_elements.dat' index 1 using 2 prefix "A" set arrow 1 from A_index_min, graph 0.1 to A_index_min, A_min fill set arrow 2 from A_index_max, graph 0.9 to A_index_max, A_max fill set label 1 at A_index_min, graph 0.1 "min" center offset 0,-1 set label 2 at A_index_max, graph 0.9 "max" center offset 0,1 plot 'orbital_elements.dat' index 1 using 0:2 title " Data" lw 2, \ A_mean title " Mean" pause -1 "Hit return to continue" set title "Use of stats command to find min/max/mean before plotting\nTwo data columns" f(x) = log(1+x) stats 'orbital_elements.dat' index 1 using (f($0)):2 prefix "B" set arrow 1 from B_pos_min_y, graph 0.1 to B_pos_min_y, B_min_y fill set arrow 2 from B_pos_max_y, graph 0.9 to B_pos_max_y, B_max_y fill set label 1 at B_pos_min_y, graph 0.1 "min" center offset 0,-1 set label 2 at B_pos_max_y, graph 0.9 "max" center offset 0,1 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" pause -1 "Hit return to continue" reset }