Mercurial > hg > jslab
comparison src/samer/mds/CorrelationTask.java @ 0:bf79fb79ee13
Initial Mercurial check in.
author | samer |
---|---|
date | Tue, 17 Jan 2012 17:50:20 +0000 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:bf79fb79ee13 |
---|---|
1 package samer.mds; | |
2 | |
3 // import samer.core.*; | |
4 import samer.maths.*; | |
5 import samer.tools.*; | |
6 | |
7 /** | |
8 Transfer correlation matrix to MDS distances. Assumes 1s | |
9 down the main diagonal and a symmetric matrix. | |
10 */ | |
11 public class CorrelationTask extends AnonymousTask | |
12 { | |
13 int N; | |
14 double [] d; // linear array of distances | |
15 double [][] _R; // matrix of correlation coefficients | |
16 | |
17 /** link each object to all the others using distances in matrix, returns a task | |
18 that can be used to refresh distances from original matrix */ | |
19 public CorrelationTask(MDSBase mds, Matrix R) { | |
20 N = R.getRowDimension(); | |
21 | |
22 d=new double[N*(N - 1)/2]; | |
23 _R=R.getArray(); | |
24 | |
25 mds.clearLinks(d); | |
26 for (int k=0, i=0; i<N; i++) | |
27 for (int j=0; j<i; j++) mds.setLink(k++,i,j); | |
28 | |
29 run(); | |
30 } | |
31 | |
32 public Vec getDistances() { return new Vec.ForArray(d); } | |
33 public void run() { | |
34 for (int k=0, i=0; i<N; i++) { | |
35 double [] Ri=_R[i]; | |
36 for (int j=0; j<i; j++) | |
37 d[k++]=Math.sqrt(-Math.log(Math.abs(Ri[j]))); | |
38 } | |
39 } | |
40 } |