r 箱线图中的下降因子组合

Drop factor combinations in r boxplot

我有一个数据框,其中有一个定量变量和两个分类变量,每个变量都有几个因素。当我用这些制作箱线图时,我得到一个图表,其中这些因素与缺失值的组合显示为空。所以,我想删除图中的这些空白区域。

我正在使用这些公式:

boxplot(anova_propagacion$crecimiento ~ localidad*transgen, las=2)

这是我的数据:

crecimiento transgen  localidad
1           19        N    YUCATAN
2            5        N    YUCATAN
3           17        N    YUCATAN
4           10        N    YUCATAN
5           10        N    YUCATAN
6            1        N    YUCATAN
7           19        N    NAYARIT
8            4        N    NAYARIT
9            7        N    NAYARIT
10          12        N    NAYARIT
11           2        N   PACIFICO
12          13        N    NAYARIT
13          19        N    NAYARIT
14           3        N        BCS
15           2        N        BCS
16           5        N        BCS
17           2        N        BCS
18           2        N        BCS
19           3        N        BCS
20          18        N   PACIFICO
21           2        N   PACIFICO
22           3        N   PACIFICO
23          17        N   PACIFICO
24           6        N   VERACRUZ
25           3        N   VERACRUZ
26           3        N   VERACRUZ
27          10        N   VERACRUZ
28           3        N        BCS
29           8        N        BCS
30           4        N     OAXACA
31           6        N     OAXACA
32           3        N        BCS
33         NaN        N    SINALOA
34         NaN        N TAMAULIPAS
35           5        N     OAXACA
36          17        Y     OAXACA
37          18        Y     OAXACA
38           3        Y TAMAULIPAS
39           6        Y TAMAULIPAS
40          19        Y    NAYARIT
41          19        Y    SINALOA
42           4        Y   PACIFICO
43          13        Y   PACIFICO
44           3        Y   PACIFICO
45          19        Y   PACIFICO
46          19        Y   PACIFICO
47          19        Y   PACIFICO
48          17        Y   VERACRUZ
49           2        Y        BCS
50          18        Y        BCS
51          19        Y        BCS
52         NaN        Y    YUCATAN

我得到这张图:

拜托,你能帮我去掉空格吗?

先谢谢了。

下图,省略空格:

boxplot(crecimiento ~ interaction(localidad, transgen, drop = TRUE),
        data = na.omit(anova_propagacion), las = 2)

与您的代码相比的两个主要变化是:

  • 为了计算公式中的交互项,我使用 interaction(localidad, transgen, drop = TRUE)。设置 drop = TRUE 意味着忽略实际没有出现在数据中的组合。

  • 仅此是不够的,因为某些组合do出现在数据中,但是crecimento的值是NaN (例如,SINALOA.N)。因此,它们没有被丢弃,但仍然没有什么可以绘制的。这可以通过删除带有 na.omit().

    的那些行来解决