Define mode limits by mouse clicks.

Description

This function allows to define limits for robust end-members by mouse clicks on a combined plot output, showing a histogram and all end-members together. Clicks must be placed in the order lower limit, upper limit - for each end-member successively.

Usage

define.limits(data, n, classunits)

Arguments

Value

Numeric matrix with limit classes. The first row contains lower limits, the second row upper limits for each end-member.

Author(s)

Michael Dietze, Elisabeth Dietze

See Also

test.robustness, robust.EM

Examples

## load example data set
data(X, envir = environment())

## Test robustness
q <- 4:7
l <- seq(from = 0, to = 0.1, by = 0.02)
TR <- test.robustness(X = X, q = q, l = l)

## define 2 limits by mouse clicks (uncomment to use).
# limits <- define.limits(data = TR, n = 2)
# limits

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(EMMAgeo)
Loading required package: GPArotation
Loading required package: limSolve
Loading required package: shape
Loading required package: shiny
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/EMMAgeo/define.limits.Rd_%03d_medium.png", width=480, height=480)
> ### Name: define.limits
> ### Title: Define mode limits by mouse clicks.
> ### Aliases: define.limits
> ### Keywords: EMMA
> 
> ### ** Examples
> 
> ## load example data set
> data(X, envir = environment())
> 
> ## Test robustness
> q <- 4:7
> l <- seq(from = 0, to = 0.1, by = 0.02)
> TR <- test.robustness(X = X, q = q, l = l)
> 
> ## define 2 limits by mouse clicks (uncomment to use).
> # limits <- define.limits(data = TR, n = 2)
> # limits
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>