我有一个非常奇怪的缺失值消息,该消息不可重现,我也找不到问题。有什么建议吗? 一旦收到警告,一次也没有,并且通过计算NA,我也找不到它。
> ggplot(jnk, aes(x=Experiment, y=Log2Intensity)) +
+ geom_boxplot(outlier.shape = NA) +
+ geom_point(aes(color=imputed),
+ position=position_jitter(width=.4),
+ alpha=.5) + #scatter dots, fill if imputed
+ scale_color_manual('Imputed', values = c(`TRUE`='blue', `FALSE`='orange')) +
+ scale_y_continuous(limits=range_intensities) +
+ facet_wrap(~ my_label, ncol = 10, nrow=5)
> ggplot(jnk, aes(x=Experiment, y=Log2Intensity)) +
+ geom_boxplot(outlier.shape = NA) +
+ geom_point(aes(color=imputed),
+ position=position_jitter(width=.4),
+ alpha=.5) + #scatter dots, fill if imputed
+ scale_color_manual('Imputed', values = c(`TRUE`='blue', `FALSE`='orange')) +
+ scale_y_continuous(limits=range_intensities) +
+ facet_wrap(~ my_label, ncol = 10, nrow=5)
Warning message:
Removed 1 rows containing missing values (geom_point).
> apply(jnk, 2, function(x) sum(is.na(x)))
my_label Protein.IDs Experiment replicate imputed Log2Intensity page
0 0 0 0 0 0 0
编辑:
谢谢你的投票。如我所说,它是不可复制的。
多亏了Markus的意见,我才发现我的range_intensities
恰好达到了jnk$Log2intensity
的最大值,这以某种方式产生了此错误,但并非总是如此。
我不知道这里的算术问题是什么。
因此,为了演示,这不是胡说八道:
> range_intensities[2] >= max(jnk$Log2Intensity)
[1] TRUE
> range_intensities[1] <= max(jnk$Log2Intensity)
[1] TRUE
> range_intensities[2] == max(jnk$Log2Intensity)
[1] TRUE
我不会收到任何警告。
答案 0 :(得分:3)
问题是抖动和限制的结合:
phpinfo()