How to perform One-Sample Wilcoxon Signed Rank Test in R?

One-Sample Wilcoxon Signed Rank Test in R?, When the data cannot be assumed to be normally distributed, the one-sample Wilcoxon signed-rank test is a non-parametric alternative to the one-sample t-test.

It’s used to see if the sample’s median is the same as a known standard value (i.e. theoretical value).

The data should be symmetrically distributed around the median. In other words, the number of numbers above and below the median should be nearly equal.

Common research questions include:

  1. whether the sample’s median (m) equals the theoretical value (m0)?

2. whether the sample’s median (m) lower than the theoretical value (m0)?

3. whether the sample’s median (m) higher than the theoretical value (m0)?

In statistics, the analogous null hypothesis (H0) is defined as follows:

H0:m=m0
H0:m≤m0
H0:m≥m0

The following are the relevant alternative hypothesis (H1):

Ha:m≠m0 (different)
Ha:m>m0 (greater)
Ha:m<m0 (less)

Note that:

Two-tailed tests are used to test hypotheses 1.

One-tailed tests are used to test hypotheses 2 and 3.

Visualize your data and do one-sample calculations. Install the ggpubr R package for data visualization to perform the Wilcoxon test in R.

R base graphs can be created as explained here: R base graphs. For easy ggplot2-based data visualization,

we’ll use the ggpubr R tool. install from CRAN as follow:

install.packages("ggpubr")

Wilcoxon one-sample test with the R function

The R function wilcox.test() can be used to do a one-sample Wilcoxon test as follows.

wilcox.test(x, mu = 0, alternative = "two.sided")

x: a numeric vector containing your data values

mu: the theoretical value of mean/median The default value is zero, but you can modify it.

alternative: the opposite hypothesis “two.sided” (default), “greater” or “less” are all valid values.

We’ll use a data collection including the weights of ten mice as an example.

If the mice’s median weight differs from 22g, we’d like to know.

set.seed(123)
data <- data.frame(
  name = paste0(rep("P_", 10), 1:10),
  weight = round(rnorm(20, 30, 2), 1))

Let’s print the first 10 rows of the data

head(data, 10)
name weight
1   P_1   28.9
2   P_2   29.5
3   P_3   33.1
4   P_4   30.1
5   P_5   30.3
6   P_6   33.4
7   P_7   30.9
8   P_8   27.5
9   P_9   28.6
10 P_10   29.1

 Statistical summaries of weight

summary(data$weight)
  Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  26.10   29.05   30.25   30.28   31.10   33.60

Minimum: the lowest value

The first quartile: Only 25% of the values are lower than this.

The median is the middle value. Half of the readings are lower, while the other half are higher.

The third quartile:  The majority of the values are higher.

Maximum: the highest value

Use box plots to visualize your data.

library(ggpubr)
ggboxplot(data$weight,
          ylab = "Weight (g)", xlab = FALSE,
         ggtheme = theme_minimal())

In R, compute a one-sample Wilcoxon signed-rank test.

If the average weight of the mice differs from 22g (two-tailed test), we want to know.

# One-sample Wilcoxon test

res <- wilcox.test(data$weight, mu = 22)
res

Wilcoxon signed-rank test with continuity correction

data:  data$weight
V = 210, p-value = 9.542e-05
alternative hypothesis: true location is not equal to 22

Now we can print only the p-value

res$p.value
9.542331e-05

The test’s p-value is less than the significance level of alpha = 0.05. We can rule out the null hypothesis and infer that the mice’s average weight differs significantly from 22g.

Keep in mind:

Type this if you wish to see if the median weight of mice is less than 22g (one-tailed test):

wilcox.test(data$weight, mu = 22,
              alternative = "less")

Wilcoxon signed-rank test with continuity correction

data:  data$weight
V = 210, p-value = 1
alternative hypothesis: true location is less than 22

Alternatively, type this to see if the median weight of mice is larger than 22g (one-tailed test).

Statistical test assumptions and requirements

wilcox.test(data$weight, mu = 22,
              alternative = "greater")

Wilcoxon signed-rank test with continuity correction

data:  data$weight
V = 210, p-value = 4.771e-05
alternative hypothesis: true location is greater than 22

You may also like...

Leave a Reply

Your email address will not be published. Required fields are marked *

3 × five =