物理学者(ポスドク)による日々の研究生活のメモ書きです ( python/emacs/html/Japascript/シェルスクリプト/TeX/Mac/C言語/Linux/git/tmux/R/ポケモンGO)
前回は塗りつぶしのプロットを書いた
■ 過去記事 : 【gnuplot 6】gsl-histogramの結果を塗りつぶしてプロットしたい
(今気づいたけどgnuplot6になってますよ・・・おいおい・・・)
今回は、with stepsでこれをプロットしたい
今まで単一のデータをプロットするときは、塗りつぶしプロットでも用は足りていたが、2種類のデータを塗りつぶしのwith boxesで重ねると
多少は透過させているとはいえ、見にくい・・・・
単純にwith stepsで2つのデータをプロットすると、一見綺麗に書けるが
set log yとしてみると、データ点が入っていないbin(y軸の値が0になるbin)の処理がlogスケールなので、想像と違うことになってしまう・・・
↑
カウントが0のbinの端の処理が嫌・・・・
with boxesでもいいかなぁ〜と思っていたが解決方法に気づいたので、ここに書いとく・・・
まずはデータの用意
gsl-histogramで0から50の間を30個のbinで分けてヒストグラムのデータを生成する
特に問題なさそう・・・・
けど、もっとデータ点が増えてくると、y軸をlogスケールにしたくなる・・・
矢印を書いたあたりの側面の処理がなんか嫌・・・
with stepsは(x_1, y_1)の点から(x_2, y_2)の点に向かって垂直な直線を引くようなプロットらしい
それでそのうちのいずれかがy=0だとy軸がログスケールだと、直線が引かれない仕様らしい
思いついた解決方法は、三項演算子でy=0になっている点はy=0.0001のようなグラフに影響が出ないような小さい数字を代わりに当てはめてプロットする方法
ツイート
■ 過去記事 : 【gnuplot 6】gsl-histogramの結果を塗りつぶしてプロットしたい
(今気づいたけどgnuplot6になってますよ・・・おいおい・・・)
今回は、with stepsでこれをプロットしたい
今まで単一のデータをプロットするときは、塗りつぶしプロットでも用は足りていたが、2種類のデータを塗りつぶしのwith boxesで重ねると
多少は透過させているとはいえ、見にくい・・・・
単純にwith stepsで2つのデータをプロットすると、一見綺麗に書けるが
set log yとしてみると、データ点が入っていないbin(y軸の値が0になるbin)の処理がlogスケールなので、想像と違うことになってしまう・・・
↑
カウントが0のbinの端の処理が嫌・・・・
with boxesでもいいかなぁ〜と思っていたが解決方法に気づいたので、ここに書いとく・・・
まずはデータの用意
gsl-randist 0 1000 rayleigh 10 | gsl-histogram 0 50 30 > hoge.txt
gsl-randistで1000個の乱数データを生成して、gsl-histogramで0から50の間を30個のbinで分けてヒストグラムのデータを生成する
plot [:][0.1:] "hoge.txt" u 1:3 w steps lw 3 lc 2
でこれ特に問題なさそう・・・・
けど、もっとデータ点が増えてくると、y軸をlogスケールにしたくなる・・・
set log y
plot [:][0.1:] "hoge.txt" u 1:3 w steps lw 3 lc 2
plot [:][0.1:] "hoge.txt" u 1:3 w steps lw 3 lc 2
矢印を書いたあたりの側面の処理がなんか嫌・・・
with stepsは(x_1, y_1)の点から(x_2, y_2)の点に向かって垂直な直線を引くようなプロットらしい
それでそのうちのいずれかがy=0だとy軸がログスケールだと、直線が引かれない仕様らしい
思いついた解決方法は、三項演算子でy=0になっている点はy=0.0001のようなグラフに影響が出ないような小さい数字を代わりに当てはめてプロットする方法
set log y
plot [:][0.1:] "hoge.txt" u 1:($3 == 0 ? 0.001 : $3) w steps lw 3 lc 2
plot [:][0.1:] "hoge.txt" u 1:($3 == 0 ? 0.001 : $3) w steps lw 3 lc 2
PR
この記事にコメントする
プロフィール
HN:coffee
職業:物理屋(自称)
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
職業:物理屋
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
カテゴリー
ブログ内検索
リンク
相互リンク募集中です
(Twitterにてお知らせください)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
最新記事
(11/20)
(03/05)
(02/29)
(02/21)
(02/21)
(02/21)
(02/21)
(01/13)
(01/05)
(01/05)