plot(data$praise_rejection, data$vignette_response, main="altruism x rejection avoidance by social norms",
xlab="Praise Rejection score ", ylab="Altruism likelihood response ", pch=19)
with(data, plot(data$praise_rejection, data$vignette_response, col = factor(data$condition)))