子集R中的回归分类协变量

Question

我正在尝试对名为“公共”的分类协变量进行子集化，该变量在我的回归模型中具有两个级别（0和1）。我想看看当模型设置为“ 1”时，通过设置子变量，该协变量如何影响回归曲线的均值和离散值。不幸的是，我收到一条错误消息：

Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels

我的数据位于名为“ erad.surv”的数据框中，其外观如下：

    Surv(Erad_Time, Event)        Dist Public (weights)
1                       12 -0.51385627      1         1
2                      39+  0.65870948      1         1
3                      39+  0.47563451      0         1
4                      39+  0.73255621      1         1
5                      39+  0.32500333      1         1
6                      39+ -0.60323102      1         1
7                       31  0.60391220      1         1

最后，我的回归模型是下面提供的模型。我似乎还尝试了“子集”命令的其他所有组合都无济于事。

model.Public <- flexsurvreg(formula = Surv(Erad_Time, Event) ~ Dist + Public, data = erad.surv, subset = (Public == "1"), dist=ig)

任何想法或帮助将不胜感激！