2016-04-13 3 views
0

У меня есть две переменные, которые я хочу сравнить: одна, содержащая переменную по возрасту (непрерывную) и одну, содержащую переменную в состоянии инсульта (имела или не имела ход, коэффициент). Я могу сделать это с удовольствием, используяИспользование sjPlot в R для выполнения wilcoxon-mann whitney на факторной переменной

wilcox.test (allcoding$age~allcoding$stroketia) 

Однако, я борюсь с делать это в sjPlot, как он хочет, и чтобы быть числовыми

sjt.mwu(mwu(allcoding$age, allcoding$stroketia)) 
Performing Mann-Whitney-U-Test... 
--------------------------------- 
showing statistics between groups (x|y) 
Error in wilcox.test.default(xsub, ysub.n, paired = TRUE) : 
'y' must be numeric 

Как сделать эквивалент "~" в sjPlot, а не "," ?

Спасибо за помощь

+0

Не могли бы вы предоставить образцы данных в своем посте? Кроме того, 'sjt.mwu' должен создавать HTML-вывод, я считаю, что проблема связана с' mwu'; «Mwu» принимает два числовых вектора. Оставив данные сбоку, второй вопрос в вашем сообщении мне не совсем понятен. 'Mwu' принимает аргументы' mwu (x = ваши номера, grp = ваши группы, ...) ',' ~ 'используется в объектах формулы, но почему и где вы хотели бы включить его? – Konrad

+0

Привет, Конрад, спасибо за вашу помощь. Его трудная вещь, например, дать примеры данных, потому что ее esentational xtab я хочу сделать. Ive уже сделал это для предыдущих данных, например. 'sjt.xtab (allcoding $ курение, allcoding $ stroketia, variableLabels = c (" Курение "," Состояние штриха "), showColPerc = TRUE, tdcol.col =" black ")' Ive использовал это для создания chi-sq test, поскольку обе переменные являются факторами. Вместо этого я хочу, по сути, сделать то же самое, но создать xtab, содержащее среднее и стандартное отклонение, и результат wilcoxon, потому что одно значение является фактором, а другое непрерывным. –

+0

Тогда вам будет сложно помочь. Не могли бы вы просто опубликовать вывод о запуске 'mwu' без' sjt.mwu'? Было бы необходимо взглянуть на ваши переменные, возможно, вы сможете решить проблему с использованием 'as.numeric', но это невозможно сказать, не глядя на данные. * Изменить: * и все эти переменные являются числовыми векторами? – Konrad

ответ

0

Я могу только воспроизвести эту ошибку, когда группировка переменного типа character.

library(sjmisc) 
library(sjPlot) 
a <- runif(50, 10, 30) 
b <- sample(LETTERS[1:3], size = 50, replace = T) 
# throws an error 
mwu(a, b) 

Однако, когда вы принуждать к фактору, он работает:

mwu(a, as.factor(b)) 
sjt.mwu(mwu(a, as.factor(b))) 

Не могли бы вы попробовать это?