R语言中,BMI的有序多分类变量的亚组分析能否做交互作用呢?

既往文章的交互作用多看到的是,二分类变量在不同亚组的交互作用,如下图:y是死亡,x为是否高膳食多样化,x和y在年龄、性别等不用亚组中的HR以及p for interaction交互作用的p值。

我想做的比如是:y是死亡,x是BMI多分类分组:过瘦、正常、超重、肥胖4个组,x和y与年龄分组的交互作用。

1. R语言代码这样对吗?图中红圈就是它的交互作用p值嘛?(我不确定多分类变量是不是也和二分类变量一样的代码做法)

glm1<-glm(mortality~BMIgroup*agegroup+age+sex+highschool+widow+income+smoke+drink+comorbidity+polypharmacy,

      data=a,

      family = binomial)

glm2<- summary(glm1);glm2

2.这里BMI group的分组值写0.1.2.3可以嘛,还是一定要1.2.3.4?

是否数值一定需要从小到大的分组过瘦1→正常2→超重3→肥胖4?如果我想把正常作为参照组,正常1→过瘦2→超重3→肥胖4,这样能行嘛?

问题有些多,辛苦老师帮忙解答下了~!万分感谢