In this case study, we use the gene expression data collected from the mammalian eyes (Scheetz et al., 2006) and investigate their correlation to Bardet-Biedl syndrome. The dataset is a matrix which contains n = 120 rows and 11 columns, Y1 I 11 T 12 I 1,10 Y2 221 I 22 12.10 yn Inl In2 In 10 with the ith row representing a vector (Yi, Ni1, X;2, ..., , /i10) from subject i, where 5 Yil, ..., 2:10: the expression levels of 10 different genes collected from subject i yi: a severity measure of Bardet-Biedl syndrome for subject i The the 1st column of the matrix y = (41, ..., Yn)T corresponds to the sever- ity measures of Bardet-Biedl syndrome from subject 1-120; the 2nd column gi = (411, ... , 1) corresponds to the expression levels of Genel from subject 1-120, and the 11th column corresponds to 910 = (21,10, . . . , In10)". First we load genetic data using the following command ## ---- Load-the-dataset data.source