For the follwoing data set
ndat=data.frame(id=c(2,2,3,3,3,3,4,4,5,5,5,5,6,6,6), vc=c("bl",\'m0\', "bl","m0",