performed LCR on chr 4 at marker class1 <- which(model.LCR@cluster == 1)