mvmesh (version 1.0)

mvhist: Multivariate histograms

Description

Tabulate and plot histograms for data, including directional histograms

Usage

histDirectional( x, k, p=2, plot.type="default", freq=TRUE, positive.only=FALSE, 
     report="summary", ... )
histRectangular( x, breaks=10, plot.type="default", freq=TRUE, report="summary", ... )
histSimplex( x, S, plot.type="default", freq=TRUE, report="summary", ... )

TallyHrep( x, H, report="summary" )  
DrawPillars( S, height, shift=rep(0.0,3), ... )

Arguments

x

data in an (n x d) matrix; rows are d-dimensional data vectors

k

number of subdivisions

p

power of p-norm

freq

TRUE for a frequency histogram, FALSE for a relative frequency histogram

breaks

specifes the subdivision of the region; see 'breaks' in RectangularMesh

plot.type

type of plot, see details below

positive.only

If TRUE, look only in the first octant

report

level of warning messages; one of "summary", "all", "none".

...

Optional arguments to plot

S

(vps x d x nS) array of simplices in V representation, see V2Hrep

H

array of simplices in H representation, see V2Hrep

height

vector of length nS giving the heights of the pillars

shift

shift of the pillars, typically (0,0,0) for 2d data or (0,0,z0) for 3d data

Value

• A plot is drawn (unless plot.type="none"). A list is returned invisibly, with fields:
  • counts - frequency count in each bin
  • nrejects - number of x values not in any bin
  • nties - number of points in more than one bin (if bins are set up to be non-overlapping, this should only occur on a shared edge between two simplices
  • nx - total number of data points in x
  • rel.freq - counts/nx
  • rel.rejects - nrejects/nx
  • mesh - object of type mvmesh, seemvmesh
  • plot.type - input value
  • report - input value

Details

Calculate and plot multivariate histograms. In all cases, the bins are simplices computed from some description. Then the number in each simplex is tallied using TallyHrep. 'plot.type' values depend on the type of plot being used. Possible values are:

• "none" - does not show a plot, just return the counts
• "index" - shows a histogram of simplex index number versus count, does not show the geometry, but works in any dimension
• "pillars" - shows a 3D plot with pillars/columns having base the shape of the simplices and height proportional to frequency counts. When the points are 2D, this works forhistRectangularandhistSimplex; when the points are 3D, this only works forhistRectangular
• "counts" - shows frequency counts as a number in the center of each simplex
• "radial" -histDirectionalonly, shows radial spikes proportional to the counts
• "grayscale" -histDirectionalonly, shows radial spikes proportional to the counts
• "orthogonal" -histDirectionalonly, shows radial spikes proportional to the counts
• "default" - type depends on the dimension of the data and type of histogram

Examples

# two dimensional, isotropic
x <- matrix( rnorm(8000), ncol=2 )
histDirectional( x, k=1 )
histRectangular( x, breaks=5 )

# three dimensional positive data
x <- matrix( abs(rnorm(9000)), ncol=3 )
histDirectional( x, k=3, positive.only=TRUE, col='blue', lwd=3 )
histRectangular( x, breaks=4 )

demo(mvhist) # shows a range of multivariate histograms