Here's my modest contribution:
smode<-function(x){
xtab<-table(x)
modes<-xtab[max(xtab)==xtab]
mag<-as.numeric(modes[1]) #in case mult. modes, this is safer
themodes<-names(modes)
mout<-list(themodes=themodes,modeval=mag)
return(mout)
}
(I called it 'smode' because there's a 'mode' func somewhere which is
related to the data type or something)
Carl