|
25
|
1 args = commandArgs(trailingOnly=TRUE)
|
|
|
2
|
|
|
3 library("amap")
|
|
|
4 library("labdsv")
|
|
|
5
|
|
|
6 data <- read.table("plink.raw",header=T)
|
|
28
|
7 rownames(data) <- data[,2]
|
|
26
|
8 distances <- Dist(data[,7:dim(data)[2]], method = "euclidean", nbproc = as.integer(args[2]))
|
|
25
|
9 pco_results <- pco(distances,k=10)
|
|
|
10 #Variance explained by first three PCOs
|
|
|
11 pco_results$eig[1]/sum(pco_results$eig[pco_results$eig>0])
|
|
|
12 #[1] 0.2565937
|
|
|
13 pco_results$eig[2]/sum(pco_results$eig[pco_results$eig>0])
|
|
|
14 #[1] 0.06878127
|
|
|
15 pco_results$eig[3]/sum(pco_results$eig[pco_results$eig>0])
|
|
|
16 #[1] 0.04340111
|
|
|
17 write.csv(pco_results$points,file=args[1])
|
|
|
18
|
