Inconsistency of matrix cov and cor in the presence of missing

[adienes]

julia> using Statistics

julia> x = [missing, 1, 2]; y = [-1, -2, -3];

julia> cor([x y])
2×2 Matrix{Union{Missing, Int64}}:
 1          missing
  missing  1

julia> cor([y x])
ERROR: MethodError:

would be fixed by JuliaLang/julia#56673