t<-c(165,158,170,166)

w<-c(56,48,57,60)

mean(t)     #平均數
## [1] 164.75
median(w)    #中位數
## [1] 56.5
as.numeric(names(table(t)))[which.max(table(t))] #眾數
## [1] 158
sd(t)   #標準差
## [1] 4.99166
var(w)  #變異數
## [1] 26.25
plot(t,w,
     pch = 17,
     col = "pink",
     main ="班上同學體重和身高",
     xlab= "體重",
     ylab= "身高")

hist(t,
     col= "pink",
     main ="班上同學體重和身高",
     xlab = "體重",
     ylab = "次數")