sigTest: Perform a null hypothesis significance test of a given...

Description Usage Arguments See Also Examples

View source: R/sig-test.R

Description

This function performs a null hypothesis significance test, for a given dataset and outputs a p-value as well as all the information needed for plotting.

Usage

1
2
sigTest(data, type = "curv", ncores = parallel::detectCores() - 1,
  nsims = 2000, conf = 0.95, prec = 0.01, range, verbose = T, ...)

Arguments

data

Data frame including latitude, azimuth and horizon altitude of sites, created using reduct.compass or reduct.theodolite, plus the uncertainty in azimuthal measurement, which should be added into a column named Azimuth.Uncertainty.

type

(Optional) Type of data visualizarion you want to conduct the test on. Current options are 'curv' for curvigram and 'hist' for histogram. Defaults to 'curv' for curvigram.

ncores

(Optional) Number of processing cores to use for parallelisation. Defaults to the number of available cores minus 1.

nsims

(Optional) Number of simulations to run. The higher this number the slower this process will be, but the lower it is the less power the method has. Defaults to 2000 as a base minimum to test for significance at the p=0.0005 level, but the recommended value is 10,000.

conf

(Optional) Confidence level for p-value calculation and output. Defaults to 0.95, i.e. a 95% confidence envelope.

prec

(Optional) Smallest possible azimuth for the random sampler, i.e. precision being considered. Defaults to 0.01º.

range

(Optional) Range of declination values to consider.

verbose

(Optional) Boolean to decide whether or not

...

Other parameters to be passed on to histogram or curvigram, as appropriate.

See Also

reduct.compass, reduct.theodolite

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## Not run: 
loc <- c(35,-7)
mag.az <- c(89.5, 105, 109.5)
data <- reduct.compass(loc, mag.az, "2016/04/02", alt=c(1,2,0))
data$Azimuth.Uncertainty <- 1  # adds the information on the preision of the azimuthal meaurement
sig <- sigTest(data)

plot(sig)

## End(Not run)

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