btvntuan1

BÀI TẬP TUẦN 2

Phân tích dữ liệu chéo của biến định lượng “students” của dataset: “CASchool” lên biến định lượng “income” thuộc package “AER” # Mô tả dữ liệu - Biến định lượng Income - Biến định lượng Students # Đánh giá dữ liệu

library(AER)

## Loading required package: car

## Loading required package: carData

## Loading required package: lmtest

## Loading required package: zoo

## 
## Attaching package: 'zoo'

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

## Loading required package: sandwich

## Loading required package: survival

data("CASchools")
### Tạo dữ liệu bản sau “S” từ “CASchools""
S <- CASchools
### Đặt lại tên cho các thuộc tính 
names(S) <- c("D","SC","C","G","ST","T","CA","L","CO","E","I","EN","R","M")

Biến income

Trung bình

mean(S$I)

## [1] 15.31659

Trung vị

median(S$I, na.rm = FALSE)

## [1] 13.7278

Phương sai

var(S$I)

## [1] 52.21348

Tính tổng

sum(S$I)

## [1] 6432.967

Độ lệch chuẩn

sd(S$I)

## [1] 7.22589

Biến student

Trung bình

mean(S$ST)

## [1] 2628.793

Trung vị

median(S$ST, na.rm = FALSE)

## [1] 950.5

Độ lệch chuẩn

sd(S$ST)

## [1] 3913.105

Phương sai

var(S$ST)

## [1] 15312391

Tính tổng

sum(S$ST)

## [1] 1104093

LẬP BẢNG TẦN SỐ

Biến Thu nhập (income)

Lập bảng tần số của biến thu nhập (income) như sau:

I <- S$I
BangI <- cut(I, breaks = c(5, 15, 25, 35, 50 ), labels = c("thấp", "trung bình" , "trung bình cao", "rất cao"), right = TRUE)
tabincome= table(BangI)
tabincome

## BangI
##           thấp     trung bình trung bình cao        rất cao 
##            255            132             17             14

Nhận xét: Dựa vào bảng tần số ta có thể thấy mức thu nhập Thấp chiếm tỷ lệ rất cao và ngược lại mức thu nhập Rất Cao lại chiếm tỷ lệ thấp nhất.

Lập bảng tần số cho biến Students

ST <- S$ST
BangST <- cut(ST, breaks = c(200,400,600,1000), labels = c("Ít", "Vừa", "Đông"), right = TRUE)
tabde=table(BangST)
tabde

## BangST
##   Ít  Vừa Đông 
##   52   61   46

Phân tích sự tác động giữa 2 biến định lượng

tabearnde = table(BangI, BangST)
tabearnde

##                 BangST
## BangI            Ít Vừa Đông
##   thấp           36  46   26
##   trung bình     15  13   15
##   trung bình cao  0   0    3
##   rất cao         1   2    1

Kiểm định chi bình phương (chi-squared):

###Trong kiểm định này, chúng ta phải kiểm tra các giá trị p-value và đặt ra bài toán kiểm định gồm có giả thuyết H0 và H1. ###Chúng ta sẽ bác bỏ giả thuyết H0 nếu giá trị p-value xuất hiện trong kết quả nhỏ hơn mức ý nghĩa xác định trước, mức ý nghĩa thường là 0,05. ####H0: Hai biến độc lập ####H1: Hai biến phụ thuộc

chisq.test(I, ST, correct=FALSE)

## Warning in chisq.test(I, ST, correct = FALSE): Chi-squared approximation may be
## incorrect

## 
##  Pearson's Chi-squared test
## 
## data:  I and ST
## X-squared = 131671, df = 131040, p-value = 0.1089

Vì giá trị p-value=0.1089 nhỏ hơn mức ý nghĩa nên bác bỏ H0.

Kết luận: Giữa 2 biến có sự liên quan đến nhau hay nói cách khác đây là hai biến phụ thuộc.

BÀI TẬP TUẦN 1

Mô tả dữ liệu

Một khung dữ liệu chứa 420 quan sát trên 14 biến.

district: Mã huyện

school: Tên trường

county: Thành phố

grades: nhân tố biểu thị nhịp độ của huyện.

