body3k <- read.table(".../body3k.txt",header=T,sep="\t",quote="")

#TeildatensŠtze vorbereiten
table(body3k$Geschlechtk)
Mbod<-subset(body3k,Geschlechtk=="m")
Fbod<-subset(body3k,Geschlechtk=="w")

#Histogramme vergleichen
dev.new(height=500)
par(mfrow=c(2,1))
hist(Mbod$Hueftumfang,col="blue")
hist(Fbod$Hueftumfang,col="red")
dev.new(height=500)
par(mfrow=c(2,1))
hist(Mbod$Hueftumfang,col="blue",xlim=c(70,140),breaks=seq(70,140,5),ylim=c(0,80))
hist(Fbod$Hueftumfang,col="red",xlim=c(70,140),breaks=seq(70,140,5),ylim=c(0,80))

#Mittelwerte vergleichen
ms<-by(body3k[,1:24],body3k$Geschlechtk,mean)
mst<-cbind(ms$m,ms$w)
mst
ratio<-ms$m/ms$w
barplot(ratio-1)

#Mittelwerte plotten
dev.new()
plot(mst)
abline(a=0,b=1)
dev.new()
plot(mst,xlim=c(0,180),ylim=c(0,180),pch=19)
abline(a=0,b=1)

#Mediane vergleichen
med<-by(body3k[,1:24],body3k$Geschlechtk,summary)
msm<-cbind(med$m[3,1:24],med$w[3,1:24])