Я хочу добавить adj.pval в корреляционную матрицу, созданную с помощью Pearson. Ниже приведены мои коды для создания парной матрицы.
from pandas import *
import numpy as np
import pandas as pd
from scipy.stats.stats import pearsonr
import itertools
# Generate sample data
data = np.random.rand(10, 6)
# Create a DataFrame with row names and column names
genes = ['Gene{}'.format(i+1) for i in range(6)]
samples = ['Sample{}'.format(i+1) for i in range(10)]
df = pd.DataFrame(data, index=samples, columns=genes)
correlations = {}
#store the gene names
columns = df.columns.tolist()
#perform pearson correlation
for col_a, col_b in itertools.combinations(columns, 2):
correlations[col_a +':' + col_b] = pearsonr(df.loc[:, col_a], df.loc[:, col_b])
#Convert the dictionary into dataframe
result = DataFrame.from_dict(correlations, orient='index')
result.columns = ['PCC', 'p-value']
#Convert index as a column
result.reset_index(level=0, inplace=True)
#rename the converted column
result.rename(columns = {'index':'Genes'}, inplace = True)
#split the column as two columns
result[['Gene-1','Gene-2']] = result.Genes.str.split(":", expand=True)
#select the required variables
result = result[['Gene-1', 'Gene-2', 'PCC', 'p-value']]
Он дал мне парную матрицу с pvalue. Я тоже хочу добавить adj.pval? Как мне добавить его?
Вы можете использовать multitest
для расчета скорректированного значения p.
from statsmodels.stats import multitest
# Calculate adjusted p-values
reject, adj_pvals, _, _ = multitest.multipletests(result['p-value'], method='bonferroni')
# Add adjusted p-value to the DataFrame
result['adj.pval'] = adj_pvals
Я не уверен насчет отрицательных голосов, но вы могли бы уточнить
adj.pval
. Вероятно, это означает «р-значения, скорректированные для множественных сравнений», но это не является явным, и существует множество возможных корректировок такого рода — Бонферрони, Беньямини-Хохберга, Беньямини-Йекутиэли и т. д., или что-то совсем другое? Также рассмотрите возможность замены таблицы скриншотов текстом.