students: tổng Số học sinh ghi danh

teachers: Số giáo viên

calworks: Phần trăm đủ điều kiện cho hỗ trợ thu nhập

lunch: phần trăm đủ điều kiện cho bữa ăn trưa giảm giá

computer: Số lượng máy tính

expenditure: Chi tiêu cho mỗi sinh viên

income: thu nhập bình quân của huyện (tính theo USD 1.000)

english: tỷ lệ người học tiếng Anh

read: Điểm đọc trung bình

math: Điểm toán trung bình

Xử lý dữ liệu

Khai báo sử dụng package “AER” từ thư viện và sử dụng tập dữ liệu “CASchool”

library(AER)

data("CASchools")
### Tạo dữ liệu bản sau “S” từ “CASchools""
S <- CASchools
### Đặt lại tên cho các thuộc tính 
names(S) <- c("D","SC","C","G","ST","T","CA","L","CO","E","I","EN","R","M")
### Xem cấu trúc dữ liệu “CASchools”
str(CASchools)

## 'data.frame':    420 obs. of  14 variables:
##  $ district   : chr  "75119" "61499" "61549" "61457" ...
##  $ school     : chr  "Sunol Glen Unified" "Manzanita Elementary" "Thermalito Union Elementary" "Golden Feather Union Elementary" ...
##  $ county     : Factor w/ 45 levels "Alameda","Butte",..: 1 2 2 2 2 6 29 11 6 25 ...
##  $ grades     : Factor w/ 2 levels "KK-06","KK-08": 2 2 2 2 2 2 2 2 2 1 ...
##  $ students   : num  195 240 1550 243 1335 ...
##  $ teachers   : num  10.9 11.1 82.9 14 71.5 ...
##  $ calworks   : num  0.51 15.42 55.03 36.48 33.11 ...
##  $ lunch      : num  2.04 47.92 76.32 77.05 78.43 ...
##  $ computer   : num  67 101 169 85 171 25 28 66 35 0 ...
##  $ expenditure: num  6385 5099 5502 7102 5236 ...
##  $ income     : num  22.69 9.82 8.98 8.98 9.08 ...
##  $ english    : num  0 4.58 30 0 13.86 ...
##  $ read       : num  692 660 636 652 642 ...
##  $ math       : num  690 662 651 644 640 ...

Xem thử 6 dòng đầu của dữ liệu “S”

head(S)

##       D                              SC       C     G   ST     T      CA
## 1 75119              Sunol Glen Unified Alameda KK-08  195 10.90  0.5102
## 2 61499            Manzanita Elementary   Butte KK-08  240 11.15 15.4167
## 3 61549     Thermalito Union Elementary   Butte KK-08 1550 82.90 55.0323
## 4 61457 Golden Feather Union Elementary   Butte KK-08  243 14.00 36.4754
## 5 61523        Palermo Union Elementary   Butte KK-08 1335 71.50 33.1086
## 6 62042         Burrel Union Elementary  Fresno KK-08  137  6.40 12.3188
##         L  CO        E         I        EN     R     M
## 1  2.0408  67 6384.911 22.690001  0.000000 691.6 690.0
## 2 47.9167 101 5099.381  9.824000  4.583333 660.5 661.9
## 3 76.3226 169 5501.955  8.978000 30.000002 636.3 650.9
## 4 77.0492  85 7101.831  8.978000  0.000000 651.9 643.5
## 5 78.4270 171 5235.988  9.080333 13.857677 641.8 639.9
## 6 86.9565  25 5580.147 10.415000 12.408759 605.7 605.4

### Lấy biến C của dữ liệu S làm dữ liệu mới
tea <- S$T
### Lấy dữ liệu có điều kiện
tea20 <- tea[tea<20]
tea20den80 <- tea[tea>20&tea<80]
### Xem dòng 3 của dữ  liệu tea
tea[3]

## [1] 82.9

### Lấy dữ liệu khi biến E>6000 hoặc biến I<10
teavv <- S[S$E>6000|S$I<10,]
### Thêm 2 biến 
S$logEN <- log(S$EN)
S$sumR <- sum(S$R)

