BÀI TẬP VỀ NHÀ SỐ 1

BÀI TẬP VỀ NHÀ SỐ 1

TravelMode - Package:AER

Mô tả: Dữ liệu về lựa chọn phương thức di chuyển cho hành trình giữa Sydney, Melbourne và Úc

Một data frame chứa 840 quan sát trên 4 chế độ cho 210 mẫu, gồm có 9 biến

individual: Hệ số biểu thị cá nhân có cấp độ từ 1 đến 210

mode: Yếu tố biểu thị phương thức di chuyển với các cấp độ “ô tô”, “hàng không”, “tàu hỏa” hoặc “xe buýt”

choice: Yếu tố thể hiện sự lựa chọn với mức độ “không” và “có”.

wait: Thời gian chờ đợi của nhà ga, số 0 thể hiện cho ô tô

vcost: Thành phần chi phí phương tiện

travel: Thời gian di chuyển trên xe.

gcost: Đo lường chi phí chung.

income: Thu nhập hộ gia đình.

size: Quy mô bữa tiệc

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

# Lấy dữ liệu của datashet TravelMode
data("TravelMode")
# Gán dữ liệu của TravelMode cho biến d
d <- TravelMode
# Mô tả chi tiết kiểu biến số của datashet TravelMode
str(d)

## 'data.frame':    840 obs. of  9 variables:
##  $ individual: Factor w/ 210 levels "1","2","3","4",..: 1 1 1 1 2 2 2 2 3 3 ...
##  $ mode      : Factor w/ 4 levels "air","train",..: 1 2 3 4 1 2 3 4 1 2 ...
##  $ choice    : Factor w/ 2 levels "no","yes": 1 1 1 2 1 1 1 2 1 1 ...
##  $ wait      : int  69 34 35 0 64 44 53 0 69 34 ...
##  $ vcost     : int  59 31 25 10 58 31 25 11 115 98 ...
##  $ travel    : int  100 372 417 180 68 354 399 255 125 892 ...
##  $ gcost     : int  70 71 70 30 68 84 85 50 129 195 ...
##  $ income    : int  35 35 35 35 30 30 30 30 40 40 ...
##  $ size      : int  1 1 1 1 2 2 2 2 1 1 ...

# Lấy dữ liệu của 6 dòng đầu của biến d
head(d)

##   individual  mode choice wait vcost travel gcost income size
## 1          1   air     no   69    59    100    70     35    1
## 2          1 train     no   34    31    372    71     35    1
## 3          1   bus     no   35    25    417    70     35    1
## 4          1   car    yes    0    10    180    30     35    1
## 5          2   air     no   64    58     68    68     30    2
## 6          2 train     no   44    31    354    84     30    2

# Gán tên viết tắt cho các biến để dễ thao tác
names(d) <- c('I','M','C','W','V','T','G','IN','S')
# Lấy dữ liệu của biến V
Vco <- d$V
# Chi phí phương tiện lớn hơn 50
Vco50 <- Vco[Vco>50]
# Chi phí phương tiện lớn hơn 30 và nhỏ hơn 100
Vco30100 <- Vco[Vco>30&Vco<100]
# Lấy dữ liệu của biến T
Tra <- d$T
# Lấy dữ liệu của biến G
Gco <- d$G
#Xem dữ liệu thứ 10 của biến G
Gco[10]

## [1] 195

# Tổng giá trị của V,T và G
tongTVG <- d$T+d$V+d$G
tongTVG

