紙張密度與樹木品種分析

場景：

紙張由 4 種不同樹種所生產，有多種相應不同紙張密度製品，列表如后：

問題：

是否紙張密度與樹種相關？

分析：

# 指定精準度 options(digits=4) # 清除殘餘變數 rm(list=ls()) # 資料整理 Wood.Type <- c('Type-I', 'Type-II', 'Type-III', 'Type-IV') Wood <- data.frame(WoodType=c(rep(Wood.Type[1], 6), rep(Wood.Type[2], 6), rep(Wood.Type[3], 6), rep(Wood.Type[4], 6)), WoodDensity=c(c(7,8,15,11,9,10), c(12,17,13,18,19,15), c(14,18,19,17,16,18), c(19,25,22,23,18,20))) colnames(Wood) <- c('Type', 'Density') # 是否紙張密度與樹種有關有關 Wood.ANOVA = aov(Density~Type, data=Wood) summary(Wood.ANOVA) head(Wood, n=10) boxplot(Density~Type, ylab='Density', data=Wood)

# 是否紙張密度與樹種有關有關 Df Sum Sq Mean Sq F value Pr(>F) Type 3 383 127.6 19.6 3.6e-06 *** Residuals 20 130 6.5 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Type Density 1 Type-I 7 2 Type-I 8 3 Type-I 15 4 Type-I 11 5 Type-I 9 6 Type-I 10 7 Type-II 12 8 Type-II 17 9 Type-II 13 10 Type-II 18