# 載入套件與讀取資料
library(haven)
pisa_data <- read_sav("PISA_tawian2022_trimmed_lab.sav")
escs <- pisa_data$ESCS

# 計算各項統計值
escs_mean <- round(mean(escs, na.rm = TRUE), 2)
escs_sd <- round(sd(escs, na.rm = TRUE), 2)
escs_max <- round(max(escs, na.rm = TRUE), 2)
escs_min <- round(min(escs, na.rm = TRUE), 2)
escs_range <- round(escs_max - escs_min, 2)
escs_median <- round(median(escs, na.rm = TRUE), 2)
escs_iqr <- round(IQR(escs, na.rm = TRUE), 2)

# 列印統計結果
cat("均值:", escs_mean, "\n")
均值: -0.27 
cat("標準差:", escs_sd, "\n")
標準差: 0.91 
cat("最大值:", escs_max, "\n")
最大值: 4.21 
cat("最小值:", escs_min, "\n")
最小值: -3.84 
cat("全距:", escs_range, "\n")
全距: 8.05 
cat("中位數:", escs_median, "\n")
中位數: -0.24 
cat("四分位距:", escs_iqr, "\n")
四分位距: 1.38