中值.default（）条件返回值的代码说明

Question

我遇到了中值median.default中的median.default我根本无法理解的代码。 当na.rm = FALSE ，它涉及返回值，并且向量确实包含至少一个NA值。

有问题的代码是

x[FALSE][NA]

在中median.default body元素4和6中median.default

> bod <- body(median.default)
> bod[[4]]
# if (na.rm) 
#     x <- x[!is.na(x)] 
# else if (any(is.na(x))) 
#     return(x[FALSE][NA])
> bod[[6]]
# if (n == 0L) 
#     return(x[FALSE][NA])

当我分解这段代码时，我就迷路了。 我只是不明白它的目的。

当您只需编写return(NA)时，为什么会这样写？ 还是使用as.numeric(NA)或其他as.*函数？

> x <- c(1, NA, 5, NA, 12)
> if(any(is.na(x))) x[FALSE][NA]
# [1] NA
> x[FALSE]
# numeric(0)
> if(any(is.na(x))) NA
# [1] NA

如果答案是将NA值转换为x的同一类，那么为什么不使用as.* ？

> Y <- c(as.integer(NA), as.double(NA), as.numeric(NA)) 
> sapply(Y, function(x) c(class(x), typeof(x), mode(x)))
#      [,1]      [,2]      [,3]     
# [1,] "numeric" "numeric" "numeric"
# [2,] "double"  "double"  "double" 
# [3,] "numeric" "numeric" "numeric" 

Answer 1

我认为这就是他们返回相同类型x的NA的方式。

来源： 中-R上的平均页面。

例：

> x <- 1:3
> y <- NA
> z <- x[FALSE][NA]
> class(x)
[1] "integer"
> class(y)
[1] "logical"
> class(z)
[1] "integer"
> x
[1] 1 2 3
> y
[1] NA
> z
[1] NA
> 

进一步的示例，说明为什么两者都需要：

> x
[1] 1 2 3
> x[FALSE]
integer(0)
> x[NA]
[1] NA NA NA
> x[FALSE][NA]
[1] NA
> 

加[FALSE]可将x减小为零长度矢量（但仍为数字）。 加[NA]使其实际上返回NA。