A_c(13.2,8.2,10.9,14.3,10.7,6.6,9.5,10.8,8.8,13.3) B_c(14,8.8,11.2,14.2,11.8,6.4,9.8,11.3,9.3,13.6) shoes_data.frame(A,B) shoes t.test(shoes$A, mu=10) attach(shoes) t.test(A, mu=10) t.test(A)$conf.int wilcox.test(A,mu=10) help(wilcox.test) var.test(A,B) help(var.test) t.test(A,B) help(t.test) t.test(A,B,paired=T) t.test(A,B,var.equal=F) wilcox.test(A,B) wilcox.test(A,B,paired=T) eda.shape_function(x) { par(mfrow=c(2,2)) hist(x) boxplot(x) iqd_summary(x)[5]-summary(x)[2] plot(density(x,width=2*iqd),xlab="x",ylab="",type="l") qqnorm(x) qqline(x) par(mfrow=c(1,1)) invisible() } eda.shape(B-A)