这个看不出来  你可以试试RStudio

> summary(fit.logit)

Call:
glm(formula = is_sp ~ X1day + X4day + X5day + X7day + X10day +
    X13day + X22day + X29day + X31day, family = binomial, data = fp.dau1.cast[,
    -1])

Deviance Residuals:
     Min        1Q    Median        3Q       Max  
-1.95538  -0.45175  -0.23178  -0.06122   2.69461  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)   
(Intercept)  -3.6036     0.4269  -8.441  < 2e-16 ***
X1day1        1.5334     0.5720   2.681  0.00735 **
X4day1        1.7753     0.6424   2.764  0.00572 **
X5day1       -1.0353     0.7622  -1.358  0.17437   
X7day1        1.7002     0.7109   2.392  0.01678 *  
X10day1      -2.6753     0.9418  -2.841  0.00450 **
X13day1       1.3726     0.7547   1.819  0.06893 .  
X22day1       1.6233     0.6382   2.543  0.01098 *  
X29day1       2.0012     0.6480   3.088  0.00201 **
X31day1       1.7310     0.8143   2.126  0.03352 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 281.20  on 251  degrees of freedom
Residual deviance: 126.73  on 242  degrees of freedom
AIC: 146.73

Number of Fisher Scoring iterations: 6

我在Rstudio中运行
options(contrasts = c("contr.treatment", "contr.poly"))
ldose <- rep(0:5, 2)
numdead <- c(1, 4, 9, 13, 18, 20, 0, 2, 6, 10, 12, 16)
sex <- factor(rep(c("M", "F"), c(6, 6)))
SF <- cbind(numdead, numalive = 20 - numdead)
budworm.lg <- glm(SF ~ sex*ldose, family = binomial)
summary(budworm.lg, cor = F)
也出现了“Error in factorial(n) : 缺少参数"n",也没有缺省值”提示。
解决办法：
点Studio菜单栏的“session”，下拉菜单点“terminate R...”，然后点“restart R”,又把命令程序载入，成功了！！！

点Studio菜单栏的“session”,点“terminate R...”，点"restart R",重新加载运行命令，OK！