Calculate the mode of a sample: both modal value(s) and the corresponding frequency
modeOf(x, na.rm = TRUE) maxFreq(x, na.rm = TRUE)
- A vector containing the observations.
- Logical value indicating whether NA values should be removed.
These two functions can be used to calculate the mode (most frequently observed value) of a sample, and the actual frequency of the modal value. The only complication is in respect to missing data. If
na.rm = FALSE, then there are multiple possibilities for how to calculate the mode. One possibility is to treat
NA as another possible value for the elements of
x, and therefore if
NA is more frequent than any other value, then
NA is the mode; and the modal frequency is equal to the number of missing values. This is the version that is currently implemented.
Another possibility is to treat
NA as meaning "true value unknown", and to the mode of
x is itself known only if the number of missing values is small enough that -- regardless of what value they have -- they cannot alter the sample mode. For instance, if
c(1,1,1,1,2,2,NA), we know that the mode of
1 regardless of what the true value is for the one missing datum; and we know that the modal frequency is between 4 and 5. This is also a valid interpretation, depending on what precisely it is the user wants, but is not currently implemented.
Because of the ambiguity of how
na.rm = FALSE should be interpreted, the default value has been set to
na.rm = TRUE, which differs from the default value that I've used elsewhere in the package.
modeOffunction returns the mode of
x. If there are ties, it returns a vector containing all values of
xthat have the modal frequency. The
maxFreqfunction returns the modal frequency as a numeric value.
This package is under development, and has been released only due to teaching constraints. Until this notice disappears from the help files, you should assume that everything in the package is subject to change. Backwards compatibility is NOT guaranteed. Functions may be deleted in future versions and new syntax may be inconsistent with earlier versions. For the moment at least, this package should be treated with extreme caution.
# simple example eyes <- c("green","green","brown","brown","blue") modeOf(eyes) maxFreq(eyes) # vector with missing data eyes <- c("green","green","brown","brown","blue",NA,NA,NA) # returns NA as the modal value. modeOf(eyes, na.rm = FALSE) maxFreq(eyes, na.rm = FALSE) # returns c("green", "brown") as the modes, as before modeOf(eyes, na.rm = TRUE) maxFreq(eyes, na.rm = TRUE)