r の相互作用項を使用してハードル回帰出力をプロットしようとしていますが、モデルのカウント部分で問題が発生しています。
Model <- hurdle(Suicide. ~ Age + gender + bullying*support, dist = "negbin", link = "logit")
各部分を個別に推定せずにこれらのデータをプロットできるアプローチ (二項ロジットと負の二項カウント) を知らないため、MASS パッケージを使用して推定値を使用してそれぞれをプロットしようとしています。これまでのところ、プロットに visreg パッケージを使用するのが最善でしたが、他の提案も受け付けています。ハードル モデルからの元のロジスティック出力を再現して正常にプロットできましたが、負の二項カウント データはできませんでした (つまり、MASS からのパラメーター推定値は、ハードル回帰出力のものと同じではありません)。
他の人が過去にハードル回帰の結果をどのようにプロットしたか、または MASS で glm.nb を使用してハードル モデルから元々得られた負の二項係数をどのように再現できるかについての洞察をいただければ幸いです。
データをプロットするために使用しているものは次のとおりです。
##Logistic
logistic<-glm(SuicideBinary ~ Age + gender + bullying*support, data = sx, family="binomial")
data("sx", package = "MASS")
##Linear scale
visreg(logistic, "bullying", by="support", xlab = "bullying", ylab = "Log odds (Suicide yes/no)")
##Logistic/probability scale
visreg(logistic, "bullying", by="support", scale = "response", xlab = "bullying", ylab = "P(Initial Attempt)")
##Count model
NegBin<-glm.nb(Suicide. ~ Age + gender + bullying*support, data = sx)
data("sx", package = "MASS")
##Linear scale
visreg(NegBin, "bullying", by="support", xlab = "bullying", ylab = "Count model (number of suicide attempts)")
##Logistic/probability scale
visreg(NegBin, "bullying", by="support", scale = "response", xlab = "bullying", ylab = "P(Subsequent Attempts)")