格拉布斯离群值检验(版本 2)(此脚本是新版本,提供By group处理。 最后! 请注意,第一个版本中报告的 p 值不再可用。)
此脚本将 Grubbs 的双尾异常值测试添加到分发平台。 打开正态分位数图和拟合优度检验以帮助评估假设样本是从正常人群中抽取的。
只需打开包含要评估的数值变量的数据表,然后打开并运行脚本即可。 选择数据列并单击是,数据。 指定所需的显着性水平(默认情况下 alpha 为 0.05)。 这个例子是基于变量高度在里面大班数据表中的样本数据文件夹。

点击好的。

正常分位数图中标记的模式似乎是线性的,并且没有标记位于点状读数曲线指定的区域之外。 在这种情况下,样本以及总体被判断为正常。 平台底部的测试在指定水平上并不显着。
现在使用经过多变的。 本示例使用性别作为分组变量。

点击好的。

探索性数据分析统计数据回复:格拉布斯离群值检验(版本 2)谢谢你的剧本。我赞同 Jans 的建议,即采用某种方法来指示检测到的异常值。
回复:格拉布斯离群值检验(版本 2)你好,谢谢你的剧本。
在输出中包含实际上被视为异常值的值会很有帮助
(或者至少提供一个指示是否涉及最低值或最高值)。
回复:格拉布斯离群值检验(版本 2)如果 by 变量是数字,则脚本会中断。解决此问题的一个简单方法是将行后的 bcol 更改为字符数据类型
dt = 当前数据表();
If( N 个项目( bCol ),
bCol[1] << 设置数据类型(字符);
/* 您还可以修复脚本以使用数字变量而不需要对它们进行修改,但这似乎更简单,您总是可以在最后将其更改回数字*/
干杯
冈特
回复:格拉布斯离群值检验(版本 2)<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />如果数据集恰好有任何缺失值,则此代码会错误地计算 Grubbs 检验的 N。如果缺失列的数量相当大,这可能会出现问题。
请参阅以下更正:
第 76 行和 110 行
n = N 行( yVal ) - N 缺失( yVal );
回复:格拉布斯离群值检验(版本 2)嗨马克
感谢您发布此内容。
如果可以通过分组变量来完成,我会很棒。
BR,玛丽安