Dear all, In SAS, we can use "proc corr data=a pearson spearman; var a; partial b c d; with e; run;" to get the partial pearson and spearman correlation coefficient and p value of a with e after adjustment b,c and d. Are there any built-in functions available in any packages in R to get the partial pearson/spearman correlation coefficient and p value now? Thank you. Xingwang Ye PhD candidate Research Group of Nutrition Related Cancers and Other Chronic Diseases Institute for Nutritional Sciences, Shanghai Institutes of Biological Sciences, Chinese Academy of Sciences P.O.Box 32 294 Taiyuan Road Shanghai 200031 P.R.CHINA