2x2 bar charts require several steps. For ease of comprehension I will list them, and then list the code below.

​

(Step 1) Make sure your IVs are set to be "factors". Load "psych" library, and then use the "describeBy" function to extract means for each condition

(Step 2) "describeBy" combines the two IVs names. Split them and assign the names you would like to use in your figure

(Step 3) Create the Figure

​

Also note: The names you give to your levels of your IV will appear on the Figure. You should make them appropriate in your data set before starting.

​

​

___Step 1_____________________

ds$IV1 <- as.factor(ds$IV1)

ds$IV2 <- as.factor(ds$IV2)

​

library(psych)

​

XYXY<-describeBy(ds$DV,

group = ds$IV1 : ds$IV2,
mat = TRUE,
type = 2,
digits = 2)

_____________________________

​

​

___Step 2_____________________

XYXY$IV1_Name <-sub("[:].*", "", XYXY$group1)   
XYXY$IV2_Name <-sub(".*[:]", "", XYXY$group1)

​

notes:

(a) "IV1_Name" is the name you want to give to the IV. It will appear on the figure

(b) group1 is the name that "describeBy" will give to the column with your IVs

(c)   "[:].*", ""    is an example of regular expression. It means "take everything after the colon and replace it with nothing." 

_____________________________

​

​

___Step 3_____________________

​library(ggplot2)
ggplot(XYXY, aes(fill=IV1_Name, y=mean, x=IV2_Name))+

geom_bar(position="dodge", stat="identity", width = .5)+
theme_classic()+
geom_errorbar(aes (ymin = mean - se, ymax = mean + se, width = .2),

position = position_dodge(0.5))+

coord_cartesian(ylim=c(####,###))+
ylab("Dependent Variable Description") 

​

notes:

(a) "mean" and "se" should be column names in XYXY. If you renamed them, then change accordingly.

(b) "coord_cartesian" allows you to set the minimum and maximum values of the y-axis. If you start at zero, replace this line with "scale_y_continuous(expand = c(0, 0), limits = c(0, #####))" to prevent white space.

(c) "ylab" allows you to add a description of your DV, in place of the column name from XYXY

_____________________________

​

​

"If you need to change your labs for x and legend you can use the following:

coord_cartesian(ylim=c(0,5))+

    ylab("Average Trial Scores") +

    xlab("Word Type")+

    labs(fill = "Word Length")"