在R中使用带有比例/百分比因变量的glmrob函数

时间:2019-06-12 07:27:48

标签: r percentage glm

我一直在利用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左右。

0 个答案:

没有答案