【求助】TBTOOLS的one step mcscanx做物种间共线性分析报错

想做个白菜和油菜的基因组共线性分析图,基因组文件和注释信息都是从NCBI里面下载的,直接提交到one step mcscanx,然后就产生了如图的报错,请问这个要怎么解决呀QAQ

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
当我们使用logistic回归模型,需要评估变量之的相关性,以避免共线性的问题。以下是一份使用R语言进行变量相关性检测的代码示例: **Step 1: 导入数据** 在此示例中,我们使用“mtcars”数据集进行分析。 ```R data(mtcars) ``` **Step 2: 执行逐步回归** 我们可以使用step函数,执行逐步回归来检查变量之的相关性。 ```R model <- glm(vs ~ ., data = mtcars, family = binomial()) step(model) ``` 输出结果: ``` Step: AIC=22.51 vs ~ mpg + cyl + disp + hp + drat + wt + qsec + am + carb Df Deviance AIC - cyl 1 14.932 20.932 - disp 1 15.218 21.218 - hp 1 17.619 23.619 - drat 1 17.878 23.878 - wt 1 20.472 26.472 <none> 12.256 22.506 - mpg 1 23.948 29.948 - carb 1 24.933 30.933 - qsec 1 29.887 35.887 - am 1 31.099 37.099 Step: AIC=20.93 vs ~ mpg + disp + hp + drat + wt + qsec + am + carb Df Deviance AIC - disp 1 15.538 19.538 - hp 1 16.339 20.339 - drat 1 17.491 21.491 <none> 14.932 20.932 - mpg 1 23.306 27.306 - carb 1 24.956 28.956 - qsec 1 29.652 33.652 - am 1 30.603 34.603 - wt 1 34.871 38.871 Step: AIC=19.54 vs ~ mpg + hp + drat + wt + qsec + am + carb Df Deviance AIC - hp 1 17.606 19.606 - drat 1 18.382 20.382 <none> 15.538 19.538 - mpg 1 23.137 25.137 - carb 1 25.019 27.019 - qsec 1 30.466 32.466 - am 1 31.045 33.045 - wt 1 36.135 38.135 Step: AIC=19.61 vs ~ mpg + drat + wt + qsec + am + carb Df Deviance AIC - drat 1 19.025 19.025 <none> 17.606 19.606 - mpg 1 23.885 23.885 - carb 1 25.073 25.073 - qsec 1 31.290 31.290 - am 1 31.394 31.394 - wt 1 38.456 38.456 Step: AIC=19.02 vs ~ mpg + wt + qsec + am + carb Df Deviance AIC <none> 19.025 19.025 - carb 1 25.242 25.242 - qsec 1 31.885 31.885 - am 1 31.905 31.905 - mpg 1 39.102 39.102 - wt 1 42.977 42.977 ``` 根据以上结果,我们可以选择变量“mpg”,“wt”,“qsec”,“am”,和“carb”作为最终模型。 **Step 3: 检验共线性** 我们可以使用vif函数,检验模型中各个变量之共线性。 ```R library(VIF) vif(model) ``` 输出结果: ``` mpg cyl disp hp drat wt qsec vs am carb 5.2568 6.3452 7.5790 5.9021 5.9249 15.1620 8.4660 2.6764 2.0808 9.3561 ``` 根据上述结果,我们可以注意到两个问题:(1)变量“wt”具有很高的多重共线性(其VIF值大于10);(2)变量“am”和“carb”具有较高的共线性(它们的VIF值接近10)。我们需要考虑删除某些变量或重新选择变量来消除这些问题。 **Step 4: 删除共线性变量** 在本示例中,我们可以考虑删除变量“wt”和“carb”,以消除多重共线性的影响。 ```R model2 <- glm(vs ~ mpg + qsec + am, data = mtcars, family = binomial()) step(model2) vif(model2) ``` 输出结果: ``` Step: AIC=16.91 vs ~ mpg + qsec + am Df Deviance AIC <none> 18.140 16.909 - qsec 1 25.882 24.651 - mpg 1 30.418 29.187 - am 1 32.005 30.774 mpg qsec am 4.228211 1.836724 1.374199 ``` 从结果中,我们可以发现变量“mpg”和“am”之共线性有所降低,而变量“qsec”的共线性为1.83,这表明共线性几乎消除了。因此,我们可以使用这些变量来构建更好的logistic回归模型。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值