R中glm函数中提示缺少参数"n",也没有缺省值怎么办

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

dau<-read.csv("G:/Rwork/《数据分析实战》书中数据和R代码/R/section7-dau.csv",header = T,stringsAsFactors = F)
fp.dau1<-dau[dau$device == "FP"&dau$region_month == "2013-01",c("region_month","region_day","user_id","device")]
fp.dau1$is_assess <- 1
fp.dau1.cast<-dcast(fp.dau1,user_id ~ region_day,value.var = "is_assess",function(x) as.character(length(x)))
names(fp.dau1.cast)[-1]<-paste0("X",1:31,"day")
fp.dau1.cast<-merge(fp.dau1.cast,fp.mau1[,c("user_id","is_sp")],by.fp.dau1.cast = "user_id")

fit.logit<- step(glm(is_sp ~ . , data = fp.dau1.cast[,-1], family = binomial))
这段代码在R中运行到最后一段逻辑回归是报错，内容：Error in factorial(n) : 缺少参数"n",也没有缺省值，一直没有找到问题所在，麻烦大家帮忙找下问题，谢谢！！
section7-dau.zip (237.42 KB) 本附件包括：

• section7-dau.csv

2018-4-16 21:36:29 上传

R中glm函数中提示缺少参数"n",也没有缺省值怎么办

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：glm函数

我运行了是正确的 你可以截图看看

Roivan 发表于 2018-5-8 18:30
我运行了是正确的 你可以截图看看

不会吧，我运行好多次都一直报错

Roivan 发表于 2018-5-8 18:30
我运行了是正确的 你可以截图看看

我把截图以附件的形式上传了，你看下，一直是这个提示

黑和白 发表于 2018-5-10 22:48
不会吧，我运行好多次都一直报错

我已经把运行后R提示的信息截图了，你看下

这个看不出来  你可以试试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！