##   [1]  229  474  512  220  194  469  509  316  369 1185 1084  844  176  459  501
##  [16]  217  286  529  569  707  229  422  488  339  423 1269 1168  982  410 1090
##  [31] 1075  965  229  474  512  267  191  457  495  247  194  443  488  250  249
##  [46]  494  534  300  222  483  525  286  224  476  516  358  361 1232 1174 1416
##  [61]  260  399  532  416  204  359  483  375  229  375  481  397  282  374  533
##  [76]  505  223  408  449  382  203  380  468  374  395  626  751  670  400  803
##  [91]  718  706  362  801  716  706  242  431  474  346  331 1249 1148 1211  272
## [106]  848  763  722  331  842  828  720  354  454  440  388  229  338  488  341
## [121]  230  367  515  355  227  424  507  399  202  392  507  380  198  389  478
## [136]  374  267  395  466  367  250  406  521  374  336  463  590  514  177  367
## [151]  480  367  273  465  560  412  405  943  859  753  372  801  716  688  409
## [166]  811  726  699  432  870  785  755  321  857  772  699  342  787  772  689
## [181]  423  874  790  744  295  825  740  691  223  466  504  373  236  428  474
## [196]  371  348  943  859  755  307  929  844  764  308  877  829  692  256  453
## [211]  526  394  220  436  530  403  174  378  481  372  404 1032 1071 1080  496
## [226] 1031 1153 1174  389 1044 1101 1029  397  994 1077 1055  313 1049 1053 1017
## [241]  427 1078 1107 1137  623 1058 1129 1137  374 1080 1172 1060  359  931  996
## [256] 1010  421  967 1102 1045  349 1202  936 1083  483 1474 1372 1158  412 1190
## [271] 1020 1063  370 1058 1152 1054  423 1031 1102 1072  437  508 1187 1087  341
## [286]  903 1061 1030  481  966 1274 1166  430  972 1024 1108  432 1143 1136 1058
## [301]  373 1113 1113 1130  344 1047 1020 1031  413 1194 1178  881  321 1184 1134
## [316] 1678  396 1181  956 1044  468 1202 1240 1089  203  454  456  217  226  377
## [331]  483  402  245  494  519  387  291  825  740  697  194  390  500  383  254
## [346]  436  441  286  202  374  493  392  254  466  465  317  205  382  500  383
## [361]  195  395  437  372  253  437  523  432  292  826  741  445  393  457  497
## [376]  367  357 1251 1149 1059  367  928  843  746  284  487  395  414  537 1270
## [391] 1169 1071  454 1251 1343 1111  418 1261 1187  932  241  479  449  435  327
## [406] 1165 1115  855  225  437  481  359  396 1188 1087 1148  431 1253 1152 1698
## [421]  406 1262  937 1082  416 1177 1076  866  324 1172 1114 1024  286  529  569
## [436]  327  358 1207 1149  853  254  486  524  290  447 1268 1166  867  308  460
## [451]  500  377  220  647  687  479  265  520  558  323  248  474  512  399  189
## [466]  455  308  381  463 1168 1017 1067  407 1184 1083 1006  235  457  495  384
## [481]  265  431  477  348  508 1219 1118  995  404 1223 1122 1127  369 1161 1146
## [496] 1036  405 1278  824 1076  402 1199  957 1078  365 1269 1196 1028  376 1172
## [511]  970 1044  236  474  512  359  259  507  518  435  391 1191 1090  764  337
## [526] 1161 1088  792  433 1254 1106 1077  371 1079  859 1052  390 1315 1213  828
## [541]  413 1308 1292  867  380 1195 1094  903  387 1194 1072 1080  316 1139  957
## [556]  985  287 1148 1075  715  390  987 1089 1054  497 1070 1088 1056  416 1093
## [571] 1205 1081  376 1006 1125 1044  402  926 1061 1039  390  949 1013 1052  383
## [586]  903 1066 1033  395  923 1151 1072  477 1116 1172 1149  487 1197 1096 1065
## [601]  403 1096 1081 1078  481 1329 1189 1162  596 1266 1164 1146  449 1178 1077
## [616] 1046  258  476  516  400  229  468  509  371  537 1310 1168 1143  380 1334
## [631] 1233 1180  275  505  545  396  246  547  589  345  240  442  480  366  293
## [646]  473  511  376  321 1018  957  563  297  884  836  661  295  837  777  603
## [661]  316  916  831  731  308  821  772  634  286  839  801  551  267  816  731
## [676]  548  309  876  815  700  347 1023  938  590  203  454  304  376  267  502
## [691]  308  417  311  850  713  701  248  504  542  234  528 1347 1023 1201  215
## [706]  461  369  357  244  493  311  396  242  837  752  726  301  463  503  383
## [721]  209  496  534  357  252  490  528  327  205  525  510  389  376  565  473
## [736]  388  201  382  496  377  279  827  778  697  183  431  474  375  214  451
## [751]  493  370  229  474  512  259  197  421  420  302  328  541  546  285  280
## [766]  470  508  367  257  447  485  375  419  879  591  752  169  336  483  381
## [781]  223  423  424  376  222  461  498  319  338  803  718  694  220  461  501
## [796]  349  218  484  526  394  410 1102 1009  955  312  852  768  657  245  451
## [811]  438  405  193  458  351  366  374 1152  929 1036  299  809  724  858  402
## [826]  801  716  704  213  429  351  359  299  740  726  619  293  880  831  646

