小果的单日小技巧DAY9:850K甲基化数据的免疫浸润分析很难吗?






小果的单日小技巧DAY9:850K甲基化数据的免疫浸润分析很难吗?

小果  生信果  2023-08-02 19:00:23

生信人R语言学习必备

立刻拥有一个Rstudio账号

开启升级模式吧

(56线程,256G内存,个人存储1T)

 



850K甲基化芯片相对于450K甲基化芯片来说,具有更高的覆盖度和更好的分辨率。它检测的位点数量是450K的两倍以上,可以检测到更多的CpG位点和其他类型的位点。另外,850K芯片还能够覆盖更多的基因,因此在免疫浸润分析中,可以更全面地评估基因和免疫细胞之间的关系。上一期小果分享了云生信平台(云生信  – 学生物信息学 (biocloudservice.com))450K甲基化数据的免疫浸润分析方法,本期跟小果一起看看850K甲基化芯片的吧。

上传数据

1、Group.850.csv(分组信息)


2、methylation850.csv(不同样本各甲基化位点得分)


输出结果


1、Wilcox.result.csv(包含了经过 Wilcoxon rank-sum test 检验的差异甲基化位点的统计信息)


2、Immune.infiltration.result.csv(包含了不同类型免疫细胞浸润得分)


3、B.pdf


env_dir <- "D:/SCP_env/"dir.create(env_dir,recursive = TRUE)setwd(env_dir)
# 使用小果曾经教过大家的renv包,在本地创建一个新的环境if(!require("renv")) install.packages(renv)renv::init(project = env_dir)



4、meth.DMP.results.csv(差异甲基化位点的分析结果)



5、Neutro.pdf





6、Eosino.pdf



7、CD4T.pdf


8、NK.pdf


9、CD8T.pdf


10、Mono.pdf


好啦,到这里全部输出文件就介绍完毕了,生信云平台同时提供根据450K和850K两种甲基化芯片的结果进行免疫浸润分析。

看了小果的介绍是不是想要迫不及待地试一下了呢?

450K甲基化数据的免疫浸润分析www.biocloudservice.com/789/789.php

850K甲基化数据的免疫浸润分析www.biocloudservice.com/111_850k/111_850k.php


在这之后还有非常多进阶且美观的代码,小果的笔记本有所限制就不一一展示啦,把代码展示给大家,后续的结果图小果会打包成压缩文件夹在文章结尾分享给大家哦!

今天的分享就到这啦~小果希望大家也可以像学习popsicleR一样,先将代码运行成功,然后有能力的同学可以可以研究一下源代码,进行深层次的学习哦!

如果小伙伴有其他数据分析需求,可以尝试使用本公司新开发的生信分析小工具云平台,零代码完成分析,非常方便奥,云平台网址为:(http://www.biocloudservice.com/home.html),其中也包括了通路表达分析(http://www.biocloudservice.com/313/313.php),单细胞的基因共表达分析(http://www.biocloudservice.com/906/906.php)等各种小工具哦~,有兴趣的小伙伴可以登录网站进行了解。


点击“阅读原文”立刻拥有

↓↓↓