# 曜日毎の月齢の割合 cat("\n[各曜日におけるそれぞれの月齢の試行回数の割合]\n") if (imgout) png(file=filename8, width=600, height=600) par(mfrow=c(3,3)) rs <- rep(0, length(dows)) for (w in 1:length(dows)) { rs[w] <- length(result[pet==dst.pet & CoP==ac & vdow==dows[w]]) } maxrs <- max(rs) r <- sqrt(0.5 * rs / maxrs) for (w in 1:length(dows)) { dat <- rep(0, times=length(moons)) for (m in 1:length(moons)) { dat[m] <- length(result[pet==dst.pet & CoP==ac & vmoon==moons[m] & vdow==dows[w]]) cat(dows[w],moons[m],dat[m],"\n") } pie(dat+1, moons, main=dows[w], radius=r[w]) } if (imgout) dev.off()