#Phân tổ dữ liệu cho biến thu nhập(income), phân tổ thành 4 tổ
Inc <- d$IN
cut(Inc,4)

##   [1] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
##   [7] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
##  [13] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
##  [19] (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
##  [25] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (1.93,19.5] (1.93,19.5]
##  [31] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
##  [37] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
##  [43] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
##  [49] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
##  [55] (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
##  [61] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
##  [67] (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
##  [73] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1]
##  [79] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
##  [85] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]  
##  [91] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
##  [97] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (54.5,72.1] (54.5,72.1]
## [103] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [109] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [115] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [121] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]  
## [127] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [133] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]  
## [139] (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [145] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (1.93,19.5] (1.93,19.5]
## [151] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [157] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [163] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [169] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
## [175] (37,54.5]   (37,54.5]   (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [181] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]  
## [187] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [193] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [199] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [205] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [211] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [217] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [223] (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [229] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [235] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [241] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [247] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [253] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [259] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [265] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [271] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [277] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
## [283] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [289] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [295] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [301] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [307] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [313] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]  
## [319] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [325] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
## [331] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [337] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [343] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [349] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [355] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [361] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5]
## [367] (1.93,19.5] (1.93,19.5] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [373] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [379] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [385] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [391] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [397] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]  
## [403] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [409] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [415] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [421] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [427] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [433] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
## [439] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [445] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]  
## [451] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [457] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
## [463] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [469] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [475] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [481] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [487] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [493] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [499] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [505] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]  
## [511] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [517] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]  
## [523] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [529] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
## [535] (1.93,19.5] (1.93,19.5] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [541] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]  
## [547] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [553] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [559] (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [565] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [571] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [577] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5]
## [583] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [589] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [595] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [601] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (37,54.5]   (37,54.5]  
## [607] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [613] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
## [619] (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [625] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [631] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [637] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1]
## [643] (54.5,72.1] (54.5,72.1] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [649] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1]
## [655] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [661] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]  
## [667] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [673] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]  
## [679] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [685] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [691] (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [697] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [703] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [709] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [715] (37,54.5]   (37,54.5]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## [721] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [727] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [733] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [739] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [745] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [751] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [757] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
## [763] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [769] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (1.93,19.5] (1.93,19.5]
## [775] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5]
## [781] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [787] (19.5,37]   (19.5,37]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [793] (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1] (37,54.5]   (37,54.5]  
## [799] (37,54.5]   (37,54.5]   (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [805] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]   (1.93,19.5] (1.93,19.5]
## [811] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]   (19.5,37]   (19.5,37]  
## [817] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [823] (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]   (37,54.5]  
## [829] (1.93,19.5] (1.93,19.5] (1.93,19.5] (1.93,19.5] (19.5,37]   (19.5,37]  
## [835] (19.5,37]   (19.5,37]   (54.5,72.1] (54.5,72.1] (54.5,72.1] (54.5,72.1]
## Levels: (1.93,19.5] (19.5,37] (37,54.5] (54.5,72.1]

table(cut(Inc,4))

## 
## (1.93,19.5]   (19.5,37]   (37,54.5] (54.5,72.1] 
##         204         284         196         156

#Phân tổ dữ liệu cho biến Choice
Choi <- d$C
table(Choi)

## Choi
##  no yes 
## 630 210