Python统计分析——参数估计

2024-02-05 07:20

本文主要是介绍Python统计分析——参数估计,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

参考资料:用python动手学统计学

       所谓参数就是总体分布的参数。

1、导入库

# 导入用于数值计算的库
import numpy as np
import pandas as pd
import scipy as sp
from scipy import stats
# 导入用于绘图的库
from matplotlib import pyplot as plt
import seaborn as sns
sns.set()

2、导入案例材料

data=np.array([4.352982,3.735304,5.944617,3.798326,4.087688,5.265985,3.272614,3.526691,4.150083,3.736104])
data

3、点估计

      直接指定总体分布的参数为某一值的估计方法叫作点估计。

      我们使用样本均值作为总体均值的点估计量,所以只需要计算出样本的均值就可以完成估计。 这看起来很简单,但要注意,正因为样本均值具有无偏性和一致性,它才可以作为总体均值的估计值。

      同理,我们使用样本的无偏方差作为总体方差的估计值。

python实现步骤如下;

mu=np.mean(data)
sigma_2=np.var(data,ddof=1)
print('总体均值的估计值为:',mu)
print('总体方差的估计值为:',sigma_2)

结果如下:

4、区间估计

      估计值具有一定范围的估计方法叫作区间估计。我们使用概率的方法计算这个范围。因为估计值是一个范围,所以可以引入估计误差。估计误差越小,区间估计的范围越小;样本容量越大,区间的范围越小。

      置信水平,是表示区间估计的区间可信度的概率。例如95%、99%都是常用的置信水平。二满足某个置信水平的区间叫作置信区间。对于同一组数据,置信水平越大,置信区间就越大。

       置信区间的计算如下:

\bar{x}-t_{1-\alpha/2}\times \sigma/\sqrt{n}<\mu<\bar{x}-t_{\alpha/2}\times \sigma/\sqrt{n}

python实现步骤如下:

# 自由度
df=len(data)-1
sigma=np.std(data,ddof=1)
se=sigma/np.sqrt(len(data))
interval=stats.t.interval(confidence=0.95,df=df,loc=mu,scale=se)
interval

结果如下:

与公式计算结果一致,如下图:

5、python函数参数介绍:

5.1 scipy.stats.t.interval()用于获取t分布的置信区间,参数介绍如下:

(1)confidence,用于设置置信水平。可以用列表的形式设置多个置信水平。如下:

(2)df为自由度,loc为样本均值,scale为样本均值的标准误。

5.2 scipy.stats.t.ppf()用于获取t分布的百分位数。

(1)q,小数形式,设置需要获取百分数对应的百分位

(2)df,设置自由度。

6、决定置信区间大小的因素

6.1 样本方差越大,置信区间越大

将样本标准差变为原来的10倍进行验证。

5.2 样本容量越大,样本均值就越可信,进而置信区间就越小

将样本容量为原来的10倍进行验证。

 5.3 置信水平越大,置信区间就会越大。

将置信水平调整为99%,进行验证。

 6、置信区间结果的解读

       如上图所示,置信水平为95%的置信区间,表示所得到的该区间包含真正的总体均值这一参数的概率为95%。

      下面用2万次的抽样结果,对置信区间的置信水平进行验证。

# 执行2万次求95%置信区间的操作
# 如果置信区间包含总体均值(本例设置为4),就为True
np.random.seed(1) # 设置随机种子,用于复现结果
# 设置数组用于存放置信区间是否包含总体均值的判断结果
be_included_array=np.zeros(20000,dtype='bool')
# 设置正态总体
pop=stats.norm(loc=4,scale=0.8)
# 完成2万次的样本抽取并对置信区间是否包含总体均值进行验证
for i in range(0,20000):sample=pop.rvs(size=10)df=len(sample)-1mu=np.mean(sample)std=np.std(sample,ddof=1)se=std/np.sqrt(len(sample))interval=stats.t.interval(0.95,df=df,loc=mu,scale=se)if(interval[0]<4 and interval[1]>4):be_included_array[i]=True# 汇总True的占比
sum(be_included_array)/len(be_included_array)

