cons.angular.dat {BMAmevt}

Angular data set generation from unit Frechet data.

Description

Builds an angular data set, retaining the points with largest radial component.

Usage

cons.angular.dat(
  coordinates = c(1, 2, 3),
  frechetDat = get("frechetdat"),
  n = 100,
  displ = TRUE,
  invisible = TRUE,
  add = FALSE,
  lab1 = "w1",
  lab2 = "w2",
  lab3 = "w3",
  npoints = 60,
  col.polygon = "white",
  ...
)

Arguments

coordinates

Index vector of the columns in frechetDat to be retained to construct the angular data set.

frechetDat

The data set. A matrix: each row is a multivariate record. May contain NA's.

n

The number of desired observations in the final angular data set. Should be less than nrow(frechetDat)

displ

logical. Should the angular data set be plotted ?

invisible

logical. Should the result be returned as invisible ?

add

logical. Only used when displ==TRUE. Should the points be added to the current plot ?

lab1

Character string: label for first component.

lab2

Character string: label for second component.

lab3

Character string: label for third component.

npoints

The number of grid nodes on the squared grid containing the desired triangle.

col.polygon

The background color outside the simplex.

...

Additional graphical parameters and arguments to be passed to function plot.window and points.

Details

The data set frechetDat is assumed to be marginally unit Frechet distributed.

Value

The angular data set: A n*length(coordinates) matrix, containing values between zero and one, which rows sum to one: Each row is thus a point on the unit simplex of dimension length(coordinates)-1. Returned as invisible if invisible==TRUE.

Examples

 ## Not run: cons.angular.dat()