aic值检验 p值_参数估计和假设检验

本文深入探讨统计推断中的参数估计和假设检验。参数估计涵盖点估计和区间估计,如矩估计法、最小二乘法、极大似然法;假设检验包括T检验、非参数检验,如秩和检验、中位数评分检验,以及分布拟合检验,如正态性检验。通过实例解析了T检验的各种类型,并介绍了非参数检验和分布拟合检验的应用。
摘要由CSDN通过智能技术生成

2b20bb6a68ad000d779f78a1f59bab84.png

前言:所谓统计推断,就是利用样本所提供的信息对总体的某些统计特征进行估计或者判断,进而认识总体。统计推断分为两大类:参数估计,假设检验。

1 .参数估计

参数估计和假设检验是统计推断的基本内容,几乎所有统计建模的PROC步都会设计参数估计以及相对应的假设检验。假设总体

的分布函数的类型已知,但其中一个或者多个参数未知,那么就需要对这些未知的参数做出合理的估计,并且对估计做出评价,这就是
参数估计。

参数估计分为 点估计和区间估计。

1.1 点估计

也称为定值估计。以一个单一的取值,去近似的作为未知的参数的估计值。

例如样本的均值

可以用来估计总体的均值
,样本的标准差
可以用来估计总体的标准差

点估计的方法:

矩估计法,最小二乘法,极大似然法

1.2 区间估计

对于未知参数

,我们除了关注点估计
外,还希望估计出一个范围,并希望知道参数
落在这个范围的可信程度,这就是区间估计。这样的区间即参数的置信区间。

1.2.1 置信区间

就是一个包含统计量值的取值范围,并且这个范围在一定的置信水平下包含参数

的真实值。

数学表达为:

若有一分布

,
为未知参数,对于给定值
,若由样本计算出两个统计量
,满足
,则随机区间
置信水平为
的置信区间,
是置信水平,
是显著性水平,
为置信上下限。

例如,95%置信区间:从总体中重复抽取100个相同样本容量的样本,并且计算了100个置信区间,那么这100个置信区间有95个包含了总体的均值。

我们都希望提高置信水平,但是在样本容量一定的情况下,提高置信水平,势必会使置信区间变宽,估计的精度降低。

1.2.2 均值的置信区间

均值标准误差:

刻画样本均值对于总体均值的变异程度。

例子:计算sashelp.fish的bream的95%置信区间。

proc means data=sashelp.fish mean std stderr clm maxdec=2 n;
where species="Bream";
var height;
title "95% confidence interval for bream";
run;

317dc349f0ed957a76ecb0dc017faa45.png
  • stderr:计算标准误差值
  • clm:计算均值的置信区间
  • 添加 alpha=选项可以计算不同置信水平的置信区间。

本个案例中,可以认为区间(14.51,15.86)包含总体均值的可能性为95%,区间的宽度很窄,可以认为样本均值是总体均值的一个比较准确的估计。

2.假设检验

在总体的分布函数只知其形式,但不知其参数的情况下,或者对总体分布完全未知的情况下,为了推断总体的某些未知特征,现提出某些关于总体的假设,然后根据样本采用适当的方法,对所提出的假设作出接受或拒绝的决策,这一过程就是假设检验。

2.1 基本原理

例子:假如我手上有一枚硬币,一般的硬币都是均匀的,如果我认为它是不均匀的,我就必须给出足够的理由来证明。为此我制定了一个原则:连续抛掷5次,5次都朝上或者朝下,则认为该硬币不均匀,否则&#x

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个利用 GARCH 类模型来分析多支股票的对数收益率数据边缘分布的 R 代码示例: ```R library(tseries) library(fGarch) library(nortest) # 读取数据 data <- read.csv("stock_returns.csv", header = TRUE) # 建立GARCH(1,1)模型 garch.spec <- ugarchspec(mean.model = list(armaOrder = c(0,0)), variance.model = list(garchOrder = c(1,1)), distribution.model = "std") garch.fit <- ugarchfit(spec = garch.spec, data = data, solver = "hybrid") # 输出参数结果、对数似然AIC summary(garch.fit) # 进行K-S检验 residuals <- residuals(garch.fit, standardize = TRUE) ks.test(residuals, "pnorm") # 绘制残差分布图 hist(residuals, freq = FALSE, main = "GARCH(1,1) Residuals", xlab = "Residuals", ylim = c(0, 0.5)) curve(dnorm(x, mean = mean(residuals), sd = sd(residuals)), col = "red", add = TRUE) ``` 在上述代码中,我们首先使用`read.csv`函数读取多支股票的对数收益率数据,并使用`ugarchspec`函数建立了一个GARCH(1,1)模型。接着,我们使用`ugarchfit`函数对模型进行参数估计,并使用`summary`函数输出了参数结果、对数似然AIC。最后,我们使用`ks.test`函数进行了K-S检验,并绘制了残差分布图。 需要注意的是,在实际应用中,我们需要对每支股票的收益率序列进行单独的建模和分析,以获取更为准确的边缘分布特征。同时,我们还需要根据实际情况选择适当的GARCH模型和Copula函数,并进行多重比较和拟合诊断,以确保模型的可靠性和准确性。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值