项目场景:
stata
问题描述
program onesample, rclass //定义程序,存储r()形式
drop _all //(删除内存已有数据 )
set obs 30//随机抽取样本容量30
gen x =runiform()//r指的是random uniformdistribution指的是均匀分布
sum x //因为是默认0-1的均匀分布,所以求和就是求均值 错,sum在stata里面是统计,求方差求均值等,
//在书上17页
return scalar mean_sample =r(mean)//scalar是1*1维标量,就是个数,mean是计算sum时存储在r里面
*的均值
end
set more off //指定stata输出结果连续翻页
simulate xbar = r(mean_sample), reps(10000) seed(101) nodots: onesample
hist xbar, normal//选择项是画正态分布做对比
/*将onesample运行10000次,xbar为
10000×1的向量*/