如何在R中有效编辑Tweedie GLM系数

时间:2018-09-25 12:20:15

标签: r glm coefficients tweedie

我正在使用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

任何想法都很棒。谢谢。

0 个答案:

没有答案