数据污染下的稳健估计_第1页
数据污染下的稳健估计_第2页
数据污染下的稳健估计_第3页
数据污染下的稳健估计_第4页
数据污染下的稳健估计_第5页
已阅读5页,还剩19页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

18/23数据污染下的稳健估计第一部分数据污染对稳健估计的影响 2第二部分污染性异常值检测方法 3第三部分污染鲁棒估计量 6第四部分稳健协方差矩阵的估计 8第五部分污染残差的自适应处理 11第六部分高维数据中的稳健估计 13第七部分稳健估计在实际应用中的案例 16第八部分稳健估计的局限性与展望 18

第一部分数据污染对稳健估计的影响数据污染对稳健估计的影响

数据污染是指存在极端值、异常值或错误值引入高方差和偏度的数据,这是统计建模和预测中的一个挑战。稳健估计旨在缓解数据污染的影响,以获得稳定可靠的结果。

稳健估计的定义

稳健估计是旨在对数据污染具有抵抗力的估计方法。它使用统计量度,如中位数和四分位数,这些量度对异常值不敏感。稳健估计器还可以基于稳健损失函数,这些函数以小于最小二乘损失函数的速度对异常值进行惩罚。

数据污染对稳健估计的影响

数据污染对稳健估计的影响取决于污染的程度和所用稳健估计器的类型。

极端值的影响:极端值会增加数据的方差,从而导致非稳健估计器的偏差。稳健估计器通过平均化或修剪异常值来减少极端值的影响。

异常值的影响:异常值会增加数据的偏度,导致非稳健估计器产生不准确的估计。稳健估计器使用中位数和四分位数等对异常值不敏感的统计量度,从而抵消了异常值的影响。

错误值的影响:错误值会引入随机噪声,导致非稳健估计器产生不可靠的估计。稳健估计器使用稳健损失函数来减轻错误值的影响,该函数以低于最小二乘函数的速度对异常值进行惩罚。

稳健估计器的类型

常用的稳健估计器包括:

*中位数:中位数是不受异常值影响的数据的中等值。

*四分位数:四分位数将数据分成四等份,从而提供分布的稳健度量。

*修剪平均值:修剪平均值通过去除一定比例的异常值来计算平均值。

*M估计器:M估计器基于对稳健损失函数的最小化,该损失函数以小于最小二乘函数的速度对异常值进行惩罚。

结论

数据污染是稳健估计面临的一个重大挑战。通过将数据污染的影响降至最低,稳健估计器提供了稳定可靠的结果,甚至在存在异常值或错误值的情况下也是如此。稳健估计器的选择取决于数据污染的性质和分析目标。第二部分污染性异常值检测方法关键词关键要点【穷举迭代法】:

1.将观测值划分为大小相近的子集,逐个计算每个子集的均值和方差。

2.迭代计算,将每个子集中的数据点与其他所有子集的数据点比较,找出显著差异的值。

3.识别差异值较大的数据点作为潜在异常值。

【基于密度的异常值检测】:

污染性异常值检测方法

在数据污染的情况下,识别和处理异常值至关重要,以确保稳健的估计。污染性异常值是指与数据生成过程显著偏离的极端值,它们可以对统计分析造成严重偏差。

为了检测污染性异常值,可以使用多种方法,包括:

1.统计方法

*协方差分析(ANOVA):ANOVA比较不同组别的均值,如果存在极端值,则会导致组内方差显著增加。

*Grubbs检验:Grubbs检验使用t分布来识别与其他数据点显著不同的极端值。

*Dixon检验:Dixon检验基于极差,用于识别数据集中最小或最大的极端值。

2.基于距离的方法

*马氏距离:马氏距离衡量数据点相对于多维中心位置的距离。远离中心位置的点可能是异常值。

*欧氏距离:欧氏距离是一种简单且常见的基于距离的测量,用于计算数据点之间的距离。

*局部异常因子(LOF):LOF计算每个数据点与其邻居的密度比。低密度点可能是异常值。

3.非参数方法

*箱形图:箱形图显示数据分布的四分位数范围。高于或低于四分位数箱体的极端值可能是异常值。

