95 パーセンタイルを超えて 5 パーセンタイルを下回る値をとる比較的大きなRデータセットのすべての値を、それぞれそれらのパーセンタイル値に置き換えたいと思います。私の目的は、これらの外れ値をデータから完全に切り取らないようにすることです。
アドバイスをいただければ幸いです。他の場所でこれを行う方法に関する情報は見つかりません。
これで十分です。
fun <- function(x){
quantiles <- quantile( x, c(.05, .95 ) )
x[ x < quantiles[1] ] <- quantiles[1]
x[ x > quantiles[2] ] <- quantiles[2]
x
}
fun( yourdata )
次を使用して、1行のコードで実行できますsquish()
。
d2 <- squish(d, quantile(d, c(.05, .95)))
はかりライブラリで、?squish
と?discard
#--------------------------------
library(scales)
pr <- .95
q <- quantile(d, c(1-pr, pr))
d2 <- squish(d, q)
#---------------------------------
# Note: depending on your needs, you may want to round off the quantile, ie:
q <- round(quantile(d, c(1-pr, pr)))
例:
d <- 1:20
d
# [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
d2 <- squish(d, round(quantile(d, c(.05, .95))))
d2
# [1] 2 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 19
私はあなたが必要とするものを得るためにこのコードを使用しました:
qn = quantile(df$value, c(0.05, 0.95), na.rm = TRUE)
df = within(df, { value = ifelse(value < qn[1], qn[1], value)
value = ifelse(value > qn[2], qn[2], value)})
df
data.frameとvalue
、データを含む列はどこにありますか。