蚌埠市网站建设_网站建设公司_全栈开发者_seo优化
2026/1/20 10:40:27 网站建设 项目流程

Python 数据统计完全指南:从入门到实战(2026 最新实用版)

Python 是当今数据统计与分析的首选语言,主要依赖以下核心库:

  • NumPy:数值计算基础
  • Pandas:数据清洗、结构化处理、描述统计
  • SciPy:高级统计函数、假设检验
  • Statsmodels:经典统计建模(回归、时间序列、ANOVA 等)
  • Seaborn / Matplotlib / Plotly:统计可视化

本文从零基础生产级实战,覆盖描述统计推断统计建模实战项目全链路。

1. 快速上手:环境与核心库安装(2026 推荐)

# 推荐使用 conda(最稳)或 pipconda create -n statspython=3.11conda activate stats condainstallpandas numpy scipy statsmodels matplotlib seaborn plotly# 或 pippipinstallpandas numpy scipy statsmodels matplotlib seaborn plotly

2. 描述统计(Descriptive Statistics)——理解数据“长什么样”

核心指标一览表
类别指标Python 实现方式含义与适用场景
中心位置均值 meandf['col'].mean()/np.mean(arr)算术平均,受极端值影响大
中位数 mediandf['col'].median()排序后中间值,抗极端值
众数 modedf['col'].mode()/scipy.stats.mode()出现频率最高值(适合分类数据)
离散程度标准差 stddf['col'].std()数据波动大小(样本标准差 ddof=1)
方差 vardf['col'].var()标准差的平方
四分位距 IQRdf['col'].quantile(0.75) - df['col'].quantile(0.25)中间 50% 数据范围,异常值检测常用
分布形状偏度 skewdf['col'].skew()/scipy.stats.skew()>0 右偏,<0 左偏,=0 对称
峰度 kurtosisdf['col'].kurt()>0 更尖,<0 更平(注意 Fisher 或 Pearson 定义)
整体概览describe()df.describe()/df.describe(include='all')一键输出数值/分类变量的统计摘要

实战小例子(Titanic 数据集)

importseabornassnsimportpandasaspd df=sns.load_dataset('titanic')print(df.describe(include='all'))# 全字段概览print(df['age'].median())# 中位数更鲁棒print("偏度:",df['fare'].skew())# 票价明显右偏print("峰度:",df['fare'].kurt())# 极端高峰

可视化描述统计(强烈推荐)

importseabornassnsimportmatplotlib.pyplotasplt# 单变量sns.histplot(df['age'],kde=True)# 直方 + 核密度plt.show()sns.boxplot(x='class',y='age',data=df)# 箱线图看分布 & 异常值plt.show()# 相关性热力图sns.heatmap(df.corr(numeric_only=True),annot=True,cmap='coolwarm')plt.show()

3. 推断统计(Inferential Statistics)——从样本推总体

核心概念与 Python 实现对照
任务统计方法Python 主要实现适用场景 & 注意事项
单样本均值检验单样本 t 检验scipy.stats.ttest_1samp(data, popmean=0)H0: μ = μ₀
两独立样本均值比较独立样本 t 检验scipy.stats.ttest_ind(group1, group2)假设方差相等(或用 Welch ttest)
配对样本均值比较配对 t 检验scipy.stats.ttest_rel(before, after)前后测量(如药物前后)
多组均值比较单因素 ANOVAscipy.stats.f_oneway(g1,g2,g3)statsmodels事后检验需 Tukey HSD
方差齐性检验Levene / Bartlettscipy.stats.levene(*groups)ANOVA 前置检验
相关性检验Pearson / Spearmanscipy.stats.pearsonr()/spearmanr()线性 / 单调关系
独立性检验卡方检验scipy.stats.chi2_contingency(table)分类变量间是否独立
正态性检验Shapiro-Wilk / KSscipy.stats.shapiro(data)小样本首选 Shapiro,大样本用 KS 或 Anderson
置信区间t 分布 / bootstrapstatsmodels.stats.api._tconfint_generic或手动bootstrap 更现代、无分布假设

经典实战:判断男女票价是否有显著差异

male_fare=df[df['sex']=='male']['fare'].dropna()female_fare=df[df['sex']=='female']['fare'].dropna()fromscipyimportstats t_stat,p_value=stats.ttest_ind(male_fare,female_fare,equal_var=False)# Welchprint(f"t ={t_stat:.3f}, p ={p_value:.4f}")# p < 0.05 → 拒绝原假设,男女票价有显著差异

4. 统计建模(Statsmodels 核心实战)

Statsmodels 是 Python 中最接近 R 的统计建模工具。

线性回归(OLS)完整流程

importstatsmodels.apiassmimportstatsmodels.formula.apiassmf# 公式法(最推荐)model=smf.ols('fare ~ age + sex + pclass + C(embarked)',data=df).fit()print(model.summary())# 一键输出:系数、p值、R²、F检验、诊断图等# 手动加常数项(当不用公式时)X=sm.add_constant(df[['age','pclass']])y=df['fare']model2=sm.OLS(y,X).fit()print(model2.summary())# 诊断图(非常重要!)sm.graphics.plot_regress_exog(model,'age')plt.show()

其他常见模型速查

  • 逻辑回归:smf.logit('survived ~ age + sex + pclass', data=df).fit()
  • 时间序列 ARIMA:sm.tsa.arima.ARIMA(data, order=(1,1,1)).fit()
  • 广义线性模型 GLM:smf.glm('... ~ ...', family=sm.families.Poisson()).fit()

5. 综合实战项目推荐(简易到进阶)

难度项目名称核心技术栈学习价值 & 简要步骤
★★泰坦尼克生存率影响因素分析Pandas + Seaborn + t检验 + 逻辑回归数据清洗 → 描述统计 → 假设检验 → 建模解释
★★★房价多因素回归分析与诊断Statsmodels OLS + 残差诊断 + VIF 多重共线性变量选择、变换、异常值处理、模型比较
★★★A/B 测试分析平台(模拟或真实数据)t检验 / Mann-Whitney / Bootstrap CI效应量 Cohen’d、功效分析 power
★★★★电商用户复购行为预测与 RFM 分析RFM + KMeans + 逻辑回归 / XGBoost商业指标 + 统计检验 + 机器学习对比
★★★★股票/加密货币波动率建模与 GARCHStatsmodels GARCH + 时间序列分解ACF/PACF、ARCH效应检验、预测

一句话总结学习路径(2026 版)

  1. 先熟练 Pandas describe() + Seaborn 画图(1 周)
  2. 掌握 SciPy 所有常见假设检验(2 周)
  3. 用 Statsmodels 跑 10 个回归模型 + 读懂 summary(2 周)
  4. 完成 1–2 个完整项目,从数据导入到结论报告(持续迭代)

如果你现在有一个具体数据集(比如 Kaggle 的某个 csv),或者想专注某个方向(金融统计、生物统计、A/B 测试、时间序列),可以告诉我,我直接给你定制代码 + 解释 + 论文/报告写法的完整方案。

统计的本质不是 p 值大小,而是用数据讲一个可信的故事。祝你学得顺、毕设/面试一把过!

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询