--- a
+++ b/clusters/scripts/ward.py
@@ -0,0 +1,46 @@
+from __future__ import print_function
+from sklearn.decomposition import TruncatedSVD
+from sklearn.feature_extraction.text import TfidfVectorizer
+from sklearn.feature_extraction.text import HashingVectorizer
+from sklearn.feature_extraction.text import TfidfTransformer
+from sklearn.preprocessing import Normalizer
+from sklearn import metrics
+
+from sklearn.cluster import KMeans, MiniBatchKMeans
+
+import logging
+import sys
+from time import time
+
+import numpy as np
+import pandas as pd
+
+data = pd.read_csv("cleanedctrpinds.csv", sep = ",", quoting = 1, quotechar = '"')
+dataset = []
+data = np.array(data)
+for x in data :
+	count = 0
+	placeholder = ""
+	for y in x :
+		if (count == 1) :
+			placeholder = y + "";
+		elif (count == 2) :
+			if y == True :
+				dataset.append(str(1) + " " + placeholder)
+			else :
+				dataset.append(str(0) + " " + placeholder)
+		count=count+1
+
+vectorizer = TfidfVectorizer(max_df=0.5, max_features=10000, min_df=2, stop_words='english', use_idf=True)
+X = vectorizer.fit_transform(dataset)
+print(X.shape)
+
+Ks = range(1, 140001, 20000)
+
+km = [KMeans(n_clusters=i, init='k-means++', max_iter=500, n_init=1,verbose=True) for i in Ks]
+score = [km[i].fit(X).score(X) for i in range(len(km))]
+
+f1 = open("./elbow9k.txt", "a")
+for x in range(len(score)) :
+	f1.write(str(Ks[x]) + " " + str(score[x]) + "\n")
+