readstata.R

library(haven)
library(sjPlot)

## Learn more about sjPlot with 'browseVignettes("sjPlot")'.

library(sjmisc)
cgss2017 <- read_dta("E:/course/cgss2017/cgss2017.dta")
#View(cgss2017)
table(cgss2017$a15)

## 
##    1    2    3    4    5   98   99 
##  593 2014 3261 4409 2300    3    2

c2017 <- sjmisc::to_label(cgss2017)
table(c2017$a15)

## 
##   很不健康 比较不健康       一般   比较健康     很健康     不知道   拒绝回答 
##        593       2014       3261       4409       2300          3          2

c2017[c2017=="拒绝回答"] <- NA
c2017[c2017=="不知道"] <- NA
c2017[c2017=="不适用"] <- NA
c2017[c2017=="无法回答"] <- NA
c2017[c2017=="无法选择"] <- NA
c2017[c2017=="不好说"] <- NA
c2017 <- droplevels(c2017)
table(c2017$a15)

## 
##   很不健康 比较不健康       一般   比较健康     很健康 
##        593       2014       3261       4409       2300