Functional principal component analysis under the Linear Mixed Models representation of smoothing splines. The method utilizes the Demmler-Reinsch basis and assumes error independence. For more details see: F. Rosales (2016) < https://ediss.uni-goettingen.de/handle/11858/00-1735-0000-0028-87F9-6>.