curvigram: Computes declination curvigram

Description Usage Arguments References See Also Examples

View source: R/data-analysis.R

Description

This function computes the curvigram of declinations, using provided measurement uncertainty and a Gaussian kernel, i.e. using the method of Silva (2017). When all measurements have the same associated uncertainty this function wraps density and accepts the same input for bw in unc.

Usage

1
curvigram(dec, unc = 2, norm = F, cut = 4, range)

Arguments

dec

Array of declination values.

unc

(Optional) Uncertainty in declination values. Either a single value or string to be applied to all measurements (see bw.nrd), or an array of values of the same length as dec. Defaults to 2 degrees.

norm

(Optional) Boolean specifying whether the resulting curvigram should be normalized to unity. Defaults to FALSE.

cut

(Optional) Number of uncertainties beyond the extremes of the data at which to trim the curvigram. Defaults to 4. See density.

range

(Optional) As an alternative to cut you can stipulate the range of declination values to output as an array of two values. See from, to in density.

References

Silva, Fabio (2017) Inferring Alignments I: Exploring the Accuracy and Precision of Two Statistical Approaches, Journal of Skyscape Archaeology 3(1), 93-111. DOI: 10.1558/jsa.31958

See Also

density

Examples

1
2
3
4
# Curvigram of Ruggles' Recumbent Stone Circle data:
data(RugglesRSC)
curv <- curvigram(RugglesRSC$Dec, 2)
plot(curv)

f-silva-archaeo/skyscapeR documentation built on Jan. 15, 2019, 3:22 a.m.