京公网安备 11010802034615号
经营许可证编号:京B2-20210330
Stata软件对截断和删失数据处理方法介绍
截断和删失是完全不同的现象,都会导致我们的样本不完整。这些现象出现在医疗科学、工程、社会科学和其他研究领域。如果忽略截断和删失,当我们分析数据时,我们的人口参数估计就会不一致。
截断和删失会出现在处理样本的过程中,那我们就从定义左截断和左删失开始:
当低于阈值的个体在样本中不存在时,我们的数据就属于左截断。比如,我们想研究某些鱼的大小,以捕鱼网为样本,鱼小于鱼网,所以在我们的样本中是不存在的。
我们的数据从K开始左删失,如果每个个体值在样本中存在并低于K,但实际值未知。例如,我们有一个测量仪器,不能检测到一定水平以下的值时,就会发生这种情况。
我们主要讨论左截断和左删失,但是我们讨论的概念可以应用到所有的截断和删失中去:右截断、右删失和区间。
当执行截断或删失数据的估计时,我们需要使用一些工具来说明这些不完整的数据。对于截断线性回归,我们可以使用truncreg命令;对于删失线性回归,我们可以使用intreg和tobit命令。
这篇文章,我们将要分析截断数据和删失数据的特征,并讨论用truncreg命令和tobit命令来说明不完整的数据。
截断数据
案例:皇家海军陆战队
Fogel et al.(1978)发布了皇家海军陆战队人员的身高的数据集,此数据可以扩展到2个世纪。它可以用来确定不同时期,英国男性的平均身高。Trussell and Bloom (1979)指出样本被截断,由于新兵最低身高的限制。数据被截断了(而不是删失),因为身高低于最低限制的个人都没有出现在样本中。考虑到这一事实,他们拟合了1800年到1809年期间皇家海军陆战队身高的截断分布。
由于Trussell和Bloom提到的问题,我们使用了人工数据集。我们假设人口数据服从正态分布μ=65和σ=3.5,并且都是左截断到64.
我们使用一个直方图来总结我们的数据:
可以看到截断点,没有小于64的数据。
如果我们忽略截断,会发生什么呢?
如果我们忽略截断,将不完整的数据视为完整的,样本均值与总体均值就会不一致,因为截断点以下的所有观测值都是缺失的。在我们的实例中,真实的均值95%都在置信区间预测平均值外。
我们可以将样本直方图与忽略截断后得出的正态分布进行比较,并且把这些值看成是人口均值和标准差的估计。
使用truncreg考虑截断
我们可以使用truncreg来估计潜在非截断分布的参数。考虑左截断64,可以使用选项ll(64)。
现在估计的值接近我们的实际模拟值μ=65,σ=3.5。
让我们将截断密度重叠到数据直方图中去。
截断分布适合我们的样本,我们分析人口分布均值等于65,标准偏差等于3.5.
删失数据
现在我们看一下删失数据的案例,看看他们和截断数据之间的区别。
案例:家庭表面尼古丁的含量情况
Matt et al.在2004年进行了一项研究,对烟草烟雾污染吸烟者家庭的整个表面进行了评估。非常有趣的一项测量是家具表面的尼古丁含量情况。每个家庭中的擦拭样本来自每件家具。然而,尼古丁污染低于一定限度的,测量仪检测不到。
数据被删失了,而不是被截断了。当尼古丁污染低于检测极限值时,样本中仍然包含了尼古丁的检测值,这个检测值就等于最低极限值。被这项研究中的这个问题启发,我随意创建了一个人工数据集。尼古丁污染水平的日志被假定为正常。在这里,lognlevel包含尼古丁含量。用于模拟日志尼古丁含量的参数,删失数据是μ=ln(5),σ=2.5,左删失数据为0.1。我们开始绘制直方图。
直方图左侧有一个尖峰,因为在检测极限以下的值被记录为等于极限值。计算样本的原始均值和标准偏差,将不会为潜在的未经审查的高斯分布提供适当的估计。
均值和标准偏差分别估计为1.68和2.4,而实际参数为ln(5) =1.61 和2.5。
使用Tobit账户审核
我们估计均值和标准偏差分布,并使用ll选项的tobit来考虑左删失值(如果审核极限值随观测值而变化,那么可以用intreg来代替)。
潜在的未经审核的分布估计的均值为1.62,标准差2.49. 我们把未经审核的分布叠加到直方图中:
潜在的未经审核的分布匹配直方图的一部分,左边尾部补偿审查点的尖峰。
总结
在抽样数据中,删失和截断是不同的两种现象。截断高斯抽样中潜在的人口参数可以用truncreg来估计。删失高斯抽样中潜在的人口参数要用intreg或tobit来估计。
结语
我们已经讨论了删失和截断的概念,也举例说明了这两个概念的意思。与本次讨论有关的要点如下:
本次讨论是基于高斯模型之上的,但是主要的概念可以扩展到任意的分布中。以上的例子在没有协变量的情况下拟合回归模型,因此,我们可以更好地可视化删失和截断分布的形状。然而,这些概念很容易扩展到协变量的回归框架中,并且特定观测值的期望值是协变量函数。
我们已经讨论过使用truncreg和tobit来处理删失和截断数据。但是这些命令也可以应用到非删失和非截断数据中,只要这些数据是特定分布中的人口抽样。
数据分析咨询请扫描二维码
若不方便扫码,搜微信号:CDAshujufenxi
主讲人简介 张冲,海归统计学硕士,CDA 认证数据分析师,前云南白药集团资深数据分析师,自媒体 Python 讲师,全网课程播放量破 ...
2026-04-10在数据可视化与业务分析中,同比分析是衡量业务发展趋势、识别周期波动的核心手段,其核心逻辑是将当前周期数据与上年同期数据进 ...
2026-04-10在机器学习模型的落地应用中,预测精度并非衡量模型可靠性的唯一标准,不确定性分析同样不可或缺。尤其是在医疗诊断、自动驾驶、 ...
2026-04-10数据本身是沉默的,唯有通过有效的呈现方式,才能让其背后的规律、趋势与价值被看见、被理解、被运用。统计制图(数据可视化)作 ...
2026-04-10在全球化深度发展的今天,跨文化传播已成为连接不同文明、促进多元共生的核心纽带,其研究核心围绕“信息传递、文化解读、意义建 ...
2026-04-09在数据可视化领域,折线图是展示时序数据、趋势变化的核心图表类型之一,其简洁的线条的能够清晰呈现数据的起伏规律。Python ECh ...
2026-04-09在数据驱动的时代,数据分析早已不是“凭经验、靠感觉”的零散操作,而是一套具备固定逻辑、标准化流程的系统方法——这就是数据 ...
2026-04-09长短期记忆网络(LSTM)作为循环神经网络(RNN)的重要改进模型,凭借其独特的门控机制(遗忘门、输入门、输出门),有效解决了 ...
2026-04-08在数据分析全流程中,数据质量是决定分析结论可靠性的核心前提,而异常值作为数据集中的“异类”,往往会干扰统计检验、模型训练 ...
2026-04-08在数字经济飞速发展的今天,数据已渗透到各行各业的核心场景,成为解读趋势、优化决策、创造价值的核心载体。而数据分析,作为挖 ...
2026-04-08在数据分析全流程中,数据处理是基础,图形可视化是核心呈现手段——前者负责将杂乱无章的原始数据转化为干净、规范、可分析的格 ...
2026-04-07在数据分析与统计推断中,p值是衡量假设检验结果显著性的核心指标,其本质是在原假设(通常为“无效应”“无差异”)成立的前提 ...
2026-04-07在数字经济深度渗透的今天,数据已成为企业生存发展的核心资产,企业的竞争本质已转变为数据利用能力的竞争。然而,大量来自生产 ...
2026-04-07Python凭借简洁的语法、丰富的生态库,成为算法开发、数据处理、机器学习等领域的首选语言。但受限于动态类型、解释性执行的特性 ...
2026-04-03在深度学习神经网络中,卷积操作是实现数据特征提取的核心引擎,更是让模型“看懂”数据、“解读”数据的关键所在。不同于传统机 ...
2026-04-03当数字化转型从企业的“战略口号”落地为“生存之战”,越来越多的企业意识到,转型的核心并非技术的堆砌,而是数据价值的深度挖 ...
2026-04-03在日常办公数据分析中,数据透视表凭借高效的汇总、分组功能,成为Excel、WPS等办公软件中最常用的数据分析工具之一。其中,“计 ...
2026-04-02在数字化交互的全场景中,用户的每一次操作都在生成动态的行为轨迹——电商用户的“浏览商品→点击详情→加入购物车”,内容APP ...
2026-04-02在数字化转型深度推进的今天,企业数据已成为驱动业务增长、构建核心竞争力的战略资产,而数据安全则是守护这份资产的“生命线” ...
2026-04-02在数据驱动决策的浪潮中,数据挖掘与数据分析是两个高频出现且极易被混淆的概念。有人将二者等同看待,认为“做数据分析就是做数 ...
2026-04-01