Rに関するmemo

グループ毎に度数分布を作成する

by(housing$VAL, housing$SERIALNO, length)

グループで平均値を計算するための式

DT[,mean(pwgtp15),by=SEX]             
sapply(split(DT$pwgtp15,DT$SEX),mean)             
mean(DT[DT$SEX==1,]$pwgtp15); mean(DT[DT$SEX==2,]$pwgtp15)             
mean(DT$pwgtp15,by=DT$SEX)             
tapply(DT$pwgtp15,DT$SEX,mean)             
rowMeans(DT)[DT$SEX==1]; rowMeans(DT)[DT$SEX==2]

NAがあるときに除外してカウントする方法

nrow(housing[housing$VAL==24 & !is.na(housing$VAL),])

NAについてはwhichを利用しても除外ができるらしい。