异常值处理

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> head(w)

  平均值 最大值 最小值

1  0.181  0.185  0.180

2  0.180  0.180  0.180

3  0.180  0.190  0.175

4  0.175  0.180  0.170

5  0.173  0.175  0.170

6  0.170  0.175  0.170    #w一共有7776行数据。

sp=boxplot(w[,2],boxwex=0.7)

> sp

$stats

      [,1]

[1,] 0.010

[2,] 0.165

[3,] 0.220

[4,] 0.270

[5,] 0.425

$n

[1] 7772

$conf

          [,1]

[1,] 0.2181182

[2,] 0.2218818

$out

[1] 1.865 0.005 0.000 0.000 0.005 0.000 0.000 0.000 0.005 0.000 0.000 0.000 0.000 0.000 0.000 0.000

[17] 0.530 0.430 0.480 0.535 0.445 1.635 0.500 0.440 0.455 0.435 0.435 0.470 0.450 0.445 0.440 0.475

[33] 0.430 0.440 0.465 0.480 0.480 0.440 0.430 2.105 0.005 0.450 0.470 0.475 0.440 0.545 0.465 0.445

[49] 0.460 0.005 0.460 0.465

$group

[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

[49] 1 1 1 1

$names

[1] ""

其中sp$out中异常值找出了，请问如何知道异常值所在的行？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：异常值 BOXPLOT Names Group stats 平均值 最大值

1. c<-list()
   
2. for(i in 1:length(sp$group)){
   
3.   d<-which(w[,2]==sp$out[i])
   
4.   c[[i]]<-d
   
5. }
   
6. unique(unlist(c))

复制代码

[1]  8478   799   836  8579  8610  8623 14458 15336 16403 22230 22235 23076 23111  8756  8778 13524 14105 21878  8779
[20] 14102 14103  8780  9623 13502 15260 11348 11650 11656 13522 14100 14104 14524 13463 13489 13490 21879 13491 14522
[39] 13501 14521 13523 14523 14101 14527 15354 14241 14526 15261 15353
> 总共都没有一万行呀，估计不对吧？

或
which(w[,2] %in% sp$out)

参考：http://stackoverflow.com/questions/11134348/how-can-i-identify-the-labels-of-outliers-in-a-r-boxplot

好的，谢谢