Bài tập về nhà

Bài tập tuần 2

Mô tả dữ liệu: Dữ liệu có nguồn gốc từ cuộc khảo sát sức khỏe SOMIPOPS cho Thụy Sĩ vào năm 1981. Dữ liệu gồm 872 quan sát và 7 biến:

Participation: Cá nhân có tham gia vào lực lượng lao động không

Income: thu nhập người không lao động

age: tuổi

education: số năm học

youngkids: số trẻ nhỏ (dưới 7 tuổi)

oldkids: số trẻ lớn (trên 7 tuổi)

foreign: có phải người nước ngoài không

# Đầu tiên chúng ta lấy dữ liệu từ gói package AER
library(AER)

## Warning: package 'AER' was built under R version 4.2.3

## Loading required package: car

## Warning: package 'car' was built under R version 4.2.3

## Loading required package: carData

## Warning: package 'carData' was built under R version 4.2.3

## Loading required package: lmtest

## Warning: package 'lmtest' was built under R version 4.2.3

## Loading required package: zoo

## Warning: package 'zoo' was built under R version 4.2.3

## 
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric

## Loading required package: sandwich

## Warning: package 'sandwich' was built under R version 4.2.3

## Loading required package: survival

# Tiếp theo lấy dữ  "SwissLabor" từ package "AER"
data("SwissLabor")
# Tạo dữ liệu bản sao "kt" từ dữ liệu gốc "SwissLabor"
kt <- SwissLabor
kt <- na.omit(kt)
# Xem cấu trúc dữ liệu của dữ liệu "kt"
str(kt)

## 'data.frame':    872 obs. of  7 variables:
##  $ participation: Factor w/ 2 levels "no","yes": 1 2 1 1 1 2 1 2 1 1 ...
##  $ income       : num  10.8 10.5 11 11.1 11.1 ...
##  $ age          : num  3 4.5 4.6 3.1 4.4 4.2 5.1 3.2 3.9 4.3 ...
##  $ education    : num  8 8 9 11 12 12 8 8 12 11 ...
##  $ youngkids    : num  1 0 0 2 0 0 0 0 0 0 ...
##  $ oldkids      : num  1 1 0 0 2 1 0 2 0 2 ...
##  $ foreign      : Factor w/ 2 levels "no","yes": 1 1 1 1 1 1 1 1 1 1 ...

# Phân tích biến "Participation":
## Tóm tắt như sau:
summary(kt$participation)

##  no yes 
## 471 401

## Tỉ lệ nam nữ như sau:
prop.table(table(kt$participation))*100

## 
##       no      yes 
## 54.01376 45.98624

Ta thấy tỉ lệ không tham gia lao động chiếm 54,01%, còn tỉ lệ tham gia lao động chiếm 45,99%

#Phân tích biến income
kt$inco <- cut(kt$income, breaks = c(0,8,10,12), labels = c("thunhapthap","thunhapvua","thunhapcao"))
summary(kt$inco)

## thunhapthap  thunhapvua  thunhapcao        NA's 
##           1          29         837           5

prop.table(table(kt$inco))

## 
## thunhapthap  thunhapvua  thunhapcao 
## 0.001153403 0.033448674 0.965397924

plot(kt$inco)

#Tỉ lệ thu nhập cao chiếm nhiều nhất khoảng 96,54%
#Tỉ lệ thu nhập vừa chiếm khoảng 3,33%
#Tỉ lệ thu nhập thấp chiếm ít nhất khoảng 0,11%

#phân tích biến age
kt$ag <- cut(kt$age, breaks = c(1,3,5,7), labels = c("thanhnien","trungnien","caonien"))
summary(kt$ag)

## thanhnien trungnien   caonien 
##       183       526       163

prop.table(table(kt$ag))*100

## 
## thanhnien trungnien   caonien 
##  20.98624  60.32110  18.69266

plot(kt$ag)

#Tỉ lệ thanh niên chiếm 20,99%
#Tỉ lệ trung nhiên chiếm nhiều nhất khoảng 60,32%
#Tỉ lệ cao niên chiếm ít nhất khoảng 18,69%

summary(kt$education)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.000   8.000   9.000   9.307  12.000  21.000

#Gọi kt$educ là nhóm bậc các năm học "education":
kt$educ <-cut(kt$education, breaks = c (12,14,16,18), labels = c("HThấp","Hvừa","Hcao"))
#Tỷ lệ của nhóm "kt$educ":
prop.table(table(kt$educ))

## 
##      HThấp       Hvừa       Hcao 
## 0.70103093 0.25773196 0.04123711

plot(kt$educ)

#Nhóm HThấp:70,1%
#Nhóm Hvừa: 25,77%
#Nhóm Hcao: 4,13%

#Phân tích biến youngkids
kt$yk <-cut(kt$youngkids, breaks = c(-1,1,3))
prop.table(table(kt$yk))

## 
##     (-1,1]      (1,3] 
## 0.93119266 0.06880734

# số gia đình có đứa trẻ dưới 7 tuổi khoảng 0 đến 1 đứa là 93,12%
# Số gia đình có đứa trẻ dưới 7 tuổi khoảng 2 đến 3 đứa là 6,88%

#Phân tích biến oldkids
kt$ok <-cut(kt$oldkids, breaks = c(-1,1,3,6))
prop.table(table(kt$ok))

## 
##    (-1,1]     (1,3]     (3,6] 
## 0.6777523 0.3016055 0.0206422

#Số gia đình có đứa trẻ trên 7 tuổi khoảng 0  đến 1 đứa chiếm cao nhất với 67,78%
#Số gia đình có đứa trẻ trên 7 tuổi khoảng 2 đến 3 đứa chiếm 30,16%
#Số gia đình có đứa trẻ trên 7 tuổi trên 3 đứa chiếm thấp nhất khoảng 2,06%

#Phân tích biến Foreign
summary(kt$foreign)

##  no yes 
## 656 216

prop.table(table(kt$foreign))

## 
##        no       yes 
## 0.7522936 0.2477064

#Tỉ lệ người nước ngoài chiếm 24,77%
#Tỉ lệ người bản xứ chiếm 75,23%

Bài tập về nhà

Ngọc Toàn

2023-06-11