Question

每个人，我都会使用SAFETY数据运行逻辑回归。这是我的代码：

library(lmtest)

data <- read.csv("C:/Users/user/Desktop/SAFETY.csv",header=TRUE)(your directory)
data$Type<-factor(data$Type)
data$Size<-factor(data$Size)
data$Region<-factor(data$Region)
mylogit<- glm(Unsafe~Size+Weight+Region, data=data,family = "binomial")
waldtest(mylogit,test="Chisq")

然后就行了。

但是如果：

df <- read.csv("C:/Users/user/Desktop/SAFETY.csv",header=TRUE)
df$Type<-factor(df$Type)
df$Size<-factor(df$Size)
df$Region<-factor(df$Region)
mylogit<- glm(Unsafe~Size+Weight+Region, data=df,family = "binomial")
waldtest(mylogit,test="Chisq")

它显示：

错误＆＃34; ＆＃39;数据＆＃39;参数类型错误＆＃34;

如果我坚持使用df的名字，我可以知道如何解决它吗？谢谢你的帮助。

Answer 1

我无法在R 3.2.2和lmtest 0.9.34上重现此问题（最后请参阅sessionInfo()输出）。

使用此代码：

library(lmtest)
dfS = read.csv("Data/SAFETY.csv")
logitS = glm(Unsafe ~ factor(Size) + Weight + 
               Region, data = dfS, family = binomial())
summary(logitS)
waldtest(logitS, test = "Chisq")

我能够制作：

> summary(logitS)

Call:
glm(formula = Unsafe ~ factor(Size) + Weight + Region, family = binomial(), 
    data = dfS)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.7970  -0.5866  -0.3119   0.7864   2.2018  

Coefficients:
                Estimate Std. Error z value Pr(>|z|)   
(Intercept)       2.7285     1.3949   1.956  0.05046 . 
factor(Size)2    -2.0200     0.6246  -3.234  0.00122 **
factor(Size)3    -2.6785     0.8810  -3.040  0.00236 **
Weight           -0.6678     0.4589  -1.455  0.14559   
RegionN America  -0.3775     0.5624  -0.671  0.50203   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 119.249  on 95  degrees of freedom
Residual deviance:  84.004  on 91  degrees of freedom
AIC: 94.004

Number of Fisher Scoring iterations: 5

> waldtest(logitS, test = "Chisq")
Wald test

Model 1: Unsafe ~ factor(Size) + Weight + Region
Model 2: Unsafe ~ 1
  Res.Df Df  Chisq Pr(>Chisq)    
1     91                         
2     95 -4 23.987  8.035e-05 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

我建议您重新启动R会话并重新加载所有必需的包。如果需要，请执行rm(list = ls())（注意：这将删除工作区中的所有对象）。

> sessionInfo()
R version 3.2.2 (2015-08-14)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] lmtest_0.9-34    zoo_1.7-12       tidyr_0.3.1      samplesize_0.2-2 ggplot2_1.0.1   
[6] lazyeval_0.1.10  RMySQL_0.10.7    DBI_0.3.1        dplyr_0.4.3     

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.1      magrittr_1.5     MASS_7.3-43      munsell_0.4.2    colorspace_1.2-6
 [6] lattice_0.20-33  R6_2.1.1         stringr_1.0.0    plyr_1.8.3       tools_3.2.2     
[11] parallel_3.2.2   grid_3.2.2       gtable_0.1.2     assertthat_0.1   digest_0.6.8    
[16] reshape2_1.4.1   labeling_0.3     stringi_0.5-5    scales_0.3.0     proto_0.3-10

在R中运行具有不同数据名称的Wald测试（逻辑回归）

1 个答案: