WinsorizedVariance[list,f]
给出把最小和最大的 f(比例)元素用剩下元素的极值替换后 list 中元素的方差.
WinsorizedVariance[list,{f1,f2}]
给出把最小的 f1(比例)元素和的最大的 f2(比例)元素用剩下元素的极值替换后列表中元素的方差.
WinsorizedVariance[list]
给出经 5% 缩尾处理后的方差 WinsorizedVariance[list,0.05].
WinsorizedVariance[dist,…]
给出单变量分布 dist 经缩尾处理后的方差.
WinsorizedVariance
WinsorizedVariance[list,f]
给出把最小和最大的 f(比例)元素用剩下元素的极值替换后 list 中元素的方差.
WinsorizedVariance[list,{f1,f2}]
给出把最小的 f1(比例)元素和的最大的 f2(比例)元素用剩下元素的极值替换后列表中元素的方差.
WinsorizedVariance[list]
给出经 5% 缩尾处理后的方差 WinsorizedVariance[list,0.05].
WinsorizedVariance[dist,…]
给出单变量分布 dist 经缩尾处理后的方差.
更多信息
- 由于偏离程度大的极值被偏离程度较小的极值所替代,WinsorizedVariance 给出对方差的稳健估计.
- 缩尾比例由参数 f1 和 f2 确定,表示把最小的 f1(比例)元素和的最大的 f2(比例)元素替换成剩下元素的极值.
- WinsorizedVariance[list,{f1,f2}] 给出 Clip[list,{z1,z2}] 的方差,其中 z1 等于 RankedMin[list,1+
],z2 等于 RankedMax[list,1+
],n 等于 list 的长度. » - 单变量 WeightedData data 的 WinsorizedVariance 给出截尾 data 的加权方差. »
- WinsorizedVariance[{{x1,y1,…},{x2,y2,…},…},f] 给出 {WinsorizedVariance[{x1,x2,…},f],WinsorizedVariance[{y1,y2,…},f],…}. »
- 对于单变量分布 dist,WinsorizedVariance[dist,{f1,f2}] 给出 Variance[CensoredDistribution[Quantile[dist,{f1,1-f2}],dist]]. »
范例
打开所有单元 关闭所有单元基本范例 (4)
WinsorizedVariance[{-10, 1, 1, 1, 1, 20}, 0.2]WinsorizedVariance[{-10, 1, 1, 1, 1, 20}, {0.2, 0}]RandomDate[6, DateGranularity -> "Month"]//SortWinsorizedVariance[%]WinsorizedVariance[ExponentialDistribution[λ]]范围 (11)
数据 (10)
WinsorizedVariance[{1, 20, 3, 4}, 2 / 5]WinsorizedVariance[{Sqrt[2], E, Pi, Pi ^ 2, 1, 2, 3}, 1 / 4]WinsorizedVariance[{5., 10., 4., 25., 2., 1.}, 0.2]WinsorizedVariance[N[{5, 10, 4, 25, 2, 1}, 30], 0.2]WinsorizedVariance[RandomReal[1, {50, 2}], .1]WinsorizedVariance[RandomReal[1, 10 ^ 6], 1 / 100]WinsorizedVariance[RandomReal[1, {10 ^ 5, 5}]]可以像对稠密数组一样使用 SparseArray 数据:
sp = SparseArray[{{i_, i_} :> i, {i_, j_} /; j == i + 1 :> i - 1}, {100, 10}, 1];WinsorizedVariance[sp]一元 WeightedData 的 WinsorizedVariance:
wd = WeightedData[RandomReal[100, 30], Range[30]];WinsorizedVariance[wd]WinsorizedVariance[wd["InputData"]]TimeSeries 的缩尾方差:
ts = TemporalData[TimeSeries, {{{3, 8, 4, 11, 9, 2}}, {{{1, 3, 5, 7, 8, 10}}}, 1, {"Continuous", 1},
{"Discrete", 1}, 1, {ResamplingMethod -> {"Interpolation", InterpolationOrder -> 1}}}, False,
10.1];WinsorizedVariance[ts, 1 / 6]//NWinsorizedVariance[ts["Values"], 1 / 6]//Ndata = Quantity[RandomReal[1, 6], "Meters"]WinsorizedVariance[data]dates = WolframLanguageData[All, "DateIntroduced"];DateHistogram[dates]WinsorizedVariance[dates]UnitConvert[%, "Years" ^ 2]RandomTime[3]WinsorizedVariance[%]{TimeObject[{12}, TimeZone -> 0], TimeObject[{12}, TimeZone -> 2], TimeObject[{12}, TimeZone -> "Asia/Tokyo"]}WinsorizedVariance[%]应用 (2)
WinsorizedVariance[{1, 5, 2, 6, 10, 10 ^ 5, 5, 4, -200, 5}, .1]//NVariance[{1, 5, 2, 6, 10, 10 ^ 5, 5, 4, -200, 5}]//Nheights = Quantity[{134, 143, 131, 140, 145, 136, 131, 136, 143, 136, 133, 145, 147,
150, 150, 146, 137, 143, 132, 142, 145, 136, 144, 135, 141}, "Centimeters"];ListPlot[heights, Filling -> Axis, AxesLabel -> Automatic]wv = WinsorizedVariance[heights]//NwinVar[f_ ? NumericQ] /; 0 ≤ f < 0.5 := WinsorizedVariance[heights, f]Plot[winVar[f], {f, 0, 0.49}, AxesLabel -> {f}]m = WinsorizedMean[heights];
n = Length[heights];
swv = Sqrt[wv];ListPlot[{heights, {{0, m}, {n, m}}, {{0, m - swv}, {n, m - swv}}, {{0, m + swv}, {n, m + swv}}}, Filling -> {1 -> 0, 3 -> {4}}, Joined -> {False, True, True, True}, PlotStyle -> {Automatic, Automatic, Gray, Gray}, PlotLegends -> {"heights", "winsorized mean", "square-root-winsorized-variance bands"}, AxesLabel -> Automatic]属性和关系 (5)
0% WinsorizedVariance 等价于 Variance:
WinsorizedVariance[Range[10], 0]Variance[Range[10]]当 f 趋于 1/2 时 WinsorizedVariance 趋于 0:
data = RandomReal[CauchyDistribution[0, 1], 1000];Plot[WinsorizedVariance[data, f], {f, 0, .499}]一个分布的 WinsorizedVariance 是它的 CensoredDistribution 的方差:
𝒟 = NormalDistribution[a, b];f1 = .2;
f2 = .3;wv = WinsorizedVariance[𝒟, {f1, f2}]适当界限下的 CensoredDistribution 的方差:
cd = CensoredDistribution[Quantile[𝒟, {f1, 1 - f2}], 𝒟]var = Variance[cd]wv == var样本的 WinsorizedVariance 给出对删失分布的方差的估计:
𝒟 = NormalDistribution[];
data = RandomVariate[𝒟, 10 ^ 6];f1 = .2;
f2 = .3;WinsorizedVariance[data, {f1, f2}]适当界限下的 CensoredDistribution 的方差:
cd = CensoredDistribution[Quantile[𝒟, {f1, 1 - f2}], 𝒟]Variance[cd]TrimmedVariance 会丢弃一定分位数之上的数据,然后再计算样本方差:
len = 100;
f = 1 / 10;
data = RandomReal[1, len];data1 = Part[Sort[data], 1 + Floor[len f] ;; len - Floor[len f]];TrimmedVariance[data, f] == Variance[data1]WinsorizedVariance 剪切一定分位数之上的数据,然后再计算样本方差:
data2 = Clip[data, {RankedMin[data, 1 + Floor[len f]], RankedMax[data, 1 + Floor[len f]]}];WinsorizedVariance[data, f] == Variance[data2]绘制排序后的数据与经过截断(移除元素)或剪切(较大极值被替换)处理的数据:
d0 = Transpose[{Range[1, len], Sort@data}];
d1 = Transpose[{Range[1 + Floor[len f], len - Floor[len f]], data1}];
d2 = Transpose[{Range[1, len], Sort@data2}];ListPlot[{d0, d1, d2}, Filling -> {2 -> 0}, PlotLegends -> {"data", "data1", "data2"}]可能存在的问题 (1)
文本
Wolfram Research (2017),WinsorizedVariance,Wolfram 语言函数,https://reference.wolfram.com/language/ref/WinsorizedVariance.html (更新于 2024 年).
CMS
Wolfram 语言. 2017. "WinsorizedVariance." Wolfram 语言与系统参考资料中心. Wolfram Research. 最新版本 2024. https://reference.wolfram.com/language/ref/WinsorizedVariance.html.
APA
Wolfram 语言. (2017). WinsorizedVariance. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/WinsorizedVariance.html 年
BibTeX
@misc{reference.wolfram_2026_winsorizedvariance, author="Wolfram Research", title="{WinsorizedVariance}", year="2024", howpublished="\url{https://reference.wolfram.com/language/ref/WinsorizedVariance.html}", note=[Accessed: 09-August-2026]}
BibLaTeX
@online{reference.wolfram_2026_winsorizedvariance, organization={Wolfram Research}, title={WinsorizedVariance}, year={2024}, url={https://reference.wolfram.com/language/ref/WinsorizedVariance.html}, note=[Accessed: 09-August-2026]}