extract_mean.R

#define variance
variance <- function(x) var(x)*(length(x)-1)/length(x)
options(digits=10)
tardir <- "./result1"
rfile <- paste(tardir,"/mlocval.txt",sep="")
dat <- read.table(rfile)
dat <- data.matrix(dat)
h <- dat[,1]
h <- unique(h)
v <- dat[,4]
boxelem <- length(h)
alltime <- length(v) / boxelem
vmean <- numeric(boxelem)
vdev  <- numeric(boxelem)
btime <- 9
mask <- (btime : alltime-1) * boxelem + 1
for (i in 1:boxelem) {
    vmean[i] <- mean(v[mask])
    vdev[i] <- variance(v[mask])
    mask <- mask + 1
}
sample <- alltime - btime
wfile <- paste(tardir,"/vxmean.txt",sep="")
vdev <- vdev / sqrt(sample)
wdat <- cbind(h,vmean,vdev)
write(t(wdat), wfile, append=F, ncolumns=3)