psvg <- data.frame(row.names=names(ausw.liste))
psvg[1:(35+8)] <- NA
pseudonull <- function (i,qrea) {
c(
mean(qrea$hl[qrea$tgl==i&qrea$tgl2==i&qrea$t>n2],trim=0.2,na.rm=T)-mean(qrea$hr[qrea$tgl==i+7&qrea$tgl2==i+7&qrea$t>n2],trim=0.2,na.rm=T),
mean(qrea$hl[qrea$tgl==i&qrea$tgl2==i&qrea$t<n2],trim=0.2,na.rm=T)-mean(qrea$hr[qrea$tgl==i+7&qrea$tgl2==i+7&qrea$t<n2],trim=0.2,na.rm=T),
mean(qrea$hl[qrea$tgl==i&qrea$tgl2==i&qrea$t>n2],trim=0.2,na.rm=T)-mean(qrea$hr[qrea$tgl==i+7&qrea$tgl2==i+7&qrea$t<n2],trim=0.2,na.rm=T),
mean(qrea$hl[qrea$tgl==i&qrea$tgl2==i&qrea$t<n2],trim=0.2,na.rm=T)-mean(qrea$hr[qrea$tgl==i+7&qrea$tgl2==i+7&qrea$t>n2],trim=0.2,na.rm=T),
mean(qrea$hl[qrea$tgl==i&qrea$tgl2==i],trim=0.2,na.rm=T)-mean(qrea$hr[qrea$tgl==i+7&qrea$tgl2==i+7],trim=0.2,na.rm=T)
)
}
fdx <- function (qrea) {
c(
mean(yrea$vg,trim=0.2,na.rm=T),
mean(qrea$vg[qrea$tgl!=0&qrea$tgl<=7&qrea$t<n/2],trim=0.2,na.rm=T),
mean(qrea$vg[qrea$tgl>7&qrea$tgl<=14&qrea$t<n/2],trim=0.2,na.rm=T),
mean(qrea$vg[qrea$tgl!=0&qrea$tgl<=7&qrea$t>n/2],trim=0.2,na.rm=T),
mean(qrea$vg[qrea$tgl>7&qrea$tgl<=14&qrea$t>n/2],trim=0.2,na.rm=T),
mean(yrea$vg[abs(yrea$hl+5)<0.5],trim=0.2,na.rm=T),
mean(yrea$vg[abs(yrea$hl  )<0.5],trim=0.2,na.rm=T),
mean(yrea$vg[abs(yrea$hl-5)<0.5],trim=0.2,na.rm=T)
)
}