2016-02-28 6 views
1

Я смущен различными здравыми методами для сравнения независимых средств. Я нашел хорошее объяснение в статистических учебниках. Например, yuen() в случае равных размеров выборки. Мои образцы довольно неравны, поэтому я хотел бы попробовать метод бутстрап-т (из книги Уилкокса: Введение в надежную оценку и тестирование гипотез, с.163). Он говорит, что yuenbt() будет возможным решением.Метод Bootstrap-t для сравнения Trimmed Средство в R

Но все учебники говорят, что я могу использовать векторы здесь:

yuenbt(x,y,tr=0.2,alpha=0.05,nboot=599,side=F) 

Если я проверяю локальное описание он говорит:

yuenbt(formula, data, tr = 0.2, nboot = 599) 

Что случилось с моим суда:

x <- c(1,2,3) 
y <- c(5,6,12,30,2,2,3,65) 
yuenbt(x,y) 

Почему я не могу использовать функцию yuenbt с моими двумя векторами? Большое спасибо

ответ

2

Глядя на помощь (для тех, кто задается вопросом, yuenbt из пакета WRS2 ...) для yuenbt, в качестве аргументов используется формула и dataframe. Мое впечатление, что он ожидает данных в длинном формате.

С вашими, например, данные, мы можем добиться того, что как так:

library(WRS2) 

x <- c(1,2,3) 
y <- c(5,6,12,30,2,2,3,65) 

dat <- data.frame(value=c(x,y),group=rep(c("x","y"), c(length(x),length(y)))) 

Затем мы можем использовать функцию:

yuenbt(value~group, data=dat) 
+0

Может там произошли изменения от пакета WRS? У меня нет глубоких знаний о документации пакетов. Но эти учебники работают с векторами ... – Mac

+0

Я понятия не имею, я не знаком с пакетом, а просто скачал его, чтобы посмотреть и попытаться ответить на ваш вопрос. Учебники с кодом могут устареть. – Heroka

+0

Большое спасибо за ваше кодирование :) Теперь я могу работать над своей интерпретацией моего бутстрапинга – Mac