For the benefit of other users that could encounter the same problem:
It was indeed connected to the fact of being the dataset so small.
In the default values of
limmaquickpca2go, I use a higher values for the number of genes to use in the PCA and to use for the highest loadings then the actual number of genes included in the dataset Reuben was using.
In your case, to avoid this (and another couple of issues related to
limma::goana), you could try forcing
rld_reu the rlog transformed object with the input Reuben was using.
To avoid the problem at all during runtime, I can recommend you to compute the
pca2go object in advance, and then pass it to the function call
# dds_reu and rld_reu are created in advance
pcaExplorer(dds_reu, rld_reu, pca2go = pca2go_reu)