Five most popular similarity measures implementation in python
The buzz term similarity distance measure has got a wide variety of definitions among the math and data mining practitioners. As a result, those terms, concepts and their usage went way beyond the head for the beginner , Who started to understand them for the very first time. So today I write this post to give more simplified and very intuitive definitions for similarity and I will drive to Five most popular similarity measures and implementation of them.
Before going to explain different similarity distance measures let me explain the effective key term similarity in datamining. This similarity is the very basic building block for activities such as Recommendation engines, clustering, classification and anomaly detection.
Similarity:
The similarity measure is the measure of how much alike two data objects are. Similarity measure in a data mining context is usually described as a distance with dimensions representing features of the objects. If this distance is small it will be the high degree of similarity where large distance will be the low degree of similarity. Similarity is subjective and is highly dependent on the domain and application. For example, two fruits are similar because of color or size or taste. Care should be taken when calculating distance across dimensions/features that are unrelated. The relative values of each feature must be normalized or one feature could end up dominating the distance calculation. Similarity are measured in the range 0 to 1 [0,1].
Two main consideration about similarity:
 Similarity = 1 if X = Y (Where X, Y are two objects)
 Similarity = 0 if X ≠ Y
That’s all about similarity let’s drive to five most popular similarity distance measures.
Euclidean distance:
Euclidean distance is the most common use of distance. In most cases when people said about distance, they will refer to Euclidean distance. Euclidean distance is also known as simply distance. When data is dense or continuous , this is the best proximity measure. The Euclidean distance between two points is the length of the path connecting them.This distance between two points is given by the Pythagorean theorem.
Euclidean distance implementation in python:
1 2 3 4 5 6 7 8 9 10 
#!/usr/bin/env python from math import* def euclidean_distance(x,y): return sqrt(sum(pow(ab,2) for a, b in zip(x, y))) print euclidean_distance([0,3,4,5],[7,6,3,1]) 
Script Output:
1 2 
9.74679434481 [Finished in 0.0s] 
Manhattan distance:
Manhattan distance is a metric in which the distance between two points is the sum of the absolute differences of their Cartesian coordinates. In a simple way of saying it is the absolute sum of the difference between the xcoordinates and ycoordinates. Suppose we have two points A and B if we want to find the manhattan distance between them, just we have, to sum up, the absolute xaxis and y – axis variation means we have to find how these to points A and B are varying in Xaxis and Y axis. In a more mathematical way of saying Manhattan distance between two points measured along axes at right angles.
In a plane with p1 at (x1, y1) and p2 at (x2, y2).
Manhattan distance = x1 – x2 + y1 – y2
This Manhattan distance metric is also known as Manhattan length,rectilinear distance, L1 distance or L1 norm ,city block distance, Minkowski’s L1 distance, taxicab metric, or city block distance.
Manhattan distance implementation in python:
1 2 3 4 5 6 7 8 9 
#!/usr/bin/env python from math import* def manhattan_distance(x,y): return sum(abs(ab) for a,b in zip(x,y)) print manhattan_distance([10,20,10],[10,20,20]) 
Script Output:
1 2 
10 [Finished in 0.0s] 
Minkowski distance:
The Minkowski distance is a generalized metric form of Euclidean distance and Manhattan distance.
In the equation, d^MKD is the Minkowski distance between the data record i and j, k the index of a variable, n the total number of variables y and λ the order of the Minkowski metric. Although it is defined for any λ > 0, it is rarely used for values other than 1, 2 and ∞.
The way distances are measured by the Minkowski metric of different orders between two objects with three variables ( In the image it displayed in a coordinate system with x, y ,zaxes).
Synonyms of Minkowski:
Different names for the Minkowski distance or Minkowski metric arise from the order:
 λ = 1 is the Manhattan distance. Synonyms are L1Norm, Taxicab or CityBlock distance. For two vectors of ranked ordinal variables, the Manhattan distance is sometimes called Footruler distance.
 λ = 2 is the Euclidean distance. Synonyms are L2Norm or Ruler distance. For two vectors of ranked ordinal variables, the Euclidean distance is sometimes called Spearman distance.
 λ = ∞ is the Chebyshev distance. Synonyms are LmaxNorm or Chessboard distance.
