This is a note to my second question. I just discovered my data containing an 'Inf' value and this causes the error I saw when I used rcorr() function from the Hmisc package. If I delete this observation, everything is fine. Best, Ming-Chung