data1<-"http://kanggc.iptime.org/book/data/chap10-2.csv" data1_dat<-as.matrix(read.csv(data1,header=T), ncol=1) var1<-data1_dat[,1] xbar=mean(var1) z<-qnorm(0.025, 0, 1, lower.tail=F) LCL1<-xbar-z*(0.6/sqrt(35)) UCL1<-xbar+z*(0.6/sqrt(35)) LCL1;UCL1 LCL2<-0.71643-z*(0.6/sqrt(35)) UCL2<-0.71643+z*(0.6/sqrt(35)) LCL2;UCL2