我目前正在进行回归分析，该回归通过几个腐败指标解释了只有一家公司参与招标的可能性。我当前的回归是：
glm(Single_Bid ~ Announcement_NoCall + Country_NL + Extr_Adv + Extr_Dec + Legal_NoCall + Procedure_Open + Type_Regionaal + Year2012 + Year2013 + Year2014 + Year2015 + Year2017 + Year2018 + Year2019, data = Dataset, family = binomial('logit) )
和我用于健壮性检查的回归是：
Rgr_robust<-glmrob(SingleBid ~ ExtrAdv + ExtrDec + Year2012 + Year2013 + Year2014 + Year2015 + Year2016 + Year2017 + Year2018 + Year2019 + Open + AnnouncNoCall + Local + LegalNoCall + CountryNL, data=Dataset_ALL, family=binomial('logit'), method="BY")
我想检查模型的健壮性。由于这是逻辑回归，因此我读到需要使用glmrob方法。 但是，当我执行此操作时，使用“ BY”方法（我读到这是否适合二进制逻辑回归？），我收到以下警告消息：
1：In（grad.BY％*％xistart）* xistart：回收长度为1的数组 在数组矢量算法中已弃用。使用c（）或as.vector（） 代替。
2：In（grad.BY％*％xistart）* xistart：回收长度为1的数组 在数组矢量算法中已弃用。使用c（）或as.vector（） 代替。
3：In（grad.BY％*％xistart）* xistart：回收长度为1的数组 在数组矢量算法中已弃用。使用c（）或as.vector（） 代替。
4：In（grad.BY％*％xistart）* xistart：回收长度为1的数组 在数组矢量算法中已弃用。使用c（）或as.vector（） 代替。
这些是我的二进制逻辑回归的估计结果：
> ESTIMATE STD. ERROR Z VALUE PR(>|Z|)
>(INTERCEPT -1.39717 0.37549 -3.721 0.000199 ***
>EXTRADV -0.13305 0.09656 -1.378 0.168242
>EXTRDEC 0.22597 0.09678 2.335 0.019545 *
>COUNTRYNL -1.02834 0.25186 -4.083 4.45e-05 ***
>YEAR2012 -1.09127 0.34454 -3.167 0.001538 **
>YEAR2013 -0.86818 0.28519 -3.044 0.002333 **
>YEAR2014 -0.92484 0.28996 -3.190 0.001425 **
>YEAR2015 -1.06884 0.29242 -3.655 0.000257 ***
>YEAR2016 -0.48629 0.29214 -1.665 0.095990 .
>YEAR2017 -0.69089 0.31622 -2.185 0.028900 *
>YEAR2018 -0.12943 0.28781 -0.450 0.652928
>YEAR2019 0.11843 0.38106 0.311 0.755960
>RESTRICTED -0.71202 0.18031 -3.949 7.85e-05 ***
>ANNOUNCNOCALL -0.98108 0.73439 -1.336 0.181576
>LOCAL -0.22801 0.09274 -2.459 0.013951 *
>LEGALNOCALL 0.34881 0.18666 1.869 0.061670