ベクトルからNA値を削除する

いくつかのNA値を持つ巨大なベクトルがあり、そのベクトルの最大値（ベクトルはすべて数値です）を見つけようとしていますが、NA値が原因でこれを行うことができません。

NA最大値を計算できるように値を削除するにはどうすればよいですか？

を試して?maxみると、実際にはna.rm =引数があり、デフォルトでに設定されていることがわかりますFALSE。（それはを含む他の多くのRの機能のための一般的なデフォルトだsum()、mean()など）

設定na.rm=TRUEはあなたが望んでいることだけを行います：

d <- c(1, 100, NA, 10)
max(d, na.rm=TRUE)

すべてのを削除する場合は、NA代わりに次のイディオムを使用します。

d <- d[!is.na(d)]

最後の注意：他の機能（例えばtable()、lm()、およびsort()）がNA別の名前を使用します（とさまざまなオプションを提供します）関連引数を。したがってNA、関数呼び出しで問題が発生する場合は、関数の引数間で組み込みのソリューションを確認することをお勧めします。通常はすでにそこにあることがわかりました。

na.omitこの関数は、回帰ルーチンの多くは内部的に使用するものです。

vec <- 1:1000
vec[runif(200, 1, 1000)] <- NA
max(vec)
#[1] NA
max( na.omit(vec) )
#[1] 1000

?maxにna.rm設定できる追加のパラメーターがあることを示していますTRUE。

それとは別に、本当にを削除したい場合はNA、次のように使用します。

myvec[!is.na(myvec)]

あなたが呼び出すことができますmax(vector, na.rm = TRUE)。より一般的には、na.omit()関数を使用できます。

Rの初心者が元の質問に対する簡単な答えを求めている場合に備えて

ベクトルからNA値を削除するにはどうすればよいですか？

ここにあります：

foo次のようなベクトルがあるとします。

foo = c(1:10, NA, 20:30)

実行length(foo)すると22になります。

nona_foo = foo[!is.na(foo)]

length(nona_foo) NA値が削除されているため、21です。

is.na(foo)ブール行列を返すことを忘れないでください。したがってfoo、この値の反対のインデックスを使用すると、NA以外のすべての要素が得られます。

purrrdiscardから使用します（リストとベクトルで機能します）。

discard(v, is.na) 

利点は、パイプを使いやすいことです。あるいは、組み込みのサブセット機能を使用します[。

v %>% discard(is.na)
v %>% `[`(!is.na(.))

ご了承ください na.omitリストは機能しないて。

> x <- list(a=1, b=2, c=NA)
> na.omit(x)
$a
[1] 1

$b
[1] 2

$c
[1] NA

私は2つのbaseアプローチを比較するクイックベンチマークを実行しましたが、それx[!is.na(x)]よりも速いことがわかりましたna.omit。ユーザーqwrは私purrr::dicardも試してみることを提案しました-これは非常に遅いことがわかりました（ただし、実装とテストについて喜んでコメントを取ります！）

microbenchmark::microbenchmark(
  purrr::map(airquality,function(x) {x[!is.na(x)]}), 
  purrr::map(airquality,na.omit),
  purrr::map(airquality, ~purrr::discard(.x, .p = is.na)),
  times = 1e6)

Unit: microseconds
                                                     expr    min     lq      mean median      uq       max neval cld
 purrr::map(airquality, function(x) {     x[!is.na(x)] })   66.8   75.9  130.5643   86.2  131.80  541125.5 1e+06 a  
                          purrr::map(airquality, na.omit)   95.7  107.4  185.5108  129.3  190.50  534795.5 1e+06  b 
  purrr::map(airquality, ~purrr::discard(.x, .p = is.na)) 3391.7 3648.6 5615.8965 4079.7 6486.45 1121975.4 1e+06   c

参考までに、x[!is.na(x)]vs の元のテストはna.omit次のとおりです。

microbenchmark::microbenchmark(
    purrr::map(airquality,function(x) {x[!is.na(x)]}), 
    purrr::map(airquality,na.omit), 
    times = 1000000)


Unit: microseconds
                                              expr  min   lq      mean median    uq      max neval cld
 map(airquality, function(x) {     x[!is.na(x)] }) 53.0 56.6  86.48231   58.1  64.8 414195.2 1e+06  a 
                          map(airquality, na.omit) 85.3 90.4 134.49964   92.5 104.9 348352.8 1e+06   b