I am using visreg to visualize a continuous by continuous interaction. The model is
model1 = DayOfYear ~ year*latitude
where year and latitude are both continuous. By default, visreg plots cross-sections at the 10th, 50th, and 90th quantiles. It is my understanding that all of the data are used to calculate the slopes at these values. I want to know how visreg breaks up the points for visualization into the separate panels.
Here is my visreg code
visreg(model1, "year", by="decimalLatitude", breaks = 3,
partial = TRUE, alpha = 0.05, overlay=FALSE,
gg = TRUE, points = list(pch = 20, alpha = .5), xlab = "",
ylab = "")