|
#田中2008、12章
x <- seq(0,5,0.2);d1 <- c(1,3,8,8);d2 <- c(10,10,10,30)
xx <- c(0,5);yy <- c(0,1.5)
par(mfrow=c(2,2))
for (j in 2:4){
y <- df(x,df1=d1[j],df2=d2[j])
ttl <- paste("F(",d1[j],",",d2[j],")")
plot(x,y,type="l",xlim=xx,ylim=yy,main=ttl)
}
#case01
x1 <- rnorm(20,5,1.5);y1 <- rnorm(30,5,1.2)
var.test(x1,y1)
#case02
x2 <- rnorm(200,5,1.5);y2 <- rnorm(300,5,1.2)
var.test(x2,y2)
vx1 <- var(x1);vy1 <- var(y1)
FF1 <- vx1/vy1;FF1
pp <- 1 -pf(FF1,df1=19,df2=29);pp
|
