【R方差分析】蛋白质表达量多组比较
2024-08-28 14:19:44
初始数据类似:
蛋白质组数据虽不是严格的正态分布,但目前最常用的检验方法还是T检验(两组比较)和方差分析(多组比较)。这个话题值得深究,这里不展开。
主要是求多个蛋白的Pvalue值或FDR,用于差异筛选。
Pvalue <- c()
type<-factor(c(rep("S01CC",3),rep("S11SC",3),rep("S12CC",3),rep("S12SC",3)))
for(i in 1:nrow(exp)){
if(sum(exp[i,1:3])==0 && sum(exp[i,4:6])==0 && sum(exp[i,7:9])==0 && sum(exp[i,9:12])==0){
Pvalue[i] <- "NA"
}else{
y=aov(as.numeric(exp[i,1:12]) ~ type)
Pvalue[i]<-summary(y)[[1]][,5][1]
# Pvalue[i]<-summary(y)[[1]][,"Pr(>F)"][1]
}
}
FDR=p.adjust(Pvalue, "BH")
out<-cbind(exp,Pvalue,FDR)
write.csv(out,"anova_out.csv",row.names = T)
最新文章
- java 堆栈 理解
- Download Free Oracle Reports Building Guide eBook
- 去除bootstrap模态框半透明阴影
- Phabricator部署手册
- vim基础命令
- 同时使用Twitter nlp 和stanford parser的解决方法
- unix PS命令和JPS命令的区别
- hibernate的懒加载问题
- Linux Shell 命令
- QT5.3.2在ARM上的移植
- Java 图片转换为字符图 CharMaps (整理)
- List、ArrayList、Vector及map、HashTable、HashMap分别的区别
- 操,escape sequence的输入方法我以前找过一次,这次又忘了,又找了一次,记下来,
- [Swift]LeetCode674. 最长连续递增序列 | Longest Continuous Increasing Subsequence
- Android 7.0 适配
- 知识点:SQL中char、varchar、text区别
- 洛谷P2054 [AHOI2005]洗牌(扩展欧几里德)
- sql server连接字符串与tcp/ip开启
- 当inline-block或者float失效的时候怎么弄
- IDEA Intellij 打开springboot项目 配置文件无法出现输入提示