odd ratio置信区间的计算,你学会了吗?
欢迎关注”生信修炼手册”!
odd ratio称之为交叉乘积比,对于如下所示的数据
Allele | A | a |
---|---|---|
Case | a | b |
Control | c | d |
其计算公式如下
这里的A表示minor alllel, a表示major allel,以major allel为参照,用minor alllel的频数去除以参照的频数。然后用case组的比值除以control组的比值就可以得到odd ratio的值了。
那么odd ratio的置信区间如何计算呢?首先将odd raio值取log, 然后用log odd raio来进行分析,计算其标准误,公式如下
对于95%的置信区间,直接套用公式进行计算
以rs4970383
为例,显性模型中基因型统计如下
genotype | AA+Aa | aa |
---|---|---|
Case | 9 | 3 |
Control | 5 | 7 |
隐性模型中基因型统计如下
genotype | AA | Aa + aa |
---|---|---|
Case | 1 | 11 |
Control | 2 | 10 |
其OR值和对应的置信区间结果如下
在R中用上述公式进行计算,代码如下
可以看到,结果完全一致。对于其他的置信区间,只需要将95%对应的1.96换成其他系数即可。对于如下所示的钟型曲线,根据置信度计算两侧的概率累计值
然后求解对应的z值即可,计算方法如下
值得强调的是,逻辑回归中的回归系数就是log odd raio,所以对比plink逻辑回归输出的OR和BETA值,可以发现,将OR值取log之后就是BETA值
在R中进行验证
这里我们可以得出结论,OR值置信区间的计算实际上就是根据逻辑回归的回归系数,即log odd ratio推导出来的。
·end·
—如果喜欢,快分享给你的朋友们吧—
往期精彩
自己动手进行逻辑回归,你也可以!
逻辑回归or线性回归,傻傻分不清楚
没想到你是这个样子的置信区间
线性回归的这些细节,你都搞明白了吗?
3分钟掌握曼哈顿图的绘制
一文搞懂Q-Q plot图的含义
Cochran-Mantel-Haenszel检验在关联分析中的应用
plink中case/control关联分析细节解析
odd ratio值在关联分析中的含义
Cochran-Armitage趋势检验在关联分析中的应用
基因型填充
GWAS中的genotype imputation简介
基因型填充中的phasing究竟是什么
基因型填充前的质控条件简介
使用shapeit进行单倍型分析
gtool:操作genotype data的利器
使用IMPUTE2进行基因型填充
使用Beagle进行基因型填充
使用Minimac进行基因型填充
使用Eagle2进行单倍型分析
X染色体的基因型填充
文献解读|不同基因型填充软件性能的比较
Haplotype Reference Consortium:最大规模的单倍型数据库
Michigan Imputation Server:基因型填充的在线工具
CNV分析
aCGH芯片简介
aCGH芯片分析简介
基于SNP芯片进行CNV分析中的基本知识点
PennCNV:利用SNP芯片检测CNV
DGV:人类基因组结构变异数据库
dbvar:染色体结构变异数据库
DGVa:染色体结构变异数据库
CNVD:疾病相关的CNV数据库
DECIPHER:疾病相关的CNV数据库
全基因组数据CNV分析简介
使用CNVnator进行CNV检测
使用lumpy进行CNV检测
CNVnator原理简介
WES的CNV分析简介
XHMM分析原理简介
使用conifer进行WES的CNV分析
使用EXCAVATOR2检测WES的CNV
靶向测序的CNV分析简介
使用CNVkit进行CNV分析
DECoN:最高分辨率的CNV检测工具
TCGA
TCGA数据库简介
使用GDC在线查看TCGA数据
使用gdc-client批量下载TCGA数据
一文搞懂TCGA中的分析结果如何来
通过GDC Legacy Archive下载TCGA原始数据
使用GDC API查看和下载TCGA的数据
使用GDC下载TCGA肿瘤患者的临床信息
使用TCGAbiolinks下载TCGA的数据
使用TCGAbiolinks进行生存分析
使用TCGAbiolinks分析TCGA中的表达谱数据
使用TCGAbiolinks进行甲基化和转录组数据的联合分析
Broad GDAC:TCGA数据分析中心
使用cBioPortal查看TCGA肿瘤数据
UCSC Xena:癌症基因组学数据分析平台
GEPIA:TCGA和GTEx表达谱数据分析平台
TANRIC:肿瘤相关lncRNA数据库
SurvNet:基于网络的肿瘤biomarker基因查找算法
TCPA:肿瘤RPPA蛋白芯片数据中心
TCGA Copy Number Portal:肿瘤拷贝数变异数据中心
生存分析
生存分析详细解读
用R语言进行KM生存分析
使用OncoLnc进行TCGA生存分析
用R语言进行Cox回归生存分析
使用kmplot在线进行生存分析
肿瘤数据库
ICGC:国际肿瘤基因组协会简介
HPA:人类蛋白图谱数据库
Oncomine:肿瘤芯片数据库
ONGene:基于文献检索的肿瘤基因数据库
oncomirdb:肿瘤相关的miRNA数据库
TSGene:肿瘤抑癌基因数据库
NCG:肿瘤驱动基因数据库
mutagene:肿瘤突变频谱数据库
mSignatureDB:肿瘤突变特征数据库
GTEx:基因型和基因表达量关联数据库
肿瘤免疫和新抗原
Cancer-Immunity Cycle:肿瘤免疫循环简介
TMB:肿瘤突变负荷简介
肿瘤微环境:Tumor microenvironment (TME)简介
肿瘤浸润免疫细胞量化分析简介
使用EPIC预测肿瘤微环境中免疫细胞构成
TIMER:肿瘤浸润免疫细胞分析的综合网站
quanTIseq:肿瘤浸润免疫细胞定量分析
The Cancer Immunome Atlas:肿瘤免疫图谱数据库
肿瘤新抗原简介
TSNAdb:肿瘤新抗原数据库
使用NetMHCpan进行肿瘤新抗原预测分析
Hi-C数据分析
chromosome-territories:染色质疆域简介
chromosome conformation capture:染色质构象捕获技术
3C的衍生技术简介
解密Hi-C数据分析中的分辨率
A/B compartment:染色质区室简介
TAD:拓扑关联结构域简介
chromatin loops:染色质环简介
Promoter Capture Hi-C:研究启动子区染色质互作的利器
使用HiCUP进行Hi-C数据预处理
Juicer:Hi-C数据处理分析的利器
Juicer软件的安装详解
Juicebox:Hi-C数据可视化利器
Juicer实战详解
HiC-Pro:灵活的Hi-C数据处理软件
HiC-Pro实战详解
3D Genome Browser:Hi-C数据可视化工具
HiCPlotter:Hi-C数据可视化工具
3CDB:基于3C技术的染色质互作信息数据库
3DIV:染色质空间互作数据库
4D nucleome project:染色质三维结构研究必不可少的参考项目
3dsnp:SNP在染色质环介导的调控网络中的分布数据库
iRegNet3D:疾病相关SNP位点在三维调控网络中的作用
HiGlass:高度定制的Hi-C数据可视化应用
Hi-C Data Browser:Hi-C数据浏览器
使用FitHiC评估染色质交互作用的显著性
使用TADbit识别拓扑关联结构域
使用pyGenomeTracks可视化hi-c数据
hi-c辅助基因组组装简介
文献解读|使用hi-C数据辅助埃及伊蚊基因组的组装
chip_seq数据分析
Chip-seq简介
chip_seq质量评估之计算样本间的相关性
chip_seq质量评估之查看抗体富集效果
chip_seq质量评估之PCA分析
chip_seq质量评估之coverage分析
chip_seq质量评估之FRiP Score
chip_seq质量评估之文库复杂度
depth, bedgraph, bigwig之间的联系与区别
bigwig归一化方式详解
使用igvtools可视化测序深度分布
使用UCSC基因组浏览器可视化测序深度分布数据
使用deeptools查看reads分布特征
使用phantompeakqualtools进行cross correlation分析
blacklist regions:NGS测序数据中的黑名单
MACS:使用最广泛的peak calling软件之一
MACS2 peak calling实战
使用SICER进行peak calling
使用HOMER进行peak calling
peak注释信息揭秘
PAVIS:对peak区域进行基因注释的在线工具
使用UPORA对peak进行注释
使用GREAT对peak进行功能注释
annoPeakR:一个peak注释的在线工具
使用ChIPpeakAnno进行peak注释
使用ChIPseeker进行peak注释
使用PeakAnalyzer进行peak注释
使用homer进行peak注释
利用bedtools预测chip_seq数据的靶基因
motif
关于motif你需要知道的事
详解motif的PFM矩阵
详解motif的PWM矩阵
使用WebLogo可视化motif
使用seqLogo可视化motif
使用ggseqlogo可视化motif
MEME:motif分析的综合性工具
使用MEME挖掘序列中的de novo motif
使用DREME挖掘序列中的de novo motif
使用MEME-ChIP挖掘序列中的de novo motif
chip_seq数据库
ENCODE project项目简介
FactorBook:人和小鼠转录因子chip_seq数据库
ReMap:人类Chip-seq数据大全
IHEC:国际人类表观基因组学联盟
Epifactors:表观因子数据库
GTRD:最全面的人和小鼠转录因子chip_seq数据库
ChIP-Atlas:基于公共chip_seq数据进行分析挖掘
Cistrome DB:人和小鼠的chip_seq数据库
chipBase:转录因子调控网络数据
unibind:human转录因子结合位点数据库
chip_seq在增强子研究中的应用
DENdb:human增强子数据库
VISTA:人和小鼠的增强子数据库
EnhancerAtlas:人和小鼠的增强子数据库
FANTOM5:人类增强子数据库
TiED:人类组织特异性增强子数据库
HEDD:增强子疾病相关数据库
HACER:human增强子数据库
SEdb:超级增强子数据库简介
dbSUPER:人和小鼠中的超级增强子数据库
dbCoRC:核心转录因子数据库
使用ROSE鉴定超级增强子
18年文章目录
2018年推文合集
扫描下方二维码,关注我们,解锁更多精彩内容!
生物信息入门
只差这一个
公众号
odd ratio置信区间的计算,你学会了吗?相关推荐
- 机器的思考故事——计算机能学会牛顿第二定律吗?
机器的思考故事--计算机能学会牛顿第二定律吗? 本项目从"牛顿第二定律"入手,介绍机器学习的思考过程,以及在过程中如何确定模型参数,模型三个关键部分(假设.评价.优化)该如何应用. ...
- R语言 总体均数的95%置信区间、几何均数(GMT)及其置信区间的计算
1.求总体均数的95%置信区间,可根据自身需要更该alpha值,计算相应的置信区间 ---------------------------------------------------------- ...
- Z检验、T检验下 P-value 和置信区间的计算
目录 1.置信区间的计算 1.1 总体方差已知 1.2 总体方差未知 2.计算 P-Value 2.1 总体方差已知 2.2 总体方差未知 1.置信区间的计算 根据总体分布(T分布或者Z分布)和规定的 ...
- python 散点图 置信区间_python计算置信区间
python 置信区间 置信区间是指由 样本统计量 所构造的总体参数的估计区间. 这句话也就是说 ( 这里统计量一般指均值 ) 利用样本均值来估计总体均值的可靠程度, 这个"可靠程度&quo ...
- 置信系数 表达了置信区间的计算_LM算法推导:阻尼法与置信域法
应用领域 LM 算法用于解决非线性最小二乘问题,问题定义如下: LM 算法有两种推导和实现,一种是算法发明者使用的 阻尼法,一种是后来学者补充的 置信域法.这里分别作出推导. 阻尼法推导 Damped ...
- 效率值计算机公式,工作效率提升怎么计算(学会十个Excel公式)
IF语法: 根据条件判断返回指定的值 COUNTIF语法: 统计满足给定条件的单元格的个数 1.查找重复内容 在C4输入公式: =IF(COUNTIF($B$4:$B$8,$B4)>1,&quo ...
- 置信系数 表达了置信区间的计算_参数估计05---置信区间, 置信限
区间估计 引言 前面,我们讨论了参数点估计.它是用样本算得的一个值去估计未知参数.但是,点估计值仅仅是未知参数的一个近似值,它没有反映出这个近似值的误差范围,使用起来把握不了.区间估计正好弥补了点估计 ...
- 物理化学 焓变的计算,学会状态的转变
- 一文掌握inbreeding coefficient近交系数的计算
欢迎关注"生信修炼手册"! inbreeding coefficient,中文翻译为近交系数,近婚系数,近亲交配系数等等,用大写字母F表示.要理解这个概念,首先要搞清楚" ...
- TCGA肿瘤数据分析专题
欢迎关注"生信修炼手册"! 癌症作为人类健康的头号杀手,其研究的意义不言而喻.目前世界范围内已经有大量的肿瘤相关数据,鉴于公共数据库的数据挖掘成为一种趋势. GEO是一个国际化的开 ...
最新文章
- sqlite在c++中的使用方法
- Jenkins部署Windows UI自动化的调度权限问题
- gdb工作原理(二)
- oracle操作字符串:拼接、替换、截取、查找
- 美酒节成就及任务攻略指引
- SQLServer2014 安装错误:等待数据库引擎恢复句柄失败
- 【机器学习】Lasso回归(L1正则,MAP+拉普拉斯先验)
- redmine 无法登录 mysql 服务器_Redmine 数据库连接错误
- /proc 虚拟文件系统(实例)
- 最大数max(x,y,z)(信息学奥赛一本通-T1152)
- 服务器cpu支持DDR4,迟来的惊喜 AMD新处理器支持DDR4内存
- 三维重建_基于RGB-D相机的三维重建总览(静态动态)
- 计算机系统管理程序,计算机系统管理操作程序
- DIY专门制作PCB的雕刻机--所需软件
- 白菜cms双端影视APP源码
- 【5G核心网】 GTP协议
- 读计算机视觉综述做的脑图(11.09更新)
- 轻松解决Tomcat启动慢的问题,只需一行代码
- 母亲节是在每年五月份的第二个星期日,给定年份,求出当年母亲节的日期
- django教程day05
热门文章
- 漏洞系列一一看我一招征服漏洞SSRF
- Prefix-Tuning: Optimizing Continuous Prompts for Generation
- C 语言面试题大汇总
- 考考你、智商题 小明借爸爸500元 又借妈妈500元 买双鞋 970元 还剩30元 还给爸爸10元...
- win10小娜助手无法搜索本地应用怎么办
- 5000在线的服务器选择,支持负载5000人10000人在线的服务器的配置及其网络要求!!...
- 【一本通】1148-连续出现的字符
- 3大奇葩排序之猴子算法
- Variable used in lambda expression should be final or effectively final报错解决方案
- windows未能启动计算机,电脑开机windows未能启动是怎么回事