MAINT.Data (version 1.0.1)

fulltle-methods: Methods for Function fulltle in Package ‘MAINT.Data’

Description

Performs maximum trimmed likelihood estimation by an exact algorithm (full enumeratiom of all k-trimmed subsets)

Usage

fulltle(Idt, alpha=0.75, reweighted=TRUE, CorrF=c("smallsmp","consistent","none"),
  outlin=c("MidPandLogR","MidP","LogR"),
  CovCase=1:4, SelCrit=c("BIC","AIC"),
  force=FALSE,otpType=c("OnlyEst","SetMD2andEst"), …)

Arguments

Idt

An IData object representing interval-valued entities.

alpha

Numeric parameter controlling the size of the subsets over which the trimmed likelihood is maximized; roughly alpha*Idt@NIVar observations are used for computing the trimmed likelihood. Allowed values are between 0.5 and 1.

reweighted

Should a (Re)weighted estimate of the covariance matrix be used in the computation of the trimmed likelihood or just a “raw” covariance estimate; default is (Re)weighting.

CorrF

Whether to use sample correction factors for the covariance estimates; options are “smallsmp” (default), for small sample correction, “consistent”, for consistency correction and “none” for no correction.

outlin

The type of outliers to be consideres. “MidPandLogR” if outliers may be present in both MidPpoints and LogRanges, “MidP” if outliers are only present in MidPpoints, or “LogR” if outliers are only present in LogRanges.

CovCase

Configuration of the variance-covariance matrix: a set of integers between 1 and 4.

SelCrit

The model selection criterion.

force

A boolean flag indicating whether, for moderate or large data sets the algorithm should proceed anyway, regardless of an expected long excution time, due to exponential explosions in the number of different subsets that need to be avaluated by fulltle

otpType

The amount of output returned by fulltle. Current options are “OnlyEst” (default) where only an ‘IdtE’ object with the fulltle estimates is returned and “SetMD2andEst” which returns a list with an ‘IdtE’ object of fulltle estimates, a vector with the final trimmed subset elements used to compute these estimates and the corresponding robust squared Mahalanobis distances.

…

Further arguments to be passed to internal functions of ‘fulltle’.

Value

If argument ‘otpType’ is set to “OnlyEst”, an object of class ‘IdtE’ with the fulltle estimates, their log-likelihood values, and the value of the comparison criterion used to select the covariance configurations.

If argument ‘otpType’ is set to “SetMD2andEst” a list with the following components:

sol

An object of class ‘IdtE’ with the fulltle estimates, their log-likelihood values, and the value of the comparison criterion used to select the covariance configurations.

Set

A vector with the final trimmed subset elements used to compute the fulltle estimates.

RobMD2

A vector with the robust squared Mahalanobis distances used to select the trimmed subset.

Methods

signature(Idt = "IData")

Performs maximum trimmed likelihood estimation for interval-valued data using an exact algorithm, and assuming a Gaussian distribution, and considering alternative variance-covariance matrix configurations

References

Brito, P., Duarte Silva, A. P. (2012), Modelling Interval Data with Normal and Skew-Normal Distributions. Journal of Applied Statistics 39(1), 3--20.

Hadi, A. S. and Luceno, A. (1997), Maximum trimmed likelihood estimators: a unified approach, examples, and algorithms. Computational Statistics and Data Analysis 25(3), 251--272.

See Also

fasttle

Examples

# Create an Interval-Data object containing the intervals for characteristics 
# of 27 cars models.

Cars <- IData(CarsData[1:8],Seq="MidPLogR_VarbyVar",
  VarNames=c("Price","EngineCapacity","TopSpeed","Acceleration"))

#Display the first and last observations

head(Cars)
tail(Cars)


# Estimate normal distributuion  parameters robustly by full maximum trimmed likelihood, 

CarsTE <- fulltle(Cars)
cat("Cars data -- normal maximum trimmed likelhiood estimation results:\n")
print(CarsTE)