作者etcsc (44)
看板NCCU03_STGAR
标题实验设计 4-6.
时间Mon Apr 12 22:43:40 2004
不嫌弃的参考一下
#4.6
#(a)
#Model=Ratio Control Algorithm Block=Time Period ACV=average cell voltage
Model Block=list(Model=paste("Ratio Control Algorithm",1:4),Block=LETTERS[1:6])
voltage.design=fac.design(c(4,6),Model Block)
ACV=c(4.93,4.85,4.83,4.89,4.86,4.91,4.88,4.77,4.75,4.79,4.90,4.94,4.95,4.85,
4.75,4.86,4.79,4.75,4.82,4.79,4.88,4.85,4.90,4.76)
ACV data=data.frame(voltage.design,ACV)
aov.ACV=aov(ACV~Block+Model,ACV data)
summary(aov.ACV)
#(b)
#SDOV=standard deviation of voltage
SDOV=c(0.05,0.04,0.09,0.03,0.04,0.02,0.13,0.04,0.05,0.03,0.11,0.05,0.06,0.05,
0.15,0.05,0.03,0.03,0.08,0.03,0.05,0.02,0.12,0.02)
SDOV data=data.frame(voltage.design,-log(SDOV))
aov.SDOV=aov(SDOV~Block+Model,SDOV data)
summary(aov.SDOV)
#(c)
qqnorm(resid(aov.SDOV))
qqline(resid(aov.SDOV))
plot(fitted(aov.SDOV),resid(aov.SDOV))
abline(h=0)
plot(Ratio Control Algorithm,resid(aov.SDOV))
abline(h=0)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 202.178.204.62