source("http://peach.l.chiba-u.ac.jp/course_folder/jn_plot.R")
protest<-read.csv("http://peach.l.chiba-u.ac.jp/course_folder/protest.csv")
lm <- lm(liking~sexism*protest,data=protest)
summary(lm)
jn_plot(lm, "protest", "sexism", alpha=.05)