foram.data=read.table("foram.txt",header=TRUE) attach(foram.data) boxplot(length~period,range=0,names=c("K","T"),ylab="Length",boxwex=.5) detach(foram.data)