Cetin Sert - 1 year ago 67

R Question

I have this example from the coin package of R:

`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"))

I want to assign alpha a different value, how can I do this??

This doesn't work!

`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"))

Recommended for you: Get network issues from **WhatsUp Gold**. **Not end users.**

Answer Source

The alpha levels are hardcoded and fixed at 0.99 If you want to change that, then you have to download the package source, change the levels and compile the package. The levels are coded in the Methods.R file. Search for binom.test or conf.level

You could ask the package author to change the package so you can set the level yourself. But bear in mind that the package author is not obliged to do that!

Recommended from our users: **Dynamic Network Monitoring from WhatsUp Gold from IPSwitch**. ** Free Download**