需要帮助将二项式（logit）输出到 R 到一个小样本

Question

我昨天坐了四个小时试图从 R 中获取输出，我知道如果将来使用 R，我需要阅读它。我使用我大学的指南下载了 Mac 的 R，但它的指南已经过时了。当我尝试下载最新版本（+ 附加功能）但不断出现错误时。

但是，我有一小部分数据需要广义线性模型 - 二项式 (logit) 输出。任何人都可以帮助我获得它吗？我将非常感谢这样的帮助。如果我也能得到 "codes" 以便我可以从中学习，我也将不胜感激。再次感谢。

sort   temp  survive total

  0     1       29    29 

  0     8       19    20 

  0    10       27    29 

  0    12       30    30 

  0    14       20    29 

  0    18       15    30 

  0    20       11    20 

  1     1       27    28

  1     8       21    28

  1    10       21    26

  1    12       20    27

  1    14        8    27

  1    18        4    30

  1    20        1    10

Answer 1

如果 survive 是 total 次试验的成功次数，temp 和 sort 是您的协变量，并且您的数据框名称是 db 那么:

result <- glm(cbind(survive, total-survive) ~ temp+sort,family=binomial(logit), data=db)

然后看看你得到了什么

summary(result)

给你（根据你的数据）

Call:
glm(formula = cbind(survive, total - survive) ~ temp + sort, 
    family = binomial(logit), data = db)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.6671  -0.6430  -0.3071   0.9203   2.5942  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  5.90102    0.64433   9.158  < 2e-16 ***
temp        -0.31415    0.03823  -8.216  < 2e-16 ***
sort        -1.83193    0.31458  -5.823 5.77e-09 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 158.445  on 13  degrees of freedom
Residual deviance:  19.682  on 11  degrees of freedom
AIC: 62.063

Number of Fisher Scoring iterations: 5

需要帮助将二项式（logit）输出到 R 到一个小样本

Need help getting binomial (logit) output to R to a small sample

r

glm

output