*直方图:直方图显示数据值的频率分布。异常值会出现为不寻常的尖峰或尾部。

*密度估计:密度估计使用平滑函数来估计数据分布。非模式值可能表明存在异常值。

4.基于机器学习的方法

*聚类分析:聚类分析将数据点分组为相似组。未被分配到任何组或被分配到小稀疏组的数据点可能是异常值。

*异常值检测算法:基于机器学习的算法,例如孤立森林和支持向量机,可以识别与训练数据集显著不同的数据点。

5.混合方法

*混合异常值检测(HAD):HAD结合多种异常值检测方法来提高准确性和鲁棒性。

*时间序列异常值检测:专用于检测时间序列数据中异常值的算法,例如滑动窗口和局部异常检测。

选择合适的方法

选择最合适的污染性异常值检测方法取决于数据类型、异常值模式和检测目标。对于小型数据集,统计方法可能是合适的。对于大型或高维数据集,基于距离或机器学习的方法更有效。

处理异常值

一旦检测到异常值,可以采取以下步骤:

*删除:删除极端异常值,但前提是它们不会包含有价值的信息。

*替换:使用替代值替换异常值,例如中位数或组内中位数。

*赢缩:调整异常值以使其更接近其他数据点。

*建模:建立一个包含异常值的统计模型,以捕获潜在的非线性或异方差。

结论

污染性异常值检测对于稳健的估计至关重要。通过识别和处理异常值,可以改善分析结果的准确性和可靠性。选择最合适的检测方法对于有效检测污染性异常值并提高统计分析的鲁棒性至关重要。第三部分污染鲁棒估计量关键词关键要点【稀疏数据估计】

1.利用稀疏性先验知识,对污染数据进行鲁棒估计。

2.开发基于ℓ0范数或正则化技术的稀疏估计方法,抑制异常值的影响。

3.研究稀疏模型的选择、正则化参数优化和算法复杂度等问题。

【非参数估计】

污染鲁棒估计

定义

污染鲁棒估计量是指在存在数据污染的情况下仍然能良好估计目标参数的估计量。

污染类型

*离群值:显著偏离总体分布的大观测值。

*缺失值:随机丢失数据的观测值。

*错误值:记录不正确或测量错误导致的观测值偏差。

*篡改值:故意改变观测值以影响结果。

污染鲁棒估计方法

*修剪法:移除一定比例的极端观测值,再使用标准估计方法。

*删帽法:在估计过程中反复移除高杠杆点观测值。

*加权法:赋予不同权重给观测值,降低离群值的影响。

*M估计:最小化一个污染鲁棒目标函数,通常涉及绝对偏差或Huber损失函数。

*MM估计:多步M估计,通过迭代过程改进估计值。

*L1惩罚法:在目标函数中加入L1范数惩罚项,以抑制离群值的影响。

*随机取样法:重复从原始数据中随机抽样,并对每个样本计算估计值,最终取平均值。

核心原则

污染鲁棒估计量通常遵循以下原则:

*高断裂点:对离群值或污染点具有较高的容忍度。

*有效率:在不存在污染的情况下接近最小方差无偏估计量。

*鲁棒性:对数据分布的轻微偏差不敏感。

主要优势

*提高估计精度:消除污染对估计结果的破坏性影响。

*增强稳定性:确保估计值不受极端观测值或数据污染的影响。

*适用广泛:可用于各种数据类型和污染模型。

主要挑战

*计算复杂性:某些鲁棒估计方法可能在计算上很密集。

*精度损失:与标准估计方法相比,鲁棒估计量可能存在一定程度的效率损失。

*模型选择:选择适当的污染鲁棒方法需要考虑污染类型和样本大小。

应用领域

污染鲁棒估计量在以下领域具有广泛的应用:

*污染检测:识别和识别污染观测值。

*异常值分析:探索离群值和数据分布中的异常情况。

*建模和预测:构建鲁棒的统计模型,即使存在污染也能准确预测。

*财务和经济:估计资产价格、风险和经济指标,不受操纵或错误数据的影响。

*医疗和生物统计:处理临床试验数据中的缺失值或测量错误。第四部分稳健协方差矩阵的估计稳健协方差矩阵的估计

引言

在数据分析中,协方差矩阵估计对于理解变量之间的相关性至关重要。然而,当数据受到污染(存在离群值或极端值)时,传统协方差矩阵估计方法可能会产生误导性的结果。稳健协方差矩阵估计技术旨在减少离群值的影响,提供更准确的协方差矩阵估计。

稳健协方差矩阵估计方法

1.M估计:

M估计通过最小化加权残差平方和来估计协方差矩阵,其中权重函数用于减少离群值的影响。常见的M估计函数包括:

*Huber函数

*Tukey的双权重函数

*Bisquare函数

2.Trimmed均值:

Trimmed均值通过剔除一定比例的极端值(例如,最高和最低10%)来估计协方差矩阵。然后,使用剩余数据的均值和协方差来计算稳健的协方差矩阵。

3.Winsorization:

Winsorization通过将极端值替换为数据的特定分位数(例如,第25或第75分位数)来减少离群值的影响。然后,使用Winsorized数据的样本均值和协方差来估计稳健的协方差矩阵。

4.协方差的稳健度量:

协方差的稳健度量通过使用对离群值不敏感的统计量来估计协方差矩阵。这些统计量包括:

*平均绝对偏差(MAD)

*中位绝对偏差(MADN)

*沿分位数回归(QRR)

稳健协方差矩阵估计的优缺点

优点:

*减少离群值的影响。

*提供更准确的协方差矩阵估计。

*对数据的非正态性和异方差性鲁棒。

缺点:

*计算复杂度较高。

*可能需要精确指定稳健度量。

*效率可能低于传统协方差矩阵估计方法。

选择稳健协方差矩阵估计方法

选择稳健协方差矩阵估计方法取决于数据污染的程度和数据的分布。以下是一些准则:

*如果数据只有轻微污染,则可以考虑M估计。

*如果数据污染严重,则可以考虑Trimmed均值或Winsorization。

*如果对数据的分布不确定,则可以考虑使用基于协方差稳健度量的估计方法。

应用

稳健协方差矩阵估计在各种应用中都很有用,包括:

*多元回归和时间序列分析。

*风险建模和金融分析。

*数据清洗和异常检测。

结论

稳健协方差矩阵估计是处理数据污染问题的强大工具。通过减少离群值的影响,它可以提供更准确的协方差矩阵估计,从而提高后续统计分析的准确性。在选择稳健协方差矩阵估计方法时,需要考虑数据污染的程度和数据的分布。第五部分污染残差的自适应处理污染残差的自适应处理

污染残差的存在会严重影响统计分析的结果,因此需要对污染残差进行处理。污染残差的自适应处理是一种稳健估计方法,可以有效地去除残差中的污染点,从而提高估计结果的准确性。

基本原理

污染残差的自适应处理基于以下基本原理:

*污染残差通常具有极端值或异常值,与正常分布的残差显著不同。

*污染残差的数量相对于正常分布的残差而言非常少。

因此,自适应处理方法通过识别并去除极端残差,从而达到去除污染残差的目的。

具体方法

污染残差的自适应处理方法有很多种,其中最常用的方法包括:

Tukey-Huber损失函数

Tukey-Huber损失函数是一种非平方的损失函数,对于小残差,它与平方损失函数类似,但对于大残差,它则增长得更慢。这使得Tukey-Huber损失函数对污染残差具有鲁棒性。

MM估计

MM估计(MaximumLikelihoodwithaMixtureDistribution)是一种混合模型估计方法。它假设残差服从混合分布,其中污染残差的权重很小。通过极大化混合分布的对数似然函数,可以得到稳健的估计结果。

LTS估计

LTS估计(LeastTrimmedSquares)是一种修剪平均方法。它选择残差绝对值最小的子集,并使用这个子集来估计模型参数。LTS估计对污染残差具有很强的鲁棒性,但它需要大量的观测数据。

自适应加权方法

