零基础入门转录组数据分析——预后模型之lasso模型

2024-08-28 04:52
文章标签 基础 入门 模型 数据分析 转录 lasso 预后

本文主要是介绍零基础入门转录组数据分析——预后模型之lasso模型，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

零基础入门转录组数据分析——预后模型之lasso模型

目录

  • 零基础入门转录组数据分析——预后模型之lasso模型
    • 1. 预后模型和lasso模型基础知识
    • 2. lasso预后模型（Rstudio）——代码实操
      • 2. 1 数据处理
      • 2. 2 构建lasso预后模型
      • 2. 3 提取Lasso预后基因
      • 2. 4 计算风险评分

1. 预后模型和lasso模型基础知识

1.1 预后模型是什么？
预后模型（Prognostic Model）是一种基于统计学和机器学习方法的预测工具，它通过分析患者的临床特征、生物标志物、治疗方案等多种因素，来预测患者未来发生某种不良事件（如疾病复发、死亡等）的概率。

1.2 预后模型的特点是什么？

  • 时间维度： 预后模型通常具有时间维度，即预测的是未来某一时刻或时间段内发生的事件，例如：在1000天的时候疾病患者的生存的概率。
  • 多因素分析： 预后模型考虑多种影响因素，包括患者的临床特征、生物学标志物、治疗方案等，以提高预测的准确性和可靠性。

1.3 构建预后模型对于临床应用有哪些帮助？

  • 指导临床决策： 预后模型可以帮助医生预测疾病的可能进展，从而为患者提供更加个性化的治疗建议。例如，在癌症治疗中，预后模型可以预测患者的复发风险和生存时间，帮助医生制定更加合理的治疗方案。
  • 疾病风险评估： 预后模型能够识别哪些患者更有可能出现不良结局，使医生可以进行早期干预，降低患者的风险。
  • 提供研究基础： 预后研究为后续的临床试验提供了有价值的数据，为疾病的进一步研究打下基础，例如：后续研究可以重点关注某些指标/特征。

1.4 风险评分是什么？
风险评分： 是一种量化评估患者疾病风险或治疗结果的方法。它通过综合考虑患者的多个临床特征和生物标志物，计算出一个数值（即风险评分），用于预测患者的预后（如生存率、复发率等）或治疗反应。

1.5 风险评分和预后模型的关系是什么？
预后模型往往通过计算风险评分来实现对患者预后的预测，风险评分是预后模型中的核心组成部分，它反映了模型对患者临床结局的量化评估。 简单理解为：风险评分就是预后模型的量化指标，后续也是通过风险评分来对预后模型进行验证

1.6 lasso是什么？
Lasso（Least Absolute Shrinkage and Selection Operator）回归，由Robert Tibshirani在1996年提出，是一种通过构造惩罚函数来实现变量选择和正则化的回归分析方法。它通过向损失函数中添加L1正则化项（即系数绝对值的和），来压缩一些回归系数，甚至使一些回归系数变为零，从而实现特征选择和降维。

1.7 lasso预后模型和之前章节介绍的lasso部分的区别是什么？
Lasso预后模型是一种结合Lasso回归方法的临床预后评估工具，它利用Lasso回归进行特征选择和系数压缩，从而构建出具有预测患者预后能力的模型。

举个栗子： 有8个基因，先通过lasso算法的特征选择思路从中筛选出4个更加重要的基因（筛选基因的目标是构建一个既简洁又准确的预测模型），之后基于这4个基因构建lasso预后模型，并计算风险评分。

综上所述： 预后模型就是一种预测工具，它可以通过多种特征/变量来预测未来某段时间患者的生存/复发概率，在预后模型中有一个量化指标叫做风险评分。在确认要输入的特征/变量之后，先通过lasso特征选择思路对输入的特征/变量进行筛选，用筛选后的特征/变量构建lasso预后模型并计算风险评分。

注意：
（1）与预后相关的分析输入的都是疾病样本，和对照样本无关！！
（2）输入的特征/变量必须是与预后相关的特征/变量，如何筛选预后特征，可见之前的章节：单因素cox筛选预后相关特征

2. lasso预后模型（Rstudio）——代码实操

本项目以TCGA——肺腺癌为例展开分析
物种：人类（Homo sapiens）
R版本：4.2.2
R包：tidyverse，survival，glmnet

废话不多说，代码如下：

2. 1 数据处理

设置工作空间：

rm(list = ls()) # 删除工作空间中所有的对象
setwd('/XX/XX/XX') # 设置工作路径
if(!dir.exists('./14_risk_model_lasso')){dir.create('./14_risk_model_lasso')
} 
setwd('./14_risk_model_lasso/')

加载包：

library(tidyverse)
library(survival)
library(glmnet)

