10 More Discussions You Might Find Interesting
1. Shell Programming and Scripting
Hi,
I have 2 files with similar structure - reference and test that I would like to BIN both and generate the comparison.
input files structure is:
a 3
b 10
c 3
d 7
e 1
f 4
g 9
h 6
I would like the output to be (lets say both reference and test are the file above - no diff)
BIN ... (3 Replies)
Discussion started by: yan1
3 Replies
2. Programming
I am new to R and would like to calculate the percentage frequency distribution of h1 and h2. How can I combine h1 and h2 in one plot? I tried the following code.
h1=c(5.18,4.61,3.30,7.58,3.00,3.80,1.95,2.67,2.77,2.73,2.33,3.36,3.50,1.91,4.25,3.87,2.86,2.26,2.00,3.86,3.33,3.59,4.00)... (0 Replies)
Discussion started by: ayyappa342
0 Replies
3. Shell Programming and Scripting
I have a single file that looks like this:
1.62816
1.62816
0.86941
0.86941
0.731465
0.731465
1.03174
1.03174
0.769444
0.769444
0.981181
0.981181
1.14681
1.14681
1.00511
1.00511
1.20385
1.20385 (2 Replies)
Discussion started by: kayak
2 Replies
4. Shell Programming and Scripting
My data looks like this:
frame phi psi
0 68.466774 -58.170494
1 75.128593 -51.646816
2 76.083946 -64.300102
3 77.578056 -76.464218
4 63.180199 -76.067680
5 77.203979 -58.560757
6 66.574913 -60.000214
7 73.218269 -70.978203
8 70.956879 -76.096558
9 65.538872 -76.716568... (19 Replies)
Discussion started by: chrisjorg
19 Replies
5. UNIX for Dummies Questions & Answers
Dear all,
I have numerous dat files (a.dat, b.dat...) containing 500 numeric values each. I would like to count them, based on their range and obtain a histogram or a counter.
INPUT:
a.dat
1.3
2.16
0.34
......
b.dat
1.54
0.94
3.13
.....
... (2 Replies)
Discussion started by: chen.xiao.po
2 Replies
6. UNIX for Dummies Questions & Answers
hello experts,
I have been trying to plot a histogram of a data like;
-54
-56
-43
-65
-67
-78
...
I have 156.000 rows of these kind of values between 0 and -100. I just want to make x axis takes values 5 spacing in between like;
-100 -95 -90 .............. 0
And I want y... (4 Replies)
Discussion started by: enes71
4 Replies
7. UNIX for Dummies Questions & Answers
Hi I am trying to fit my histogram data with a gaussian model and am encountering two problems:
1. I can't seem to fit the histogram data with a model
fit y(x) 'bin.txt' using 2:xtic(1) via a,b,c (error: need 2 to 7 using specs)
2. Even when I manually guess the correct parameters for my fit and... (1 Reply)
Discussion started by: goffinj
1 Replies
8. UNIX for Dummies Questions & Answers
Dear all,
Im hoping you may be able to help with a little problem Im having using xmgrace.
I want to plot a dataset as a line curve e.g 1,2,3,4,5,6,15,16,17,18,19,20...
what i would like to do is not plot the missing bits of the sequence, i.e 1-6 would be a straight line then there would... (0 Replies)
Discussion started by: Mish_99
0 Replies
9. UNIX for Dummies Questions & Answers
Dear All,
I have just discovered Xmgrace these days. I find it quite useful for my work so I would like to use it further. Nevertheless there are few things which I did not understood how can be done.
1) Let's say that I have one file with 4 columns and I want to plot column 1 versus the... (1 Reply)
Discussion started by: Eduard
1 Replies
10. Linux
I want to install grace (plotting package; formerly called "xmgr"; also called xmgrace) on a standard PC running Red Hat Linux. I was hoping to find someprecompiled executables, but I only found the complete source files with Makefiles, ".rpm" files, and ".so" files, but I am not experienced enough... (2 Replies)
Discussion started by: karman
2 Replies
funhist(1) SAORD Documentation funhist(1)
NAME
funhist - create a 1D histogram of a column (from a FITS binary table or raw event file) or an image
SYNOPSIS
funhist [-n|-w|-T] <iname> [column] [[lo:hi:]bins]
OPTIONS
-n # normalize bin value by the width of each bin
-w # specify bin width instead of number of bins in arg3
-T # output in rdb/starbase format (tab separators)
DESCRIPTION
funhist creates a one-dimensional histogram from the specified columns of a FITS Extension binary table of a FITS file (or from a non-FITS
raw event file), or from a FITS image or array, and writes that histogram as an ASCII table. Alternatively, the program can perform a 1D
projection of one of the image axes.
The first argument to the program is required, and specifies the Funtools file: FITS table or image, raw event file, or array. If "stdin"
is specified, data are read from the standard input. Use Funtools Bracket Notation to specify FITS extensions, and filters.
For a table, the second argument also is required. It specifies the column to use in generating the histogram. If the data file is of type
image (or array), the column is optional: if "x" (or "X"), "y" (or "Y") is specified, then a projection is performed over the x (dim1) or y
(dim2) axes, respectively. (That is, this projection will give the same results as a histogram performed on a table containing the equiva-
lent x,y event rows.) If no column name is specified or "xy" (or "XY") is specified for the image, then a histogram is performed on the
values contained in the image pixels.
The argument that follows is optional and specifies the number of bins to use in creating the histogram and, if desired, the range of bin
values. For image and table histograms, the range should specify the min and max data values. For image histograms on the x and y axes,
the range should specify the min and max image bin values. If this argument is omitted, the number of output bins for a table is calcu-
lated either from the TLMIN/TLMAX headers values (if these exist in the table FITS header for the specified column) or by going through the
data to calculate the min and max value. For an image, the number of output bins is calculated either from the DATAMIN/DATAMAX header val-
ues, or by going through the data to calculate min and max value. (Note that this latter calculation might fail if the image cannot be fit
in memory.) If the data are floating point (table or image) and the number of bins is not specified, an arbitrary default of 128 is used.
For binary table processing, the -w (bin width) switch can be used to specify the width of each bin rather than the number of bins. Thus:
funhist test.ev pha 1:100:5
means that 5 bins of width 20 are used in the histogram, while:
funhist -w test.ev pha 1:100:5
means that 20 bins of width 5 are used in the histogram.
The data are divvied up into the specified number of bins and the resulting 1D histogram (or projection) is output in ASCII table format.
For a table, the output displays the low_edge (inclusive) and hi_edge (exclusive) values for the data. For example, a 15-row table con-
taining a "pha" column whose values range from -7.5 to 7.5 can be processed thus:
[sh] funhist test.ev pha
# data file: /home/eric/data/test.ev
# column: pha
# min,max,bins: -7.5 7.5 15
bin value lo_edge hi_edge
------ --------- --------------------- ---------------------
1 22 -7.50000000 -6.50000000
2 21 -6.50000000 -5.50000000
3 20 -5.50000000 -4.50000000
4 19 -4.50000000 -3.50000000
5 18 -3.50000000 -2.50000000
6 17 -2.50000000 -1.50000000
7 16 -1.50000000 -0.50000000
8 30 -0.50000000 0.50000000
9 16 0.50000000 1.50000000
10 17 1.50000000 2.50000000
11 18 2.50000000 3.50000000
12 19 3.50000000 4.50000000
13 20 4.50000000 5.50000000
14 21 5.50000000 6.50000000
15 22 6.50000000 7.50000000
[sh] funhist test.ev pha 1:6
# data file: /home/eric/data/test.ev
# column: pha
# min,max,bins: 0.5 6.5 6
bin value lo_edge hi_edge
------ --------- --------------------- ---------------------
1 16 0.50000000 1.50000000
2 17 1.50000000 2.50000000
3 18 2.50000000 3.50000000
4 19 3.50000000 4.50000000
5 20 4.50000000 5.50000000
6 21 5.50000000 6.50000000
[sh] funhist test.ev pha 1:6:3
# data file: /home/eric/data/test.ev
# column: pha
# min,max,bins: 0.5 6.5 3
bin value lo_edge hi_edge
------ --------- --------------------- ---------------------
1 33 0.50000000 2.50000000
2 37 2.50000000 4.50000000
3 41 4.50000000 6.50000000
For a table histogram, the -n(normalize) switch can be used to normalize the bin value by the width of the bin (i.e., hi_edge-lo_edge):
[sh] funhist -n test.ev pha 1:6:3
# data file: test.ev
# column: pha
# min,max,bins: 0.5 6.5 3
# width normalization (val/(hi_edge-lo_edge)) is applied
bin value lo_edge hi_edge
------ --------------------- --------------------- ---------------------
1 16.50000000 0.50000000 2.50000000
2 6.16666667 2.50000000 4.50000000
3 4.10000000 4.50000000 6.50000000
This could used, for example, to produce a light curve with values having units of counts/second instead of counts.
For an image histogram, the output displays the low and high image values (both inclusive) used to generate the histogram. For example, in
the following example, 184 pixels had a value of 1, 31 had a value of 2, while only 2 had a value of 3,4,5,6, or 7:
[sh] funhist test.fits
# data file: /home/eric/data/test.fits
# min,max,bins: 1 7 7
bin value lo_val hi_val
------ --------------------- --------------------- ---------------------
1 184.00000000 1.00000000 1.00000000
2 31.00000000 2.00000000 2.00000000
3 2.00000000 3.00000000 3.00000000
4 2.00000000 4.00000000 4.00000000
5 2.00000000 5.00000000 5.00000000
6 2.00000000 6.00000000 6.00000000
7 2.00000000 7.00000000 7.00000000
For the axis projection of an image, the output displays the low and high image bins (both inclusive) used to generate the projection. For
example, in the following example, 21 counts had their X bin value of 2, etc.:
[sh] funhist test.fits x 2:7
# data file: /home/eric/data/test.fits
# column: X
# min,max,bins: 2 7 6
bin value lo_bin hi_bin
------ --------------------- --------------------- ---------------------
1 21.00000000 2.00000000 2.00000000
2 20.00000000 3.00000000 3.00000000
3 19.00000000 4.00000000 4.00000000
4 18.00000000 5.00000000 5.00000000
5 17.00000000 6.00000000 6.00000000
6 16.00000000 7.00000000 7.00000000
[sh] funhist test.fits x 2:7:2
# data file: /home/eric/data/test.fits
# column: X
# min,max,bins: 2 7 2
bin value lo_bin hi_bin
------ --------------------- --------------------- ---------------------
1 60.00000000 2.00000000 4.00000000
2 51.00000000 5.00000000 7.00000000
You can use gnuplot or other plotting programs to graph the results, using a script such as:
#!/bin/sh
sed -e '1,/---- .*/d
/^$/,$d' |
awk '
BEGIN{print "set nokey; set title "funhist"; set xlabel "bin"; set ylabel "counts"; plot "-" with boxes"}
{print $3, $2, $4-$3}' |
gnuplot -persist - 1>/dev/null 2>&1
Similar plot commands are supplied in the script funhist.plot:
funhist test.ev pha ... | funhist.plot gnuplot
SEE ALSO
See funtools(7) for a list of Funtools help pages
version 1.4.2 January 2, 2008 funhist(1)