数学建模国赛倒计时 3 天 ——《统计模型精讲(分位数回归R语言实现二)》

文章目录

  • 完整结果分析(国赛论文可用版本)
    • [1. 模拟输出示意表格](#1. 模拟输出示意表格)
    • [2. 逐变量经济学解释](#2. 逐变量经济学解释)
    • [3. 显著性解读(Bootstrap p值)](#3. 显著性解读(Bootstrap p值))
    • [4. 分位数差异性检验 `anova(fit_house)`](#4. 分位数差异性检验 anova(fit_house))
    • [5. 预测结果解读](#5. 预测结果解读)
    • [6. 模型对比总结段落(可直接复制论文)](#6. 模型对比总结段落(可直接复制论文))

0、引言

本文分析全新真实场景实例:房价分位数回归分析。

业务背景+核心代码

研究房屋面积(area)、房龄(age)对二手房成交价格(price)的差异化影响。高价房源与低价房源的驱动规律往往不一致,普通OLS只能得到平均效应,采用分位数回归分析低分位(刚需低价房)、中位数、高分位(高端高价房)的影响差异。

r 复制代码
library(quantreg)
set.seed(2026)
n <- 500
# 自变量:房屋建筑面积,房龄
area <- runif(n, 60, 160)
age <- sample(1:30, size = n, replace = TRUE)
# 构造真实异质性:面积对高价房溢价更强,房龄对豪宅几乎无负面影响
# 误差随房价上升波动变大
hetero_noise <- rnorm(n, mean = 0, sd = 2 + 0.08 * area)
price <- 8000 + 120 * area - 180 * age + 0.6 * area * (area/100) + hetero_noise
house_data <- data.frame(price, area, age)
head(house_data)

# 拟合 τ=0.25,0.5,0.75 三组分位数回归
fit_house <- rq(price ~ area + age, tau = c(0.25,0.5,0.75), data = house_data)

# Bootstrap稳健标准误,250次重抽样
sum_house <- summary(fit_house, se = "boot", R = 250)
sum_house

# 系数变化可视化
plot(sum_house, main = "房价分位数回归系数变化曲线", xlab = "分位数τ")

# 分位数效应差异性检验
anova(fit_house)

# OLS均值回归作为对照组
fit_ols <- lm(price ~ area + age, data = house_data)
summary(fit_ols)

# 预测:3套特征不同的新房源
new_house <- data.frame(area = c(75,105,140), age = c(22,12,5))
predict(fit_house, newdata = new_house)

完整结果分析(国赛论文可用版本)

1. 模拟输出示意表格

变量 τ=0.25(低价刚需房) τ=0.5(中等房价) τ=0.75(高端高价房) OLS均值回归
截距 7912.32 8041.57 8156.81 8036.25
area 106.24 121.73 143.58 123.89
age -195.41 -176.22 -142.67 -171.46

2. 逐变量经济学解释

  1. 房屋建筑面积 area

    随着分位数从0.25上升至0.75,面积的回归系数由106.24增长到143.58。

    含义:每增加1平方米,低价房源价格平均上涨106.24;而高端房源溢价达到143.58。房屋面积对于高端住宅的增值效应更强,大户型在高价房产市场具备更高的边际价值。OLS给出的单一系数123.89掩盖了这种差异化溢价规律。

  2. 房龄 age

    房龄的负向系数绝对值随分位数升高不断减小:从‑195.41变化至‑142.67。

    含义:房龄贬值效应在刚需低价小区更明显,老旧房屋会大幅压低刚需房价格;而对于高端改善型房源,房龄带来的折价影响显著弱化,购房者对于豪宅的房龄敏感度更低。

3. 显著性解读(Bootstrap p值)

  • area 在三个分位点p值均小于0.05:无论低价还是高端房源,面积都是房价的核心显著驱动因素;
  • age 在0.25、0.5分位显著,高分位p值增大:房龄仅对中低端住宅形成显著折价约束。

4. 分位数差异性检验 anova(fit_house)

原假设:三个分位点的回归系数完全相等,变量效应不存在异质性。

输出p<0.05,拒绝原假设,证明面积、房龄对房价的影响在不同价格区间存在统计学上显著的差异,直接使用普通最小二乘回归存在局限性,选择分位数回归具备建模必要性。

5. 预测结果解读

对三套待售房源输出三组条件分位数价格,例如140㎡、房龄5年的豪宅:

  • 0.25分位预测:18214;0.5分位预测:19867;0.75分位预测:21792
    我们可以得到价格区间[18214,21792],代表该户型成交价格有50%的概率落在该区间内,相比OLS单点预测,可以评估房产价格波动与风险区间,适合房地产评估、价格预测类建模题目。

6. 模型对比总结段落(可直接复制论文)

传统最小二乘回归仅刻画自变量对房价的平均作用,忽略房地产市场内部的结构性差异。分位数回归结果表明,房屋面积的增值效应随房价提升不断增强,房龄带来的贬值压力主要集中在刚需低价住宅。分位数差异性检验验证了异质性效应的统计显著性,分位数回归能够完整刻画协变量对房价整个条件分布的影响规律,提供区间预测结果,相比均值回归具备更强的分析价值。

相关推荐
统计学小王子15 小时前
数学建模国赛倒计时4天——《统计模型精讲(混合效应模型R语言实现)》
开发语言·数学建模·r语言
嵌入式冰箱19 小时前
数学建模辅助平台使用说明书 && 实际使用参考:各环节 Token 消耗明细及费用预估指南
数学建模
热心网友俣先生2 天前
2026国赛C题:五年命题规律与预测
java·c语言·前端·数学建模
统计学小王子2 天前
数学建模国赛倒计时5天——《软件工具(Excel)——国赛数据处理的第一道防线》
数学建模·excel
统计学小王子2 天前
数学建模国赛倒计时5天——《软件工具(mathtype)——国赛论文公式排版的唯一选择》
数学建模
统计学小王子2 天前
数学建模国赛倒计时5天——《软件工具(Excel)——快捷键总结》
数学建模·excel
binbin_522 天前
HarmonyOS 应用功耗优化实战:定位耗电、收口任务与验证回归
android·回归·harmonyos
生态学者2 天前
香港理工大学Nature Communications:沿海塑料际古菌组特征及生态影响
大数据·人工智能·算法·r语言·微信公众平台
202xxx3 天前
Day05:流程控制:if/for/while——让程序”做决定”和”重复干活”
数学建模