R与F分布(3) 分布检验

我们依然用Kolmogorov-Smirnov连续分布检验法来检验一个连续分布是否是服从F分布。

原假设为H0:数据集符合F分布
研究假设H1:样本所来自的总体分布不符合F分布。
令F0(x)表示预先假设的理论分布,Fn(x)表示随机样本的累计概率(频率)函数.

统计量D为: D=max|F0(x) - Fn(x)|

D值越小,越接近0,表示样本数据越接近F分布
p值,如果p-value小于显著性水平α(0.05),则拒绝H0

set.seed(1)
data<-rf(1000,1,1,2)
ks.test(data, "pf", 1,1,2)

One-sample Kolmogorov-Smirnov test
data: data
D = 0.0113, p-value = 0.9996
alternative hypothesis: two-sided

D值很小, p-value>0.05,不能拒绝原假设,所以数据集data符合df1=1, df2=1, ncp=2的F分布
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在Python,可以使用fitter库来进行分布拟合检验。首先,你需要安装fitter库,可以使用pip install fitter命令进行安装。接下来,你可以使用scipy库的stats模块生成一段模拟数据,例如使用stats.norm.rvs函数生成一个服从正态分布数据样本。然后,你可以使用fitter库的Fitter类来拟合数据样本的分布,可以通过指定distributions参数来选择要尝试的分布类型。最后,可以使用Fitter类的summary方法来查看拟合结果的摘要信息。 例如,下面的代码演示了如何使用fitter库进行分布拟合检验: ```python # 导入所需库 from scipy import stats import numpy as np from fitter import Fitter # 生成模拟数据 data1 = list(stats.norm.rvs(loc=0, scale=2, size=70000)) data2 = list(stats.norm.rvs(loc=0, scale=20, size=30000)) data = np.array(data1 + data2) # 使用fitter拟合数据样本的分布 f = Fitter(data, distributions=\['norm', 't', 'laplace'\]) f.fit() f.summary() ``` 在上述代码,我们生成了一个模拟数据样本,其包含了两个不同的正态分布。然后,我们使用Fitter类来拟合数据样本的分布,指定了要尝试的分布类型为正态分布(norm)、t分布(t)和拉普拉斯分布(laplace)。最后,我们使用summary方法来查看拟合结果的摘要信息。 希望这个例子能够帮助你理解如何使用fitter库进行分布拟合检验。 #### 引用[.reference_title] - *1* [Python拟合数据样本的分布](https://blog.csdn.net/sinat_23971513/article/details/111874082)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [概率统计Python计算:假设检验应用——分布拟合检验](https://blog.csdn.net/u012958850/article/details/118366266)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *3* [利用python通过拟合优度检验判别数据是否为正态分布](https://blog.csdn.net/qq_45019507/article/details/105842372)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值