自适应加权方法给不同的残差分配不同的权重。污染残差的权重较小,而正常分布残差的权重较大。通过迭代更新残差的权重,可以得到稳健的估计结果。

优点

污染残差的自适应处理方法具有以下优点:

*鲁棒性强:可以有效地去除残差中的污染点,提高估计结果的准确性。

*适用性广:可以适用于各种线性回归模型和非线性回归模型。

*易于实现:大多数统计软件都提供了污染残差的自适应处理功能,方便使用。

局限性

污染残差的自适应处理方法也存在一些局限性:

*计算量大:某些方法,如LTS估计,需要大量的计算。

*可能去除有价值的信息:在某些情况下,自适应处理方法可能会去除有价值的信息,导致估计结果的偏差。

*对污染类型敏感:不同的自适应处理方法对不同类型的污染残差具有不同的鲁棒性。

应用

污染残差的自适应处理方法在实际应用中非常广泛,例如:

*金融数据分析:去除异常交易数据的影响,提高财务模型的准确性。

*图像处理:去除图像中的噪声和杂质,提高图像质量。

*医疗数据分析:去除极端值病例的影响,提高医疗研究结果的可靠性。

*气象数据分析:去除异常天气事件的影响,提高气候预测的准确性。

总之,污染残差的自适应处理是一种有效且实用的稳健估计方法,可以有效地提高估计结果的准确性,但在使用时需要注意其局限性,并根据具体情况选择最合适的处理方法。第六部分高维数据中的稳健估计关键词关键要点高维数据中的稳健估计

主题名称:高维数据中的维度灾难

1.在高维数据中,观测值的样本量往往远小于变量的维度,导致数据过稀疏,难以估计参数。

2.维度灾难会导致传统参数估计方法(如最小二乘法)失效,产生不稳定和有偏差的估计结果。

3.需要采用专门针对高维数据设计的稳健估计技术来解决维度灾难问题。

主题名称:降维技术

高维数据中的稳健估计

高维数据是具有大量特征或维度的数据集,在现代数据分析中越来越普遍。高维数据给稳健估计带来了独特的挑战,因为传统方法可能对异常值和噪声数据高度敏感。稳健估计旨在抵御这些异常值的影响,从而产生可靠且可信的结果。

稳健估计的挑战

在高维数据中,异常值可能对传统估计方法产生过度影响。这是因为高维空间中数据点之间的距离更近,异常值更有可能被认为是局部极值。此外,随着维数的增加,数据分布变得更加复杂,这使得识别异常值变得困难。

稳健估计的方法

有几种方法可以实现稳健估计,包括:

*重加权方法:将较小的权重分配给异常值,从而降低其对估计的影响。

*截断方法:删除超过一定阈值的极端值。

*中位数方法:使用数据集的中位数作为估计量,因为中位数不受异常值的影响。

*M估计量:通过最大化一个稳健的目标函数来获得估计量,该函数对异常值不那么敏感。

具体方法

以下是一些具体的高维数据稳健估计方法:

*重加权最小二乘法(RWLS):通过将由异常值产生的残差赋予较小的权重,对最小二乘法进行稳健化。

*L1正则化:向损失函数中添加L1范数惩罚项,这会惩罚异常值。

*凸优化方法:利用凸优化框架对估计问题进行公式化,这允许使用快速和有效的求解器。

*几何中位数:通过求解一组几何方程来计算数据点的中位数,这些方程对异常值不敏感。

应用

高维数据中的稳健估计在许多领域具有应用,包括:

*异常值检测:识别数据集中的异常值,这些异常值可能表明数据损坏或欺诈。

*数据降维:通过仅选择与目标变量相关的特征来对高维数据进行降维。

*机器学习:开发对异常值不敏感的机器学习模型,从而提高预测性能。

局限性

尽管稳健估计在高维数据中很有用,但它也有一些局限性:

*效率:稳健估计方法通常比非稳健方法计算效率低。

*偏差:稳健估计量可能比非稳健估计量有更大的偏差,尤其是在数据中存在极端值的情况下。

结论

