#
# Use an array to hold the sum of values in each columnn
#
array Sum[20]
do for [i=5:13] {
stats 'immigration.dat' skip 2 using (column(i)) nooutput
Sum[i] = STATS_sum
print sprintf("Sum[%2d] = %d",i,Sum[i])
}
#
# Use the stored sum to normalize values in each column
#
set style fill transparent solid 0.2 border -1
set style data filledcurve y=0
set tics scale 0
set xrange [] noextend
set xlabel "Decade"
set title "Use of an array to aid normalization and to plot summed values"
set multiplot
plot for [i=5:13] 'immigration.dat' using 1:(column(i)/Sum[i]) title columnhead
#
# Plot the array values themselves in a separate inset
#
set size 0.4, 0.3
set origin 0.05, 0.6
unset xtics
unset ytics
set border 0
set xlabel "Total"
set xrange [4:15]
set yrange [0:*]
unset title
set boxwidth 1.0 absolute
plot Sum using 1:2:($1-4) with boxes lc variable notitle
unset multiplot
|
reset
set title "Fit function to values stored in an array"
set xtics nomirror
array A[100]
do for [i=1:100] { A[i] = sin(2*pi*i/100.) + 0.1*rand(0) }
a = b = c = 0.01
set key title "before fit"
plot A with points title "Array A", c+cos(a+b*x) with lines
|
set fit nolog fit c+cos(a+b*x) A using 1:2 via a,b,c set key title "after fit" plot A with points title "Array A", c+cos(a+b*x) with lines |
reset
#
# Illustrate loading an array from a data file
# This could be done using a preliminary 'plot' command, perhaps together
# with 'set table' if you don't actually want generate a plot, but here
# we use 'stats'. The first run is just to pick up the array size needed.
# The second run stores one value into each array element.
#
set title "Illustrate loading an array from a column in a data file\n".\
"Note that first data point in the file is 'line 0'\n".\
"but it goes into array element Data[1]"
set yrange [0:20]
set border 3
set tics nomirror
set key left center nobox
FILE = 'orbital_elements.dat'
stats FILE index 2 using 2 prefix "A"
array Data[A_records]
stats FILE index 2 using (Data[int($0+1)] = $2) prefix "B"
plot FILE index 2 using 2 title "Data from file" with lines lw 2, \
Data using 0:2:(sprintf("Data[%d]",$1)) with labels point pt 6 ps 2 center offset 0,1 title "Data[i]"
|
reset |