To produce a plot with 90% CIs for the mean and PIs for the regression of species on ph levels in the "acidfish" dataset, you would type the following
into the command line (after loading the function as above):
> regplot.confbands.fun(acidfish$ph,acidfish$species,.90,T,T,F)