warehouse <- read.table("http://www.stat.ufl.edu/~winner/sta4210/data/CH04TA02.txt", header=F,col.names=c("X","Y")) attach(warehouse) # pdf("F:\\sta4210\\CH04TA02.pdf") (n <- length(Y)) (sum_XY <- sum(X*Y)) (sum_X2 <- sum(X^2)) (b1 <- sum_XY/sum_X2) (Y_hat <- b1*X) (e <- Y-Y_hat) (MSE <- sum(e^2)/(n-1)) (s2_b1 <- MSE/sum_X2) (s_b1 <- sqrt(s2_b1)) (b1_95CI <- c(b1-qt(.975,n-1)*s_b1,b1+qt(.975,n-1)*s_b1)) plot(X,Y,main="Warehouse regression Thru Origin",xlim=c(0,200), ylim=c(0,1000),xlab="Work Units Performed",ylab="Variable Labor Cost") abline(a=0,b=b1) warehouse.reg1 <- lm(Y ~ X -1) summary(warehouse.reg1) # dev.off()