To all: I am trying to understand if there is a way to develop my effect plots further. I was wondering if there is a way to get the data points used to create the effect plots onto the effect plot itself, so it sort of looks like a scatter diagram.

I am using the packages: (ggplot2), (stargazer), (effects), (gridExtra), (plyr)

Code Used:

setwd("F:/Sussex PhD/data/constresult19")

my_gen <- read.csv("genhis.csv") #data set name

#linear models

const.qualsvc <- lm(Conpv1019 ~ + Noqual19 + Alevel19 + Degplus19 + Const5ac + Gdschlsrate + FSMKS2 + FSMAlevel, data = my_gen)

const.qualsvl <- lm(Conpv1019 ~ + Noqual19 + Alevel19 + Degplus19 + Const5ac + Gdschlsrate + FSMKS2 + FSMAlevel data = my_gen)

stargazer(const.qualsvc, const.qualsvl, type = "text")

#effect plots

eddeg.cona <- effect(term = "Degplus19", mod = const.qualsvc)

eddeg.cona <- plot(eddeg.cona, ylim=c(-20,30), xlab = "Percent Degree Qualification Plus", ylab = "Percentage Change in Con vote 2010-19", main="Effect of Degrees on Change in Con % Vote", col.line = "blue")

eddeg.laba <- effect(term = "Degplus19", mod = const.qualsvl)

eddeg.laba <- plot(eddeg.laba, ylim=c(-20,35), xlab = "Percent Degree Qualification Plus", ylab = "Percentage Change in Lab vote 2010-19", main="Effect of Degrees on Change in Lab % Vote",col.line = "red")

grid.arrange(eddeg.cona, eddeg.laba, ncol=2) #effect plots side by side

Kind Regards, James Prentice (Only a fairly basic level user of R)