data <- read.csv("C:/Users/wenzao/Downloads/table1_1.csv")




names(data)  # 確保有 "Japanese" 欄位
## [1] "Name"       "Statistic"  "Math"       "Japanese"   "Management"
## [6] "Accounting"
japanese <- data$Japanese


mean_jp <- mean(japanese)

# 中位數
median_jp <- median(japanese)


get_mode <- function(x) {
  uniq_vals <- unique(x)
  uniq_vals[which.max(tabulate(match(x, uniq_vals)))]
}
mode_jp <- get_mode(japanese)


sd_jp <- sd(japanese)


var_jp <- var(japanese)


q1_jp <- quantile(japanese, 0.25)


q3_jp <- quantile(japanese, 0.75)


result <- data.frame(
  統計量 = c("平均數", "中位數", "眾數", "標準差", "變異數", "第一四分位數", "第三四分位數"),
  數值 = c(mean_jp, median_jp, mode_jp, sd_jp, var_jp, q1_jp, q3_jp)
)

print(result)
##         統計量      數值
## 1       平均數  67.90000
## 2       中位數  62.00000
## 3         眾數  84.00000
## 4       標準差  16.25115
## 5       變異數 264.10000
## 6 第一四分位數  54.50000
## 7 第三四分位數  82.75000