Sorted the factor levels in the order required in the legend to avoid problems later. Reorganized ggplot calls.

added na.rm, apply_to and alpha parameters and their checks. Changed na.rm=T for na.rm=na.rm everywhere. Added documentation.