我真的需要這方面的幫助。我想為我的 glm quasipoisson 制作一個預測模型。我有一個問題，因為我錯誤地用我的資料集制作了 glm 模型。我曾經根據我的glm quasipoisson對我的所有引數進行預測模型，但我最終對每個引數進行了預測，結果與glm quasipoisson資料不同。
這是我的資料集。我對所有資料集使用 csv 檔案。我知道如何在這篇文章中上傳這個 csv 資料，請原諒我。
Richness = as.matrix(dat1[,14])
Richness
8
3
3
4
3
5
4
3
7
8
Parameter = as.matrix(dat1[,15:22])
Parameter
JE Temp Hmdt Sond HE WE L MH
1 31.3 93 63.3 3.89 4.32 80 7.82
2 26.9 92 63.5 9.48 8.85 60 8.32
1 27.3 93 67.4 1.23 2.37 60 10.10
3 31.6 99 108.0 1.90 3.32 80 4.60
1 29.3 99 86.8 2.42 7.83 460 12.20
2 29.4 85 86.1 4.71 15.04 200 10.10
1 29.4 87 93.5 3.65 14.70 200 12.20
1 29.5 97 87.5 1.42 3.17 80 4.07
1 25.9 95 62.3 5.23 16.89 140 10.03
1 29.5 95 63.5 1.85 6.50 120 6.97
Rich = glm(Richness ~ Parameter, family=quasipoisson, data = dat1)
summary(Rich)
Call:
glm(formula = Richness ~ Parameter, family = quasipoisson, data = dat1)
Deviance Residuals:
1 2 3 4 5
-0.017139 0.016769 -0.008652 0.002194 -0.003153
6 7 8 9 10
-0.016828 0.022914 -0.013823 -0.012597 0.030219
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -7.4197959 0.5061733 -14.659 0.0434 *
ParameterJE 0.1833651 0.0224198 8.179 0.0775 .
ParameterTemp 0.2441301 0.0073380 33.269 0.0191 *
ParameterHmdt 0.0393258 0.0032176 12.222 0.0520 .
ParameterSond -0.0319313 0.0009662 -33.050 0.0193 *
ParameterHE -0.0982213 0.0060587 -16.212 0.0392 *
ParameterWE 0.1001758 0.0027575 36.329 0.0175 *
ParameterL -0.0014170 0.0001554 -9.117 0.0695 .
ParameterMH 0.0137196 0.0073704 1.861 0.3138
---
Signif. codes:
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for quasipoisson family taken to be 0.002739787)
Null deviance: 7.8395271 on 9 degrees of freedom
Residual deviance: 0.0027358 on 1 degrees of freedom
AIC: NA
Number of Fisher Scoring iterations: 3
這是我嘗試用 ggplot 制作的模型
ggplot(dat1, aes(Temp, Richness))
geom_point()
geom_smooth(method = "glm", method.args = list(family = quasipoisson),
fill = "grey", color = "black", linetype = 2)``
這就是結果。
我為每個引數制作，但我只知道這個結果是錯誤的，因為它為每個引數使用了準泊松資料，我想要的是基于準泊松資料的預測模型，如上面的摘要。
我嘗試使用帶有 95% CIs 的多個解釋的結果 glm 中的代碼，但我真的很困惑如何像那里的示例一樣設定我的資料。但是那個例子中的結果幾乎就像我想要的那樣。誰能幫我這個？如何使用 ggplot 將所有引數的 glm 預測模型放在一幀中？希望任何人都可以幫助我解決這個問題。太感謝了！
uj5u.com熱心網友回復：
您是否嘗試過 sjplot 包中的 plot_model 函式？我是用手機寫的，但代碼是這樣的。
庫（sjPlot）
plot_model（glm_model）
更多資訊：http : //www.strengejacke.de/sjPlot/reference/plot_model.html
代碼：
data("mtcars")
glm_model<-glm(am~.,data = mtcars)
glm_model
library(sjPlot)
plot_model(glm_model, vline.color = "red")
plot_model(glm_model, show.values = TRUE, value.offset = .3)
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/403332.html
標籤：