критерий Вилкоксона с использованием стратификации данных - PullRequest
0 голосов
/ 19 декабря 2018

У меня действительно основная проблема.У меня есть концентрации одного химического вещества, хранящиеся в одном столбце, и пол участника исследования во втором столбце.Каков код для проведения теста Уилкоксона, чтобы увидеть, есть ли разница между концентрациями, обнаруженными у мальчиков, и концентрациями, обнаруженными у девочек?Некоторое объяснение кода также было бы полезно для меня, чтобы понять, как он работает.Спасибо!

Я получил этот код для теста ANOVA, что тоже хорошо.Может кто-нибудь сказать мне, если он делает то, что мне нужно?

av <- aov(UC_MEHP ~ BQF05C1, data=data)
av
summary(av)

вывод выглядит так

> av <- aov(UC_MEHP ~ BQF05C1, data=data)
> av
Call:
   aov(formula = UC_MEHP ~ BQF05C1, data = data)

Terms:
                  BQF05C1 Residuals
Sum of Squares     0.3445 2917.4564
Deg. of Freedom         1       151

Residual standard error: 4.395555
Estimated effects may be unbalanced
21 observations deleted due to missingness
> summary(av)
             Df Sum Sq Mean Sq F value Pr(>F)
BQF05C1       1    0.3   0.344   0.018  0.894
Residuals   151 2917.5  19.321               
21 observations deleted due to missingness

Извините, я знаю, что это не очень сложный вопрос...

1 Ответ

0 голосов
/ 19 декабря 2018

С ?wilcox.test:

## S3 метод для класса 'формула' wilcox.test (формула, данные, подмножество, na.action, ...)
...
формула: формула вида 'lhs ~ rhs', где 'lhs' - это числовая переменная, задающая значения данных, а 'rhs' - коэффициент с двумя уровнями, дающими соответствующие группы.

Так что wilcox.test(UC_MEHP ~ BQF05C1, data=data) должно работать (при условии, что BQF05C1 - это столбец с указанием пола, а UC_MEHP - это концентрация).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...