annotate src/edu/unc/genomics/visualization/KMeansRow.java @ 2:e16016635b2a

Uploaded
author timpalpant
date Mon, 13 Feb 2012 22:12:06 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
2
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
1 package edu.unc.genomics.visualization;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
2
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
3 import java.io.Serializable;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
4 import java.util.Collection;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
5
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
6 import org.apache.commons.math.stat.clustering.Clusterable;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
7
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
8 /**
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
9 * @author timpalpant
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
10 * Holds a row of data for the KMeans program
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
11 */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
12 class KMeansRow implements Clusterable<KMeansRow>, Serializable {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
13
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
14 private static final long serialVersionUID = -323598431692368500L;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
15
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
16 private final String id;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
17 /** Point coordinates. */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
18 private final float[] point;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
19
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
20 /**
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
21 * Build an instance wrapping an float array.
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
22 * <p>The wrapped array is referenced, it is <em>not</em> copied.</p>
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
23 * @param point the n-dimensional point in integer space
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
24 */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
25 public KMeansRow(final String id, final float[] point) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
26 this.point = point;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
27 this.id = id;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
28 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
29
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
30 /**
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
31 * Get the n-dimensional point in float space.
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
32 * @return a reference (not a copy!) to the wrapped array
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
33 */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
34 public float[] getPoint() {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
35 return point;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
36 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
37
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
38 /** {@inheritDoc} */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
39 public double distanceFrom(final KMeansRow p) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
40 double sumSquares = 0;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
41 float[] otherPoint = p.getPoint();
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
42 for (int i = 0; i < point.length; i++) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
43 sumSquares += Math.pow(point[i]-otherPoint[i], 2);
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
44 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
45 return Math.sqrt(sumSquares);
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
46 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
47
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
48 /** {@inheritDoc} */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
49 public KMeansRow centroidOf(final Collection<KMeansRow> points) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
50 float[] centroid = new float[getPoint().length];
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
51 for (KMeansRow p : points) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
52 for (int i = 0; i < centroid.length; i++) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
53 centroid[i] += p.getPoint()[i];
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
54 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
55 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
56 for (int i = 0; i < centroid.length; i++) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
57 centroid[i] /= points.size();
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
58 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
59 return new KMeansRow(id, centroid);
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
60 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
61
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
62 /** {@inheritDoc} */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
63 @Override
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
64 public boolean equals(final Object other) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
65 if (!(other instanceof KMeansRow)) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
66 return false;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
67 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
68 final float[] otherPoint = ((KMeansRow) other).getPoint();
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
69 if (point.length != otherPoint.length) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
70 return false;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
71 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
72 for (int i = 0; i < point.length; i++) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
73 if (point[i] != otherPoint[i]) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
74 return false;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
75 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
76 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
77 return true;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
78 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
79
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
80 /** {@inheritDoc} */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
81 @Override
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
82 public int hashCode() {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
83 int hashCode = 0;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
84 for (Float i : point) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
85 hashCode += i.hashCode() * 13 + 7;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
86 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
87 return hashCode;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
88 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
89
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
90 /**
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
91 * {@inheritDoc}
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
92 * @since 2.1
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
93 */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
94 @Override
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
95 public String toString() {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
96 final StringBuilder buff = new StringBuilder(id);
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
97 for (float value : getPoint()) {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
98 buff.append("\t").append(value);
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
99 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
100 return buff.toString();
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
101 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
102
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
103 /**
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
104 * @return the id
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
105 */
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
106 public String getId() {
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
107 return id;
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
108 }
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
109
e16016635b2a Uploaded
timpalpant
parents:
diff changeset
110 }