# HG changeset patch
# User cropgeeks
# Date 1523970642 14400
# Node ID c9686a670e6037f84f6c9bfe604e51536bc01581
# Parent  dde8fec1b31ef6998a3807bb63b0db9aaaae6398
Uploaded

diff -r dde8fec1b31e -r c9686a670e60 dartseq_seeduk_2.R
--- 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