cerardat_estim_nf {SPARTAAS}

Number of axes to keep

Description

Estimate the correct number of axes to keep in the regression model.

Usage

cerardat_estim_nf(df, row.sup, date)

Arguments

df

The data (data.frame) is a contingency table with the observations in rows and the technical groups in columns.

row.sup

Index of supplementary rows in df (vector).

date

The dates of each observation or NA (vector).

Value

nf

Number of axes to keep (minimal PRESS value)

MSE

plot of the Mean Squared Error.

PRESS

plot of the PRediction Error Sum Of Squares.

adj.R_sq

plot of the Coefficient of determination R².

data

data frame of MSE, PRESS and R_sq values.

Author(s)

A. COULON

L. BELLANGER

P. HUSI

References

Bellanger L. and Husi P. (2012) Statistical tool for dating and interpreting archaeological contexts using pottery. Journal of Archaeological Science, Elsevier, 39 (4), pp.777-790. doi:10.1016/j.jas.2011.06.031.

Examples

library(SPARTAAS)
data(datacerardat)
res = cerardat_estim_nf(datacerardat$df, datacerardat$row.sup, datacerardat$date)

#Number of axes to keep (minimal PRESS value)
res$nf

#the plots
res$MSE
res$PRESS
res$adj.R_sq