R

パイプ

R

library(dplyr) # ←この時点でわかりにくい。これで合ってるかな。なんの略だこれ。 x <- rnorm (100) y <- rnorm (100) lm (y ~ x) %>% summary わー便利。さて、lm (y ~ x) の中身はと・・・中身はと・・・ 結局コードしなおしじゃん。 パイプの中をごちゃ…

はじめてのべいずせいこう

R

かけんひのさいたくりつ(p)の事後分布 (by Rstan)。Stanでは事前分布は一様分布になるそうだ。ここでは デフォルトの0 to 1。 うけるわ。なんやこれ。一生無理やん(右端は1じゃないからね)。生々しいなー、こうするとなんか。 ある意味ベイズ統計のちから…

 覚え:切片だけのロジスティック回帰

R

(1つの飼育容器に10匹入れて選好性試験:AとBのどっちに行くか)×20反復 母比率=0.5を検定したい 切片だけのロジスティック回帰?(もちろん各個体は他の個体と独立に選ぶという前提で) # ロジスティック関数・切片だけモデル・シミュレーション# 関数定…

R R

昨日の会議が5時間半もあったのでかなり疲れが残っている. 日を超えるのはなんとか免れた。 そういうわけで今日の午前はゆるゆるやる。夜の部もあるし。 図はこんなもんか。 x軸をあと書きするとき plot(..., xaxis = F); axis(1,at = 1:6, labels = 文字ベ…

 FAQ

R

はい、そうです。 たった2群でも分散分析できるんです。 等分散を仮定した場合の2群の分散分析はシンプルなt検定と同じ結果になるんです。 卒論時期には必ず一回は来る質問。 「2群の時はt検定、3群以上のときは分散分析」 という教条化した、硬直化し…

R R

条件にあう要素だけを置換する alpha[which(alpha > 8)]

R R

lev lev levels(dat$host) 以下にエラー `levelsvalue = c("hati", "kene", "meda", : number of levels differs 消えないな。

 削除して、もはや存在しない level を消す

R

除去したはずのlevelがまだ有効・・・不思議な仕様だ。 horai, susukiが生き残っている。なぜだろう?追記: http://www.okada.jp.org/RWiki/?Tips%2F%BB%B3%DC%A5Tips%BD%B8#a8c9e343factor(x)で消えるとな

R R

NA行を除去dat hostがAでもBでもないものを抜く dat[dat$host != "horai" & dat$host != "susuki",]

 軸目盛の変更lab

R

なんかいいところだけど、整数軸にしたい・・・こういうちょっとしたことを直すのは簡単ではない。 plot( ... ,lab=c(4,7,7)) # x, y, サイズ むずい!追記: plot(..., yaxp = c(0, 3, 3)) # 下限、上限、目盛数

R R

ラベル軸の大きさとか変える方法とか忘却のかなた。 なんでデフォがこんなに小さいのか、馬鹿なんじゃないのか。 cexがラベルに効かない。おかしいね。

R R

プロットの並べ換えってどうやるんだっけ・・・?もうね、毎度毎度。 っていうか、この程度のことも簡単にできないRはやはり市民権を得ることはできないのではないだろうか。

 防備録

R

NAを除いたベクトルの要素を取り出す x[!is.na(x)] #横方向へのデータから差分をとる xx for (i in 3:15){ xx } xx # タテにならべてデータフレームを作る n host ind day fec dat2 平均と標準偏差プロット library(gplots) plotmeans(y~x)

R R

一応できた。おやつのあんみつがよかったか。 trt d n m rvector(t(m)) trt1 response dat 071003追記(こっちのほうがいいか) # 別解 trt d n rvector(t(cbind(d,n))) trt1 response dat

R

Rで二値応答データ treatment diapause nondiapause A 0 3 B 1 3 C 3 1 を trt response A 1 A 1 A 1 B 0 B 1 B 1 B 1 C 0 C 0 C 0 C 1 みたいにする方法。前に自分でやったか人に教えてもらったことがあるような気がするが completely に忘れてしまった。ご…