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