### Lập bảng tần số của biến Cou thuộc dữ liệu S
Cou <- S$CO
cut(Cou,5)

##   [1] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##   [4] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##   [7] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [10] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [13] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [16] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [19] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
##  [22] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [25] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
##  [28] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [31] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
##  [34] (-3.32,665]         (1.33e+03,1.99e+03] (-3.32,665]        
##  [37] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [40] (-3.32,665]         (1.33e+03,1.99e+03] (-3.32,665]        
##  [43] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [46] (-3.32,665]         (1.33e+03,1.99e+03] (-3.32,665]        
##  [49] (2.66e+03,3.33e+03] (665,1.33e+03]      (665,1.33e+03]     
##  [52] (-3.32,665]         (665,1.33e+03]      (665,1.33e+03]     
##  [55] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [58] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [61] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
##  [64] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [67] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [70] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [73] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
##  [76] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [79] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
##  [82] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [85] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [88] (665,1.33e+03]      (665,1.33e+03]      (-3.32,665]        
##  [91] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [94] (-3.32,665]         (-3.32,665]         (-3.32,665]        
##  [97] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [100] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [103] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [106] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [109] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [112] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
## [115] (-3.32,665]         (-3.32,665]         (1.33e+03,1.99e+03]
## [118] (1.33e+03,1.99e+03] (-3.32,665]         (-3.32,665]        
## [121] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [124] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [127] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [130] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [133] (665,1.33e+03]      (-3.32,665]         (1.33e+03,1.99e+03]
## [136] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [139] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [142] (-3.32,665]         (1.33e+03,1.99e+03] (-3.32,665]        
## [145] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [148] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [151] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [154] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [157] (665,1.33e+03]      (-3.32,665]         (665,1.33e+03]     
## [160] (2.66e+03,3.33e+03] (-3.32,665]         (-3.32,665]        
## [163] (-3.32,665]         (1.99e+03,2.66e+03] (665,1.33e+03]     
## [166] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [169] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
## [172] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [175] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
## [178] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [181] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [184] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [187] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [190] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [193] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [196] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [199] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [202] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [205] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [208] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [211] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [214] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [217] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [220] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
## [223] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [226] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [229] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [232] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [235] (1.99e+03,2.66e+03] (665,1.33e+03]      (-3.32,665]        
## [238] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [241] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [244] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [247] (-3.32,665]         (-3.32,665]         (1.33e+03,1.99e+03]
## [250] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [253] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [256] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [259] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [262] (-3.32,665]         (-3.32,665]         (1.33e+03,1.99e+03]
## [265] (-3.32,665]         (1.33e+03,1.99e+03] (-3.32,665]        
## [268] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [271] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [274] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [277] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [280] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [283] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [286] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [289] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [292] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [295] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [298] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
## [301] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [304] (1.33e+03,1.99e+03] (665,1.33e+03]      (-3.32,665]        
## [307] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [310] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [313] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [316] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [319] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [322] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [325] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [328] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [331] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [334] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [337] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [340] (665,1.33e+03]      (665,1.33e+03]      (-3.32,665]        
## [343] (665,1.33e+03]      (665,1.33e+03]      (-3.32,665]        
## [346] (-3.32,665]         (665,1.33e+03]      (-3.32,665]        
## [349] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [352] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [355] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [358] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [361] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [364] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [367] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [370] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [373] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [376] (665,1.33e+03]      (-3.32,665]         (-3.32,665]        
## [379] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [382] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [385] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [388] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [391] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [394] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [397] (-3.32,665]         (-3.32,665]         (1.99e+03,2.66e+03]
## [400] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [403] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [406] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [409] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [412] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## [415] (-3.32,665]         (-3.32,665]         (665,1.33e+03]     
## [418] (-3.32,665]         (-3.32,665]         (-3.32,665]        
## 5 Levels: (-3.32,665] (665,1.33e+03] ... (2.66e+03,3.33e+03]

table(cut(Cou,5))

## 
##         (-3.32,665]      (665,1.33e+03] (1.33e+03,1.99e+03] (1.99e+03,2.66e+03] 
##                 365                  39                  11                   3 
## (2.66e+03,3.33e+03] 
##                   2