昔々，『プログラム書法』と『ソフトウェア作法』という本があったそうな。異国の著者 Brian W.Kernighan と P.J.Plauger の著作を 木村 泉 氏が翻訳したものであったそうな。「作法」は「作り方の法則」が本来の意図であろうが，木村氏はそれを「お作法（さ…

「人柄ではなく『コード柄』あなたは大丈夫ですか？」に書いてあることだけど， > この人は美しいコードを書くとか，行き届いたコードを書くとか。逆に，乱雑なコードを書く人もいれば，センスのないコードを書く人もいる…と。 > 「コード柄」は見る人が見れ…

http://gihyo.jp/dev/serial/01/codeiq/0007「第7回 t検定による問題解決，Rで実践できますか？～データサイエンティストの統計学─倉橋一成からの問題」に書かれていることなんだけど，> A群とB群⇒ 平均値が170cmの集団> C群⇒ 平均値が175cmの集団> > # A群…

「二群の平均値の差の検定」に関するシミュレーションを行う shiny アプリ http://glimmer.rstudio.com/rder/difference-of-two-means/ のように指定する（指定されたディレクトリに server.R, ui.R を置く） 曲線は母集団の密度関数 信頼区間は，母平均の 9…

「Rでlmコマンドの結果からStd. Errorを取得する」に書いてあることについて 回帰分析の係数の区間推定を求めたいときがあり，そのためには，回帰係数の標準誤差を求める必要があるということですが。 その後，自分で計算式を書いて求めているのだけど， 回…

「MeanDecreaseGiniが大きい変数の箱ひげ図を描く」の記事についてだけど... ベクトル要素の大きい順に 1, 2, ... と順序づけを行いたかったときに，「order 関数を二重に使えばよい」と twitter で教えてもらったということのようだ。 > x <- c(5,2,1,6,4,7…

「reshape2 でデータ整形」に書いてあることについてなのだが library(reshape2)aq.m <- melt(airquality, id=c("Month", "Day"), na.rm=TRUE)dcast(aq.m, Month ~ variable, mean) によって Month Ozone Solar.R Wind Temp1 5 23.61538 181.2963 11.622581 …

「全ての変数の分布を箱ひげ図にする」に書かれていることについて 提示されているプログラムは，作者本人も認めていること以外にもいろいろ不都合がある 道具は，とにかく「汎用的でなくてはならない」 そのためには，特定の状況を仮定したものではダメ。例…