facetのstripのtextを変更する

R

ggplot2でfacet_gridやfacet_wrapを使ってグラフを分割描画した際に、各グラフについているラベル(strip)のテキストは、何もしなければ分割に用いたグループ名が自動的に入る。 例: ggplot(iris, aes(x = Petal.Length, y = Petal.Width)) + geom_point()…

飽和水蒸気圧の求め方

何年か前に飽和水蒸気圧の求め方に関して、Goff-Gratchの式やTetensの式やMurrayの式について少し触れたことがあった*1のだが、最近Wikipediaを見ていたらTetensの式だと思っていたものが Tetens(1930)のパラメータ値によるAugust他の式 飽和水蒸気量 - Wiki…

ThinkPad X61 SSD化+その他いろいろ

2万ほど投資したらまだまだ使えそうな感じになった。 用意したもの SSD(A-DATA ASP900S3-256GM-C-7MM ADATA 2.5"SSD 256GB SATA6G A-DATA ASP900S3-256GM-C-7MM) 2.5インチHDケース(【日本正規代理店】 ORICO 2.5インチ HDD/SSD 外付け ドライブ ケース S…

ggplot2で任意の線分

R

baseライブラリ中のlines()に対応するggplot2のレイヤー関数はgeom_line()であると説明される場合が多く、大抵の場合はそれで事足りるのだが、任意の位置に線分を描き込みたいという場合はあまり使い勝手が良くない。 例えば、「統計学:Rを用いた入門書(Mi…

分布の集中度を調べる

R

ランダムか、集中か 水田の調査圃場において、正方形に区切った調査区画毎にニカメイガの卵塊が次のように見つかったとする。 egg <- c(2,4,1,1, 4,0,3,3, 4,3,0,1, 3,2,1,1, 2,1,4,3, 2,2,1,0, 1,0,5,1, 2,1,3,2, 2,2,2,2, 4,1,1,2, 3,2,0,1, 1,2,1,0) もし…

Rで逆推定

R

JMPには逆推定というコマンドがある 次のようなデータを用意して、線形回帰を行うとする。 temp <- rep(c(15, 20, 25), c(5, 5, 5)) day <- c(34, 33, 36, 37, 35, 30, 28, 29, 25, 28, 23, 25, 22, 24, 26) こんなかんじで適当に。 > result <- lm(day~temp…

不審なファイルの存在確認等

事の発端 【注意喚起】潜伏しているかもしれないウイルスの感染検査を今すぐ!:IPA 独立行政法人 情報処理推進機構www.ipa.go.jp 「バッチでの配布でも可能です」とか書くならバッチファイル用意しといてほしい…。 バッチファイル 上からやれといわれたけど…

ggplot2で縦に並べたグラフの横幅を揃える

R

(※2015/06/03 20:27追記あり) 下記の記事に基づく発言に関連してどうも某所で勘違いが発生しているようなので。R Recipe: Aligning Axes in ggplot2 | Exegetic Analyticswww.exegetic.biz 使用データ 下記のように生成したものを用いる。 x <- 0:100 x <-…

RStudioの折りたたみ機能

R

今更ながらコードの折りたたみ機能が付いている事に気付いた。 詳細はCode Folding and Sections – RStudio Supportを参照。 折りたたみ設定されるもの ブレース(波括弧{})で括られた領域。 R SweaveやR Markdownドキュメントにおけるコードチャンク。 R M…

S3

R

RStudioではじめるRプログラミング入門を読んだ。RStudioについて詳しい情報があるかと思っていたのだが、これについてはそこそこの記述だった。 ただ、全体的に見て、Rとプログラミングをゼロから学ぼうと思うのであれば良い情報量だと感じた。リゲス本なん…

S-Insertでペーストできない

.vimperatorrc中でペースト操作をしたいところに今までS-Insertと書いていたのだけど、気付いたら動作しなくなっていた。 今まではVimperator側でペースト操作に上書きしていたが(vimperatorでのコマンドラインへの貼り付け方法が分からなかった, S-Fマガジ…

キュウリに対して、おそらく同時期に放飼したであろうスワルスキーカブリダニが容易に見つかるほ場と30分探しても見つからないほ場があった。 スワルスキーカブリダニはナスやキュウリの害虫を食べるダニで、生き物だが農薬として登録されている。生き物だろ…

3週連続、洗車の翌日に雨が降っている。雨乞いに使えるのではないかというレベル。

紙の日記、全く長続きしないのでやはりここに書こうと思う。 朝、今日薬剤散布をするのであれば今か夕方かという相談があった。この時期の薬剤散布は難しい。だいたい、農薬というのは散布後に乾かなかったり、散布時に温度が高すぎたりすると薬害の危険性が…

gap.plot()で省略される部分のラインを途中まで描く

R

plotrixパッケージに含まれるgap.plot()関数は、データの途中を省略したプロットを行う関数で、以前に少し紹介した。 ギャップのあるデータを途中を省略してプロットする - もうカツ丼はいいよな このとき、typeに"o"、"l"、"b"などを指定した場合、省略範囲…

Majestouch MINILA買った&また改造した

職場で使う用にMINILAの黒軸を購入したので、前回(Majestouch MINILA買った&Mac用にとりあえず改造した - もうカツ丼はいいよな)と同じような感じで少々いじった。 基本的には前回と同様で、Fnキーへの変換・無変換割り当てとCapsLockへのFnキー割り当てを…

Rでinstall.packages()しようとしたらエラー

R

すごく久しぶりにRでinstall.packages()してみたら、 > install.packages("ggplot2") --- Please select a CRAN mirror for use in this session --- Error: .onLoad failed in loadNamespace() for 'tcltk', details: call: dyn.load(file, DLLpath = DLLpa…

TypeError:['@@iterator'] is not a function

Firefox34.0.5にした(された)ら見出しのようなエラーが出てVimperatorがうんともすんとも言わなくなってしまって大分焦ったけど、_libly.jsを新しいものにしたら直った。 fix "TypeError: ['@@iterator'] is not a function" on Firefox 34 · 0d4ec6d · vi…

ギャップのあるデータを途中を省略してプロットする

R

ギャップのあるデータ 例えばこのようなデータが有るとする。 testdata <- data.frame(x = 1:10, y = c(1,3,2,4,3,21,23,22,24,23)) 前半と後半の間にデータの大きな飛躍があるため、そのままプロットすると前半部、後半部の変化が相対的に小さくなって見難…

Majestouch MINILA買った&Mac用にとりあえず改造した

なんかMacでも使えるらしいという曖昧な情報を元にMINILAを買ってしまった。スペースキーの両隣にFnキーが付いてるのが若干気になったけど、まあKarabinerとかでなんとかなるだろうと思っていた。 … どうにもならなかった…。なんか単独ではイベントが発生し…

Excelでエラー値を表示させなくする方法

Excelで印刷用の資料を作成していると(そもそもするなという話だが)、VLOOKUPなんかでエラーが帰ってきた場合にそれを表示したくないという場面がしばしばある。その解決方法3つ。 1. 条件付き書式を使う Excel for Mac 2011の場合だがスタイルを「クラシ…

害虫の薬剤抵抗性

害虫の薬剤抵抗性 ある害虫に対し特定の農薬を使い続けていると農薬の効果が低下してくる現象があり、薬剤抵抗性と呼ばれている。その歴史は古く、1908年にアメリカで柑橘の害虫であるナシマルカイガラムシの石灰硫黄合剤に対する感受性が低下したのが最初の…

Elpy使用時のFlymakeエラー

いつからか定かでないがElpy使用時に Flymake: Failed to launch syntax check process 'pyflakes' with args (/path/to/file/filename_flymake.py): Searching for program: no such file or directory, pyflakes. Flymake will be switch OFFとかいうエラ…

OS X Yosemiteアップデート後にやったこと

XCode6.1インストール まだApp Storeにアップデートが来てないが、brew doctorすると怒られるので、Xcode - ダウンロード - Apple Developerの「その他のツール」からDownloads for Apple Developersに入って6.1を落としてきてインストール。一度起動して色…

最も出現頻度の高い項目を表示する

例えば、 {a, b, b, c, c, c, d, d, e, e, e, e} というようなデータがあったとして、eを表示したい。これをExcelでやりたい。一時計算用のセルなどを使わずに。 最初は簡単そうに思えたのだが、ちょっと試行錯誤した程度では思いつかなかったので検索。結果…

Excelの空白セルと空文字列

「空のセル」 Excelにおいて空のセルをどう処理すれば良いのかというのは難しい問題である。そもそも「空のセル」と言った場合に、それがそもそも何を意味しているのかという点からして問題がある。 見た目が空白になるセルの例として、代表的なものは以下の…

トマト・ナス・キュウリの世界における単収推移

ちょっと気になったので調べてみた。データは主にFAOSTATから。 トマト 国は日本の他にオランダとスペインと中国を選んだ。深い理由は無い。浅い理由ならある。 オランダ:施設園芸がスゴイ(単収的な意味で) スペイン:施設園芸がスゴイ(施設一箇所に集ま…

ユークリッドの互除法

tex記法のテスト。 割られる数を、割る数をとする()。 商を、余りをとする()。 また、としておく。 このとき、 …(1) であり、 である。また、、に対し、、はただ一つ定まる。 、の共通の約数のうち最大のものを最大公約数と呼び、 であらわす。なお、と…

海陸風

昨日は大学の後輩の温室を見せてもらいに田原まで出かけていた。 ひと通り色々案内してもらって最後にサーファー見学に海を見ていた。 珍しく休日に出歩いた結果の軽い疲労感もあり海からの風がやけに涼しく感じ……。海風だ(゚д゚)! 海陸風(かいりくふう) 海…

ベンズイミダゾール系殺菌剤と薬剤耐性菌

ベンズイミダゾール系殺菌剤 ベンズイミダゾール系殺菌剤という言葉をそのまま解釈するとベンズイミダゾール環をふくむ殺菌剤ということで、植物に対する病害を防除する農薬のうち、該当するのはベノミル(商品名:ベンレート)なのだが、実際にはチオファネ…