adding nearPD to compute the nearest Positive definate of the matrix phi
matrix phi is not positive definate only on one value with h2 = 1
nearPD is used in get_m1_Vbeta_j_mu