python计算spearman_如何在python中使用scipy计算spearman相关性?

这篇博客介绍了如何在Python脚本中利用Scipy库计算基因表达数据的Spearman相关性。作者询问是否可以在一个命令中使用Scipy来比较所有复制品的Spearman等级相关性,并分享了处理数据和计算Z分数的代码片段,期望得到社区的帮助。
摘要由CSDN通过智能技术生成

我有一个python脚本,它从基因表达式数据生成z分数。我还希望脚本创建一个spearman相关输出。我想知道复制品之间的相互关系。所以在输出中,行和列都是复制名。在

我还没有找到任何关于如何使用scipy计算两个以上列表的spearman等级相关性的例子。我能用西皮的吗统计斯皮尔曼在一个命令中比较所有的复制?我把数据存储在一个字符串中。以下是我的一些代码:t_vals=[]

z_outline = []

z_outlines = []

header = None

for line in input_file:

line = formatLine(line)

if line.startswith("#"):

header = line

header += "\n"

z_output.write(header)

sp_output.write(header)

continue

lineList = line.split("\t")

vals = []

in_vals =[]

for item in lineList[1:]:

val = float(item)

in_vals.append(val)

vals.append(val)

median = robjects.r['mean'](robjects.FloatVector(vals))[0]

mad = robjects.r['sd'](robjects.FloatVector(vals))[0]

z_vals

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值