GIẢI THÍCH Dữ liệu Medicaid1986 của package AER.visits: Số bác sĩ đến thăm. exposure: Thời gian quan sát đối với dịch vụ chăm sóc ngoại trú. children: số đứ trẻ có trong gia đình. age: số tuổi. income Annual household: thu nhập thường niên trong gia đình. health1: Thành phần chính đầu tiên của ba biến tình trạng sức khỏe: chức năng hạn chế, tình trạng cấp tính và tình trạng mãn tính. health2: Thành phần chính thứ 2 của ba biến tình trạng sức khỏe: chức năng hạn chế, tình trạng cấp tính và tình trạng mãn tính. access: Sự sẵn có của các dịch vụ y tế (0 = khả năng tiếp cận thấp, 1 = khả năng tiếp cận cao). married: có kết hôn hay không. gender: giới tính. ethnicity: dân tộc. school: số năm học ở trường. enroll; có phải cá nhân đăng ký nhập viện. program: chương trình chăm sóc
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#Lấy dữ liệu "Medicaid1986" từ package "AER"
data("Medicaid1986")
#Tạo dữ liệu bản sao "a" từ dữ liệu gốc "Medicaid1986"
a <- Medicaid1986
#Xem cấu trúc dữ liệu của dữ liệu "a"
str(a)## 'data.frame': 996 obs. of 14 variables:
## $ visits : int 0 1 0 0 11 3 0 6 1 0 ...
## $ exposure : int 100 90 106 114 115 102 92 92 117 101 ...
## $ children : int 1 3 4 2 1 1 2 1 1 1 ...
## $ age : int 24 19 17 29 26 22 24 21 21 24 ...
## $ income : num 14.5 6 8.38 6 8.5 ...
## $ health1 : num 0.495 0.52 -1.227 -1.524 0.173 ...
## $ health2 : num -0.854 -0.969 0.317 0.457 -0.599 0.062 0.202 -0.981 0.317 -0.562 ...
## $ access : num 0.5 0.17 0.42 0.33 0.67 0.25 0.5 0.67 0.25 0.67 ...
## $ married : Factor w/ 2 levels "no","yes": 1 1 1 1 1 1 1 2 1 2 ...
## $ gender : Factor w/ 2 levels "male","female": 2 2 2 2 2 2 2 2 2 2 ...
## ..- attr(*, "contrasts")= num [1:2, 1] 1 0
## .. ..- attr(*, "dimnames")=List of 2
## .. .. ..$ : chr [1:2] "male" "female"
## .. .. ..$ : chr "male"
## $ ethnicity: Factor w/ 2 levels "cauc","other": 1 1 1 1 1 2 1 1 1 1 ...
## ..- attr(*, "contrasts")= num [1:2, 1] 1 0
## .. ..- attr(*, "dimnames")=List of 2
## .. .. ..$ : chr [1:2] "caucasian" "other"
## .. .. ..$ : chr "caucasian"
## $ school : int 13 11 12 12 16 12 11 11 12 15 ...
## $ enroll : Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 2 2 2 2 ...
## $ program : Factor w/ 2 levels "afdc","ssi": 1 1 1 1 1 1 1 1 1 1 ...#Hiển thị 3 dòng đâu của dữ liệu "a"
head(a,3)## visits exposure children age income health1 health2 access married gender
## 1 0 100 1 24 14.500 0.495 -0.854 0.50 no female
## 2 1 90 3 19 6.000 0.520 -0.969 0.17 no female
## 3 0 106 4 17 8.377 -1.227 0.317 0.42 no female
## ethnicity school enroll program
## 1 cauc 13 yes afdc
## 2 cauc 11 yes afdc
## 3 cauc 12 yes afdc#Gán tên viết tắt cho các biến thuộc dữ liệu "g" cho tiện thao tác
names(a) <- c("vi","ex","ch","ag","in","he1","he2","ac","ma","ge","et","sc","en","pr")
#Chỉ lấy biến "vi" của dữ liệu "a"
visi <- a$vi
#Chỉ lấy biến "ag" của dữ liệu "a"
ag <- a$ag
#Xem dữ liệu thứ 6 trong dữ liệu "visi"
visi[6]## [1] 3#Liệt kê các biến có trong dữ liệu "a"
list=ls(a)#Tạo các dữ liệu mới bằng các phép toán trên dữ liệu "visi" và "ag"
phepcong <- visi + ag
pheptru <- visi - ag
phepnhan <- visi * ag
phepchia <- visi / ag
#Chỉ lấy những giá trị trên 2 của biến "visi"
visi2 <- visi[visi>2]
#Chỉ lấy dữ liệu có giá trị từ 2 đến 4 của biến "visi"
visi2den4 <- visi[visi>1&visi<5]
#Thêm vào 2 biến ""sqrtex" và "sqrtch" cho dữ liệu "a"
a$sqrtex <- sqrt(a$ex)
a$sqrtch <- sqrt(a$ch)#Tạo bảng tần số cho dữ liệu định tính "school" của dữ liệu "a"
school <- a$sc
table(school)## school
## 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
## 75 13 24 32 44 27 52 22 101 53 77 94 213 64 50 24 16 4 11#Tạo bảng tần số cho biến định lượng "in" của dữ liệu "a"
inc <- a$"in"
table(cut(inc,3))##
## (0.483,6.17] (6.17,11.8] (11.8,17.5]
## 465 425 106