Laser-induced breakdown spectroscopy (LIBS) is being proposed more and more as a high-throughput technology to assess the elemental composition of materials. When a specific element is of interest, semiquantification is possible by building a calibration model using the emission line intensity of this element for known samples. However, a unique element has usually more than one emission line, and there are many examples where several emission lines used in combination give dramatically better results than any of the individual variables used alone. With a multivariate approach, models can be constructed that take into account all the emission lines related to a specific element; therefore more robust models can be developed. In this work, chemometric methods such as principal component analysis and partial least squares are proposed to resolve and extract useful information from the LIBS spectral data collected on biological materials.