R Graphical Manual

Browse All

Last data update: 2014.03.03

R: S3 methods for manipulating eDist objects.

eDist

R Documentation

S3 methods for manipulating eDist objects.

Description

S3 methods for manipulating eDist objects

Usage

## S3 method for class 'eDist'
logLik(object, ...)

## S3 method for class 'eDist'
AIC(object, ..., k = 2)

AICc(object)

## S3 method for class 'eDist'
AICc(object, ...)

## S3 method for class 'eDist'
vcov(object, ..., corr = FALSE)

BIC(object)

## S3 method for class 'eDist'
BIC(object, ...)

MDL(object)

## S3 method for class 'eDist'
MDL(object, ...)

## S3 method for class 'eDist'
print(x, ...)

## S3 method for class 'eDist'
plot(x, ...)

Arguments

`object,`	x An object of class eDist, usually the output of a parameter estimation function.
`...`	Additional parameters
`k`	numeric, The penalty per parameter to be used; the default k = 2 is the classical AIC.
`corr`	logical; should vcov() return correlation matrix (instead of variance-covariance matrix).
`x,`	A list to be returned as class eDist.
`plot`	logical; if TRUE histogram, P-P and Q-Q plot of the distribution returned else only parameter estimation is returned.

Note

The MDL only works for parameter estimation by numerical maximum likelihood.

Author(s)

A. Jonathan R. Godfrey, Sarah Pirikahu, and Haizhen Wu.

References

Myung, I. (2000). The Importance of Complexity in Model Selection. Journal of mathematical psychology, 44(1), 190-204.

Examples

X <- rnorm(20)
est.par <- eNormal(X, method ="numerical.MLE")
logLik(est.par)
AIC(est.par)
AICc(est.par)
BIC(est.par)
MDL(est.par)
vcov(est.par)
vcov(est.par,corr=TRUE)
print(est.par)
plot(est.par)

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(ExtDist)

Attaching package: 'ExtDist'

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

    BIC

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ExtDist/eDist.Rd_%03d_medium.png", width=480, height=480)
> ### Name: eDist
> ### Title: S3 methods for manipulating eDist objects.
> ### Aliases: AIC.eDist AICc AICc.eDist BIC BIC.eDist MDL MDL.eDist eDist
> ###   logLik.eDist plot.eDist print.eDist vcov.eDist
> 
> ### ** Examples
> 
> X <- rnorm(20)
> est.par <- eNormal(X, method ="numerical.MLE")
> logLik(est.par)
[1] -29.12681
> AIC(est.par)
> AICc(est.par)
[1] 62.9595
> BIC(est.par)
[1] 64.24509
> MDL(est.par)
[1] 32.39431
> vcov(est.par)
              mean            sd
mean  5.388367e-02 -1.457494e-08
sd   -1.457494e-08  2.694185e-02
> vcov(est.par,corr=TRUE)
              mean            sd
mean  1.000000e+00 -3.825291e-07
sd   -3.825291e-07  1.000000e+00
> print(est.par)

Parameters for the Normal distribution. 
(found using the  numerical.MLE method.)

 Parameter     Type   Estimate      S.E.
      mean location -0.1924116 0.2321286
        sd    scale  1.0381105 0.1641397


> plot(est.par)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>