Meditation_Summary.R

y1111<-y1000 <-y0000<-y1011<-y1101<-y1110<-y1001<-y1100<-y1010<-NULL

i<-0
for(i in 1:1500){
  ii <- i+500
  y1111[i] <- mean(Y_1111[[ii]])
  y1000[i] <- mean(Y_1000[[ii]])
  y1011[i] <- mean(Y_1011[[ii]])
  y1101[i] <- mean(Y_1101[[ii]])
  y1110[i] <- mean(Y_1110[[ii]])
  y1001[i] <- mean(Y_1001[[ii]])
  y1100[i] <- mean(Y_1100[[ii]])
  y1010[i] <- mean(Y_1010[[ii]])
  y0000[i] <- mean(Y_0000[[ii]])
  
}

quantile(y1111-y0000, c(0.025,0.975))
quantile(y1000-y0000, c(0.025,0.975))
quantile(y1111-y1000, c(0.025,0.975))
quantile(y1111-y1011, c(0.025,0.975))
quantile(y1111-y1101, c(0.025,0.975))
quantile(y1111-y1110, c(0.025,0.975))
quantile(y1111-y1001, c(0.025,0.975))
quantile(y1111-y1100, c(0.025,0.975))
quantile(y1111-y1010, c(0.025,0.975))

mean(y1111-y0000)
mean(y1000-y0000)
mean(y1111-y1000)
mean(y1111-y1011)
mean(y1111-y1101)
mean(y1111-y1110)
mean(y1111-y1001)
mean(y1111-y1100)
mean(y1111-y1010)


mean(y1111-y1011+y1111-y1101-y1111+y1001)
quantile(y1111-y1011+y1111-y1101-y1111+y1001, c(0.025,0.975))

mean(y1111-y1011+y1111-y1110-y1111+y1010)
quantile(y1111-y1011+y1111-y1110-y1111+y1010, c(0.025,0.975))

mean(y1111-y1101+y1111-y1110-y1111+y1100)
quantile(y1111-y1101+y1111-y1110-y1111+y1100, c(0.025,0.975))