如何正确实施clogit的对比？

Question

我想使用变量测试的不同基类来运行条件logictic回归。我总是收到消息：参数对比不匹配

我试图制作一个例子。我的实际数据是一个小组。但我无法用相应的ID生成时间间隔。但我希望这是有道理的：

                         a<-seq(1,1000, 1)
                         b<-rnorm(1000, mean = 0, sd = 1)
                         c<-rnorm(1000, mean = 0, sd = 1)
                         d<-rbinom(1000,size=3, prob=.5)
                         df<-data.frame(a,b,c,d)

创造因素

                      df<- within(df,{

                             test<-NA
                             test[d==0 ]<-"very low"
                             test[d==1 ]<-"low"
                             test[d==2]<-"middle"
                             test[d==3]<-"high"
                            })

               df$test<-as.factor(df$test)

               contrasts(df$test)

                             low middle very low
                   high       0      0        0
                   low        1      0        0
                   middle     0      1        0
                   very low   0      0        1

           library(survival)
           clogit<-clogit(b~c+test+strata(a),
           contrasts=list(test = contr.treatment(4,base= 1)),data=df)