#使用法: 將此檔案 (命名為 example7.1.txt) 置於 R 的工作目錄中, 執行指令 source('example7.1.txt')

set.seed(1129)		#選定種子
d=c()			#宣告物件
for(i in 1:4000) d[i]=sum(sample(1:6, 3, replace=TRUE))		#由 1 至 6 中重覆取樣 3 次, 
								#求其和, 並放入 d[i] 中, 共 4000 次
hist(d/3, breaks=seq(2/3, 6, by=1/3), col=gray(.9), 
	main='Number of rolls = 4000', xlab='mean of 3 dice')	#繪直方圖

ybar=mean(d/3)
sigma2=var(d/3)
sigma=sd(d/3)

for(k in 3:18){
	cat('Prob(ybar =', k, '/3)=', sum(d==k)/4000, '\n')
}