table=read.delim("labor_force.txt") boxplot(table[,2:3],col=c("red","yellow"),ylab="Women labor participation rate",main="Women in the Labor Force") attach(table) stem(X1972) stem(X1968) #Pooled t.test t.test(X1972,X1968,alternative=c("greater"),var.equal=TRUE) #Matched pairs test t.test(X1972,X1968,alternative=c("greater"),paired=T)