fdth (version 1.1-8)

plot.fdt: Plot fdt.default and fdt.multiple objects

Description

S3 methods for fdt.default and fdt.multiple objects. It is possible to plot histograms and poligons (absolute, relative and cumulative).

Usage

## S3 method for class 'fdt.default':
plot(x,
     type=c('fh', 'fp', 'rfh', 'rfp', 'rfph', 'rfpp', 'd',
            'cdh', 'cdp', 'cfh', 'cfp', 'cfph', 'cfpp'),
     v=FALSE,
     v.round=2,
     v.pos=3,
     xlab="Class limits",
     ylab=NULL,
     col="gray",
     xlim=NULL,
     ylim=NULL,
     main=NULL,
     x.round=2,
     x.las=1, ...)

## S3 method for class 'fdt.multiple':
plot(x,
     type=c('fh', 'fp', 'rfh', 'rfp', 'rfph', 'rfpp', 'd',
            'cdh', 'cdp', 'cfh', 'cfp', 'cfph', 'cfpp'),
     v=FALSE,
     v.round=2,
     v.pos=3,
     xlab="Class limits",
     ylab=NULL, 
     col="gray",
     xlim=NULL,
     ylim=NULL,
     main=NULL,
     main.vars=TRUE,
     x.round=2, 
     x.las=1, ...)

Arguments

x

A fdt object.

type

The type of the plot: fh: Absolut frequency histogram, fp: Absolut frequency poligon, rfh: Relative frequency histogram, rfp: Relative frequency poligon, rfph: Relative frequency (%) hi

v

Logical flag: should the values be added to the plot?

v.round

If v=TRUE, it rounds the values to the specified number of decimal places (default 0).

v.pos

If v=TRUE, a position specifier for the text. Values of 1, 2, 3 and 4, respectively indicate positions below, to the left of, above and to the right of the coordinates (default 3).

xlab

A label for the x axis.

ylab

A label for the y axis.

col

A vector of colors.

xlim

The x limits of the plot.

ylim

The y limits of the plot.

main

Title of the plot(s). This option has priority over main.vars, i.e, if any value is informed, the variable names will not be used as title of the plot(s). For fdt.multiple, the value should be a vector of characters, in t

main.vars

Logical flag: should the variables names be added as title of each plot (default TRUE)?

x.round

A numeric value to round the x ticks.

x.las

An integer which controls the orientation of the x axis labels: 0: parallel to the axes, 1: horizontal, 2: perpendicular to the axes, 3: vertical.

...

Optional plotting parameters.

Details

The result is a single histogram or poligon (absolute, relative or cummulative) for fdt.default or a set of histograms or poligon (absolute, relative or cummulative) for fdt.multiple objects. Both default and multiple try to compute the maximum number of histograms that will fit on one page, then it draws a matrix of histograms. More than one graphical device may be opened to show all histograms.

See Also

hist.data.frame provided by Hmisc package.

Examples

library(fdth)

#======================
# Vectors: univariated
#======================
set.seed(1)

x <- rnorm(n=1e3,
           mean=5,
           sd=1)

(d <- fdt(x))

# Histograms
plot(d)  # Absolut frequency histogram

plot(d,
     main='My title')

plot(d,
     x.round=3,
     col='darkgreen')

plot(d,
     x.las=2)

plot(d,
     x.round=3,
     x.las=2,
     xlab=NULL)

plot(d, 
     v=TRUE, 
     cex=.8,
     x.round=3,
     x.las=2,
     xlab=NULL,
     col=rainbow(11))

plot(d,
     type='fh')    # Absolut frequency histogram

plot(d,
     type='rfh')   # Relative frequency histogram

plot(d,
     type='rfph')  # Relative frequency (%) histogram

plot(d,
     type='cdh')   # Cumulative density histogram

plot(d,
     type='cfh')   # Cumulative frequency histogram

plot(d,
     type='cfph')  # Cumulative frequency (%) histogram

# Poligons
plot(d,
     type='fp')    # Absolut frequency poligon

plot(d,
     type='rfp')   # Relative frequency poligon

plot(d,
     type='rfpp')  # Relative frequency (%) poligon

plot(d,
     type='cdp')   # Cumulative density poligon

plot(d,
     type='cfp')   # Cumulative frequency poligon

plot(d,
     type='cfpp')  # Cumulative frequency (%) poligon

# Density
plot(d,
     type='d')     # Density

# Theoretical curve and fdt
x <- rnorm(1e5,
           mean=5,
           sd=1)

plot(fdt(x,
         k=100), 
     type='d', 
     col=heat.colors(100))

curve(dnorm(x,
            mean=5,
            sd=1),
      col='darkgreen',
      add=TRUE,
      lwd=2)

#=============================================
# Data.frames: multivariated with categorical
#=============================================
mdf <- data.frame(X1=rep(LETTERS[1:4], 25),
                  X2=as.factor(rep(1:10, 10)),
                  Y1=c(NA, NA, rnorm(96, 10, 1), NA, NA),
                  Y2=rnorm(100, 60, 4),
                  Y3=rnorm(100, 50, 4),
                  Y4=rnorm(100, 40, 4))

# Histograms
(d <- fdt(mdf))

plot(d,
     v=TRUE,
     cex=.8)

plot(d,
     col='darkgreen',
     ylim=c(0, 40))

plot(d,
     col=rainbow(8), 
     ylim=c(0, 40),
     main=LETTERS[1:4])

plot(d,
     type='fh')

plot(d,
     type='rfh')

plot(d,
     type='rfph')

plot(d,
     type='cdh')

plot(d,
     type='cfh')

plot(d,
     type='cfph')

# Poligons
plot(d,
     v=TRUE,
     type='fp')

plot(d,
     type='rfp')

plot(d,
     type='rfpp')

plot(d,
     type='cdp')

plot(d,
     type='cfp')

plot(d,
     type='cfpp') 

# Density
plot(d,
     type='d')     

levels(mdf$X1)

plot(fdt(mdf,
         k=5,
         by='X1'),
     ylim=c(0, 12))

levels(mdf$X2)

plot(fdt(mdf,
         breaks='FD', 
         by='X2'))

plot(fdt(mdf,
         k=5,
         by='X2'))  # It is dificult to compare

plot(fdt(mdf,
         k=5,
         by='X2'),
     ylim=c(0, 8))  # Easy

plot(fdt(iris,
         k=5))

plot(fdt(iris,
         k=5),
     col=rainbow(5))

plot(fdt(iris, 
         k=5,
         by='Species'),
     v=TRUE)

d <- fdt(iris,
         k=10)

plot(d)

plot(d, 
     type='d')

#=========================
# Matrices: multivariated
#=========================
plot(fdt(state.x77))

plot(fdt(volcano))