Python实现R包brainwaver中的compute.FDR函数
2024-09-07 23:57:48
FDR(false discovery rate),是统计学中常见的一个名词,翻译为伪发现率,其意义为是 错误拒绝(拒绝真的(原)假设)的个数占所有被拒绝的原假设个数的比例的期望值。
compute.FDR {brainwaver}使用图论进行可视化和参数化的多变量时间序列的基本小波分析
def computeFDR(parray, q=0.05):
if min(parray) < 0 or max(parray) > 1: raise ValueError("请检查P值是否在0~1范围内!\n") probs = sorted(parray)
l = len(probs) correct = sum([1/i for i in range(1, l + 1)]) fdr = [i/l*(q/correct) for i in range(1, l + 1)] sig = [0] * l for i in range(l):
if probs[i] <= fdr[i]:
sig[i] = 1 maxsig = max([i*j for i, j in zip(sig, range(1, l))]) return 0 if maxsig==0 else probs[maxsig]
参考资料
http://www.gipsa-lab.grenoble-inp.fr/~sophie.achard/brainwaver/compute.FDR.html
最新文章
- 移动端(h5)开发笔记
- [poj1679]The Unique MST(最小生成树)
- ArcGIS Server 增加缓存路径
- WCF初识
- HDOJ 1856
- css3常用动画效果集合01
- EBS R12中重新enable失效用户之后,丢失职责
- 深入学习netty系列(1)
- Android记录6--ViewPage+Fragment的使用例子
- 点击listview 的列头对其item进行自动排序
- Phonegap 安卓的自动升级与更新。当版本为4.0的时候
- SKTexture类
- JSP 9 大内置对象详解
- Oracle 导入、导出DMP(备份)文件
- Swift基础之:新的访问控制fileprivate和open
- MinGW安装与使用简介
- Jodd
- 微软必应地图加载错误:Uncaught TypeError: Microsoft.Maps.Location is not a constructor
- LNMP 下使用命令导出导入 MySQL 数据库
- p 最多两行 多的显示省略号
热门文章
- vim 的复制粘贴命令,以及使用寄存器来存放要复制的内容;
- 2.微信小程序-B站:需要先知道这些
- 开发人员需要具备的DBA技术
- python基础-循环语句for\嵌套循环
- Firefox 不知道如何打开此地址,因为协议 (javascrpit) 未和任何程序关联.
- Javascript-- jQuery动画篇(1)
- codeforces 651A A. Joysticks (模拟)
- Struts2 - 通过实现 Aware 接口访问 Web 资源
- 【遍历二叉树】01二叉树的前序遍历【Binary Tree Preorder Traversal】
- VC用MCI播放mp3等音乐文件