Proj.test: Test accuracy of projection method.

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

Usage

Proj.test(x, fun = Proj)

Arguments

x
fun

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 (x, fun = Proj) 
{
    help = "\n     # sample test:\n     zz <- matrix( rnorm(120), ncol = 3) %*% diag(c(10000,1,.000001))\n     Proj.test( zz, fun = Proj)\n     Proj.test( zz, fun = Proj.1)\n"
    pp <- fun(x)
    ret <- pp %*% pp - pp
    attr(ret, "maxabs") <- max(abs(ret))
    attr(ret, "maxabs")
  }

gmonette/spida documentation built on May 17, 2019, 7:25 a.m.