Pyspark计算RDD中所有向量之间的自定义距离

时间:2017-03-07 22:10:14

标签: python pyspark rdd similarity

我有一个由密集向量组成的RDD,其中包含如下概率分布

[DenseVector([0.0806, 0.0751, 0.0786, 0.0753, 0.077, 0.0753, 0.0753, 0.0777, 0.0801, 0.0748, 0.0768, 0.0764, 0.0773]),
 DenseVector([0.2252, 0.0422, 0.0864, 0.0441, 0.0592, 0.0439, 0.0433, 0.071, 0.1644, 0.0405, 0.0581, 0.0528, 0.0691]),
 DenseVector([0.0806, 0.0751, 0.0786, 0.0753, 0.077, 0.0753, 0.0753, 0.0777, 0.0801, 0.0748, 0.0768, 0.0764, 0.0773]),
 DenseVector([0.0924, 0.0699, 0.083, 0.0706, 0.0766, 0.0708, 0.0705, 0.0793, 0.09, 0.0689, 0.0758, 0.0743, 0.0779]),
 DenseVector([0.0806, 0.0751, 0.0785, 0.0753, 0.077, 0.0753, 0.0753, 0.0777, 0.0801, 0.0748, 0.0768, 0.0764, 0.0773]),
 DenseVector([0.0806, 0.0751, 0.0786, 0.0753, 0.077, 0.0753, 0.0753, 0.0777, 0.0801, 0.0748, 0.0768, 0.0764, 0.0773])

我想计算向量和所有其他向量之间的相似性,并将结果存储在矩阵中。

我可以将完整的RDD转换为矩阵,然后获取每一行并计算与所有其他行的距离。我想知道是否有更有效的方法来使用pyspark RDD方法。

1 个答案:

答案 0 :(得分:2)

据我所知,没有用于在行之间进行余弦相似的功能。所以你必须要有点狡猾才能到达你想要的地方。

首先使用rdd.cartesian(rdd)以列格式创建行对,这将使所有行成对匹配。接下来,您需要定义余弦相似度函数并将其映射到rdd上。最后,将结果转换为np.array并重新整形为6x6。

示例:

def cos_sim(row):
    dot_product = row[0].dot(row[1])
    norm_a = np.sqrt(np.sum(row[0] * row[0]))
    norm_b = np.sqrt(np.sum(row[1] * row[1]))
    sim = dot_product / (norm_a * norm_b)
    return sim

rdd2 = rdd.cartesian(rdd)
cosine_similarities = rdd2.map(lambda x: cos_sim(x)).collect()
cosine_similariteis = np.array(cosine_similarities).reshape((6,6))