I'm using package pls to perform partial least squares regression. Data is four spectra (13178 intensity measurements each) that correspond to 4 concentrations. out <- plsr(concentration ~ intensity, ncomp=2, data=libs, validation = "LOO") Does the package have a command to calculate limits of detection? I have seen conflicting descriptions on how to do this in the literature. Thanks, Alonso. [[alternative HTML version deleted]]