員工上班搭乘工具需時分析

場景

    • 公司抽樣調查員工搭乘大眾運輸工具與自用車,到達公司所需要的時間資料如附件檔。

問題

    • 求算平均值(Mean) ,中數(Median),眾數(Mode),四分位數(Quartile)

    • 求算變異數(Variance),標準差(Std. Deviation)

GNU R

# 員工上班搭乘交通工具需時分析 source("Statistics.R") ######################################################## # 3-66 # 已知 myData <- read.table(file="Descriptive-Statistics-B66.dat", header=TRUE) print(myData) print(summary(myData)) # 求算 myMean <- c(Mean2(myData$大眾運輸), Mean2(myData$自用車)) myMedian <- c(Median2(myData$大眾運輸), Median2(myData$自用車)) myVariance <- c(Variance(myData$大眾運輸, FALSE), Variance(myData$自用車, FALSE)) myStdDeviation <- c(StdDeviation2(myData$大眾運輸, FALSE), StdDeviation2(myData$自用車, FALSE)) print(sprintf("%s: Mean=%.4f Median=%.4f Variance=%.4f StdDeviation=%.4f", "大眾運輸", myMean[1], myMedian[1], myVariance[1], myStdDeviation[1])) print(sprintf("%s: Mean=%.4f Median=%.4f Variance=%.4f StdDeviation=%.4f", "自用車", myMean[2], myMedian[2], myVariance[2], myStdDeviation[2])) # 驗證 par(mfrow=c(4,2)) barplot(myData$大眾運輸 - myMean[1], main="大眾運輸二極長條圖", ylab="需時均值差(分鐘)") barplot(myData$自用車 - myMean[2], main="自用車二極長條圖", ylab="需時均值差(分鐘)") hist(myData$大眾運輸, main="大眾運輸直方圖", xlab="需時(分鐘)", ylab="人數級距", nclass=6) hist(myData$自用車, main="自用車直方圖", xlab="需時(分鐘)", ylab="人數級距", nclass=6) boxplot(myData$大眾運輸, main="大眾運輸盒狀圖", xlab="需時(分鐘)", horizontal=TRUE) boxplot(myData$自用車, main="自用車盒狀圖", xlab="需時(分鐘)", horizontal=TRUE) dotchart(myData$大眾運輸, main="大眾運輸點圖", xlab="需時(分鐘)") dotchart(myData$自用車, main="自用車點圖", xlab="需時(分鐘)")

解答

大眾運輸 自用車 1 28 29 2 29 31 3 32 33 4 37 32 5 33 34 6 25 30 7 29 31 8 32 32 9 41 35 10 34 33 大眾運輸 自用車 Min. :25.00 Min. :29 1st Qu.:29.00 1st Qu.:31 Median :32.00 Median :32 Mean :32.00 Mean :32 3rd Qu.:33.75 3rd Qu.:33 Max. :41.00 Max. :35 大眾運輸: Mean=32.0000 Median=29.0000 Variance=19.4000 StdDeviation=4.4045 自用車: Mean=32.0000 Median=32.0000 Variance=3.0000 StdDeviation=1.7321

1