reference.
Minkowski distance implementation in python:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
#!/usr/bin/env python from math import* from decimal import Decimal def nth_root(value, n_root): root_value = 1/float(n_root) return round (Decimal(value) ** Decimal(root_value),3) def minkowski_distance(x,y,p_value): return nth_root(sum(pow(abs(ab),p_value) for a,b in zip(x, y)),p_value) print minkowski_distance([0,3,4,5],[7,6,3,1],3) 
Script Output:
1 2 
8.373 [Finished in 0.0s] 
Cosine similarity:
Cosine similarity metric finds the normalized dot product of the two attributes. By determining the cosine similarity, we would effectively try to find the cosine of the angle between the two objects. The cosine of 0° is 1, and it is less than 1 for any other angle. It is thus a judgement of orientation and not magnitude: two vectors with the same orientation have a cosine similarity of 1, two vectors at 90° have a similarity of 0, and two vectors diametrically opposed have a similarity of 1, independent of their magnitude. Cosine similarity is particularly used in positive space, where the outcome is neatly bounded in [0,1]. One of the reasons for the popularity of cosine similarity is that it is very efficient to evaluate, especially for sparse vectors.
Cosine similarity implementation in python:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
#!/usr/bin/env python from math import* def square_rooted(x): return round(sqrt(sum([a*a for a in x])),3) def cosine_similarity(x,y): numerator = sum(a*b for a,b in zip(x,y)) denominator = square_rooted(x)*square_rooted(y) return round(numerator/float(denominator),3) print cosine_similarity([3, 45, 7, 2], [2, 54, 13, 15]) 
Script Output:
1 2 
0.972 [Finished in 0.1s] 
Jaccard similarity:
We so far discussed some metrics to find the similarity between objects. where the objects are points or vectors .When we consider about jaccard similarity this objects will be sets. So first let’s learn some very basic about sets.
Sets:
A set is (unordered) collection of objects {a,b,c}. we use the notation as elements separated by commas inside curly brackets { }. They are unordered so {a,b} = { b,a }.
Cardinality:
The cardinality of A denoted by A which counts how many elements are in A.
Intersection:
The intersection between two sets A and B is denoted A ∩ B and reveals all items which are in both sets A,B.
Union:
Union between two sets A and B is denoted A ∪ B and reveals all items which are in either set.
Now going back to Jaccard similarity. The Jaccard similarity measures the similarity between finite sample sets and is defined as the cardinality of the intersection of sets divided by the cardinality of the union of the sample sets. Suppose you want to find Jaccard similarity between two sets A and B it is the ration of cardinality of A ∩ B and A ∪ B
Jaccard similarity implementation:
1 2 3 4 5 6 7 8 9 10 11 12 
#!/usr/bin/env python from math import* def jaccard_similarity(x,y): intersection_cardinality = len(set.intersection(*[set(x), set(y)])) union_cardinality = len(set.union(*[set(x), set(y)])) return intersection_cardinality/float(union_cardinality) print jaccard_similarity([0,1,2,5,6],[0,2,3,5,7,9]) 
Script Output:
1 2 
0.375 [Finished in 0.0s] 
Implementaion of all 5 similarity measure into one Similarity class:
file_name : similaritymeasures.py
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 
#!/usr/bin/env python from math import* from decimal import Decimal class Similarity(): """ Five similarity measures function """ def euclidean_distance(self,x,y): """ return euclidean distance between two lists """ return sqrt(sum(pow(ab,2) for a, b in zip(x, y))) def manhattan_distance(self,x,y): """ return manhattan distance between two lists """ return sum(abs(ab) for a,b in zip(x,y)) def minkowski_distance(self,x,y,p_value): """ return minkowski distance between two lists """ return self.nth_root(sum(pow(abs(ab),p_value) for a,b in zip(x, y)), p_value) def nth_root(self,value, n_root): """ returns the n_root of an value """ root_value = 1/float(n_root) return round (Decimal(value) ** Decimal(root_value),3) def cosine_similarity(self,x,y): """ return cosine similarity between two lists """ numerator = sum(a*b for a,b in zip(x,y)) denominator = self.square_rooted(x)*self.square_rooted(y) return round(numerator/float(denominator),3) def square_rooted(self,x): """ return 3 rounded square rooted value """ return round(sqrt(sum([a*a for a in x])),3) def jaccard_similarity(self,x,y): """ returns the jaccard similarity between two lists """ intersection_cardinality = len(set.intersection(*[set(x), set(y)])) union_cardinality = len(set.union(*[set(x), set(y)])) return intersection_cardinality/float(union_cardinality) 
Using Similarity class:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
#!/usr/bin/env python from similaritymeasures import Similarity def main(): """ the main function to create Similarity class instance and get used to it """ measures = Similarity() print measures.euclidean_distance([0,3,4,5],[7,6,3,1]) print measures.jaccard_similarity([0,1,2,5,6],[0,2,3,5,7,9]) if __name__ == "__main__": main() 
You can get all the complete codes of dataaspirant at dataaspirant data science codes
Related Courses:
Do check out unlimited data science courses
Course  Course Link  What You will Learn 
Introduction to natural language processing 
Introduction to Natural Language Processing


Applied Text mining in python 
Applied Text Mining in Python


Easy Natural Language Processing (NLP) in Python 
Easy Natural Language Processing (NLP) in Python 

Follow us:
FACEBOOK QUORA TWITTER GOOGLE+  LINKEDIN REDDIT  FLIPBOARD  MEDIUM  GITHUB
I hope you like this post. If you have any questions then feel free to comment below. If you want me to write on one specific topic then do tell it to me in the comments below.
A few points:
1. A measure of similarity need not be symmetrical
2. A measure of similarity is not a metric space
3. Information theoretic measures, like KL and Mutual Information tend to be the most powerful, but the most difficult to manipulate mathematically.
Thought you might cover Mahalanobis distance.
I don’t think there is no need to write your own implementation. All of them and a lot more are already available in scipy.spatial.distance module in python.
Hi Jitesh Khandelwal!
You said true, but we have to explain how we can implement them.
Reblogged this on Random and commented:
A good blog, explaining some important similarity metrics.
Your post seems to cover just one similarity measure: Jaccard. The remaining four are distance metrics; they must be transformed to provide similarity.
I actually found Jaccard’s metric to work nicely for weighted sets as well: if an item occurs in both A and B, its weight in the intersection is the minimum of the two weights, and its weight in the union is the maximum of the two weights.
fantastic images
Thank you 🙂
Nice Post It is easily understood with list of x and y (two lists). so please I want to know more how to implement for large documents especially for cosine similarity in IR
Hi Misgu,
Thanks for your compliment.
Nice Post.Thanks a lot.:)
Excellent work bro. I love you.
Thanks 🙂
Good post
However multiple grammatical errors
Would love to correct them and contribute towards the site
Hi Chaitanya Bapat,
I am glad for your suggestion. IF you would like to contribute content please join us.