グループ毎に度数分布を作成する
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を利用しても除外ができるらしい。