PCA图显示分组无差异,怎么办?
最近接到粉丝提问,感兴趣的数据集做差异分析,发现很勉强,不好把握。因为我以前在生信技能树写过教程: PCA都分不开的两个组强行找差异是为何 ,所以征求我的意见。但是我很忙啊,所以就把这个数据集安排给了实习生和学徒。我一直强调,做表达矩阵分析一定要有三张图,见: 你确定你的差异基因找对了吗? ,基本上看到这3张图,就明白问题出在哪里了。 粉丝求助的数据集介绍如下: 小白菊内酯处理人胆管癌细胞的基因变化 药物处理:小白菊内酯 Platform:GPL6102 Series:GSE22633 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE22633 第一步是下载数据 rm(list = ls()) ## 魔幻操作,一键清空~ options(stringsAsFactors = F)#在调用as.data.frame的时,将stringsAsFactors设置为FALSE可以避免character类型自动转化为factor类型 # 注意查看下载文件的大小,检查数据 f='GSE22633_eSet.Rdata' library(GEOquery) # 这个包需要注意两个配置,一般来说自动化的配置是足够的。 #Setting options('download.file.method.GEOquery'='auto')