r 语言 怎么计算fdr

Posted 2023-04-18

命令如下：
p<-c(0.001315146,0.001236789,0.001229388,0.000889006,0.000876515,0.000578648,0.000565415,0.000536447,0.000517434,0.000487215,0.000364518,0.000364518,0.000247193,7.93E-05);
p.adjust(p,method=”fdr”,n=length(p));
这样得到的修正后的fdr值为：
0.0013151460,0.0013151460,0.0013151460,0.0011314622,0.0011314622,0.0009001191,0.0009001191,
0.0009001191,0.0009001191,0.0009001191,0.0009001191,0.0009001191,0.0009001191,0.0009001191
可以看到，修正前后面几个P值并不相等，但是修正后的fdr后面几个都变一样的值了。
而且根据fdr的定义，用命令length(p)*p/rank(p)计算出来的结果也和用命令p.adjust(“fdr”)计算出来的结果不一样。 参考技术A p.adjust(p,method="fdr",length(p))

r语言绘制核密度图怎么计算重叠

核密度图是一种用于可视化数据分布的统计图形。它通过在数据点周围绘制一个曲线来估计数据分布的密度。重叠是指多个核密度图中，不同密度区域之间的重叠度。计算重叠最常用的方法是Kolmogorov-Smirnov检验（K-S检验）。K-S检验是一种统计检验，用于比较两个样本的分布是否相同。K-S检验的基本思想是比较两个样本的累积分布函数（CDF），如果两个样本的CDF之间的最大差值越大，则两个样本之间的差异越大。所以K-S检验可以用来估计两个核密度图之间的重叠度。 参考技术A 使用函数dnorm()来计算每个点的核密度值。绘制重叠核密度图，首先需要确定每个点的核密度值，然后以相同的方式绘制多个核密度图，最后将这些图叠加在一起。