vec <- scan('stdin',quiet=TRUE) vec <- vec[-1] modenum <- as.numeric(names(which(table(vec) == max(table(vec))))) if (length(modenum) > 1) { cat(min(modenum)) } else{ cat(modenum) }