taicrim <- rep(c(160,162.6,165.1,167.4,170.2,172.7,177.8,185.5,195.6), c(2,2,1,3,2,2,2,1,1))
DPB <- c(3.7,3.2,3.1,4.3,3.8,3.4,3.8,3.6,3.9,2.7,2.9,3.1)
DRE <- c(3.9,5,4.2,4.2,4,3.4,4.8,4.4,4.2,5,4.3,4.5,4.8,3.4)