CODEWITHSUNDEEP
rggplot2histogramdensity-plot
rf7
rf7

Reputation: 2191

ggplot - histogram - density plot by group : Only one group appears

I have a small dataset and I want to plot a histogram/density plot using ggplot by group. My dataset is the following:

> data_Test_augm
     mpg cyl4 cyl6 cyl8  disp  hp drat    wt  qsec vs0 vs1 gear3 gear4 gear5 carb1 carb2 carb3 carb4 carb6 carb8 am  predictions
 1: 21.0    0    1    0 160.0 110 3.90 2.620 16.46   1   0     0     1     0     0     0     0     1     0     0  1 1.000000e+00
 2: 21.4    0    1    0 258.0 110 3.08 3.215 19.44   0   1     1     0     0     1     0     0     0     0     0  0 7.884922e-12
 3: 17.8    0    1    0 167.6 123 3.92 3.440 18.90   0   1     0     1     0     0     0     0     1     0     0  0 7.884924e-12
 4: 32.4    1    0    0  78.7  66 4.08 2.200 19.47   0   1     0     1     0     1     0     0     0     0     0  1 1.000000e+00
 5: 30.4    1    0    0  75.7  52 4.93 1.615 18.52   0   1     0     1     0     0     1     0     0     0     0  1 7.884886e-12
 6: 19.2    0    0    1 400.0 175 3.08 3.845 17.05   1   0     1     0     0     0     1     0     0     0     0  0 7.884923e-12
 7: 26.0    1    0    0 120.3  91 4.43 2.140 16.70   1   0     0     0     1     0     1     0     0     0     0  1 1.000000e+00
 8: 30.4    1    0    0  95.1 113 3.77 1.513 16.90   0   1     0     0     1     0     1     0     0     0     0  1 7.884916e-12
 9: 19.7    0    1    0 145.0 175 3.62 2.770 15.50   1   0     0     0     1     0     0     0     0     1     0  1 1.000000e+00
10: 21.4    1    0    0 121.0 109 4.11 2.780 18.60   0   1     0     1     0     0     1     0     0     0     0  1 7.884918e-12

My code is the following:

data_Test_augm$am <- factor(data_Test_augm$am, levels = c("1" , "0"))  #Sets the faactor levels in the desired order

ggplot(data_Test_augm, aes(x = predictions, fill = am , color = am )) +
  geom_histogram(aes(y=..density..), position="identity",alpha = 0.4) + guides(color = FALSE) +
  geom_density (alpha = 0.5)+ 
  labs(title = "Predicted Probabilities per am in the Test Dataset", 
       x = "Predicted Probability of being in am 1", y = "Count") + 
  scale_fill_manual(limits=c('1', '0'),              # Defines the mapping between factor levels, labels and colors
                    labels = c("Positive", "Negative"), 
                    values = c("red", 'blue')) + 
  labs(fill="am")+  # Sets the title of the legend
  guides(color=FALSE)  # Hides the legend for Color

And this is the output: enter image description here

I don't understand why only one of the levels of the grouping variable am appears in the plot.

Upvotes: 1

Views: 792

Answers (1)

aschmsu
aschmsu

Reputation: 312

I think it's more a comment, but I can't write it because of my low reputation.

Are you sure you want a density? Without it you get a nice histogram:

dt <- data.frame(am = c(1,0,0,1,1,0,1,1,1,1), 
predictions = c(1,7.884922e-12,7.884924e-12,1,7.884886e-12,7.884923e-12,1,7.884916e-12,1,7.884918e-12))
dt$am <- factor(dt$am, levels = c("1" , "0"))

g <- ggplot(dt, aes(x = predictions, fill = am))
g+geom_histogram(aes(y=..density..), position="identity",alpha = 0.4) + 
labs(title = "Predicted Probabilities per am in the Test Dataset",
x = "Predicted Probability of being in am 1", y = "Count") +
scale_fill_manual(limits=c('1', '0'),              
labels = c("Positive", "Negative"),
values = c("red", 'blue')) +
labs(fill="am")

I just deleted geom_density (and color mapping, you already have fill!). And that's what I got: Histogram

The problem with density may occur because the data is wierd. For am = 0 you have a 100% probability to get in an interval between 7.88492e-12 and 7.88493e-12. Though it doesn't even use your colors...

Upvotes: 2

Related Questions