dnaMethyAge包学习笔记

 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.xdnf.cn/news/18154.html

如若内容造成侵权/违法违规/事实不符,请联系一条长河网进行投诉反馈,一经查实,立即删除!

相关文章

详解八大排序(四)------(归并排序)

文章目录 前言&#xff1a;1 递归版本&#xff08;MergeSort&#xff09;1.1 核心思路1.2 实现代码 2 非递归版本&#xff08;MergeSortNonR&#xff09;2.1 核心思路2.2 实现代码 3.完整代码 前言&#xff1a; 归并排序的核心思路是把数组里面的数两两分成一组&#xff0c;组内…

商城小程序的流程渠道拓展

传统印象里&#xff0c;小程序的开发制作似乎很难&#xff0c;尤其是商城类型且功能体系完善的&#xff0c;事实也确实如此&#xff0c;没有较高的技术和成本投入或团队各个流程的专业人员合作&#xff0c;很难开发出来成品&#xff0c;或者质量较低。 当然对于大公司来说&…

小程序-基于java+SpringBoot+Vue的超市购物系统设计与实现

项目运行 1.运行环境&#xff1a;最好是java jdk 1.8&#xff0c;我们在这个平台上运行的。其他版本理论上也可以。 2.IDE环境&#xff1a;IDEA&#xff0c;Eclipse,Myeclipse都可以。推荐IDEA; 3.tomcat环境&#xff1a;Tomcat 7.x,8.x,9.x版本均可 4.硬件环境&#xff1a…

Android开发-Pokémon界面设计

实现效果图&#xff0c;还没更新完 二、功能说明&#xff1a; 在上次实验的基础之上把recycleviewb列表完善并且增加点击效果&#xff0c;点击之后可以跳转到另外一个activity上&#xff0c;并且添加返回按钮&#xff0c;可以放回原列表页面&#xff0c;列表中的每一行都对应的…

jenkins的安装(War包安装)

‌Jenkins是一个开源的持续集成工具&#xff0c;基于Java开发&#xff0c;主要用于监控持续的软件版本发布和测试项目。‌ 它提供了一个开放易用的平台&#xff0c;使软件项目能够实现持续集成。Jenkins的功能包括持续的软件版本发布和测试项目&#xff0c;以及监控外部调用执行…

HopToDesk 安全加密、免费开源,远程桌面新选择!

远程桌面工具越来越成为现代工作生活的刚需。你是否还在为寻找一个既安全又免费的工具而苦恼&#xff1f;HopToDesk&#xff0c;一款支持安全加密、免费开源的远程桌面软件&#xff0c;或许正是你的不二之 HopToDesk与传统的远程桌面工具相比有哪些独特优势&#xff1f;它如何…

yum工具的学习

Linux下安装软件的方法 1.源代码安装 2.rmp包安装 3.包管理器进行安装 --- yum/apt Linux下载软件的过程 操作系统的好坏评估 -- 生态问题 yum具体操作 Linux软件安装所有人都能用&#xff0c;是以other的身份去执行可执行程序 文件拷贝&#xff08;sudo&#xff09;-- &g…

react 如何修改弹出的modal的标题

原来标题的样子&#xff1a; 修改为&#xff1a; 实现方式&#xff1a; <Modal title<span>股价趋势/{this.state.pccode}</span> visible{this.state.isPriceModalOpen} style{{ top: 20 }} width{1320} height{400} footer{null} onCancel{()>this.hideMo…

计算机网络-理论部分(一):概览

重点 计算机网络的重点是协议&#xff0c;各种协议&#xff0c;每种协议都有自己对应的应用场景以及对应功能&#xff0c;学好协议&#xff0c;就学好了计算机努网络。 协议分层 协议分层围绕数据的传递展开。数据的传递需要包括&#xff1a;打包数据&#xff0c;控制传递&a…

开源科学工程技术软件介绍 – EDA工具KLayout

link 今天向各位知友介绍的 KLayout是一款由德国团队开发的开源EDA工具。 KLayout是使用C开发的&#xff0c;用户界面基于Qt。它支持Windows、MacOS和Linux操作系统。安装程序可以从下面的网址下载&#xff1a; https://www.klayout.de/build.html KLayout图形用户界面&…

