diff --git a/.classpath b/.classpath new file mode 100644 index 0000000..d5cae95 --- /dev/null +++ b/.classpath @@ -0,0 +1,27 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/.project b/.project new file mode 100644 index 0000000..c012ce4 --- /dev/null +++ b/.project @@ -0,0 +1,17 @@ + + + japsa + + + + + + org.eclipse.jdt.core.javabuilder + + + + + + org.eclipse.jdt.core.javanature + + diff --git a/bin/.gitignore b/bin/.gitignore new file mode 100644 index 0000000..efaee1f --- /dev/null +++ b/bin/.gitignore @@ -0,0 +1,2 @@ +/japsa/ +/japsadev/ diff --git a/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering.java b/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering.java index ddd55f9..d7767b8 100644 --- a/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering.java +++ b/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering.java @@ -15,7 +15,7 @@ public class KmeanClustering { public static void Clustering() throws Exception{ - FileInputStream file1 = new FileInputStream("TRfile.fasta"); + FileInputStream file1 = new FileInputStream("TRfiletemp.fasta"); BufferedReader br = new BufferedReader(new InputStreamReader(file1)); String line = null; while((line = br.readLine())!=null){ diff --git a/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering1.java b/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering1.java deleted file mode 100644 index 6e39a5a..0000000 --- a/src/dev/java/japsadev/bio/hts/clustering/KmeanClustering1.java +++ /dev/null @@ -1,300 +0,0 @@ -package japsadev.bio.hts.clustering; - -import java.util.ArrayList; -import java.io.*; - -import japsa.seq.SequenceOutputStream; -import japsadev.bio.hts.clustering.PairDistance; -import japsadev.bio.hts.clustering.GettingTreadsFromFasta; - -/** - * @author buvan.suji - * - */ - -public class KmeanClustering1 { - - public static void Clustering() throws Exception{ - FileInputStream file1 = new FileInputStream("TRfile1.fasta"); - BufferedReader br = new BufferedReader(new InputStreamReader(file1)); - String line = null; - while((line = br.readLine())!=null){ - ArrayList descript = new ArrayList(); - ArrayList reads = new ArrayList(); - //ArrayList MaxReadLn = new ArrayList(); - int MaxReadLn = 0; - int Nreads; - String FileName; - int NumberElements; - int ClustElements; - int n; - int d[]; - int k[][]; - final int p =2;//number of clusters - int tempk[][]; - int m[]; - int Nclusters = 0; - - double max = 0; - double temp1 = 0; - int index1=0; - int index2=0; - int count1=0,count2=0; - long startTime = System.nanoTime(); - - - GettingTreadsFromFasta.DestReads(line); - //GettingTreadsFromFasta.DestReads(); - descript = GettingTreadsFromFasta.GetRname(); - reads = GettingTreadsFromFasta.GetTReads(); - Nreads = GettingTreadsFromFasta.NumberReads(); - FileName = GettingTreadsFromFasta.GetFileName(); - NumberElements = GettingTreadsFromFasta.SeqLength(); - ClustElements = GettingTreadsFromFasta.NelementsClustering(); - - double[][] table = new double[Nreads][Nreads]; - - for (int i=0; im[i]){ - diff[i] = list[m[i]][a]; - } - else{ - diff[i]=0; - } - } - - int val=0; - double temp=diff[0]; - for(int i=0;i<2;++i){ - if(diff[i] t = new ArrayList() ; - - for(int i=0;i<2;++i){ - m[i]=0; // initializing means to 0 - } - - //int cnt=0; - for(int i=0;i<2;++i){ - - for(int j=0;j s = new ArrayList(); - - for(int x1 = 0; x1t.get(x2)){ - sum = sum+(list[t.get(x2)][t.get(x1)]*list[t.get(x2)][t.get(x1)]); - } - else{ - sum = sum+0; - } - } - //System.out.println(sum); - s.add(sum/t.size()); - } - //System.out.println(s.size()); - double min=s.get(0); - int d2 = 0; - for(int x3=1;x3s.get(x3)){ - min = s.get(x3); - d2 = x3; - } - } - - - m[i]=t.get(d2); - s.clear(); - t.clear(); - - } - } - -//This checks if previous k ie. tempk and current k are same.Used as terminating case. - static int VerifyEqual(int n, int k[][], int tempk[][]){ - for(int i=0;i<2;++i){ - for(int j=0;j