来自geom_point

时间:2018-08-20 07:59:03

标签: r ggplot2

我有一个非常奇怪的缺失值消息,该消息不可重现,我也找不到问题。有什么建议吗? 一旦收到警告,一次也没有,并且通过计算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

我不会收到任何警告。

1 个答案:

答案 0 :(得分:3)

问题是抖动和限制的结合:

phpinfo()