eeptools (version 0.3.1)

statamode: to mimic the mode function in Stata.

Description

This function mimics the functionality of the mode function in Stata. It does this by calculating the modal category of a vector and replacing tied categories with a "." to represent a single mode does not exist.

Usage

statamode(x,method)

Arguments

x

a numeric vector, missing values are allowed

method

a character vector of length 1 specifying the way to break ties in cases where more than one mode exists; either "stata", "sample", or "last". "stata" provides a "." if more than one mode exists. "sample" randomly samples from among the tied values for a

Value

• The modal value of a vector if a unique mode exists, else output determined by method

Details

Specifying method="stata" will result in ties for the mode being replaced with a "." character. Specifying "sample" will result in the function randomly sampling among the tied values and picking a single value. Finally, specifying "last" will result in the function picking the value that appears last in the dataset. The default behavior is stata.

See Also

table which this function uses

Examples

# for vectors
a<-c(month.name,month.name)
statamode(a,method="stata") # returns "." to show no unique mode; useful for ddply
statamode(a,method="sample") # randomly pick one
a<-c(LETTERS,"A","A")
statamode(a)