Some minor rework:
using StatsPlots
function boxplot_missing(m)
bp=boxplot()
for i in 1:size(m,2)
bp=boxplot!(filter(!ismissing,m[:,i]))
end
display(bp)
end
a=[ 0.000163187 0.000121895 0.000229244 6.98374e-6
missing missing missing missing
0.000347736 0.000385244 0.00025956 missing]
boxplot_missing(a)
Works in Julia 1.3.1, 1.4.2 and 1.5.0.