模型：y~a/(1+exp(b-c*x))
现想引入不同区域作为虚拟变量，目前已经完成虚拟变量的设置。
请问该如何在R语言操作中引入？