高维数据中的稳健估计是应对异常值和噪声挑战的至关重要的工具。通过使用适当的稳健估计方法,数据分析人员可以获得可靠且可信的结果,即使在具有挑战性的高维数据集的情况下也是如此。第七部分稳健估计在实际应用中的案例稳健估计在实际应用中的案例

1.计量经济学模型

*OLS回归:OLS回归假设误差项正态分布,但若数据分布非正态,则OLS系数估计量可能存在偏差。稳健估计方法,如M-估计,可以降低异常值对估计结果的影响,得到更准确的系数估计。

*二元选择模型:二元选择模型(如Logit、Probit)假设误差项服从特定分布(例如正态分布或逻辑分布)。然而,当数据分布偏离这些假设时,稳健估计方法可以提高估计的稳定性。

2.金融风险建模

*价值atRisk(VaR):VaR衡量金融投资组合在特定置信水平下可能遭受的最大损失。稳健估计方法可以处理极端值,从而产生更可靠的VaR估计,从而降低金融风险管理的风险。

*尾部风险建模:尾部风险是指极端事件发生的可能性。稳健估计方法可以捕获尾部分布的特征,从而提高尾部风险模型的准确性,并为风险管理提供更好的指导。

3.生物统计学

*生存分析:生存分析研究个体经历特定事件(如死亡或复发)的时间。稳健估计方法可以降低异常值对生存函数估计结果的影响,从而得到更可靠的生存率估计。

*药物疗效评估:药物疗效评估需要比较不同治疗组的疗效。稳健估计方法可以处理异常值和偏斜数据,从而得到更准确的治疗效果比较结果。

4.环境科学

*空气污染建模:空气污染建模需要预测空气中污染物的浓度。稳健估计方法可以处理异常值和极端观测值,从而提高预测模型的准确性,并更准确地评估空气污染的影响。

*水质监测:水质监测需要测量水体中的特定参数,如溶解氧和pH值。稳健估计方法可以降低异常值对水质特征估计的影响,从而得到更可靠的水质评估结果。

5.社会科学

*问卷调查:问卷调查数据通常包含缺失值和异常值。稳健估计方法可以处理这些数据问题,从而得到更准确的调查结果。

*心理学研究:心理学研究需要测量个体的认知和情感特征。稳健估计方法可以降低异常值对测量结果的影响,从而提高心理测量工具的准确性。

案例研究:金融风险建模

背景:一家金融机构需要对投资组合进行VaR估计,以评估其在特定置信水平下的最大潜在损失。

方法:

*使用传统OLS方法估计VaR模型。

*使用稳健M-估计方法估计VaR模型。

结果:

*OLS模型估计的VaR值为1000万美元。

*M估计模型估计的VaR值为1500万美元。

解释:

M估计模型估计出的VaR值更高,表明该模型可以捕获投资组合分布中的极端值。这对于风险管理至关重要,因为它可以更准确地评估投资组合面临的潜在损失风险。

结论:

稳健估计方法在实际应用中具有广泛的应用,因为它可以降低异常值和偏斜数据的影响,从而提高估计结果的准确性和稳定性。在金融风险建模、生物统计学、环境科学和社会科学等领域,稳健估计方法为可靠的决策制定提供了基础。第八部分稳健估计的局限性与展望关键词关键要点1.稳健估计的计算复杂性

1.稳健估计方法通常具有较高的计算复杂性,尤其是当数据量较大时。

2.随着数据集的增大,鲁棒协方差矩阵的计算时间呈指数级增长。

3.这对实时应用和大型数据集的分析构成了挑战,需要开发更有效的算法和优化技术。

2.有限样本性能

稳健估计的局限性

尽管稳健估计方法具有应对数据污染的优势,但它们也存在局限性:

*效率损失:稳健估计器通常会牺牲一些效率来提高稳健性,这意味着它们在无污染数据的情况下可能不如经典估计器准确。

*有限的分布适用性:大多数稳健估计方法针对特定分布族(例如正态分布)进行优化。在非正态分布的情况下,稳健估计器的性能可能会下降。

*数据类型限制:稳健估计方法通常仅适用于连续数据。对于分类或有序数据,适用性可能会受到限制。

