drawnorm x1-x20,n(1000) clear
forv i=1/20 {
replace x`i'=x`i'^2
}
egen y8=rsum(x1-x8)
egen y20=rsum(x1-x20)
keep y8 y20
drawnorm x1-x10,n(1000)
forv i=1/10 {
replace x`i'=x`i'^2
}
egen z2=rsum(x1-x2)
egen z6=rsum(x1-x6)
keep y* z2 z6
g f28=(z2/2)/(y8/8)
g f620=(z6/6)/(y20/20)
tw hist f28,bin(100) bfcolor(red) || hist f620,bin(30) blcolor(yellow) ||
function y=Fden(2,8,x), range(0 20) || function y=Fden(6,20,x), range(0 20)
这个是用stata构造F检验,我不是很明白你的意思,F检验一般是用来检验总体回归方程的显著性,一般输出结果的时候会直接给出的。本人菜鸟,希望能帮到你。