r - multinom() はデフォルトで NA 値をどのように扱いますか?

Question

を実行しているときに、特定の行が1 つ (multinom()つまり、欠落している) であるが、すべてに値がある場合、この行全体が (SAS のように) 破棄されますか? で欠損値はどのように扱われますか?Y ~ X1 + X2 + X3X1NAYX2X3multinom()

Answer 1

以下は、?multinom異なる.nnetna.action

> library(nnet)
> library(MASS)
> example(birthwt)
> (bwt.mu <- multinom(low ~ ., bwt))

意図的にNA値を作成します。

> bwt[1,"age"]<-NA # Intentionally create NA value
> nrow(bwt)
[1] 189

4 つの異なるテストna.action:

> predict(multinom(low ~ ., bwt, na.action=na.exclude)) # Note length is 189
# weights:  12 (11 variable)
initial  value 130.311670
iter  10 value 97.622035
final  value 97.359978
converged
  [1] <NA> 0    0    0    0    0    0    0    0    0    0    0    1    1    0
 [16] 0    0    0    0    0    0    0    0    0    0    0    0    0    0    0
 ....

> predict(multinom(low ~ ., bwt, na.action=na.omit)) # Note length is 188
# weights:  12 (11 variable)
initial  value 130.311670
iter  10 value 97.622035
final  value 97.359978
converged
  [1] 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0
 [38] 0 0 0 0 0 0 1 1 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 1 1 0 0 0 0
 .....

> predict(multinom(low ~ ., bwt, na.action=na.fail))    # Generates error
Error in na.fail.default(list(low = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,  :
  missing values in object

> predict(multinom(low ~ ., bwt, na.action=na.pass))    # Generates error
Error in qr.default(X) : NA/NaN/Inf in foreign function call (arg 1)

そのため、予測でna.excludea が生成されますが、完全に省略されます。モデルを作成しません。が指定されていない場合、これはデフォルトを示します。NAna.omitna.passna.failna.action

> getOption("na.action")
[1] "na.omit"