vplot: Collection of functions to help teach matrix geometry in 2...

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

vplot(mat, type = "p", new = F, pch = 16, pop = 0, ...)

Arguments

mat
type
new
pch
pop
...

Details

vplot - plots the columns of a 2 x n matrix or a vector of length 2 - vplot adds to the current plot resizing it to include all plotted objects in a 'euclidean' frame - to start a new plot, use 'new = T' - to remove the last element added use 'vplot(pop=1)' Associated functions: - vell( mean, var) generates an ellipse, default = unit circle - vbox() generates a box - vobj() generates a circle in a box - orthog(theta) generates an orthog matrix rotating through angle theta - orthog.proj generates the matrix of an orthog. projection into span (x) - vmat( .... ) generates a 2 by n matrix Examples: vplot( new = T ) vplot( vell(), 'l' ) vplot( cbind(c(3,1),c(1,4)) %*% vell()) vplot( pop = 1) vplot( cbind(c(3,1),c(1,4)) %*% vell(), type = 'l', col = 'red')

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (mat, type = "p", new = F, pch = 16, pop = 0, ...) 
{


    if (new || !exists(".vplot")) 
        assign(".vplot", list(list(x = 0, y = 0, type = "n")), 
            pos = 1)
    a <- .vplot
    if (!missing(mat)) {
        mat <- cbind(mat)
        if (type == "v") {
            zz <- rbind(0 * mat, mat, mat/0)
            mat <- matrix(c(zz), nrow = 2)
            type = "b"
        }
        d <- dim(mat)
        if (d[1] != 2 && d[2] == 2) {
            mat <- t(mat)
            warning("mat is n x 2 and has been transposed")
        }
        a <- c(a, list(list(x = mat[1, ], y = mat[2, ], type = type, 
            pch = pch, ...)))
    }
    dat <- NULL
    for (i in seq(along = a)) {
        dat <- c(dat, a[[i]]$x, a[[i]]$y)
    }
    par(pty = "s")
    plot(range(na.omit(dat)), range(na.omit(dat)), type = "n", 
        xlab = "", ylab = "")
    if (pop > 0) {
        keep <- 1:max(1, (length(a) - (pop + 1)))
        a <- a[keep]
    }
    abline(h = 0, v = 0)
    for (i in seq(along = a)) do.call("points", a[[i]])
    assign(".vplot", a, pos = 1)
    invisible(a)
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.