Hw5

dta <- read.table("/Users/Tjlee/Desktop/weather/readingtimes.txt", header = T)
dta

##   Snt Sp Wrds New   S01   S02   S03   S04   S05    S06    S07   S08   S09   S10
## 1   1  1   13   1 3.429 2.795 4.161 3.071 3.625  3.161  3.232 7.161 1.536 4.063
## 2   2  2   16   3 6.482 5.411 4.491 5.063 9.295  5.643  8.357 4.313 2.946 6.652
## 3   3  3    9   2 1.714 2.339 3.018 2.464 6.045  2.455  4.920 3.366 1.375 2.179
## 4   4  4    9   2 3.679 3.714 2.866 2.732 4.205  6.241  3.723 6.330 1.152 3.661
## 5   5  5   10   3 4.000 2.902 2.991 2.670 3.884  3.223  3.143 6.143 2.759 3.330
## 6   6  6   18   4 6.973 8.018 6.625 7.571 8.795 13.188 11.170 6.071 7.964 7.866
## 7   7  7    6   1 2.634 1.750 2.268 2.884 3.491  3.688  2.054 1.696 1.455 3.705

##擷取S1~S10
data <- dta[,5:14]
data

##     S01   S02   S03   S04   S05    S06    S07   S08   S09   S10
## 1 3.429 2.795 4.161 3.071 3.625  3.161  3.232 7.161 1.536 4.063
## 2 6.482 5.411 4.491 5.063 9.295  5.643  8.357 4.313 2.946 6.652
## 3 1.714 2.339 3.018 2.464 6.045  2.455  4.920 3.366 1.375 2.179
## 4 3.679 3.714 2.866 2.732 4.205  6.241  3.723 6.330 1.152 3.661
## 5 4.000 2.902 2.991 2.670 3.884  3.223  3.143 6.143 2.759 3.330
## 6 6.973 8.018 6.625 7.571 8.795 13.188 11.170 6.071 7.964 7.866
## 7 2.634 1.750 2.268 2.884 3.491  3.688  2.054 1.696 1.455 3.705

#算平均
dtanew <- apply(data, 2, mean)
dtanew

##      S01      S02      S03      S04      S05      S06      S07      S08 
## 4.130143 3.847000 3.774286 3.779286 5.620000 5.371286 5.228429 5.011429 
##      S09      S10 
## 2.741000 4.493714

#排名
order(dtanew)

##  [1]  9  3  4  2  1 10  8  7  6  5

#平均估計閱讀一個字需要多長時間。
dtaword <- data/dta[,3]
dtaw <- apply(dtaword, 2, mean)
dtaw

##       S01       S02       S03       S04       S05       S06       S07       S08 
## 0.3563579 0.3218649 0.3285283 0.3283256 0.4939310 0.4616710 0.4297786 0.4474266 
##       S09       S10 
## 0.2205573 0.3949539