(a)
# 讀進本題十個受試者的資料,觀察資料型態
dta <- read.table("ex3data.txt",header=T)
head(dta)
## Snt Sp Wrds New S01 S02 S03 S04 S05 S06 S07 S08 S09
## 1 1 1 13 1 3.429 2.795 4.161 3.071 3.625 3.161 3.232 7.161 1.536
## 2 2 2 16 3 6.482 5.411 4.491 5.063 9.295 5.643 8.357 4.313 2.946
## 3 3 3 9 2 1.714 2.339 3.018 2.464 6.045 2.455 4.920 3.366 1.375
## 4 4 4 9 2 3.679 3.714 2.866 2.732 4.205 6.241 3.723 6.330 1.152
## 5 5 5 10 3 4.000 2.902 2.991 2.670 3.884 3.223 3.143 6.143 2.759
## 6 6 6 18 4 6.973 8.018 6.625 7.571 8.795 13.188 11.170 6.071 7.964
## S10
## 1 4.063
## 2 6.652
## 3 2.179
## 4 3.661
## 5 3.330
## 6 7.866
# 判斷S01~S10為閱讀所需之時間,進行不同個體閱讀所需時間之平均與排序
rank(colMeans(dta[,5:14]))
## S01 S02 S03 S04 S05 S06 S07 S08 S09 S10
## 5 4 2 3 10 9 8 7 1 6
(b)
# 估計平均來說,十位學生平均會花多少時間閱讀一個字詞
colMeans(dta[,5:14])/mean(dta$Wrds)
## S01 S02 S03 S04 S05 S06 S07
## 0.3569259 0.3324568 0.3261728 0.3266049 0.4856790 0.4641852 0.4518395
## S08 S09 S10
## 0.4330864 0.2368765 0.3883457
# 而倘若不考量個別差異,再對其作整體平均
mean(colMeans(dta[,5:14])/mean(dta$Wrds))
## [1] 0.3802173