n_list<-c(2,5,7,9) # Ç¥º»¼ö(n)
df_list<-n_list-1 # ÀÚÀ¯µµ

curve(dchisq(x, 1, ncp=0), col="blue", xlim=c(0, 16), ylim=c(0, 0.8), xlab="chisq", ylab="f(chisq)")
curve(dchisq(x, 4, ncp=0), add=T, col="red", xlim=c(0, 16), ylim=c(0, 0.8), xlab="chisq", ylab="f(chisq)")
curve(dchisq(x, 6, ncp=0), add=T, col="green", xlim=c(0, 16), ylim=c(0, 0.8), xlab="chisq", ylab="f(chisq)")
curve(dchisq(x, 8, ncp=0), add=T, col="black", xlim=c(0, 16), ylim=c(0, 0.8), xlab="chisq", ylab="f(chisq)")

par(mfrow=c(2,2))

for (i in 1:length(df_list)) {
  
  curve(dchisq(x, df_list[i], ncp=0), add=F, xlim=c(0, 16), ylim=c(0, 0.8), xlab="chisq", ylab="f(chisq)", main=paste("df=", df_list[i]))
  
}