导入要分析的表达矩阵train_data ，并对train_data 的列名进行处理（这是因为在读入的时候系统会默认把样本id中的“-”替换成“.”，所以要给替换回去）

train_data <- read.csv("./data_fpkm.csv", row.names = 1, check.names = F)  # 行名为全部基因名，每列为样本名
colnames(train_data) <- gsub('.', '-', colnames(train_data), fixed = T)

train_data 如下图所示，行为基因名（symbol），列为样本名
在这里插入图片描述
导入分组信息表group 

group <- read.csv("./data_group.csv", row.names = 1) # 为每个样本的分组信息（tumor和control）
colnames(group) <- c('sample', 'group')

group 如下图所示，第一列sample为样本名，第二列为样本对应的分组 （分组为二分类变量：disease和control）
在这里插入图片描述
导入样本对应的生存信息表survival_data

survival_data <- read.csv('./data_survival.csv', row.names = 1, check.names = F) # 样本生存信息（生存时间，生存状态）

survival_data 如下图所示，第一列sample为样本名，第二列为样本对应的生存状态（0表示存活，1表示死亡） ，第三列为样本对应的生存时间（单位是天）
在这里插入图片描述

接下来从group，train_data 和 survival_data中筛选出疾病样本

注：涉及到复发/生存的分析点基本都是用的疾病样本，不要附加对照样本！！！！

group <- group[group$group == 'disease', ] ## 筛选疾病样本
train_data <- train_data[, group$sample] ## 从全部表达矩阵中同样筛选出疾病样本
survival_data <- survival_data[survival_data$sample %in% group$sample, ] ## 从生存信息表中提取出疾病样本的生存信息

导入要筛选的基因hub_gene (12个基因)

hub_gene <- data.frame(symbol = gene <- c('VPS13D', 'MFF', 'ACSL1', 'VDAC1', 'PRELID1', 'BAK1','CYCS', 'PPIF', 'GLS2', 'BCL2L10', 'MPV17L', 'PHB'))
colnames(hub_gene) <- "symbol"

hub_gene 如下图所示，只有一列：12个基因的基因名
在这里插入图片描述

从train_data中取出这12个基因对应的表达矩阵，并且与之前准备的生存信息表从survival_data进行合并

dat <- train_data[rownames(train_data) %in% hub_gene$symbol, ] %>%t() %>% as.data.frame() # 整理后行为样本名，列为基因名
dat$sample <- rownames(dat)
dat <- merge(survival_data, dat, var = 'sample')
dat <- column_to_rownames(dat, var = "sample") %>% as.data.frame()

dat 如下图所示，行为样本名，第一列OS为生存状态，OS.time为生存时间，后面的列均为基因的表达量。
在这里插入图片描述

2. 2 构建lasso预后模型

先准备要输入的数据x_all 和y_all
x_all 就是dat中去除了OS和OS.time这两列剩余的基因表达矩阵
y_all 就是dat中仅选择OS和OS.time这两列

## 准备输入数据
x_all <- subset(dat, select = -c(OS, OS.time))
y_all <- subset(dat, select = c(OS, OS.time))

接下来构建lasso预后模型：使用cv.glmnet函数通过调整正则化参数（lambda）来寻找最优的lasso模型，以平衡模型的复杂度和拟合度，cv.glmnet函数常用参数介绍如下：

  • x参数——是要输入的基因表达矩阵（也称为特征或自变量），as.matrix函数将其转换为矩阵格式，因为glmnet函数需要输入数据为矩阵形式
  • y 参数——这部分使用survival包中的Surv函数创建生存对象。OS.time是生存时间，OS是一个指示变量，表示每个观测在观测结束时是否发生了事件（比如：死亡）。Surv函数将这两个变量结合，生成一个包含生存时间和事件状态的对象，供Cox模型使用。
  • family参数——这个参数指定了模型的类型，"cox"表示指定了模型家族为Cox比例风险模型
  • nfolds 参数——这个参数指定了交叉验证的折数（folds）。在这个例子中，数据集被分为15个大致相等的部分，模型在其中的14部分上进行训练，并在剩下的1部分上进行测试。这个过程重复15次，每次选择不同的部分作为测试集。通过这种方式，可以对整个数据集的不同部分评估模型的性能，并得到一个平均的性能指标，从而帮助选择最佳的lambda值

（cv.glmnet函数中比较关注的参数就是上述的这些，当然还有其他参数，如果想深入了解可自行查看官方说明文档）

set.seed(754) ## 设置种子
cvfit <- cv.glmnet(x = as.matrix(x_all), y = Surv(y_all$OS.time, y_all$OS), nfolds=15,family = "cox")

注：在构建模型的时候切记要设置种子（设置随机种子是为了确保结果的可重复性。由于交叉验证涉及随机分割数据，因此设置种子可以确保每次运行代码时，数据的分割方式都是相同的，从而得到相同的模型结果）

