前言

此篇文章对于贝叶斯假设检验不做定义以及公式的推导，仅仅是想通过实例R语言操作来理解贝叶斯假设检验。

一、贝叶斯假设检验的优点

简单，无需抽样分布和显著性水平
可以推广至多重检验

二、R语言实例操作

1.题目叙述

Measurements~ N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2), μ \mu μ is true weight and σ \sigma σ = 3 pounds.
Interested in accessing if his true weight >175 pounds, test : H 0 : μ ≤ 175 H_0: \mu \leq175 H0:μ≤175, H 1 : μ > 175 H_1: \mu >175 H1:μ>175.
Normal prior N(170,25)
Measurement:182, 172, 173,176,176,180,173,174,179,175

2. 先验概率

代码如下（示例）：

pmean=170;pvar=25
probH=pnorm(175,pmean,sqrt(pvar))
probA=1-probH
prior.odds=probH/probA
prior.odds#5.302974

2.后验概率

weights=c(182, 172, 173, 176, 176, 180, 173, 174, 179, 175)
xbar=mean(weights)
sigma2=3^2/length(weights)
post.precision=1/sigma2+1/pvar
post.var=1/post.precision
post.mean=(xbar/sigma2+pmean/pvar)/post.precision
c(post.mean,sqrt(post.var))#175.7915058   0.9320546
post.odds=pnorm(175,post.mean,sqrt(post.var))/
(1-pnorm(175,post.mean,sqrt(post.var)))
post.odds# 0.2467017

3. 贝叶斯因子

BF = post.odds/prior.odds
BF#0.04652139

4. H 0 H_0 H0的后验概率

postH=probH*BF/(probH*BF+probA)
postH#0.1978835

5. 结论

Conclude: unlikely that his weight is at most 175 pounds.

贝叶斯统计R语言操作3——假设检验相关推荐

贝叶斯统计R语言操作1——共轭先验分布
文章目录一.共轭先验分布定义二.正态均值的共轭先验分布 1.公式理论 2.Example: Normal mean 3.R code 一.共轭先验分布定义设 π(θ)\pi(\theta)π(θ ...
二元置信椭圆r语言_一般加性模型的简介、应用举例及R语言操作
一般加性模型的简介.应用及R语言操作举例前文在"平滑回归举例"中,提到当未知自变量和响应变量间的关系,难以选择合适的参数模型描述二者间的响应状态,或者期望探索二者可能的响应曲线形式 ...
R语言操作excel文件的数据
简述令人惊讶的是,网上的教程关于r语言处理excel数据的都是在说如何读.. 下面以用xlsx这个包为例安装这个包,并读取数据 [两行命令]R语言读取excel数据操作数据其实,导入进来之后, ...
R语言操作PDF文件的包
背景有过文章投稿经验的朋友们知道,文章中的图片排版是个很繁琐的过程.涉及到多个图片的组合跟字体字号的调整.刚开始自己没有这方面经验,对输出的结果图格式比较随意:*.jpg,*.png等等.直到后面需 ...
时间序列R语言操作2——白噪声和随机游走模型
文章目录一.白噪声 1.白噪声是什么? 2.白噪声的性质 3.样本自相关函数 4.白噪声检验二.随机游走模型 1.简单随机游走序列 2.R语言 3.实例说明 1.醉汉回家 2.久赌必输一.白噪声 ...
R语言操作pdf文档
今天给大家介绍一个可以读取pdf文件信息的R包pdftools.此包基于Poppler库(https://poppler.freedesktop.org/)进行解析pdf文件.在这里我们就不深入剖析这 ...
R语言数据分析笔记——t检验（含正态性检验和方差齐性检验在SPSS和R语言中的操作t检验（单样本、双独立样本、配对样本）在Excel、SPSS、R语言中的操作）
前言:本文为个人学习笔记,为各大网站上的教学内容之综合整理,综合整理了①假设分析的基础知识.②正态性检验和方差齐性检验在SPSS和R语言中的操作.③t检验(单样本.双独立样本.配对样本)在Excel. ...
R语言-单个总体均值的t检验
#3.2.1 单个总体均值的t检验 # (1)什么是检验? # 检验(test)是统计学中最重要的概念之一,在科学研究和实际业务中都有着广泛的应用.用一句话来概括就是:人们希望通过掌握的数据和其他背景 ...
在r中弄方差分析表_医学统计与R语言： qvalue
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发 (FalseDiscoveryRate(FDR)=Expected(FalsePositive/(FalsePositive+TruePos ...

贝叶斯统计R语言操作3——假设检验

文章目录

前言

一、贝叶斯假设检验的优点

二、R语言实例操作

1.题目叙述

2. 先验概率

2.后验概率

3. 贝叶斯因子

4. H 0 H_0 H0的后验概率

5. 结论

贝叶斯统计R语言操作3——假设检验相关推荐

最新文章

热门文章

贝叶斯统计R语言操作3——假设检验

文章目录

前言

一、贝叶斯假设检验的优点

二、R语言实例操作

1.题目叙述

2. 先验概率

2.后验概率

3. 贝叶斯因子

4. H 0 H_0 H0​的后验概率

5. 结论

贝叶斯统计R语言操作3——假设检验相关推荐

最新文章

热门文章

4. H 0 H_0 H0的后验概率