有一回參加科學素養新式評量的研習，講師帶大家用R作試題分析，那次的分析理論用試題反應理論，我自己也沒有很了解，所以就是學完了就沒再繼續用了。
後來查看看有哪些用R作試題分析的案例，看到有一個爸爸用R幫小孩的班級成績做了成績分析，他也把一些code詳細地寫出來，我只要照著改一改就可以用了。
正好第一次段考試我出題的，所以我就剛好可以拿來應用了。
準備好成績檔案之後，要先安裝R和RStudio，然後貼code進去就可以了，全程不用10秒鐘就可以完成了，這就是文字界面的威力啊。以後再試試看用R來寫SP Table，這樣我就不用為了作試題分析，還特地去開windows來用了啊。
code在此！
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
setwd("/media/pancala/3T/試題分析/0test")
scores <-read.table("bio.txt",header=TRUE,row.names="學號")
View(scores)
attach(scores)
summary(成績)
aver <-tapply(成績,班級,mean)
png(file="成績直方圖20格.png", width=1500, height=850, res=120)
hist(成績,breaks=20,main="成績分佈圖",ylab="人數")
dev.off()
png(file="成績直方圖10格.png", width=1500, height=850, res=120)
hist(成績,breaks=10,main="成績分佈圖",ylab="人數")
dev.off()
png(file="成績頻率圖.png", width=1500, height=850, res=120)
plot(density(成績),main="成績頻率圖",ylab="頻率", col='blue', lwd=2)
dev.off()
png(file="各班平均成績.png", width=1500, height=800, res=120)
boxplot(成績)
#boxplot.stats(成績)
boxplot(成績~班級,data=scores,main="各班生物成績盒狀圖", xlab="班級", ylab="成績")
lines(tapply(成績,班級,mean),col='blue',type='b')
dev.off()
png(file="全年級成績盒狀圖.png", width=500, height=800, res=120)
boxplot(成績,main="成績盒狀圖",xlab="全年級", ylab="成績")
dev.off()
shapiro.test(成績) #檢測是否常態分佈，p小於0.05，表示為非常態分佈
