Mercurial > repos > cropgeeks > ukseed

diff dartseq_seeduk_2.R @ 4:c9686a670e60 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author cropgeeks
date Tue, 17 Apr 2018 09:10:42 -0400
parents a8f30e5859c9
children
line wrap: on
line diff
--- a/dartseq_seeduk_2.R	Tue Apr 17 08:28:05 2018 -0400
+++ b/dartseq_seeduk_2.R	Tue Apr 17 09:10:42 2018 -0400
@@ -5,7 +5,7 @@
 
 data <- read.table("plink.raw")
 distances <- Dist(data[2:338,7:75189], method = "euclidean", nbproc = 144)
-pco_results <- cmdscale(distances,k=10)
+pco_results <- pco(distances,k=10)
 #Variance explained by first three PCOs
 pco_results$eig[1]/sum(pco_results$eig[pco_results$eig>0])
 #[1] 0.2565937