SpringMVC的视图

目录 一.视图分类 &#xff08;1&#xff09;转发视图&#xff08;Forward View&#xff09;&#xff1a; &#xff08;2&#xff09;重定向视图&#xff08;Redirect View&#xff09;&#xff1a; &#xff08;3&#xff09;其他视图技术 二.转发视图 三.重定向视图 四…

Spring IOCDI

1. 什么是Spring 前面介绍了Spring Boot&#xff0c;Spring MVC&#xff0c;那么Spring和他们之间有什么关系呢&#xff1f; Spring简单一句话总结就是&#xff1a;它是一个包含众多工具方法的IOC容器。前面我们也接触过容器&#xff0c;比如List/Map&#xff0c;他俩是数据存…

螺旋矩阵II(leetcode 59)

转圈过程&#xff08;边界处理&#xff09;遵循循环不变量的原则&#xff0c;坚持一个原则处理每一条边&#xff0c;左闭右开处理 class Solution { public:vector<vector<int>> generateMatrix(int n) {vector<vector<int>> num(n, vector<int>…

Vue 中的透传,插槽,依赖注入

1. 透传attributes 在组件上使用透传attribute&#xff1a; 当你在父组件中使用子组件时&#xff0c;你可以添加一些attribute到子组件上&#xff0c;即使这些attribute没有在子组件的props中声明。 父组件&#xff1a; <!-- 父组件&#xff0c;例如 ParentComponent.vue…

基于K8S1.28.2实验rook部署ceph

基于K8S1.28.2实验rook部署ceph 原文链接&#xff1a; 基于K8S1.28.2实验rook部署ceph | 严千屹博客 Rook 支持 Kubernetes v1.22 或更高版本。 rook版本1.12.8 K8S版本1.28.2 部署出来ceph版本是quincy版本 主机名ip1(NAT)系统新硬盘磁盘内存master1192.168.48.101Centos7.910…

node.js中express的基本了解

定义 Express是基于Node.js平台&#xff0c;快速、开放、极简的Web开发框架。 本质 Express是一个npm上的第三方包&#xff0c;提供了快速创建Web服务器的便捷方法。 作用 与Node.js内置的http模块类似&#xff0c;Express也是专门用来创建Web服务器的&#xff0c;但它极大地简…

容器运行时 AND Docker

容器运行时 and Docker 什么是Docker Docker 使用 Google 公司推出的 Go 语言 进行开发实现&#xff0c;基于 Linux 内核的 cgroup&#xff0c;namespace&#xff0c;以及 AUFS 类的 Union FS 等技术&#xff0c;对进程进行封装隔离&#xff0c;属于 操作系统层面的虚拟化技术…

基于Python深度学习的【垃圾识别系统】实现~TensorFlow+人工智能+算法网络

一、介绍 垃圾识别分类系统。本系统采用Python作为主要编程语言&#xff0c;通过收集了5种常见的垃圾数据集&#xff08;‘塑料’, ‘玻璃’, ‘纸张’, ‘纸板’, ‘金属’&#xff09;&#xff0c;然后基于TensorFlow搭建卷积神经网络算法模型&#xff0c;通过对图像数据集进…

Scala-数据类型-概述(Scala 3.x 类型层次结构)

Scala Scala-数据类型 Scala1. Any — 顶级类型2. Matchable — 匹配类型3. AnyVal — 值类型的父类4. AnyRef — 引用类型的父类5. Null - 引用类型的子类型Tips: 为什么 null 不推荐使用&#xff1f; 6. Nothing - 底层类型 (Bottom Type)整理不易&#xff0c;对您有帮助的话…

Linux:权限相关知识详解

1.shell命令以及运行原理 1.1初步理解认识shell Linux严格意义上说的是一个操作系统&#xff0c;我们称之为“核心&#xff08;kernel&#xff09;“ &#xff0c;但我们一般用户&#xff0c;不能直接使用kernel。而是通过kernel的“外壳”程序&#xff0c;也就是所谓的shell&…