#################################################
#豊田秀樹(2006) 購買心理を読み解く統計学 ---実例で見る心理・調査データ 東京図書
#5-4 一般化可能性理論 p.215--222
#################################################
source("http://blue.zero.jp/yokumura/R/testtheory/ctt.txt")
dat <- data.frame(
            value = c(5,5,3,5,5, 5,5,5,5,5, 4,4,4,4,5,
                      1,3,2,5,1, 2,3,3,4,2, 2,3,3,4,2,
                      5,2,2,4,1, 5,4,3,5,3, 4,3,3,4,3, 
                      5,5,3,5,4, 5,4,4,5,4, 4,5,3,3,3, 
                      1,1,2,3,1, 3,2,3,3,1, 1,1,1,2,1,
                      
                      3,3,2,4,2, 3,3,3,4,2, 2,2,2,3,1, 
                      3,5,1,3,1, 1,5,1,4,1, 2,4,3,4,1,
                      4,3,3,5,3, 5,3,4,5,3, 4,3,4,5,2,
                      4,3,3,5,4, 4,3,4,5,2, 3,3,3,4,3,
                      5,2,2,4,3, 5,2,3,4,3, 4,1,2,2,2),
            judge = factor(rep(paste("judge", 1:5, sep=""), times=30)),
            subj  = factor(rep(paste("subj", 01:10, sep=""), each=15), levels=paste("subj", 01:10, sep="")),
            item =  factor(rep(paste("item", 1:3, sep=""), each=5, times=10))
        )
output <- DG.2facets(data=dat, sim.judge = seq(1, 20, by=1),sim.item = seq(1, 20, by=1))
graphics.off()

##図5-4-2,  図5-4-3
output$sim1
par(mfrow=c(2, 1))
plot(output$sim1[, "item3"], xlab="面接官の人数", ylab="一般化可能性係数", type="b", ylim=c(.5, 1))
plot(output$sim1["judge5", ], xlab="評価観点の数", ylab="一般化可能性係数", type="b", ylim=c(.5, 1))