denpro (version 0.9.2)

profhist: Calculates the level set tree of a histogram

Description

Given a data matrix, returns the level set tree of a histogram, which is constructed from bins of equal size.

Usage

profhist(dendat, binlkm, cvol=TRUE, ccen=TRUE, cfre=FALSE)

Arguments

dendat

n*d data matrix

binlkm

positive integer; number of bins of the histogram estimate in one direction. We use the same number of bins for every direction, thus the total number of bins is binlkm to the power of d

cvol

TRUE if one wants that the volumes of the separated parts of the level sets are returned, note that one needs volumes for drawing the volume plot

ccen

TRUE if one wants that the barycenters of the separated parts of the level sets are returned

cfre

TRUE if one wants that the frequencies of the separated parts of the level sets are returned

Value

An augmented level set tree. The level set tree is a list of vectors. The elements of the vectors supply information for each node of the tree. Below we denote with "nodenum" the number of nodes of the tree.

parent

"nodenum"-vector of integers in range 0,..., nodenum-1; links to the parent of each node. Root nodes are marked with 0.

level

"nodenum"-vector of positive real numbers; level of the level set from which the set corresponding to the node is a part of.

invalue

"nodenum"-vector of positive integers; level of the level set in terms of the original frequencies (these values are not normalized so that estimate would integrate to one

volume

"nodenum"-vector of positive real numbers; gives volumes for sets corresponding to each node

center

d*nodenum-matrix; gives the barycenter for sets corresponding to each node

nodefrek

"nodenum"-vector of positive integers; number of observations in the set corresponding to the node. This is useful in cluster analysis applications.

recs

atomnum*(2*d)-matrix, where atomnum is the number of non-empty bins. Matrix defines the non-empty bins, these are the "atoms" of level sets

hisfrek

atomnum-vector of positive integers; number of observations in non-empty bins

lsets

nodenum*atomnum-matrix: describes the sets associated with the nodes. We have 1 in column "c" if the atom described in c:th row of recs is a part of the set associated with this node

See Also

plotvolu, plotbary, plottree, profgene

Examples

set.seed(1)
dendat<-matrix(rnorm(20),10)   #10*2 data-matrix
ph<-profhist(dendat,binlkm=3,cfre=TRUE)