接下来从构建的最优模型中提取每个基因对应的系数（简单说就是提取基因）

## 提取指定lambda时特征的系数
coef.min <- coef(cvfit, s = "lambda.min") ## lambda.min & lambda.1se 取一个

注：lambda.min 或 lambda.1se 取一个就行，看个人需求，我这边选择lambda.min旨在考虑最佳预测性能，不考虑稳健性

  • lambda.min——是交叉验证过程中使得平均偏差（或其他指定的损失函数）最小的lambda值。换句话说，它是根据交叉验证结果，直接选择出的最优lambda值，因为它在平均意义上给出了最好的预测性能
  • lambda.1se——lambda.1se是lambda.min之后，使得平均偏差增加不超过一倍标准误差的最大的lambda值。这个选择标准旨在在模型的预测性能和复杂度之间找到一个折中。通过选择一个稍微更大的lambda值（即更强的正则化），lambda.1se可以帮助避免过拟合，在接受一些预测性能的损失的同时来换取更好的稳健性

coef.min中就对应着每个基因的系数，如下图所示，这些系数在后面就参与了风险评分的计算
在这里插入图片描述

2. 3 提取Lasso预后基因

提取出最优Lasso模型中筛选出来的那些基因及其对应的系数并保存 （筛选出来的基因被认为相对更加重要一些）

# 找出那些回归系数没有被惩罚为0的
active.min <- which(coef.min@i != 0)# 提取基因名称
lasso_geneids <- coef.min@Dimnames[[1]][coef.min@i+1]df.coef <- cbind(gene = rownames(coef.min), coefficient = coef.min[,1]) %>% as.data.frame()
df.coef <- subset(df.coef, coefficient != 0) %>% as.data.framewrite.csv(df.coef, "Lasso_Coefficients.csv")
write.csv(lasso_geneids, './lasso_genes.csv')

df.coef如下图所示，这里面就是筛选出来的8个预后基因，可以看到相较于一开始输入的12个基因，剔除了4个与预后不相干的基因。
在这里插入图片描述

2. 4 计算风险评分

接下来就用这8个基因构成的lasso模型计算风险评分，这个风险评分就是模型对患者未来健康状况或疾病进展风险的量化评估。

使用glmnet包的predict函数来根据之前通过交叉验证得到的最佳模型（cvfit）来预测风险评分

riskScore <- predict(cvfit, newx = as.matrix(x_all), s = cvfit$lambda.min) %>% as.data.frame() # 行为每个样本名，第一列为其对应的风险评分
colnames(riskScore) <- 'riskscore'

riskScore 如下图所示，行名为样本名，第一列riskscore就是对应样本的风险评分——就是给每个样本分配一个具体的分数，以反映其未来发生不良事件（如疾病复发、死亡等）的风险大小。
在这里插入图片描述

将风险评分和生存信息，基因表达量等整理到一张表格中，后续验证会用到

riskScore$sample <- rownames(riskScore)lasso_data <- dplyr::select(dat, OS, OS.time, lasso_geneids)
lasso_data$sample <- rownames(lasso_data)risk <- merge(lasso_data, riskScore, by = 'sample')
risk$risk <- ifelse(risk$riskscore > median(risk$riskscore), 1, 0)write.csv(risk, './risk.csv')

risk 如下图所示，第一列sample为样本名，第二列OS就是样本对应的生存状态，OS.time就是样本对应的生存时间，中间的列就是那8个预后基因对应的表达量，倒数第二列riskscore就是前一步计算出来的风险评分，最后一列risk就是对应的风险分组，1表示高风险组，0表示低风险组，区分度是风险评分的中位值。
在这里插入图片描述

计算出来的风险评分，就是模型对患者未来健康状况或疾病进展风险的量化评估，后续将针对计算出来的风险评分去研究如何对预后模型进行验证，在本章中不做介绍，会在全部预后模型介绍完之后统一进行说明。

注：lasso结果的可视化在此不做过多介绍，在之前的章节中有过说明

结语：

以上就是构建预后模型之lasso模型的所有过程，如果有什么需要补充或不懂的地方，大家可以私聊我或者在下方评论。

如果觉得本教程对你有所帮助，希望广大学习者能够花点自己的小钱支持一下（点赞旁的打赏按钮）作者创作（可以的话一杯蜜雪奶茶即可），感谢大家的支持~~~~~~ ^_^ !!!

祝大家能够开心学习，轻松学习，在学习的路上少一些坎坷~~~

请添加图片描述

  • 目录部分跳转链接：零基础入门生信数据分析——导读

这篇关于零基础入门转录组数据分析——预后模型之lasso模型的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/1113779。

相关文章

Python基于火山引擎豆包大模型搭建QQ机器人详细教程(2024年最新)

Python基于火山引擎豆包大模型搭建QQ机器人详细教程(2024年最新)

