是用来筛选不同组之间高表达基因的代码,好像不是用R实现的,或许有人知道怎么用吗
**r语言貌似也能做,下面是我用py实现的一个例子:
import pandas as pd
df = pd.read_csv('你的文件读取路径')
groups = df.groupby('group')
for name, group in groups:
mean_expr = group.mean(axis=0)
diff_expr = mean_expr - groups.mean().mean()
threshold = 1.5
high_expr_genes = diff_expr[diff_expr > threshold].index.tolist()
print(f"High expression genes in group {name}: {', '.join(high_expr_genes)}")