我正在使用Tweedie GLM模型。在训练集上使用分类和连续预测变量建立模型后,我得到了系数,然后使用它们来预测测试集上的数据。
NonPenalizedTweedie <-glm(ModelingVars1,
family = tweedie(var.power=1.5, link.power = 0),
weight = Units, data = data1)
但是,在行业中,由于业务原因,我们经常需要手动更改一些系数,然后在原始截距/其他系数保持不变的情况下重新运行更新的模型。
data1$Pred <-predict(NonPenalizedTweedie, data=data1, type="response")
我的问题是:有时我只想更改1个类别变量的级别,然后在其余系数保持不变的情况下重新运行模型。我对此进行了无休止的研究,但找不到有效的方式来进行此操作的好方法。我希望能够为多个分类级别手动更改任何可变系数,然后连续更改,然后使用新系数重新运行分析以查看拟合图和提升图。
过去,我只是从R中获取输出,然后将其放到Excel中,在这里我可以动态更改电平系数,然后拟合图/提升图将自动更新,但是我想找到一种很好的动态方法在R中执行此操作如果可能的话。
我尝试过的一些早期想法:
coeffs <- exp(NonPenalizedTweedie$coefficients)
coeffs["data1$Var1Level1"] <- 1.2
coeffs["data1$Var1Level2"] <- 1.2
coeffs["data1$Var2Level3"] <- 1.2
任何想法都很棒。谢谢。