install.packages("https://cran.rstudio.com/bin/windows/contrib/4.1/faraway_1.0.7.zip", repos = NULL)
## 將程式套件安載入 'C:/Users/user/Documents/R/win-library/4.1'
## (因為 'lib' 沒有被指定)
## package 'faraway' successfully unpacked and MD5 sums checked
library(faraway)
## Warning: 套件 'faraway' 是用 R 版本 4.1.2 來建造的
data(hsb)
summary(hsb)
## id gender race ses schtyp
## Min. : 1.00 female:109 african-amer: 20 high :58 private: 32
## 1st Qu.: 50.75 male : 91 asian : 11 low :47 public :168
## Median :100.50 hispanic : 24 middle:95
## Mean :100.50 white :145
## 3rd Qu.:150.25
## Max. :200.00
## prog read write math science
## academic:105 Min. :28.00 Min. :31.00 Min. :33.00 Min. :26.00
## general : 45 1st Qu.:44.00 1st Qu.:45.75 1st Qu.:45.00 1st Qu.:44.00
## vocation: 50 Median :50.00 Median :54.00 Median :52.00 Median :53.00
## Mean :52.23 Mean :52.77 Mean :52.65 Mean :51.85
## 3rd Qu.:60.00 3rd Qu.:60.00 3rd Qu.:59.00 3rd Qu.:58.00
## Max. :76.00 Max. :67.00 Max. :75.00 Max. :74.00
## socst
## Min. :26.00
## 1st Qu.:46.00
## Median :52.00
## Mean :52.41
## 3rd Qu.:61.00
## Max. :71.00
names(hsb)
## [1] "id" "gender" "race" "ses" "schtyp" "prog" "read"
## [8] "write" "math" "science" "socst"
hsb[,1:5] <- lapply(hsb[, 1:5], factor)
aggregate(cbind(read, write, math, science, socst) ~ gender, data = hsb, mean)
## gender read write math science socst
## 1 female 51.73394 54.99083 52.39450 50.69725 52.91743
## 2 male 52.82418 50.12088 52.94505 53.23077 51.79121
aggregate(cbind(read, write, math, science, socst) ~ race, data = hsb, mean)
## race read write math science socst
## 1 african-amer 46.80000 48.20000 46.75000 42.80000 49.45000
## 2 asian 51.90909 58.00000 57.27273 51.45455 51.00000
## 3 hispanic 46.66667 46.45833 47.41667 45.37500 47.79167
## 4 white 53.92414 54.05517 53.97241 54.20000 53.68276
aggregate(cbind(read, write, math, science, socst) ~ gender, data = hsb, sd)
## gender read write math science socst
## 1 female 10.05783 8.133715 9.151015 9.038503 10.23441
## 2 male 10.50671 10.305161 9.664784 10.732171 11.33384
aggregate(cbind(read, write, math, science, socst) ~ race, data = hsb, sd)
## race read write math science socst
## 1 african-amer 7.120024 9.322299 6.487843 9.445690 10.850540
## 2 asian 7.660999 7.899367 10.120187 9.490665 9.746794
## 3 hispanic 10.239169 8.272422 6.983936 8.218815 9.250049
## 4 white 10.276783 9.172558 9.383011 9.094870 10.813253
install.packages("moments", repos='https://cran.rstudio.com/')
## 將程式套件安載入 'C:/Users/user/Documents/R/win-library/4.1'
## (因為 'lib' 沒有被指定)
## package 'moments' successfully unpacked and MD5 sums checked
##
## The downloaded binary packages are in
## C:\Users\user\AppData\Local\Temp\RtmpcRVHig\downloaded_packages
library(moments)
## Warning: 套件 'moments' 是用 R 版本 4.1.1 來建造的
aggregate(cbind(read, write, math, science, socst) ~ gender+race, data = hsb, skewness)
## gender race read write math science socst
## 1 female african-amer 1.0908453 0.2163986 1.48488184 0.73940428 -0.30003839
## 2 male african-amer -0.5014614 0.2805474 0.21617314 -0.98118254 -0.64522406
## 3 female asian -0.0720213 -0.8936840 -0.18751759 -0.15405549 0.13198759
## 4 male asian 0.7071068 -0.6993411 -0.45563859 -0.59517006 0.38180177
## 5 female hispanic 0.8336115 0.4399599 -0.07735766 0.56561271 0.43962390
## 6 male hispanic 0.3505098 0.6740362 0.38302113 0.10751219 0.01372901
## 7 female white 0.2652688 -0.8495136 0.05270082 0.01337789 -0.49165481
## 8 male white -0.1477733 -0.3818725 0.14096486 -0.50797997 -0.58469114
aggregate(cbind(read, write, math, science, socst) ~ gender+race, data = hsb, kurtosis)
## gender race read write math science socst
## 1 female african-amer 5.293553 2.441087 4.698841 2.597492 2.554017
## 2 male african-amer 2.571704 1.213593 1.944268 2.552533 2.315460
## 3 female asian 1.788177 2.911449 1.718401 1.825103 2.708239
## 4 male asian 1.500000 1.500000 1.500000 1.500000 1.500000
## 5 female hispanic 3.454724 1.920618 2.362541 2.095000 1.880657
## 6 male hispanic 2.303890 3.191498 1.951063 2.659911 1.923078
## 7 female white 2.254607 3.211956 2.336165 2.573160 2.609000
## 8 male white 2.193766 2.045456 2.284359 2.425587 2.560890