Rにアルファを設定するには？

Question

私はRのコインパッケージからthis exampleあります

library(coin) 
    library(multcomp) 
    ### Length of YOY Gizzard Shad from Kokosing Lake, Ohio, 
    ### sampled in Summer 1984, Hollander & Wolfe (1999), Table 6.3, page 200 
    YOY <- data.frame(length = c(46, 28, 46, 37, 32, 41, 42, 45, 38, 44, 
           42, 60, 32, 42, 45, 58, 27, 51, 42, 52, 
           38, 33, 26, 25, 28, 28, 26, 27, 27, 27, 
           31, 30, 27, 29, 30, 25, 25, 24, 27, 30), 
        site = factor(c(rep("I", 10), rep("II", 10), 
            rep("III", 10), rep("IV", 10)))) 

    ### Nemenyi-Damico-Wolfe-Dunn test (joint ranking) 
    ### Hollander & Wolfe (1999), page 244 
    ### (where Steel-Dwass results are given) 
    NDWD <- oneway_test(length ~ site, data = YOY, 
     ytrafo = function(data) trafo(data, numeric_trafo = rank), 
     xtrafo = function(data) trafo(data, factor_trafo = function(x) 
      model.matrix(~x - 1) %*% t(contrMat(table(x), "Tukey"))), 
     teststat = "max", distribution = approximate(B = 90000)) 

    ### global p-value 
    print(pvalue(NDWD)) 

    ### sites (I = II) != (III = IV) at alpha = 0.01 (page 244) 
    print(pvalue(NDWD, method = "single-step")) 

私はアルファに異なる値を割り当てたい、私はこれをどのように行うことができます？

これは機能しません。

library(coin) 
    library(multcomp) 
    ### Length of YOY Gizzard Shad from Kokosing Lake, Ohio, 
    ### sampled in Summer 1984, Hollander & Wolfe (1999), Table 6.3, page 200 
    YOY <- data.frame(length = c(46, 28, 46, 37, 32, 41, 42, 45, 38, 44, 
           42, 60, 32, 42, 45, 58, 27, 51, 42, 52, 
           38, 33, 26, 25, 28, 28, 26, 27, 27, 27, 
           31, 30, 27, 29, 30, 25, 25, 24, 27, 30), 
        site = factor(c(rep("I", 10), rep("II", 10), 
            rep("III", 10), rep("IV", 10)))) 

    ### Nemenyi-Damico-Wolfe-Dunn test (joint ranking) 
    ### Hollander & Wolfe (1999), page 244 
    ### (where Steel-Dwass results are given) 
    NDWD <- oneway_test(length ~ site, data = YOY, 
     ytrafo = function(data) trafo(data, numeric_trafo = rank), 
     xtrafo = function(data) trafo(data, factor_trafo = function(x) 
      model.matrix(~x - 1) %*% t(contrMat(table(x), "Tukey"))), 
     teststat = "max", distribution = approximate(B = 90000), 
     alpha = 0.05) 

    ### global p-value 
    print(pvalue(NDWD)) 

    ### sites (I = II) != (III = IV) at alpha = 0.05 (default was 0.01) (page 244) 
    print(pvalue(NDWD, method = "single-step")) 

Answer 1

アルファレベルはハードコードされており、0.99に固定されています。変更したい場合は、パッケージソースをダウンロードし、レベルを変更してパッケージをコンパイルする必要があります。レベルはMethods.Rファイルにコード化されています。 binom.testまたはconf.levelを検索

自分でレベルを設定できるようにパッケージの変更をパッケージ作成者に依頼することができます。しかし、パッケージの作者はそれをする義務がないことに注意してください！

Answer 2

あなたができないように思われる：oneway_test()はwilcox_testとnormal_test行うのに対し、引数conf.levelを持っていません。これはすべて文書化されています。help(oneway_test)を参照してください。