I am doing a microRNA-mRNA paired expression analysis and because of collinearities between the miRNAs covariates and large number of microRNAs (in comparison to the sample size), ordinary linear regression model is not applicable. Therefor I thought about partial least square (PLS) regression methods. Has anybody performed such an analysis on miR-mRNA data? Would be grateful if anybody recommend some R/Python scrip?