題名のまんまの話です。 環境省大気汚染物質広域監視システム （通称：そらまめ君）では、日本全国の様々な地域の大気汚染物質について 情報の公開を行っていますが、過去のものが1週間くらいしか遡れない。 大抵こういうのって、しばらくしてから欲しくなっ…

データの９５%信頼区間が 場合によっては真値を含まないケースを図示するので 以前、どこかを参考に書いたはずの絵が見つからないので ggplot2で描いてみた。 もっと良い方法があるのかもしれないけど ggplot2の練習もかねて、これで一応、良しとする。 data…

GMTで書いたものをPSファイル以外にしたいと 思っていたら、ありました。 単なる勉強不足でしたなぁ。 psconvert map.ps -Tg -E250 -A psconvertを使えばいいとか。 -TgはPNGにするオプション。 -E250は解像度。 -Aは余白を極力削るオプション。 まだまだオ…

Rでのちょっと戸惑ったこと。 データからsubset()で検索値をとりだして その値を使ってさらに足し算したいとき。 検索結果がヒットしなかったら ０を返してほしかったのだけど そうじゃなくてlogical(0)なんて出てきて これが煮ても焼いても食えなかった。 i…

後方流跡線を描くとき つまづいたのが、geom_line と geom_pathの違い ずっとlineで描いていたけど、ある線が、おかしくなった。 （例のデータがどれだったのか、解らないので 例図がありません。そのうち見つけたら、更新） なんだか、ギザギザになって明ら…

後方流跡線のデータを取り出す方法も分かった。 地図も描いた。 なので、あとはこれを くっつけるだけ。 いろいろと捗るといいなぁと思って データのファイル名を引数に関数にしてみた。 ちなみに、後方流跡線はNOAAのものです。 気象データはGDASの０．５度…

後方流跡線のデータを読み込んだら あとはこっちのもの（なにが？） 地図を描いて、後方流跡線を乗っけるだけ。 今回、地図はNature Earthのものを使用しました。 荒めの国境データですが、今回にはこれで十分かと。 shapeファイルを読み込んで、ggplot2で描…

オフラインで得た後方流跡線のテキストデータ。 これまではExcelマクロ＆GMTで描いていたけど 様々な処理をRにまとめたくて Rで描くことに挑戦。 で、まず最初に戸惑ったのがデータの読み込み。 今まではExcelで読み込んでいたから問題なかったけど Rで行ご…