PART3 HCAPTER7のシミュレーション。
130人の体温をシミュレーション作成してプロットしてみる。
tm<- rnorm(130,36.5,0.8) plot(rep(1,length(tm)),tm,cex=0.5,col="red",ylim=c(30,40)) par(new=TRUE) boxplot(tm,ylim=c(30,40)) par(mfrow=c(2,2)) hist(tm,main="automatically") hist(tm,breaks=max(tm),main="breaks=max(tm)") hist(tm,breaks=1,main="breaks=1") hist(tm,breaks=5,main="breaks=5") par(mfrow=c(1,2)) tm_sort<- sort(tm) plot(tm_sort,(1:length(tm_sort))/length(tm_sort),ylab="Frequency",type="o",cex=0.5) plot(tm_sort,(1:length(tm_sort))/length(tm_sort),ylab="Frequency",type="o",cex=0.5,log="y")