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