# \donttest{
require(scam)
require(emmeans)
set.seed(2)
n <- 100
x <- sort(runif(n,-.5,.5))
f <- -20*x^2
fac <- as.factor(sample(1:3,n,replace=TRUE))
y <- ifelse(fac == 1, 4, ifelse(fac == 2, 1, -2)) +
f + rnorm(n, 0, 1.5)
dat <- data.frame(y=y,x=x,fac=fac)
m <- scam(y ~ fac+s(x,bs="cv"), data=dat)
em <- emmeans(m, ~ fac)
em
pairs(em)
# }
Run the code above in your browser using DataLab