##(9)创建区域(东中西)
data<-data[data$provcd20>=0,] #剔除选项为负数的(1831),注:ear201a不行(负数太多)
dim(data)
table(data$provcd20)
data$east=rep(NA,nrow(data)) #东部地区
data$west=rep(NA,nrow(data)) #西部地区
data$central=rep(NA,nrow(data)) #中部地区
for(i in 1:nrow(data)){
if(data$provcd20[i]==11|data$provcd20[i]==12|data$provcd20[i]==13|data$provcd20[i]==21|data$provcd20[i]==31|data$provcd20[i]==32|data$provcd20[i]==33|data$provcd20[i]==35|data$provcd20[i]==37|data$provcd20[i]==44|data$provcd20[i]==46) data$east[i]=1
else data$east[i]=0
if(data$provcd20[i]==15|data$provcd20[i]==45|data$provcd20[i]==50|data$provcd20[i]==51|data$provcd20[i]==52|data$provcd20[i]==53|data$provcd20[i]==54|data$provcd20[i]==55|data$provcd20[i]==61|data$provcd20[i]==62|data$provcd20[i]==63|data$provcd20[i]==64|data$provcd20[i]==65) data$west[i]=1
else data$west[i]=0
if(data$provcd20[i]==14|data$provcd20[i]==22|data$provcd20[i]==23|data$provcd20[i]==34|data$provcd20[i]==36|data$provcd20[i]==41|data$provcd20[i]==42|data$provcd20[i]==43) data$central[i]=1
else data$central[i]=0
}
#