グラフが簡単に描けます。

[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 ~]$

2018y08m25d_010308514.png

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'

2018y08m25d_022707074.png


トップ   編集 添付 複製 名前変更     ヘルプ   最終更新のRSS
Last-modified: 2018-08-25 (土) 02:35:46 (56d)