全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
1764 3
2021-05-19
背景:利用环境指标解释同种动物在不同环境下的性状(因变量)差异。采用的方法:利用广义线性模型构建模型,通过AIC准则筛选最优模型,通过对△AICc<2的模型进行模型平均。筛选出来的最优模型中包含作用不显著的环境因子(Pr(>|t|)=0.09753 .)作用不显著的环境因子也能包含在最优模型里面吗?下面是我的模型构建代码和最优模型结果
##############################线性模型构建######################################
global.model <- lm(PC2~ LD+FZZ+GM+DLWH, data=Ant.aic)             #模型构建
ant.model <- glmulti(global.model,level = 1,crit="aicc")              #选用AICc进行评判模型  #glmulti()函数在所有可能的模型(由用户指定的候选集)中找出n个最佳模型(模型的置信度集)
summary(ant.model)
#summary(ant.model)$icvalue                                               #查看模型的AICc结果
weightable(ant.model)

f <- lm(PC2 ~ 1 + LD + FZZ + GM + DLWH, data=Ant.aic)  #手动将筛选的bestmodel添加到函数中。
summary(f)                                                                     #查看函数结果。


结果:
Call:
lm(formula = PC2 ~ 1 + LD + FZZ + GM + DLWH, data = Ant.aic)

Residuals:
    Min      1Q  Median      3Q     Max
-8.4534 -0.3036  0.0841  0.3745  1.6864

Coefficients:
             Estimate Std. Error t value Pr(>|t|)   
(Intercept) -0.466486   0.214704  -2.173  0.03125 *
LD           0.017273   0.010365   1.667  0.09753 .               #该环境因子影响不显著,可以纳入最优模型吗??????
FZZ          0.005514   0.002955   1.866  0.06381 .
GM           0.009212   0.003455   2.667  0.00844 **
DLWH        -0.800702   0.297627  -2.690  0.00788 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.9671 on 163 degrees of freedom
Multiple R-squared:  0.08703,        Adjusted R-squared:  0.06463
F-statistic: 3.885 on 4 and 163 DF,  p-value: 0.004844

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2021-5-21 08:58:25
不显著就不用加了,加了会使模型不必要的复杂。除了看显著程度,你还可以看看它和其他因子之间的correlation,如果是相关太大,也可以不加进去模型里面。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-5-21 14:51:22
jcww918 发表于 2021-5-21 08:58
不显著就不用加了,加了会使模型不必要的复杂。除了看显著程度,你还可以看看它和其他因子之间的correlatio ...
谢谢回答。模型中预测因子是通过相关性分析后,排除相关性高的预测因子而获得的。筛选出来的最优模型(bestmodel)中就包含了不显著的预测因子。函数已经筛选出来解释因子的组合,在最优的组合中可以直接筛除掉吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2021-5-22 08:56:21
15182952842 发表于 2021-5-21 14:51
谢谢回答。模型中预测因子是通过相关性分析后,排除相关性高的预测因子而获得的。筛选出来的最优模型(be ...
首先你现在的最优模型的R^2只有0.06,会不会在你选的过程中出现了问题?其次这是不知道你可不可以加上interaction term,如能,你可以试试把环境和其他的因子做两项,甚至三项的interaction,说不定AIC, R^2会有进步。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群