Scribd
Upload
40 views4 pages

Welfare

Uploaded by

dowemo9017
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
40 views4 pages

Welfare

Uploaded by

dowemo9017
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

welfare <- welfare %>%

mutate(ageg = ifelse(age < 30,"young",


ifelse(age <= 59,"middle","old")))
table(welfare$ageg)
qplot(welfare$ageg)
ageg_income <- welfare %>%
filter(!is.na(income)) %>%
group_by(ageg) %>%
summarise(mean_income = mean(income))
ageg_income

ggplot(data=ageg_income,aes(x=ageg,y=mean_income))+geom_col()
ggplot(data=ageg_income,aes(x=ageg,y=mean_income))+
geom_col()+
scale_x_discrete(limits=c("young","middle","old"))

sex_income <- welfare %>%


filter(!is.na(income)) %>%
group_by(ageg,sex) %>%
summarise(mean_income=mean(income))
sex_income

ggplot(data = sex_income,aes(x=ageg,y=mean_income,fill=sex))+
geom_col()+
scale_x_discrete(limits=c("young","middle","old"))

ggplot(data = sex_income,aes(x=ageg,y=mean_income,fill=sex))+
geom_col(position = "dodge")+
scale_x_discrete(limits=c("young","middle","old"))

sex_age <- welfare %>%


filter(!is.na(income)) %>%
group_by(age,sex) %>%
summarise(mean_income=mean(income))
head(sex_age)

ggplot(data=sex_age,aes(x=age,y=mean_income,col=sex))+geom_line()

class(welfare$code_job)
table(welfare$code_job)
library(readxl)
list_job <- read_excel("Koweps_Codebook.xlsx",col_names = T,sheet=2)
head(list_job)
dim(list_job)
welfare <- left_join(welfare,list_job,id="code_job")
welfare %>%
filter(!is.na(code_job)) %>%
select(code_job,job) %>%
head(10)

job_income <- welfare %>%


filter(!is.na(job)&!is.na(income)) %>%
group_by(job) %>%
summarise(mean_income=mean(income))
head(job_income)

top10 <- job_income %>%


arrange(desc(mean_income)) %>%
head(10)
top10

ggplot(data=top10,aes(x=reorder(job,mean_income),y=mean_income))+
geom_col()+
coord_flip()

bottom10<- job_income %>%


arrange(mean_income) %>%
head(10)
bottom10

ggplot(data=bottom10, aes(x=reorder(job,-mean_income),
y=mean_income))+
geom_col()+
coord_flip()+
ylim(0,850)

job_male<-welfare %>%
filter(!is.na(job)&sex=="male") %>%
group_by(job) %>%
summarise(n=n()) %>%
arrange(desc(n)) %>%
head(10)
job_male

job_female<-welfare %>%
filter(!is.na(job)&sex=="female") %>%
group_by(job) %>%
summarise(n=n()) %>%
arrange(desc(n)) %>%
head(10)
job_female

ggplot(data=job_male,aes(x=reorder(job,n),y=n))+
geom_col()+
coord_flip()

ggplot(data=job_female,aes(x=reorder(job,n),y=n))+
geom_col()+
coord_flip()

class(welfare$religion)
table(welfare$religion)
welfare$religion <- ifelse(welfare$religion==1, "yes", "no")
table(welfare$religion)
qplot(welfare$religion)
class(welfare$marriage)
table(welfare$marriage)
welfare$group_marriage<-ifelse(welfare$marriage==1,"marriage",
ifelse(welfare$marriage==3,"divorce",NA))
table(welfare$group_marriage)
table(is.na(welfare$group_marriage))
qplot(welfare$group_marriage)
religion_marriage<-welfare %>%
filter(!is.na(group_marriage)) %>%
group_by(religion,group_marriage) %>%
summarise(n=n()) %>%
mutate(tot_group=sum(n)) %>%
mutate(pct=round(n/tot_group*100,1))
religion_marriage

religion_marriage<-welfare %>%
filter(!is.na(group_marriage)) %>%
count(religion,group_marriage) %>%
group_by(religion) %>%
mutate(pct=round(n/sum(n)*100,1))

divorce <- religion_marriage %>%


filter(group_marriage == "divorce") %>%
select(religion,pct)
divorce
ggplot(data=divorce,aes(x=religion,y=pct))+geom_col()

ageg_marriage<- welfare %>%


filter(!is.na(group_marriage)) %>%
group_by(ageg,group_marriage) %>%
summarise(n=n()) %>%
mutate(tot_group=sum(n)) %>%
mutate(pct=round(n/tot_group*100,1))
ageg_marriage
ageg_marriage<- welfare %>%
filter(!is.na(group_marriage)) %>%
count(ageg,group_marriage) %>%
group_by(ageg) %>%
mutate(pct=round(n/sum(n)*100,1))
ageg_marriage
ageg_divorce<-ageg_marriage %>%
filter(ageg!="young"&group_marriage=="divorce") %>%
select(ageg,pct)
ageg_divorce
ggplot(data=ageg_divorce,aes(x=ageg,y=pct))+geom_col()

You might also like