I have a matrix of reads counts for each gene for each strain of an animal. I am using DESeq2 to analyse the clustering of the strains. I have performed a Variance Stabilising Transform on my DESeq2 object and can do a PCA:
p <- DESeq2::plotPCA(vsd, intgroup=c('Condition'))
p <- p + labs(title = 'PCA on gene expression levels')
Now, I would like to know the coefficients/weights of each gene for each principal component. Does anyone know how to do this?