比较麻烦,在Stata10.1下测试能用:
clear
local i=1
set obs 100
gen R=.
save R, replace //先生成保存相关系数的文件
clear
set obs 150
gen x1=int(100*uniform()) //生成虚拟数据
gen x2=int(100*ln(1/x1))
forvalues i=1/100 { //循环100次重复抽样
preserve
bsample 100
corr x1 x2
matrix rc=r(C)
use R, clear
replace R=rc[1,2] in `i' //保存相关系数
save R, replace
restore
local i=`i'+1
}
use R, clear
histogram R, percent
呵呵,也许还能从bootstrap得出更简单的解决办法