#屏東縣精進計畫資料分析 #2022/11/15 屏東大學教育學系陳新豐 #chensf@nptu.edu.tw library(readr) #讀入資料 pdata <- read_csv("D:/data/p20221114_2.csv") print(pdata) head(pdata) tail(pdata) #描述性統計 #次數分配 tp011 <- table(pdata$'P01') tp012 <- prop.table(tp011) tp021 <- table(pdata$'P02') tp022 <- prop.table(tp021) tp031 <- table(pdata$'P03') tp032 <- prop.table(tp031) tp041 <- table(pdata$'P04') tp042 <- prop.table(tp041) tp051 <- table(pdata$'P05') tp052 <- prop.table(tp051) print(tp011) print(tp012) print(tp021) print(tp022) print(tp031) print(tp032) print(tp041) print(tp042) print(tp051) print(tp052) library(moments) my_stats <- function(x){ funs <- c(mean, sd, skewness, kurtosis) sapply(funs, function(f) f(x, na.rm=TRUE)) } pdata_1 <- pdata[,c('A00','B01','B02','B03','B04','B00','C01','C02','C03','C00','D01','D02','D03','D00','T00')] pdata_2 <- apply(pdata_1, 2, my_stats) #print(pdata_2) rownames(pdata_2) <- c("平均數","標準差","偏態","峰度") pdata_3 <- as.data.frame(t(pdata_2)) print(pdata_3) #假設考驗 #不同性別P01其分數是否有所不同(2) #1.男2.女 #性別轉換為類別變項 pdata$P01 <- factor(pdata$P01) #未來建議加入同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P01, data=pdata, center=mean) library(Rmisc) #A00 sp01 <- summarySE(data=pdata,groupvars="P01",measurevar="A00") print(sp01) mp01 <- t.test(A00 ~ P01, data=pdata, var.equal=TRUE) print(mp01) #B00 sp01 <- summarySE(data=pdata,groupvars="P01",measurevar="B00") print(sp01) mp01 <- t.test(B00 ~ P01, data=pdata, var.equal=TRUE) print(mp01) #C00 sp01 <- summarySE(data=pdata,groupvars="P01",measurevar="C00") print(sp01) mp01 <- t.test(C00 ~ P01, data=pdata, var.equal=TRUE) print(mp01) #D00 sp01 <- summarySE(data=pdata,groupvars="P01",measurevar="D00") print(sp01) mp01 <- t.test(D00 ~ P01, data=pdata, var.equal=TRUE) print(mp01) #T00 sp01 <- summarySE(data=pdata,groupvars="P01",measurevar="T00") print(sp01) mp01 <- t.test(T00 ~ P01, data=pdata, var.equal=TRUE) print(mp01) #不同最高學歷P02其分數是否有所不同(3) #1.學士2.碩士3.博士 #P02轉換為類別變項 pdata$P02 <- factor(pdata$P02) library(Rmisc) #P02A00 sp01 <- summarySE(data=pdata,groupvars="P02",measurevar="A00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P02, data=pdata, center=mean) #ANOVA mp01 <- aov(A00 ~ P02, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P02B00 sp01 <- summarySE(data=pdata,groupvars="P02",measurevar="B00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(B00 ~ P02, data=pdata, center=mean) #ANOVA mp02 <- aov(B00 ~ P02, data=pdata) mp02_a <- anova(mp02) print(mp02_a) #若有顯著需要進行事後比較 library(DescTools) mp02_s <- PostHocTest(mp02, method='scheffe') print(mp02_s) #P02C00 sp01 <- summarySE(data=pdata,groupvars="P02",measurevar="C00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(C00 ~ P02, data=pdata, center=mean) #ANOVA mp02 <- aov(C00 ~ P02, data=pdata) mp02_a <- anova(mp02) print(mp02_a) #若有顯著需要進行事後比較 library(DescTools) mp02_s <- PostHocTest(mp02, method='scheffe') print(mp02_s) #P02D00 sp01 <- summarySE(data=pdata,groupvars="P02",measurevar="D00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(D00 ~ P02, data=pdata, center=mean) #ANOVA mp02 <- aov(D00 ~ P02, data=pdata) mp02_a <- anova(mp02) print(mp02_a) #若有顯著需要進行事後比較 library(DescTools) mp02_s <- PostHocTest(mp02, method='scheffe') print(mp02_s) #P02T00 sp01 <- summarySE(data=pdata,groupvars="P02",measurevar="T00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(T00 ~ P02, data=pdata, center=mean) #ANOVA mp02 <- aov(T00 ~ P02, data=pdata) mp02_a <- anova(mp02) print(mp02_a) #若有顯著需要進行事後比較 library(DescTools) mp02_s <- PostHocTest(mp02, method='scheffe') print(mp02_s) #不同教學年資P03其分數是否有所不同(8) pdata$P03 <- factor(pdata$P03) library(Rmisc) #P03A00 sp01 <- summarySE(data=pdata,groupvars="P03",measurevar="A00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P03, data=pdata, center=mean) #ANOVA mp01 <- aov(A00 ~ P03, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P03B00 sp01 <- summarySE(data=pdata,groupvars="P03",measurevar="B00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P03, data=pdata, center=mean) #ANOVA mp01 <- aov(B00 ~ P03, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P03C00 sp01 <- summarySE(data=pdata,groupvars="P03",measurevar="C00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P03, data=pdata, center=mean) #ANOVA mp01 <- aov(C00 ~ P03, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P03D00 sp01 <- summarySE(data=pdata,groupvars="P03",measurevar="D00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P03, data=pdata, center=mean) #ANOVA mp01 <- aov(D00 ~ P03, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P03T00 sp01 <- summarySE(data=pdata,groupvars="P03",measurevar="T00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P03, data=pdata, center=mean) #ANOVA mp01 <- aov(T00 ~ P03, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #不同擔任職務P04其分數是否有所不同(4) #1.級任導師2.科任老師3.教師兼組長4.教師兼主任5.校長 pdata$P04 <- factor(pdata$P04) library(Rmisc) #P04A00 sp01 <- summarySE(data=pdata,groupvars="P04",measurevar="A00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P04, data=pdata, center=mean) #ANOVA mp01 <- aov(A00 ~ P04, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P04B00 sp01 <- summarySE(data=pdata,groupvars="P04",measurevar="B00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P04, data=pdata, center=mean) #ANOVA mp01 <- aov(B00 ~ P04, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P04C00 sp01 <- summarySE(data=pdata,groupvars="P04",measurevar="C00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P04, data=pdata, center=mean) #ANOVA mp01 <- aov(C00 ~ P04, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P04D00 sp01 <- summarySE(data=pdata,groupvars="P04",measurevar="D00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P04, data=pdata, center=mean) #ANOVA mp01 <- aov(D00 ~ P04, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P04T00 sp01 <- summarySE(data=pdata,groupvars="P04",measurevar="T00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P04, data=pdata, center=mean) #ANOVA mp01 <- aov(T00 ~ P04, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #不同學校規模P05其分數是否有所不同(4) #1.6班以下2.7-12班3.13-24班4.24班以上 pdata$P05 <- factor(pdata$P05) library(Rmisc) #P05A00 sp01 <- summarySE(data=pdata,groupvars="P05",measurevar="A00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P05, data=pdata, center=mean) #ANOVA mp01 <- aov(A00 ~ P05, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P05B00 sp01 <- summarySE(data=pdata,groupvars="P05",measurevar="B00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P05, data=pdata, center=mean) #ANOVA mp01 <- aov(B00 ~ P05, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P05C00 sp01 <- summarySE(data=pdata,groupvars="P05",measurevar="C00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P05, data=pdata, center=mean) #ANOVA mp01 <- aov(C00 ~ P05, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P05D00 sp01 <- summarySE(data=pdata,groupvars="P05",measurevar="D00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P05, data=pdata, center=mean) #ANOVA mp01 <- aov(D00 ~ P05, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s) #P05T00 sp01 <- summarySE(data=pdata,groupvars="P05",measurevar="T00") print(sp01) #同質性檢定 #library(DescTools) #LeveneTest(A00 ~ P05, data=pdata, center=mean) #ANOVA mp01 <- aov(T00 ~ P05, data=pdata) mp01_a <- anova(mp01) print(mp01_a) #若有顯著需要進行事後比較 library(DescTools) mp01_s <- PostHocTest(mp01, method='scheffe') print(mp01_s)