Beaton's sweep function

Usage

1
sweeper(a, indi, rev = FALSE)

Arguments

a
indi
rev

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
##---- 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( a, indi, rev=FALSE ){
    d = dim( a )
    if( d[1]!=d[2] ) stop( "Only works on square matrix" )
    if( rev ){ rv = -1 } else { rv = 1 }
    z <- .Call( "SWEEP",
                as.integer( indi-1 ),
                as.double( a ),
                as.integer( d[1] ),
                as.double( rv )
              );
    return( array( z, d ) );
  }