A numeric vector of genome-wide mean ancestry (or any independent variable on the unit interval).

A numeric vector of sample sizes for each value in x. E.g., for individual diploid data, n=2.

A numeric vector containing the dependent variable: usually an allele count for each x

A vector of starting values: u and v for the logit-logistic, μ and ν for the Beta cline, a and b for the Barton cline, and U, L, m, and b for the Richards cline.

Character string indicating which cline function to fit:"logit-logit", "Beta", "Barton", "multinom", "logistic" or "Richards"

Character string indicating which algorithm to use to find the MLE: "L-BFGS-B" and "SANN" are used by the native optimization function optim, "mcmc" is a Markov Chain Monte Carlo using Metropolis-Hastings sampling. If "mcmc" is used, the following four items are also used.

The desired number of MCMC generations. The larger this number is, the greater the chance that the chain will find the maximum likelihood.

Dispersion parameters for the "mcmc" and "SANN" methods. In these methods, new parameter values are proposed by drawing values from normal distributions centered on the current value and with standard deviations from SD.

Logical: if TRUE and method="mcmc" or "SANN", starting values will be found by first using optim with "L-BFGS-B".

Logical: if TRUE and method="mcmc" and model="Beta", starting values for the Markov Chain will be found by finding the highest likelihood on a 100 x 100 grid made by mu <- seq(from=0.02,to=0.90,length.out=10); nu <- 2^(0:9)/10.

The cline function used

The optimization method used

Maximum likelihood parameter estimates

The maximum likelihood and the likelihood of the data given the naive null model E(y)=x

The number of fitted coefficients

Akaike's information criterion with sample size correction

From optim: 0 means the algorithm thinks it did a good job, otherwise the MLE may be incorrect. If mcmc was used, this will be the full Markov Chain, which could be analyzed for convergence (e.g., see coda).