異常值處理

是 否 +2 論壇幣

k人 參與回答

經管之家送您一份

應屆畢業(yè)生專屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經管之家聯(lián)合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行數(shù)據(jù)。

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

好的，謝謝