view fosm_cluster/kmeans.R @ 0:68a3648c7d91 draft default tip

Uploaded
author matteoc
date Thu, 22 Dec 2016 04:45:31 -0500
parents
children
line wrap: on
line source

#!/usr/bin/env Rscript
args = commandArgs(trailingOnly=TRUE);
data<-read.table(args[1],sep="\t",row.names=1);
data<-data[,1:258];
data[,2]<-data[,2]/sum(data[,2])
file<-(args[3])
num=as.integer(args[2])
K<-kmeans(data,num,nstart=1500)
clust=K$cluster

write.table(clust,file);