*超参数依赖性:一些稳健估计器需要手动设置超参数,例如权值函数或内核带宽。优化这些超参数对于稳健估计的性能至关重要,但可能是一项耗时的过程。

稳健估计的展望

尽管存在局限性,稳健估计仍然是应对数据污染和提高统计模型鲁棒性的宝贵工具。未来的研究重点可能包括:

*分布无关的方法:开发对分布假设不敏感的稳健估计器。

*优化超参数选择:探索自动化超参数优化技术,以最大化稳健估计器的性能。

*新稳健统计量:开发新的稳健统计量,以更全面地捕获数据的鲁棒性特征。

*与机器学习技术的整合:探索稳健估计与机器学习技术的融合,以创建对异常值和噪声更鲁棒的机器学习模型。

*行业特定应用:探索稳健估计在特定行业中的应用,例如金融、医疗保健和工业。

此外,稳健估计的教育和推广对于提高人们对数据污染和稳健统计方法重要性的认识至关重要。通过加强教育和提供易于使用的工具,我们可以提高人们有效处理和分析污染数据的技能,从而做出更可靠的决策。关键词关键要点【数据污染对稳健估计的影响】

关键词关键要点主题名称:稳健协方差矩阵估计方法

关键要点:

1.最小覆盖椭球估计(MVE):使用数据中所有点估计协方差矩阵,但对极端值的影响较小。

2.最小协方差行列式估计(MCDE):最小化协方差矩阵的行列式,对数据中异常值具有鲁棒性。

3.最小中位差估计(MMD):使用数据点之间的中位差计算协方差矩阵,对异常值和厚尾分布具有鲁棒性。

主题名称:稳健协方差矩阵估计的应用

关键要点:

1.回归分析:在存在异常值或数据影响点时,计算稳健的回归系数。

2.主成分分析(PCA):对受异常值影响的主成分进行鲁棒估计。

3.多元检验:在多元假设检验中,估计具有鲁棒性的协方差矩阵,以避免异常值的影响。

主题名称:稳健协方差矩阵估计的局限性

关键要点:

1.效率:稳健协方差矩阵估计器通常比非稳健估计器效率较低。

2.计算成本:一些稳健协方差矩阵估计器,例如MVE,可能计算成本很高,尤其是在高维数据集中。

3.解读性:稳健协方差矩阵估计器可能难以解释,因为它们可能使用非标准的估计程序。

主题名称:稳健协方差矩阵估计的未来趋势

关键要点:

1.分布自由方法:开发不需要假设数据分布就能提供稳健估计的方法。

2.稀疏估计:针对高维稀疏数据开发稳健协方差矩阵估计器。

3.机器学习技术:探索结合机器学习技术来增强稳健协方差矩阵估计的可能性。关键词关键要点主题名称:稳健估计中的污染残差自适应处理

关键要点:

1.污染残差的识别:利用统计量或模型选择标准(如残差标准误差,AIC或BIC)识别是否存在污染残差。

2.污染残差的分类:将污染残差分为可观测的异常值和不可观测的异常值,前者可以通过数据清理或变换修复,而后者需要使用稳健估计方法。

主题名称:稳健估计的分类

关键要点:

1.M估计:使用最大似然或最小化加权平方差函数,其中权值函数对污染残差具有较低敏感性。

2.MM估计:在M估计的基础上,进一步引入一个权值函数以减少权值分配对结果的影响。

3.其他稳健估计方法:包括最小绝对偏差(LAD)估计、加权最小绝对偏差(WLS)估计和最小相对偏差估计。

主题名称:稳健估计的性能评估

关键要点:

1.效率:将稳健估计方法与经典估计方法的方差进行比较,以评估其效率损失。

2.稳健性:通过模拟受污染的数据集,评估稳健估计方法对污染残差的鲁棒性。

3.鲁棒性:评估稳健估计方法在不同分布或模型假设下的稳定性。

主题名称:稳健估计的应用

关键要点:

1.回归模型:在存在污染数据的回归模型中,稳健估计方法可以提供更准确和稳定的估计。

2.时序模型:在

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论