Dữ liệu được sử dụng ở đây là từ tất cả 420 học khu K-6 và K-8 ở California với dữ liệu có sẵn cho năm 1998 và 1999. Điểm kiểm tra là bài kiểm tra tiêu chuẩn Stanford 9 dành cho học sinh lớp 5. Các đặc điểm của trường (tính trung bình trên toàn học khu) bao gồm số học sinh ghi danh, số lượng giáo viên (được đo bằng “tương đương toàn thời gian”, số lượng máy tính trên mỗi lớp học và chi tiêu cho mỗi học sinh. Các biến số nhân khẩu học cho học sinh được tính trung bình trên toàn học khu. Các biến số nhân khẩu học bao gồm tỷ lệ học sinh trong chương trình hỗ trợ công cộng CalWorks (trước đây là AFDC), tỷ lệ học sinh đủ điều kiện nhận bữa trưa giảm giá và tỷ lệ học sinh là người học tiếng Anh (nghĩa là học sinh sử dụng tiếng Anh là ngôn ngữ thứ hai).
Một khung dữ liệu chứa 420 quan sát trên 14 biến.
district: tính cách. Mã huyện.
school: tính cách. Tên trường.
county: yếu tố chỉ quận.
grades: nhân tố biểu thị nhịp độ của huyện.
students: tổng Số Ghi Danh.
teachers: Số giáo viên.
calworks: Phần trăm đủ điều kiện cho CalWorks (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.
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: survivalXỬ LÝ DỮ LIỆU
Lấy dữ liệu CASchools từ “AER”
data("CASchools")
#Tạo dữ liệu bản sau “S” từ “CASchools""
S <- CASchools
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 CO 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