在contr.sum下重置参考水平

Question

当对比度设置为contr.sum contr.poly时，如何更改参考级别？

玩具数据

df <- data.frame(id = rep(1:5, 4),
                 group = rep(c("A","B"), 10),
                 time = factor(rep(paste0("time", 1:4), each = 5)),
                 score = c(rnorm(5, 0 ,1),
                           rnorm(5, 4, 1),
                           rnorm(5, 8, 1),
                           rnorm(5, 12, 1)))
df

现在，如果您检查time变量的对比度，您会看到每个级别的时间都与参考级别进行比较，即time1

contrasts(df$time)

如果我们更改为contr.sum以通过

获取效果编码

options(contrasts = c("contr.sum", "contr.poly"))

现在将每个其他级别与time4

进行比较

contrasts(df$time)

如何更改此对比度矩阵，以便将2-4级与1级进行比较，同时仍然保留编码之和为零的效果？