# metabolic turnover data
data(turnover)
X <- turnover$X
D <- turnover$D
ospca <- os_pca(X,D,0.999)
ospca <- ospca_loading(ospca)
ospca$loading$R
ospca$loading$p.value
# metabolome data
data(greentea)
X <- greentea$X
D <- greentea$D
M <- greentea$M
ospca <- os_pca(X,D,0.1,M)
ospca <- ospca_loading(ospca)
ospca$loading$R
ospca$loading$p.value
Run the code above in your browser using DataLab