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: survivaldata("CASchools")
S <- CASchools
names(S) <- c("D","SC","C","G","ST","T","CA","L","CO","E","I","EN","R","M")
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