"R var.test()"의 두 판 사이의 차이

2020년 6월 22일 (월) 21:00 판

options(echo=T)
set.seed(42)
(x <- rnorm(18, mean = 0, sd = 2))
(y <- rnorm(18, mean = 1, sd = 1))
var.test(x, y)
var.test(lm(x ~ 1), lm(y ~ 1))

@@ 4번째 줄: / 4번째 줄: @@
 <source lang='r'>
+options(echo=T)
 set.seed(42)
 (x <- rnorm(18, mean = 0, sd = 2))
-##  [1]  2.7419169 -1.1293963  0.7262568  1.2657252  0.8085366 -0.2122490
-##  [7]  3.0230440 -0.1893181  4.0368474 -0.1254282  2.6097393  4.5732908
-## [13] -2.7777214 -0.5575775 -0.2666427  1.2719008 -0.5685058 -5.3129108
 (y <- rnorm(18, mean = 1, sd = 1))
-##  [1] -1.4404669  2.3201133  0.6933614 -0.7813084  0.8280826  2.2146747
-##  [7]  2.8951935  0.5695309  0.7427306 -0.7631631  1.4600974  0.3600051
-## [13]  1.4554501  1.7048373  2.0351035  0.3910736  1.5049551 -0.7170087
 var.test(x, y)
-##
-## 	F test to compare two variances
-##
-## data:  x and y
-## F = 3.8994, num df = 17, denom df = 17, p-value = 0.007605
-## alternative hypothesis: true ratio of variances is not equal to 1
-## 95 percent confidence interval:
-##   1.458635 10.424189
-## sample estimates:
-## ratio of variances
-##           3.899371
-##
 var.test(lm(x ~ 1), lm(y ~ 1))
-##
-## 	F test to compare two variances
-##
-## data:  lm(x ~ 1) and lm(y ~ 1)
-## F = 3.8994, num df = 17, denom df = 17, p-value = 0.007605
-## alternative hypothesis: true ratio of variances is not equal to 1
-## 95 percent confidence interval:
-##   1.458635 10.424189
-## sample estimates:
-## ratio of variances
-##           3.899371
 </source>