R语言复数类型NA的行为不一致性问询及优化提案

我来深入拆解你发现的R语言复数类型NA处理的直觉冲突问题，先回顾你的测试代码：

z <- complex(real = c( 1, 1, NA, NA, NA), imaginary = c(NA, NA, 1, 1, NA))
print(z) #> [1] NA NA NA NA NA
Re(z) #> [1] 1 1 NA NA NA
Im(z) #> [1] NA NA 1 1 NA
unique(z) #> [1] NA
Re(unique(z)) #> [1] 1
identical(z[1], z[3]) #> FALSE

现有行为的合理性判断 #

我们逐个看这些行为是否属于设计预期，哪些存在问题：

• print(z)输出全NA：这是符合R设计逻辑的——只要复数的实部或虚部包含NA，整个复数就被认定为NA值，所以打印时统一显示NA是预期行为。
• Re(z)和Im(z)返回原始实/虚部值：这是当前R的实现漏洞，属于非预期的Bug。从用户视角，一个被标记为NA的复数，它的实部和虚部都应该被视为NA，而不是返回存储的原始数值。
• unique(z)返回单个NA：这和print(z)的逻辑一致，R从概念上认为所有复数NA都是等价的，所以去重后只剩一个，这是预期行为。
• Re(unique(z))返回1：这是Re()函数现有Bug的延续——unique(z)返回的NA复数实际存储的是第一个元素的实部（1）和虚部（NA），但从语义上这个复数是NA，它的实部也应该是NA，所以这属于Bug。
• identical(z[1], z[3])返回FALSE：identical()默认会严格检查对象的二进制表示，z[1]是1+NAi，z[3]是NA+1i，二进制存储不同，所以返回FALSE。但从复数NA的语义角度，这两个都应该被视为同一个NA值，所以这个行为不符合直觉，需要优化。

一致性修复提案 #

为了维持R对复数NA“等价表示”的概念（即只要实部或虚部含NA，就视为同一个复数NA），同时实现行为一致性，我提出以下调整方案：

核心规则：如果两个复数是NA的不同表示形式，那么所有函数对它们的处理结果必须等价。

调整后的预期行为如下：

# Re()/Im()：只要复数是NA，就返回NA
Re(z) #> [1] NA NA NA NA NA
Im(z) #> [1] NA NA NA NA NA

# 提取NA复数的实部：语义上NA的实部也是NA
Re(unique(z)) #> [1] NA

# 默认将所有复数NA视为等价
identical(z[1], z[3]) #> TRUE

# 保留底层二进制检查的入口：通过single.NA参数切换
identical(z[1], z[3], single.NA = FALSE) #> FALSE

这个方案既保留了R对复数NA的核心定义，又修复了不符合直觉的行为，同时通过single.NA参数给需要底层操作的用户保留了灵活度，兼顾了语义一致性和实用性。

内容的提问来源于stack exchange，提问作者Patrick Perry