cos_bray计算区别

import numpy as np
from sklearn.metrics.pairwise import cosine_similarity

A = np.array([10, 8, 4, 1, 1]).reshape(1,-1)
B = np.array([7, 3, 8, 4, 0]).reshape(1,-1)
A = np.array([10, 8, 4, 1, 1]).reshape(1,-1)
B = np.array([10, 8, 4, 1, 1]).reshape(1,-1)
print(A,B)
[[10 8 4 1 1]] [[10 8 4 1 1]]
#函数进行cos以及Bray(包含丰度信息)计算,两个方法相互验证
sim1=cosine_similarity(A,B)
print(sim1)
def cosine_similarity_define(x,y):#相似性
num = x.dot(y.T)
denom = np.linalg.norm(x) * np.linalg.norm(y)
return num / denom
sim2=cosine_similarity_define(A,B)
print(sim2)
def Bray_Curtis(x, y):#群落之间组成差异
import numpy as np
x = np.array(x)
y = np.array(y)
return np.sum(np.abs(x-y))/(np.sum(x)+np.sum(y))
sim3=Bray_Curtis(A,B)
print(sim3)

结果:
[[0.82029132]]
[[0.82029132]]
0.34782608695652173

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值