我一直在利用R中的glmrob函数来描述属于某种策略(3类)是否增加成功报告的百分比(0-1)。但是,我的结果显示系数为2.1+,我认为这是不可能的,因为DV只能在0-1的范围内。
我运行的函数的tcc = 1.345和tcc = 3.5,结果没有差异。我不确定我的代码是否有问题,因为我一直遵循我在堆栈交换和堆栈溢出中可以找到的内容。
rob5 <- glmrob(Num_ofSuccess/Total_No_Reports ~ RGamDesStrat3,
weights = Total_No_Reports, family = binomial, data = MyMasterThesisDataRSTUDIOV3,
method = "Mqle", control = glmrobMqle.control(tcc = 3.5))
我的代码示例。
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.5108 0.3266 1.564 0.118
RGamDesStrat3(0,1] 2.2088 0.4744 4.656 3.22e-06 ***
RGamDesStrat3(1,2] 2.1638 0.4346 4.979 6.40e-07 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Robustness weights w.r * w.x:
All 40 weights are ~= 1.
但是,我预计两个预测因子的系数都将低于0.3或0.4左右。