Last data update: 2014.03.03

R: Value of the Oja Criterion
ojaMedianFnR Documentation

Value of the Oja Criterion

Description

The function returns the value of the Oja criterion function for a given point.

Usage

ojaMedianFn(X, x)

Arguments

X

a numeric data set or data matrix.

x

a numeric vector with the coordinates of interest.

Value

a numeric value containing the ciriterion evaluated at x for the data X.

Author(s)

ported to R by Daniel Fischer. Original C++ code by Tommi Ronkainen.

References

Oja, H. (1983), Descriptive statistics for multivariate distributions, Statistics and Probability Letters, 1, 327–332.

Ronkainen, T., Oja, H. and Orponen, P. (2002), Computation of the multivariate Oja median, in Dutter R., Filzmoser P.,Gather U. and Rousseeuw, P. J.: Developments in Robust Statistics, Heidelberg: Springer, 344–359.

See Also

ojaMedian

Examples

data(biochem)
X <- as.matrix(biochem[,1:2])
x <- ojaMedian(X)
ojaMedianFn(X, x)
ojaMedianFn(X, c(1.1, 0.4))

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(OjaNP)
Loading required package: ICS
Loading required package: mvtnorm
Loading required package: ICSNP
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/OjaNP/ojaMedianFn.Rd_%03d_medium.png", width=480, height=480)
> ### Name: ojaMedianFn
> ### Title: Value of the Oja Criterion
> ### Aliases: ojaMedianFn
> ### Keywords: multivariate nonparametric
> 
> ### ** Examples
> 
> data(biochem)
> X <- as.matrix(biochem[,1:2])
> x <- ojaMedian(X)
> ojaMedianFn(X, x)
[1] 1.5151
> ojaMedianFn(X, c(1.1, 0.4))
[1] 1.63245
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>