fmean {StratigrapheR}

Fischer mean

Description

Fischer mean

Usage

fmean(dec = NA, inc = NA, int = 1, x = NA, y = NA, z = NA, id = NULL, cart = F)

Arguments

dec

declination of the data; it is the angle from the north taken on an horizontal plane. It is measured clockwise from North and ranges from 0 to 360° (Tauxe 2010). Values outside this range are corrected by incfix().

inc

inclination of the data; it is the angle from the horizontal, is positive downward, and ranges from +90° for straight down to -90° for straight up (Tauxe, 2010). Values outside this range are corrected by incfix().

int

intensity of the data. Defaults to one (unit sphere).

x, y, z

cartesian coordinates. x is the North, y the East, and z straight down. If dec and inc are not provided they are used to be converted back in dec, inc and int data. Output is corrected by incfix().

id

a name for each point, identifying each group of points you would like to treat separately

cart

whether to output as cartesian coordinates, defaults to F

Value

a list of coordinates for the fischer mean, in cartesian form or dec, inc, int form

See Also

fmod, dipfix and incfix

Examples

dec <- c(rnorm(10, mean = 45, sd = 5), rnorm(10, mean = 20, sd = 5))
inc <- c(rnorm(10, mean = 45, sd = 5), rnorm(10, mean = 20, sd = 5))
id  <- c(rep(1, 10), rep(2, 10))

earnet()
earpoints(dec, inc)

fm <- fmean(dec, inc, id = id)

earpoints(fm $dec, fm$inc, l = list(bg = "red"))