R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Detection of Outliers in Circular-Circular Regression

DMCEE

R Documentation

Detection of Outliers in Circular-Circular Regression

Description

This function calculates the absolute values of the difference between the values of MCE and MCe statistic. Then, it draws the scatter plot of them and estimates the concentration parameter of k. Given the sample size and the estimated value of K, cut-off point from the table DMCE in the significance level of 0.05 or 0.1 will be found. Outliers are locatedon the top of teh line corresponding to the cut-off point.

Usage

DMCEE(x, y, b)

Arguments

`x`	independent variable on model y_i=α+β x_i+ε_i (mod 2π) (i=1,2,...,n)
`y`	the response variable on model y_i=α+β x_i+ε_i (mod 2π) (i=1,2,...,n)
`b`	the level of significance (0.05 or 0.1)

Details

The ith observation is identified as an outlier if the absolute values of the difference between the values of MCE and MCe statistic exceeds a pre-specified cut-off point.

Author(s)

Azade Ghazanfarihesari, Majid Sarmad

References

A. H. Abuzaid, A. G. Hussin & I. B. Mohamed (2013) Detection of outliers in simple circular regression models using the mean circular error statistics

Examples

data(wind2)
DMCEE(wind2[,1], wind2[,2], .9)

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(CircOutlier)
Loading required package: CircStats
Loading required package: MASS
Loading required package: boot
Loading required package: circular

Attaching package: 'circular'

The following objects are masked from 'package:CircStats':

    A1, A1inv, I.0, I.1, I.p, deg, plot.edf, pp.plot, rad, rose.diag,
    rstable

The following objects are masked from 'package:stats':

    sd, var


Attaching package: 'CircOutlier'

The following object is masked from 'package:circular':

    wind

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/CircOutlier/DMCEE.Rd_%03d_medium.png", width=480, height=480)
> ### Name: DMCEE
> ### Title: Detection of Outliers in Circular-Circular Regression
> ### Aliases: DMCEE
> 
> ### ** Examples
> 
> data(wind2)
> DMCEE(wind2[,1], wind2[,2], .9)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>