指定高杠杆点的ID

时间:2015-08-20 04:12:54

标签: r

看一下模型

     > summary(left_int3)
     Call:
     lm(formula = Left ~ Cursor + PostCursor + CTLE + I(Cursor^2) + 
                  I(PostCursor^2), data = QPI)

      Residuals:
      Min       1Q   Median       3Q      Max 
     -3585033  -845980   -58093   718849  4289610 

      Coefficients:
                        Estimate Std. Error t value Pr(>|t|)    
      (Intercept)      1.455e+07  6.651e+05  21.880  < 2e-16 ***
      Cursor           2.299e-06  1.302e-06   1.766  0.07754 .  
      PostCursor       1.150e-06  2.147e-06   0.536  0.59231    
      CTLE            -4.772e+00  4.548e-01 -10.493  < 2e-16 ***
      I(Cursor^2)     -2.162e-18  6.854e-19  -3.154  0.00163 ** 
      I(PostCursor^2) -2.775e-19  5.977e-19  -0.464  0.64253    
      ---
      Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

      Residual standard error: 1175000 on 2794 degrees of freedom
      Multiple R-squared:  0.3942,  Adjusted R-squared:  0.3932 
      F-statistic: 363.7 on 5 and 2794 DF,  p-value: < 2.2e-16

我想知道有任何关于高帽值的观察,我使用下面的命令,我用0.04作为标准

     >hatvalues(left_int3)>0.04
     1     2     3     4     5     6     7     8        
     FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE         

我用命令获取了True和False的数字:

      > hat<-hatvalues(left_int3)>0.04
      > summary(hat)
      Mode      FALSE    TRUE    NA's 
      logical    2780      20       0 

所以我想知道是否有任何命令知道这些观察的ID号是什么&#34; TRUE&#34;。

0 个答案:

没有答案