R ggplot - Boxplot with fill displaying different values than bargraph

Question

I have a boxplot generated using the following code, and after checking the dataset all the values are correct here.

myplot <- inDATA %>% filter(PARAMCD=="param1") %>% 
            ggplot(aes(x=ACTARMCD,y=AVAL,fill=ACTARMCD))+
            geom_boxplot()+
            stat_summary(fun.y=mean,na.rm=TRUE,shape=25,col='black',geom='point')

I want to generate a second boxplot where I split the x variable into different groups by applying a different variable as a fill. I use the following code, but the values present in the graph are incorrect.

myplot <- inDATA %>% filter(PARAMCD=="param1") %>% 
            group_by(ACTARMCD, RESPFL) %>%  
            ggplot(aes(x=ACTARMCD,y=AVAL))+
            geom_boxplot(aes(fill=RESPFL))

However when I generate a bargraph using this code, the numbers are correct.

myplot <- inDATA %>%  
            filter(PARAMCD=="param1") %>% 
            group_by(ACTARMCD,RESPFL) %>%
            dplyr::mutate(AVAL = mean(AVAL, na.rm=TRUE)) %>% 
            ggplot(aes(x=ACTARMCD,y=AVAL,fill=RESPFL))+
            geom_bar(stat="identity",position="dodge")

Can anyone please help me understand what I am doing incorrectly with the second boxplot?

Answer 1

I ended up solving the issue by using plotly instead of ggplot. The code that worked is:

myplot <- inDATA %>% filter(PARAMCD=="param1") %>% 
plot_ly(x = ~ACTARMCD, y = ~AVAL, color = ~RESPFL, type = "box",boxmean=TRUE) %>% layout(boxmode = "group")