以GSE27284为例，提取出U值和M值后

GSE27284是一个基因表达谱数据集，可以使用R语言中的GEOquery包或Python中的GEOparse包来下载和处理数据。经过下载和预处理后，可以使用以下代码来提取U值和M值：

library(limma)
library(GEOquery)

gse <- getGEO("GSE27284", destdir = "./")
gse <- gse[[1]]

exprs <- exprs(gse)
design <- model.matrix(~factor(pData(gse)$tissue))
fit <- lmFit(exprs, design)
fit <- eBayes(fit)

U <- fit$coefficients[,2]
M <- fit$Amean

# U值和M值的前10个数
head(U, 10)
head(M, 10)

import GEOparse
import numpy as np
import pandas as pd
from statsmodels.formula.api import ols
from statsmodels.stats.anova import anova_lm

gse = GEOparse.get_GEO("GSE27284", destdir="./")
df = pd.DataFrame(gse.table)

df.set_index("ID_REF", inplace=True)
df.drop(columns=["GeneSymbol"], inplace=True)
df = np.log2(df.astype(float) + 1)

metadata = gse.phenotype_data
metadata["Intercept"] = 1
model = ols("tissue ~ 0 + Intercept", metadata).fit()
U = model.params["Intercept"]
M = df.mean(axis=1)

# U值和M值的前10个数
print(U.head(10))
print(M.head(10))