plotPost is one of John Krushke's functions for his 2011 textbook

'pphi' computes the lower or upper tail probabilities of a phi distribution.  The phi distribution is the transformation of a N(mean, sd) distribution to the [0, 1] interval by means of the N(0, 1) cumulative distribution.  

'pp' = the probability on the [0, 1] scale whose lower or upper tail probability is to be computed.

'mean', 'sd' are the mean and sd of the normal distribution that produced this phi distribution.

'lower.tail = TRUE' (default) if the lower tail probability is to be returned.

Prelec (1998) function (Eq. 3.2) is W(p)  =  exp(-beta*(-log(p))^alpha).
ARGUMENTS:
'p' = a vector of probability masses (not cumulative probabilities)
'beta' & 'alpha' are the parameters shown above.
If 'cum == TRUE' (default), the p-vector is converted to a vector of cumulative
probabilities which are then converted via the Prelec function.  If
'cum == FALSE', each component of 'p' is converted separately as if
are not part of a cumulative probability distribution.
If 'mat.out == FALSE' (default), then only the transformed cumulative
probabilities are returned; if 'mat.out == TRUE', the initial 'p' vector,
the cumulative probabilities, and the transformed cumulative
probabilities are returned.

'qphi' computes the number in the phi scale (the [0, 1] interval) which has a cumulative or decumulative probability of p.  The phi distribution is the transformation of a N(mean, sd) distribution to the [0, 1] interval by means of the N(0, 1) cumulative distribution.  

'p' = the probability whose quantile on the phi scale is to be computed.  

'mean', 'sd' are the mean and sd of the normal distribution that produced this phi distribution.

'lower.tail = TRUE' (default) if the lower tail probability is to be returned.

qt.jm computes the p-th quantile of the noncentral t distribution with noncentrality parameter ncp.

'R2' takes one or more 'lm' outputs as input and produces a
regression anova table and a table that displays R-squared.
If there is more than one 'lm' input, it also produces a table
of SS change and R-squared changes.