0

したがって、 と の SAT スコアである 2 つの変数がVerbal(SATV)ありQuantitative(SATQ)ます。500 行あります。に7つありNA'sますSATQ。私の目標は、lm()and gvlma()with SATVand SATQasを実行することIVsです。しかし、変数の長さが異なるため
、R がコードを実行しないというエラーが表示されます。変数が同じ長さになるように再コーディングするにはどうすればよいですか。 非正規データと仮定違反を無視します。（また、Rで何をしているのかわかりません。アドバイスを提供できる場合は、Rをまったく理解していない人と話しているふりをしてください）NAsSATQNA's

> summary(SATQ)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
  200.0   525.0   610.0   604.5   700.0   800.0       7 

> SATQ2<-na.omit(SATQ)

> summary(SATQ2)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  200.0   525.0   610.0   604.5   700.0   800.0 

> summary(SATV)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  200.0   537.5   600.0   604.4   690.0   800.0

> summary(ms)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  765.6  1844.0  2133.0  2093.0  2395.0  2877.0 

> #ms= monthly salary
> m1 = lm(ms~SATV+SATQ2)
Error in model.frame.default(formula = ms ~ SATV + SATQ2, drop.unused.levels = TRUE) : 
  variable lengths differ (found for 'SATQ2')

> m1 = lm(ms~SATV+SATQ2)
Error in model.frame.default(formula = ms ~ SATV + SATQ2, drop.unused.levels = TRUE) : 
  variable lengths differ (found for 'SATQ2')

> summary(m1)

Call:
lm(formula = ms ~ SATV + SATQ2)

Residuals:
     Min       1Q   Median       3Q      Max 
-1551.58   -12.48    45.32    99.77   168.46 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  57.5656    55.1658   1.044    0.297    
SATV          1.4313     0.1030  13.890   <2e-16 ***
SATQ2         1.9350     0.1025  18.871   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 206.3 on 490 degrees of freedom
  (7 observations deleted due to missingness)
Multiple R-squared:  0.7419,    Adjusted R-squared:  0.7409 
F-statistic: 704.3 on 2 and 490 DF,  p-value: < 2.2e-16

> gvlma(m1)

Call:
lm(formula = ms ~ SATV + SATQ2)

Coefficients:
(Intercept)         SATV        SATQ2  
     57.566        1.431        1.935  


ASSESSMENT OF THE LINEAR MODEL ASSUMPTIONS
USING THE GLOBAL TEST ON 4 DEGREES-OF-FREEDOM:
Level of Significance =  0.05 

Call:
 gvlma(x = m1) 

                       Value  p-value                   Decision
Global Stat        7.904e+03 0.00e+00 Assumptions NOT satisfied!
Skewness           1.261e+03 0.00e+00 Assumptions NOT satisfied!
Kurtosis           6.593e+03 0.00e+00 Assumptions NOT satisfied!
Link Function      2.317e-02 8.79e-01    Assumptions acceptable.
Heteroscedasticity 5.036e+01 1.28e-12 Assumptions NOT satisfied!
4

1 に答える 1

1

おそらく最も簡単なオプションは、次のようにすることです。

dta = data.frame(SATV=SATV, SATQ=SATQ, ms = ms)
lm(ms ~ SATV + SATQ, data = na.omit(dta))

NAを行ごとに削除します。

于 2014-01-17T01:23:30.910 に答える