beta1的95%置信区间

3 Ansichten(30天)的

克里斯我2017年11月14日

0
Verknupfen

这位Frage祖茂堂Direkter链接

https://de.mathworks.com/matlabcentral/answers/367041 - 95 -间隔的beta1信心

Beantwortet: 约阿希姆Sogm我2017年11月18日

我应该模拟n线性回归和用我的估计贝塔和SE构建一个95%置信区间为了找到真正的测试的覆盖率。我试图建立一个for循环,使用我估计贝塔和SEs新循环产生许多置信区间。然而,数据不会改变在我的for循环和置信区间是一样的我跑的回归。

这是我的代码:

                         Nsim = 1000;
                        
                         Nsteps = 50;
                        
                         x = normrnd (0, 1, N, 1);
                        
                         b0 = 2;
                        
                         b1 = 2;
                        
                         Nsteps succesRate = 0 (1);
                        
                         nIndex = (1: Nsteps) * 5;
                        
                         为j = (1: Nsteps)
                        
                         n = nIndex (j);
                        
                         lm (i) = fitlm (y、x);
                        
                         上= 0 (1、Nsim);
                        
                         低= 0 (1、Nsim);
                        
                         betahat (i) = lm.Coefficients.Estimate (2);
                        
                         SEhom (i) = lm.Coefficients.SE (2);
                        
                         为我= (1:Nsim)
                        
                         上(我)= betahat(我)+ 1.96 * SEhom(我);
                        
                         低(我)= betahat(我)-1.96 * SEhom(我);
                        
                         结束
                        
                         印第安纳州=低(< = b1) & (b1 < =上);
                        
                         succesRate (j) =(印第安纳州)/ Nsim总和;
                        
                         结束