43597


# RのF検定について

以下の二つの標本に対するF検定を考えます。

[code]

    data1 = rnorm(5, mean=0, sd=1)
    data2 = rnorm(10, mean=0, sd=2)
    
[/code]

一般的には、以下のようにtest関数でF検定を行います。

[code]

    var.test(data1, data2)
    
[/code]

出力結果の例を示します。

[code]

    [1] "var.test:"
    
            F test to compare two variances
    
    data:  data1 and data2
    F = 3.4344, num df = 4, denom df = 9, p-value = 0.1149
    alternative hypothesis: true ratio of variances is not equal to 1
    95 percent confidence interval:
      0.7279282 30.5824393
    sample estimates:
    ratio of variances
              3.434423
    
[/code]

現在、このtest関数の内部で行っていることを実装しているところです。  
しかし、出力結果中の

  * num df
  * denom df
  * p-value

を求める式がわからずに困っています。  
よろしくお願いします。


CC BY-SA 3.0
2018-04-27T03:44:57.320
43595
2018-04-27T05:12:47.313
28333
question
r
R言語でF検定を行うコードを実装したい（p値の算出など）


クラシカルな方法としては、

[code]

    result = var.test(data1, data2)
    print(result$statistic)
    print(result$parameter["num df"])
    print(result$parameter["denom df"])
    print(result$p.value)
    
[/code]

今後のおすすめの方法はbroomを使用することです。

[code]

    require("broom")    # インストール未ならinstall.package("broom")
    result = tidy(var.test(data1, data2))
    print(result["statistic"])
    print(result["num.df"])
    print(result["denom.df"])
    print(result["p.value"])
    
[/code]


CC BY-SA 3.0
2018-04-27T05:12:47.313
43597
2018-04-27T05:12:47.313
6092
43595
answer