2016-05-06 3 views
1

Я не продвинутый пользователь R.Anyway Я просто хочу сделать тест Wilcox на двух наборах данных:R: получить значение слота от объекта S4 (ScalarIndependenceTest)

library(coin) 
GroupA = c(2,4,3,1,2,3,3,2,3,1) 
GroupB = c(3,5,4,2,4,3,5,5,3,2) 
g = factor(c(rep("GroupA", length(GroupA)), rep("GroupB", length(GroupB)))) 
v = c(GroupA, GroupB) 
ss = wilcox_test(v ~ g, distribution="exact") 

То, что я хочу, получить Zvalue и вычислить размер эффекта. ss - объект S4, когда я его печатаю, я вижу, что он сообщает значение Z, но я не могу найти его в ss @ statistic. Единственный слот, который имеет значение, которое я ожидаю «тестовая статистика», но даже если это г-значение теста, когда я делаю:

slot(ss, 'teststatistic') 

Я получаю ошибку:

no slot of name "teststatistic" for this object of class "ScalarIndependenceTest" 

может кто-нибудь пожалуйста, дайте мне подсказку? Спасибо

+0

В дальнейшем укажите любые необорудованные пакеты. Является 'wilcox_test' из пакета' coin'? – nrussell

+0

Да, я отредактировал вопрос. благодаря – Alex

ответ

1

print эквивалент классов S4 является show метод, который может быть проверен с помощью getMethod. В этом случае

ss 
# 
# Exact Wilcoxon-Mann-Whitney Test 
# 
#data: v by g (GroupA, GroupB) 
#Z = -2.1095, p-value = 0.0385 
#alternative hypothesis: true mu is not equal to 0 

getMethod("show","ScalarIndependenceTest") 
#Method Definition: 
# 
# function (object) 
# { 
#  distname <- switch(class([email protected]), AsymptNullDistribution = "Asymptotic", 
#   ApproxNullDistribution = "Approximative", ExactNullDistribution = "Exact") 
#  RET <- list(statistic = setNames([email protected]@teststatistic, 
#   nm = "Z"), p.value = [email protected]@pvalue([email protected]@teststatistic), 
#   alternative = [email protected]@alternative, data.name = varnames([email protected]), 
#   method = paste(distname, [email protected])) 
#... 
#... 
# } 

[email protected]@teststatistic 
# GroupA 
#-2.109531 

[email protected]@pvalue([email protected]@teststatistic) 
#[1] 0.03850484 
Смежные вопросы