various changes, coef<- generic, moved set contrast to inside getX
added simulate method to simulate parameters from their MV normal distribution