グラフが簡単に描けます。
[root@c ~]# yum install gnuplot
描く †
[illya@c ~]$ gnuplot
G N U P L O T
Version 4.6 patchlevel 2 last modified 2013-03-14
Build System: Linux x86_64
Copyright (C) 1986-1993, 1998, 2004, 2007-2013
Thomas Williams, Colin Kelley and many others
gnuplot home: http://www.gnuplot.info
faq, bugs, etc: type "help FAQ"
immediate help: type "help" (plot window: hit 'h')
Terminal type set to 'x11'
gnuplot>
gnuplot> plot sin(x)
gnuplot>
gnuplot> quit
[illya@c ~]$

2軸 †
「micrographs_ctf.star」から「_rlnCtfMaxResolution」を表示したい。
[saber@c CtfFind]$ less job003/micrographs_ctf.star
data_
loop_
_rlnMicrographNameNoDW #1
_rlnMicrographName #2
_rlnCtfImage #3
_rlnDefocusU #4
_rlnDefocusV #5
_rlnDefocusAngle #6
_rlnVoltage #7
_rlnSphericalAberration #8
_rlnAmplitudeContrast #9
_rlnMagnification #10
_rlnDetectorPixelSize #11
_rlnCtfFigureOfMerit #12
_rlnCtfMaxResolution #13
MotionCorr/job002/Micrographs/s2_100001_noDW.mrc MotionCor
MotionCorr/job002/Micrographs/s2_100002_noDW.mrc MotionCor
MotionCorr/job002/Micrographs/s2_100003_noDW.mrc MotionCor
MotionCorr/job002/Micrographs/s2_100004_noDW.mrc MotionCor
MotionCorr/job002/Micrographs/s2_100005_noDW.mrc MotionCor
MotionCorr/job002/Micrographs/s2_100006_noDW.mrc MotionCor
:
[saber@c CtfFind]$
[saber@c CtfFind]$ awk '{if($2 ~/#/)N=NR;}END{print N}' job003/micrographs_ctf.star
17
[saber@c CtfFind]$ awk '/_rlnCtfMaxResolution/{print $2}' job003/micrographs_ctf.star | cut -c 2-
13
[saber@c CtfFind]$ tail -n +18 job003/micrographs_ctf.star | awk '{print NR,$13}' | tee job003.list
1 3.174000
2 4.131000
3 3.476000
4 3.830000
5 4.037000
6 3.699000
:
[saber@c CtfFind]$
[saber@c CtfFind]$ awk '{if($2 ~/#/)N=NR;}END{print N}' job004/micrographs_ctf.star
17
[saber@c CtfFind]$ tail -n +18 job004/micrographs_ctf.star | awk '{print $13}' | tee job004.list
[saber@c CtfFind]$ paste job003.list job004.list
1 6.357333 3.174000
2 7.553267 4.131000
3 9.418272 3.476000
4 6.410756 3.830000
5 7.479216 4.037000
6 4.768000 3.699000
;
:
[saber@c CtfFind]$ gnuplot
gnuplot> plot "CtfFind.list" using 1:2 with line title 'ctffind4', "CtfFind.list" using 1:3 with line title 'Gctf'
