dglm パッケージを使用して、R に double glm を適合させようとしています。これは、tweedie モデルを使用するために statmod パッケージと組み合わせて使用されます。問題の再現は次のとおりです。
library(dglm)
library(statmod)
p <- 1.5
y <- runif(10)
x <- runif(10)
dglm(y~x,~x,family=tweedie(link.power=0, var.power=p))
#doesnt work
dglm(y~x,~x,family=tweedie(link.power=0, var.power=1.5))
#works
変数のすべてのエントリで dglm が実行されるループを使用したいので、var.power を変数で定義する必要があります。