dnaMethyAge包学习笔记

news/2024/11/19 16:18:09/

 1.introduction

      许多对甲基化年龄进行计算的文章都是采用网站实现计算的,能够实现对甲基化年龄的计算的R包相对比较少,其中应用最广的是dnaMethyAge包。作者本想寻找能够计算Grimage和Grimage2的R包,奈何没有寻找到,因此只能记录一下能够计算其他许多甲基化年龄的R包dnaMethyAge的学习笔记。

2.example

#安装和加载R包dnaMethyAge
devtools::install_github("yiluyucheng/dnaMethyAge")
library('dnaMethyAge')## prepare betas dataframe
data('subGSE174422') ## load example betas,>=0.6,<=0.2
#                 GSM5310260_3999979009_R02C02 GSM5310261_3999979017_R05C01
#cg00000029                               0.29                         0.31
#cg00000108                               0.84                         0.83
#cg00000109                               0.74                         0.70
#cg00000165                               0.09                         0.09
#cg00000236                               0.55                         0.55
#cg00000289                               0.53                         0.43
#cg00000292                               0.70                         0.67
#cg00000321                               0.12                         0.14
#cg00000363                               0.24                         0.20
#cg00000622                               0.03                         0.03
availableClock()#List all supported clocks
# [1] "HannumG2013"    "HorvathS2013"   "LevineM2018"    "ZhangQ2019"    
# [5] "ShirebyG2020"   "YangZ2016"      "ZhangY2017"     "LuA2019"       
# [9] "HorvathS2018"   "DunedinPACE"    "McEwenL2019"    "CBL_specific"  
#[13] "PCGrimAge"      "PCHorvathS2013" "PCHannumG2013"  "PCHorvathS2018"
#[17] "PCPhenoAge"     "CBL_common"     "Cortex_common"  "epiTOC2"       
#[21] "BernabeuE2023c" "LuA2023p1"      "LuA2023p2"      "LuA2023p3"  
#设置想计算的甲基化时钟
clock_name <- 'HorvathS2013'
#计算Horvath2013
horvath_age <- methyAge(betas, clock=clock_name)
print(horvath_age)
#                        Sample     mAge
#1 GSM5310260_3999979009_R02C02 74.88139
#2 GSM5310261_3999979017_R05C01 62.36400
#3 GSM5310262_3999979018_R02C02 68.04759
#4 GSM5310263_3999979022_R02C01 61.62691
#5 GSM5310264_3999979027_R02C01 59.65161
#6 GSM5310265_3999979028_R01C01 60.95991
#7 GSM5310266_3999979029_R04C02 52.48954
#8 GSM5310267_3999979031_R06C02 64.29711

        还能通过纳入age年龄计算得到age acceleration。

print(info)
#                        Sample  Age    Sex
#1 GSM5310260_3999979009_R02C02 68.8 Female
#2 GSM5310261_3999979017_R05C01 45.6 Female
#3 GSM5310262_3999979018_R02C02 67.4 Female
#4 GSM5310263_3999979022_R02C01 45.6 Female
#5 GSM5310264_3999979027_R02C01 62.5 Female
#6 GSM5310265_3999979028_R01C01 45.1 Female
#7 GSM5310266_3999979029_R04C02 53.2 Female
#8 GSM5310267_3999979031_R06C02 63.8 Female
horvath_age <- methyAge(betas, clock=clock_name, age_info=info, fit_method='Linear', do_plot=TRUE)
print(horvath_age)
#                        Sample  Age    Sex     mAge Age_Acceleration
#1 GSM5310260_3999979009_R02C02 68.8 Female 74.88139         7.334461
#2 GSM5310261_3999979017_R05C01 45.6 Female 62.36400         3.318402
#3 GSM5310262_3999979018_R02C02 67.4 Female 68.04759         1.013670
#4 GSM5310263_3999979022_R02C01 45.6 Female 61.62691         2.581311
#5 GSM5310264_3999979027_R02C01 62.5 Female 59.65161        -5.586763
#6 GSM5310265_3999979028_R01C01 45.1 Female 60.95991         2.097534
#7 GSM5310266_3999979029_R04C02 53.2 Female 52.48954        -9.340977
#8 GSM5310267_3999979031_R06C02 63.8 Female 64.29711        -1.417638

 参数"do_plot=TRUE"能生成methyage和age的散点图


http://www.ppmy.cn/news/1548272.html

相关文章

css uniapp背景图宽度固定高度自适应可以重复

page {height: 100%;background-image: url(https://onlinekc.a.hlidc.cn/uploads/20241115/350f94aaf493d05625a7ddbc86c7804e.png);background-repeat: repeat;background-size: contain;} 如果不要重复 把background-repeat: repeat;替换background-repeat: no-repeat;

《七大排序算法:原理、实现及性能剖析》

目录 一、七大排序的实现、原理及性能剖析1. 插入排序2. 希尔排序3. 选择排序4. 堆排序5. 冒泡排序6. 快速排序7. 快速排序&#xff08;非递归&#xff09;8. 归并排序9. 归并排序&#xff08;非递归&#xff09; 二、七大排序的复杂度和稳定性分析 一、七大排序的实现、原理及…

报错java: java.lang.NoSuchFieldError: Class com.sun.tools.javac.tree.JCTree$JCImport does not ...解决方法

在运行项目时出现java: java.lang.NoSuchFieldError: Class com.sun.tools.javac.tree.JCTree$JCImport does not have member field com.sun.tools.javac.tree.JCTree qualidzz这样的报错 解决方法 1.第一步&#xff1a;在pom文件中将lombok的版本改成最新的 此时1.18.34是新…

【网络安全】SSL(一):为什么需要 Keyless SSL?

未经许可,不得转载。 文章目录 背景正文背景 随着网站和应用程序向云端迁移,使用 HTTPS(SSL/TLS)加密流量已成为行业标准。然而,传统的 HTTPS 配置要求服务器持有网站的私钥,这在云计算环境中引发了一系列安全性和合规性问题。一旦云服务器遭到攻击,私钥泄露可能带来不…

集群聊天服务器(9)一对一聊天功能

目录 一对一聊天离线消息服务器异常处理 一对一聊天 先新添一个消息码 在业务层增加该业务 没有绑定事件处理器的话消息会派发不出去 聊天其实是服务器做一个中转 现在同时登录两个账号 收到了聊天信息 再回复一下 离线消息 声明中提供接口和方法 张三对离线的李…

springboot 获取spring上下文

有一个工具类可以直接获取&#xff0c;不知道为啥所有的文章基本都写了获取spring上下文的4种方法&#xff0c;抄来抄去的&#xff0c;然后获取的还特别麻烦 SpringContextUtil.getApplicationContext()直接就可以获取了 抄博客的真恶心&#xff01;&#xff01;&#xff01;…

十一:HTTP 状态码详解:解读每一个响应背后的意义

HTTP(超文本传输协议)是网络通信的基石之一,主要用于客户端(例如浏览器)和服务器之间的通信。为了让服务器能准确地向客户端反馈请求的处理状态,HTTP设计了一套标准的状态码。每一个状态码代表了特定的含义,指示了请求的状态、潜在的问题或成功的信息。 1. 信息响应 (1…

Spring Boot中集成Redis与MySQL

1. 环境准备与依赖配置 1.1 Maven 依赖管理 为了在 Spring Boot 项目中使用 Redis 和 MySQL&#xff0c;我们需要在 pom.xml 中添加必要的依赖。主要包括以下几个依赖&#xff1a; Spring Data Redis&#xff1a;用于在 Spring Boot 中集成 Redis&#xff0c;提供 RedisTem…