RSD (R Stochastic Dominance) is designed and developed to calculate Stochastic Dominance (SD) and Almost Stochastic Dominance (ASD) in general. In more details, given two probability mass functions (PMF), this package helps with:

• Calculate the first and second order SD values (FSD and SSD, respectively),
• Calculate the first and second order ASD values (AFSD and ASSD, respectively),
• Compare two probability distributions by these methods.

Install the released version of RSD from CRAN by:

install.packages("RSD")

Or from GitHub by:

install.packages("pak")
pak::pkg_install("ShayanTohidi/RSD")

RSD provides a function, createStochasticDominance() to create the SD object to be used in all other functions of this package. This function requires two discrete distributions. Here, the example data set data_ex will be used for creating the object:

library(RSD)
outcome1 = data_ex$yield[data_ex$gen == 'B73/PHM49']
outcome2 = data_ex$yield[data_ex$gen == 'LH74/PHN82']
pr = rep(1/29,29)
sd.obj = createStochasticDominance(outcome1, outcome2, pr, pr)

Using this code, we can compare the distributions of the performance (yield) of two cultivars. The output of this code snippet, sd.obj contains all information for performing SD and ASD comparisons. This is the main argument in the other functions of this package.