fit.joel: Function fit.joel

In GrimR: Calculate Optical Parameters from Spindle Stage Measurements

Description

Calculate the angle between the optical axes 2V, the optical axes in cartesian and polar coordinates and the principal axes of the dielectric tensor in cartesian and polar coordinates.

Usage

fit.joel(Data, MR = NULL, cw = c("ccw", "cw"),optimMR=FALSE)

Arguments

Data	(data frame) containing the spindle angles S and the extinction angles ES
MR	(numeric) The reference azimuth; If numeric and optimMR==TRUE, this value will be used as a starting value for further optimization. If NULL, a starting value will be guessed.
cw	(character) string "cw" for a clockwise graduated table, "ccw" for a counter-clockwise graduated table (default)
optimMR	(logical) If FALSE, the provided MR will be used without further refinement, if TRUE, the MR will be refined so as to minimize the deviance

Value

(list) with elements:

coeffs	list of the fitted parameters
covmat	matrix of covariances of the parameters
delta2V	list of estimate of 2V, its standard deviation and upper and lower confidence limits
kart	data frame with cartesian coordinates of the axes, sd, and confidence intervals
sphaer	data frame with S and ES values of the axes, sd, and confidence intervals
principal	data frame with S and MS angles to bring axes into extinction
Extinctions	data frame with S, MS, ES, calculated ES and ES-ES calculated
Wulffdat	data necessary to create a plot on the Wulff stereonet

Author(s)

Florian Dufey <[e-mail]>

Examples

# With 360 deg. data:
res<-fit.joel(Carman,MR=NULL,cw="ccw",optimMR=TRUE)
Wulffplot(res) #Plot data on a Wulff net
#with 180 degree data:
res<-fit.joel(Gunter,MR=-0.89,cw="cw",optimMR=FALSE)
Wulffplot(res) #Plot data on a Wulff net

GrimR documentation built on May 1, 2019, 9:12 p.m.