Как проверить U-тест Манна-Уитни, когда у меня есть два столбца, которые нужно сравнить

Question

Вот ссылка на мой набор данных: https://people.ucsc.edu/~mclapham/eart125/data/estuary.csv.

Вопрос, который я пытаюсь решить, заключается в следующем: соответствует ли отношение c углерод / сера к организму лимана отличаются от современных осадков и доиндустриальных образцов? Введите W statisti c ниже.

Я новичок в кодировании r и ищу небольшое руководство, я могу найти индивидуальную W statisti c для органического c углерода и для серы, но я не уверен, как найти его вместе.

Вот мой текущий код:

estuary <- read.csv("https://people.ucsc.edu/~mclapham/eart125/data/estuary.csv")

preind <- estuary$organic_carbon[estuary$timeint == "preindustrial"]
modern <- estuary$organic_carbon[estuary$timeint == "modern"]
wilcox.test(preind, modern)
preind_sul <- estuary$sulfur[estuary$timeint == "preindustrial"]
modern_sul <- estuary$sulfur[estuary$timeint == "modern"]
wilcox.test(preind_sul, modern_sul)

Answer 1

В вашей задаче кажется, что вам сначала нужно вычислить соотношение углерода / серы и сравнить это соотношение между timeint. Таким образом, вы можете сделать:

estuary$ratio <- estuary$organic_carbon / estuary$sulfur
wilcox.test(ratio ~ timeint, data = estuary)