《Python基于火山引擎豆包大模型搭建QQ机器人详细教程(2024年最新)》：本文主要介绍Python基于火山引擎豆包大模型搭建QQ机器人详细的相关资料,包括开通模型、配置APIKEY鉴权和SD... 目录豆包大模型概述开通模型付费安装 SDK 环境配置 API KEY 鉴权Ark 模型接口Prompt
阅读更多...
Spring Security 从入门到进阶系列教程

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门（一）：登录与退出》 《Spring-Security-入门（二）：基于数据库验证》 《Spring-Security-入门（三）：密码加密》 《Spring-Security-入门（四）：自定义-Filter》 《Spring-Security-入门（五）：在 Sprin
阅读更多...
大模型研发全揭秘：客服工单数据标注的完整攻略

大模型研发全揭秘：客服工单数据标注的完整攻略

在人工智能（AI）领域，数据标注是模型训练过程中至关重要的一步。无论你是新手还是有经验的从业者，掌握数据标注的技术细节和常见问题的解决方案都能为你的AI项目增添不少价值。在电信运营商的客服系统中，工单数据是客户问题和解决方案的重要记录。通过对这些工单数据进行有效标注，不仅能够帮助提升客服自动化系统的智能化水平，还能优化客户服务流程，提高客户满意度。本文将详细介绍如何在电信运营商客服工单的背景下进行
阅读更多...
Andrej Karpathy最新采访：认知核心模型10亿参数就够了，AI会打破教育不公的僵局

Andrej Karpathy最新采访：认知核心模型10亿参数就够了，AI会打破教育不公的僵局

夕小瑶科技说 原创  作者 | 海野 AI圈子的红人，AI大神Andrej Karpathy，曾是OpenAI联合创始人之一，特斯拉AI总监。上一次的动态是官宣创办一家名为 Eureka Labs 的人工智能+教育公司 ，宣布将长期致力于AI原生教育。 近日，Andrej Karpathy接受了No Priors（投资博客）的采访，与硅谷知名投资人 Sara Guo 和 Elad G
阅读更多...
零基础学习Redis(10) -- zset类型命令使用

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合，内部除了存储元素外，还会存储一个score，存储在zset中的元素会按照score的大小升序排列，不同元素的score可以重复，score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]
阅读更多...
Retrieval-based-Voice-Conversion-WebUI模型构建指南

Retrieval-based-Voice-Conversion-WebUI模型构建指南

一、模型介绍 Retrieval-based-Voice-Conversion-WebUI（简称 RVC）模型是一个基于 VITS（Variational Inference with adversarial learning for end-to-end Text-to-Speech）的简单易用的语音转换框架。 具有以下特点 简单易用：RVC 模型通过简单易用的网页界面，使得用户无需深入了
阅读更多...
透彻！驯服大型语言模型（LLMs）的五种方法，及具体方法选择思路

透彻！驯服大型语言模型（LLMs）的五种方法，及具体方法选择思路

引言 随着时间的发展，大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用，随着人们期望的不断增加，目标也发生了巨大的变化。在短短的几个月的时间里，人们对大模型的认识已经从对其zero-shot能力感到惊讶，转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构（例如Transformer）对海量的、多种多样的数据分布进行建模得到，它包含了大量的先验
阅读更多...
图神经网络模型介绍（1）

图神经网络模型介绍（1）

我们将图神经网络分为基于谱域的模型和基于空域的模型，并按照发展顺序详解每个类别中的重要模型。 1.1基于谱域的图神经网络         谱域上的图卷积在图学习迈向深度学习的发展历程中起到了关键的作用。本节主要介绍三个具有代表性的谱域图神经网络：谱图卷积网络、切比雪夫网络和图卷积网络。 （1）谱图卷积网络 卷积定理：函数卷积的傅里叶变换是函数傅里叶变换的乘积，即F{f*g}
阅读更多...
秋招最新大模型算法面试，熬夜都要肝完它

秋招最新大模型算法面试，熬夜都要肝完它

💥大家在面试大模型LLM这个板块的时候，不知道面试完会不会复盘、总结，做笔记的习惯，这份大模型算法岗面试八股笔记也帮助不少人拿到过offer ✨对于面试大模型算法工程师会有一定的帮助，都附有完整答案，熬夜也要看完，祝大家一臂之力 这份《大模型算法工程师面试题》已经上传CSDN，还有完整版的大模型 AI 学习资料，朋友们如果需要可以微信扫描下方CSDN官方认证二维码免费领取【保证100%免费
阅读更多...
数论入门整理（updating）

数论入门整理（updating）

一、gcd lcm 基础中的基础，一般用来处理计算第一步什么的，分数化简之类。 LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; } <pre name="code" class="cpp">LL lcm(LL a, LL b){LL c = gcd(a, b);return a / c * b;} 例题：
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2