由模拟结果可以看出,总体均值包含在置信区间的比例约为95%,与置信水平基本一致。

这篇关于Python统计分析——参数估计的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/680140

相关文章

使用Python和OpenCV库实现实时颜色识别系统

《使用Python和OpenCV库实现实时颜色识别系统》:本文主要介绍使用Python和OpenCV库实现的实时颜色识别系统,这个系统能够通过摄像头捕捉视频流,并在视频中指定区域内识别主要颜色(红... 目录一、引言二、系统概述三、代码解析1. 导入库2. 颜色识别函数3. 主程序循环四、HSV色彩空间详解

一文深入详解Python的secrets模块

《一文深入详解Python的secrets模块》在构建涉及用户身份认证、权限管理、加密通信等系统时,开发者最不能忽视的一个问题就是“安全性”,Python在3.6版本中引入了专门面向安全用途的secr... 目录引言一、背景与动机:为什么需要 secrets 模块?二、secrets 模块的核心功能1. 基

python常见环境管理工具超全解析

《python常见环境管理工具超全解析》在Python开发中,管理多个项目及其依赖项通常是一个挑战,下面:本文主要介绍python常见环境管理工具的相关资料,文中通过代码介绍的非常详细,需要的朋友... 目录1. conda2. pip3. uvuv 工具自动创建和管理环境的特点4. setup.py5.

Python常用命令提示符使用方法详解

《Python常用命令提示符使用方法详解》在学习python的过程中,我们需要用到命令提示符(CMD)进行环境的配置,:本文主要介绍Python常用命令提示符使用方法的相关资料,文中通过代码介绍的... 目录一、python环境基础命令【Windows】1、检查Python是否安装2、 查看Python的安

Python UV安装、升级、卸载详细步骤记录

《PythonUV安装、升级、卸载详细步骤记录》:本文主要介绍PythonUV安装、升级、卸载的详细步骤,uv是Astral推出的下一代Python包与项目管理器,主打单一可执行文件、极致性能... 目录安装检查升级设置自动补全卸载UV 命令总结 官方文档详见:https://docs.astral.sh/

Python并行处理实战之如何使用ProcessPoolExecutor加速计算

《Python并行处理实战之如何使用ProcessPoolExecutor加速计算》Python提供了多种并行处理的方式,其中concurrent.futures模块的ProcessPoolExecu... 目录简介完整代码示例代码解释1. 导入必要的模块2. 定义处理函数3. 主函数4. 生成数字列表5.

Python中help()和dir()函数的使用

《Python中help()和dir()函数的使用》我们经常需要查看某个对象(如模块、类、函数等)的属性和方法,Python提供了两个内置函数help()和dir(),它们可以帮助我们快速了解代... 目录1. 引言2. help() 函数2.1 作用2.2 使用方法2.3 示例(1) 查看内置函数的帮助(

Python虚拟环境与Conda使用指南分享

《Python虚拟环境与Conda使用指南分享》:本文主要介绍Python虚拟环境与Conda使用指南,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、python 虚拟环境概述1.1 什么是虚拟环境1.2 为什么需要虚拟环境二、Python 内置的虚拟环境工具

Python实例题之pygame开发打飞机游戏实例代码

《Python实例题之pygame开发打飞机游戏实例代码》对于python的学习者,能够写出一个飞机大战的程序代码,是不是感觉到非常的开心,:本文主要介绍Python实例题之pygame开发打飞机... 目录题目pygame-aircraft-game使用 Pygame 开发的打飞机游戏脚本代码解释初始化部

Python pip下载包及所有依赖到指定文件夹的步骤说明

《Pythonpip下载包及所有依赖到指定文件夹的步骤说明》为了方便开发和部署,我们常常需要将Python项目所依赖的第三方包导出到本地文件夹中,:本文主要介绍Pythonpip下载包及所有依... 目录步骤说明命令格式示例参数说明离线安装方法注意事项总结要使用pip下载包及其所有依赖到指定文件夹,请按照以