# Chapter 5: Making Controlled Comparisons
library(poliscidata) # Loads R companion package in session
# MEAN COMPARISON ANALYSIS WITH A CONTROL VARIABLE
cat("Controlled Mean Comparison Table:\n")
imeansC(~modsex_scale, ~link_wom_scale + dem_raceeth2, nesD)
# Next commands create plot with legend
iplotC(~modsex_scale, ~link_wom_scale+dem_raceeth2, nesD,
modsex_scale~dem_raceeth2+link_wom_scale,
xlab="Linked Fate Scale", ylab="Modern Sexism Scale",
main=list("Modern Sexism Opinions by \n Linked Fate and Race", cex=1.2))
legend("topleft", legend=c("Black","White"), lty=c(2,1), lwd=2,
title="Race", inset=0.1, bty="n")
