Getting loadings of components from pcaMethods

I want to get the loadings of components from pcaMethods.

To illustrate the problem, I have a matrix x

    a   b   c
1   ... ... ...
2   ... ... ...
3   ... ... ...

By loadings(pca(x, method="nipals", nPcs=2)), I get

    PC1 PC2
a   ... ...
b   ... ...
c   ... ...

However, I want to get

    PC1 PC2
1   ... ...
2   ... ...
3   ... ...

Anyone knows how to do that?

Thanks a lot!

pcamethods pca

• 1,477 views • link

updated 7.8 years ago by Philipp Bayer

• written 7.8 years ago by johnnytam100

1 answer

Wouldn't you just need to transpose your input table while keeping names intact?

tx = setNames(data.frame(t(x[,-1])), x[,1])
loadings(pca(tx, method="nipals", nPcs=2))

• 1 views • link

written 7.8 years ago by Philipp Bayer

Good idea! I had a similar thought but does it work for a PCA analysis that we just swap the columns and rows?

• 1 views • link

written 7.8 years ago by johnnytam100

More posts like this

remove the letter a from a biplot in ggplot
written 1.9 years ago by Assa Yeroslaviz

I'm trying to create a biplot in ggplot2 and manage to do it after a lot of trail and errors, But I can't get rid …
How can I make all my PCA plots the same shape using ggplots?
written 3.3 years ago by mgranada3 •

I am using DESeq2 for gene expression analysis. I have 3 different conditions which I am normalizing to my control so I have produced 3 …
Extraction of high variance or contributing genes using pcatools
written 4.2 years ago by mohammedtoufiq91

Hi, I am working with the qRT-PCR log2FC data (96 gene on rows * samples on columns) containing healthy controls and patients treated with different …
Loading PCA - different value depending on the package used
written 5.3 years ago by camillab. •

Hi, I want to find the variable contribution (loadings) after I compute my PCA but I noticed that depending on which package I used they …
Plot biplot (pca loadings) onto PCA plot from DESeq2?
written 5.7 years ago by NBR •

Very curious how I can modify my code below (which I got from this forum :) ), such that I can get the loadings biplot …
Differences between PCA plots in DESeq and FactoMineR
written 5.9 years ago by caitlin.mcandry •

Hello! I'm creating PCA plots with a whole transcriptome dataset, which includes 26,828 'genes', to assess similarity between three groups (E, M and L). I'm …
change axis scale on ggplot2
written 6.0 years ago by niyomiw88 •

Hi, I am trying to plot my pca in ggplot and I've manage to get something (picture attached). My question is, How do I spread …
How to check how much loading each SNP had on the first or second principal component
written 6.7 years ago by Hann

Hello, I did PCA analysis using pca function from the LEA r package It produces these files: > eigenvalue file: genotypes_cultivated.eigenvalues > eigenvector file: genotypes_cultivated.eigenvectors …
Adding Percentage in PCA
written 8.7 years ago by 1769mkc

This is my code for PCA using SVD , i get a neat plot , I want to add percatge to the axis I used …
color 2 hclust groups on pc1, pc2 plot
written 9.3 years ago by mms140130

Hello I have performed a PCA analysis on gene expression data 1096 patients and 20248 genes and plotted the pc1 vs. pc2 then using hclust …