かっこよくなった

わーい。乱かい法とか知らないで、二元分散分析やったったけど、気にしないぞー。

```a1 <- c(12,12.7,13.3,11.3,14.7,15.0)
a2 <- c(12.0,14.3,14.7,13.7,16.3,16.0)
a3 <- c(15.3,16.0,18.3,15.7,20.0,20.3)
a4 <- c(11.3,13.7,14.3,11.3,14.7,15.7)

A <- factor(c(rep("a1",6),rep("a2",6),rep("a3",6),rep("a4",6)))
B <- rep(c("b1","b2","b3"),8)
y <- c(a1,a2,a3,a4)
bunsan <- data.frame(A,B,y)
summary(aov(y~A+B,data=bunsan))

ab <- c()
for(i in 1:3){
ab[4*(i-1)+1] <- (a1[i]+a1[i+3])/2
ab[4*(i-1)+2] <- (a2[i]+a2[i+3])/2
ab[4*(i-1)+3] <- (a3[i]+a3[i+3])/2
ab[4*(i-1)+4] <- (a4[i]+a4[i+3])/2
}

mat <- matrix(ab,3,4,byrow=T)
colnames <- c("A1","A2","A3","A4")
colnames(mat) <- colnames
rownames <- c("B1","B2","B3")
rownames(mat) <- rownames

par(mar=c(6,6,6,3))
matplot(t(mat),type="b",pch=21,xlab="接着温度",ylab="プレス圧",axes=F,col=c("#1E5692", "#3E9A3B", "#F9A015"),lwd=5,cex=2,lty=1)

my.col=c("darkorchid1","cyan3","aquamarine3")

par(xpd=TRUE)
title(main="データのグラフ化")
rownames
legend(1, max(mat), legend=rownames,
col=my.col,
lty=c(1:3), pch=c(3,2,1)
)

t(cbind(a1,a2,a3,a4))

matpoints(t(cbind(a1,a2,a3,a4)),pch=c(3,2,1),col=my.col)
for(i in 1:3){
segments(1,a1[i],1,a1[i+3],lty=2)
segments(2,a2[i],2,a2[i+3],lty=2)
segments(3,a3[i],3,a3[i+3],lty=2)
segments(4,a4[i],4,a4[i+3],lty=2)
}

axis(1,1:4,colnames,)
axis(2)
```