data=scan("E:/FSU-FALL07-School/Fall07-school/STA5166/BHH2-Data/datahw3.dat") data Mean1wout = mean(c(8,6)) Mean2with = mean(c(3,4)) diff_means = Mean2with-Mean1wout y = c(rep(NA, (109))) x = c(rep(NA, (109))) n1=0 mean(data[1:2]) for(i in 1:109){ y[i] = mean(data[i:i+1]) } y for(j in 1:109){ x[j] = y[j+2] - y[j]; if(diff_means>= x[j])n1=n1+1 } x sort(x) diff_means n1 n1/109