[Distr-commits] r1346 - in branches/distr-2.9/pkg/distrMod: man tests/Examples
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Wed May 1 14:21:36 CEST 2019
Author: ruckdeschel
Date: 2019-05-01 14:21:36 +0200 (Wed, 01 May 2019)
New Revision: 1346
Modified:
branches/distr-2.9/pkg/distrMod/man/0distrMod-package.Rd
branches/distr-2.9/pkg/distrMod/man/BetaFamily.Rd
branches/distr-2.9/pkg/distrMod/man/ExpScaleFamily.Rd
branches/distr-2.9/pkg/distrMod/man/GammaFamily.Rd
branches/distr-2.9/pkg/distrMod/man/InfoNorm.Rd
branches/distr-2.9/pkg/distrMod/man/MCEstimator.Rd
branches/distr-2.9/pkg/distrMod/man/NormType-class.Rd
branches/distr-2.9/pkg/distrMod/man/NormType.Rd
branches/distr-2.9/pkg/distrMod/man/QFNorm.Rd
branches/distr-2.9/pkg/distrMod/man/SelfNorm.Rd
branches/distr-2.9/pkg/distrMod/man/returnlevelplot.Rd
branches/distr-2.9/pkg/distrMod/tests/Examples/distrMod-Ex.Rout.save
Log:
[distrMod] merged back changes from trunk to branch distr-2.9
Modified: branches/distr-2.9/pkg/distrMod/man/0distrMod-package.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/0distrMod-package.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/0distrMod-package.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -2,7 +2,6 @@
\alias{distrMod-package}
\alias{distrMod}
\docType{package}
-
\title{
distrMod -- Object Oriented Implementation of Probability Models
}
@@ -11,14 +10,12 @@
provides a flexible framework which allows computation of estimators like
maximum likelihood or minimum distance estimators for probability models.
}
-
\details{
\tabular{ll}{
Package: \tab distrMod \cr
-Version: \tab 2.9.0 \cr
-Date: \tab 2019-03-13 \cr
-Depends: \tab R(>= 3.4), distr(>= 2.8.0), distrEx(>= 2.8.0), RandVar(>= 1.2.0),
- MASS, stats4,methods \cr
+Version: \tab 2.8.2 \cr
+Date: \tab 2019-05-01 \cr
+Depends: \tab R(>= 3.4), distr(>= 2.8.0), distrEx(>= 2.8.0), RandVar(>= 1.2.0), MASS, stats4,methods \cr
Imports: \tab startupmsg, sfsmisc, graphics, stats, grDevices \cr
Suggests: \tab ismev, evd, \cr
Enhances: RobExtremes\cr
@@ -25,13 +22,11 @@
ByteCompile: \tab yes \cr
License: \tab LGPL-3 \cr
URL: \tab http://distr.r-forge.r-project.org/\cr
-VCS/SVNRevision: \tab 1324 \cr
+VCS/SVNRevision: \tab 1344 \cr
}}
-
\section{Classes}{
\preformatted{
[*]: there is a generating function with the same name
-
##########################
ProbFamily classes
##########################
@@ -38,7 +33,6 @@
slots: [<name>(<class>)]
name(character), distribution(Distribution),
distrSym(DistributionSymmetry), props(character)
-
"ProbFamily"
|>"ParamFamily" [*]
additional slots:
@@ -74,10 +68,8 @@
|>|>|>|>|>|>"NormLocationScaleFamily" [*]
|>|>|>|>|>|>"CauchyLocationScaleFamily" [*]
|>|>|>|>|>|>"LogisticLocationScaleFamily" [*]
-
and a (virtual) class union "L2ScaleUnion" between
"L2LocationScaleUnion" and "L2ScaleShapeUnion"
-
##########################
ParamFamParameter
##########################
@@ -85,7 +77,6 @@
Additional slots:
main(numeric), nuisance(OptionalNumeric), fixed(OptionalNumeric),
trafo(MatrixorFunction)
-
##########################
Class unions
##########################
@@ -93,24 +84,18 @@
"PrintDetails" = union("Estimate", "Confint",
"PosSemDefSymmMatrix",
"ParamFamParameter", "ParamFamily")
-
-
##########################
Symmetry classes (other classes moved to package "distr")
##########################
slots:
type(character), SymmCenter(ANY)
-
"Symmetry" (from package "distr")
|>"FunctionSymmetry"
|>|>"NonSymmetric" [*]
|>|>"EvenSymmetric" [*]
|>|>"OddSymmetric" [*]
-
list thereof
"FunSymmList" [*]
-
-
##########################
Matrix classes (moved to package "distr")
##########################
@@ -118,14 +103,11 @@
none
"PosSemDefSymmMatrix" [*] is subclass of class "matrix" of package "base".
|>"PosDefSymmMatrix" [*]
-
-
##########################
Norm Classes
##########################
slots:
name(character), fct(function)
-
"NormType" [*]
|>"QFNorm" [*]
Additional slots:
@@ -132,14 +114,11 @@
QuadForm(PosSemDefSymmMatrix)
|>|>"InfoNorm" [*]
|>|>"SelfNorm" [*]
-
-
##########################
Bias Classes
##########################
slots:
name(character)
-
"BiasType"
|>"symmetricBias" [*]
|>"onesidedBias"
@@ -148,14 +127,11 @@
|>"asymmetricBias" [*]
Additional slots:
nu(numeric)
-
-
##########################
Risk Classes
##########################
slots:
type(character)
-
"RiskType"
|>"asRisk"
|>|>"asCov" [*]
@@ -171,7 +147,6 @@
|>|>"fiUnOvShoot" [*]
Additional slots:
width(numeric)
-
Risk with Bias:
"asRiskwithBias"
slots: biastype(BiasType), normtype(NormType),
@@ -185,8 +160,6 @@
Additional slots:
width(numeric)
|>|>"asSemivar" [*]
-
-
##########################
Estimate Classes
##########################
@@ -199,12 +172,10 @@
untransformed.estimate(ANY),
untransformed.asvar(OptionalMatrix)
criterion.fct(function), method(character),
-
"Estimate"
|>"MCEstimate",
Additional slots:
criterion(numeric)
-
##########################
Confidence interval class
##########################
@@ -214,20 +185,16 @@
trafo.estimate(list[of function, matrix]),
nuisance.estimate(OptionalNumeric)
"Confint"
-
}
}
\section{Methods}{
besides accessor and replacement functions, we have
methods
-
\code{solve}, \code{sqrt} for matrices
\code{checkL2deriv}, \code{existsPIC} for class \code{L2ParamFamily}
\code{LogDeriv} for class \code{L2GroupParamFamily}
-
\code{validParameter} for classes \code{ParamFamily}, \code{L2ScaleFamily},
\code{L2LocationFamily}, and \code{L2LocationScaleFamily}
-
\code{modifyModel} for the pairs of classes
\code{L2ParamFamily} and \code{ParamFamParameter},
\code{L2LocationFamily} and \code{ParamFamParameter},
@@ -235,9 +202,7 @@
\code{L2LocationScaleFamily} and \code{ParamFamParameter},
\code{GammaFamily} and \code{ParamFamParameter}, and
\code{ExpScaleFamily} and \code{ParamFamParameter}
-
\code{mceCalc} for the pair of classes \code{numeric} and \code{ParamFamily}
-
\code{mleCalc} for the pairs of classes
\code{numeric} and \code{ParamFamily},
\code{numeric} and \code{BinomFamily},
@@ -245,38 +210,27 @@
\code{numeric} and \code{NormLocationFamily},
\code{numeric} and \code{NormScaleFamily}, and
\code{numeric} and \code{NormLocationScaleFamily}
-
\code{coerce} from class \code{MCEstimate} to class \code{mle}
-
\code{confint} for class \code{Estimate}
\code{profile} for class \code{MCEstimate}
-
}
\section{Functions}{
-
\preformatted{
Management of global options:
"distrModOptions", "distrModoptions", "getdistrModOption",
-
check for ker of matrix: "isKerAinKerB"
-
particular norms: "EuclideanNorm", "QuadFormNorm"
onesided bias: "positiveBias", "negativeBias",
-
Estimators:
"Estimator", "MCEstimator", "MLEstimator", "MDEstimator"
-
special location/scale models:
"L2LocationUnknownScaleFamily", "L2ScaleUnknownLocationFamily"
-
some special normal models:
"NormScaleUnknownLocationFamily", "NormLocationUnknownScaleFamily",
}}
-
\note{
Some functions of packages \pkg{stats}, \pkg{base} have intentionally been masked,
but completely retain their functionality --- see \code{distrModMASK()}.
-
If any of the packages \pkg{stats4}, \pkg{fBasics} is to be used
together with \pkg{distrMod}, the latter must be attached \emph{after} any of the
first mentioned. Otherwise \code{confint()} defined as \emph{method}
@@ -284,13 +238,11 @@
may use \code{confint <- distrMod::confint}.
See also \code{distrModMASK()}
}
-
\author{
Peter Ruckdeschel \email{peter.ruckdeschel at uni-oldenburg.de},\cr
Matthias Kohl \email{Matthias.Kohl at stamats.de}\cr
\emph{Maintainer:} Peter Ruckdeschel \email{peter.ruckdeschel at uni-oldenburg.de}
}
-
\references{
M. Kohl and P. Ruckdeschel (2010):
R Package distrMod: S4 Classes and Methods for Probability Models.
@@ -297,42 +249,33 @@
Journal of Statistical Software, 35(10), 1-27.
\url{https://www.jstatsoft.org/v35/i10/}
(see also \code{vignette("distrMod")})
-
-
P. Ruckdeschel, M. Kohl, T. Stabla, F. Camphausen (2006):
S4 Classes for Distributions, \emph{R News}, \emph{6}(2), 2-6.
\url{https://CRAN.R-project.org/doc/Rnews/Rnews_2006-2.pdf}
-
A vignette for packages \pkg{distr}, \pkg{distrSim}, \pkg{distrTEst}, and
\pkg{distrEx} is included into the mere documentation package \pkg{distrDoc}
and may be called by
\code{require("distrDoc");vignette("distr")}
}
-
\section{Start-up-Banner}{
You may suppress the start-up banner/message completely by setting
\code{options("StartupBanner"="off")} somewhere before loading this
package by \code{library} or \code{require} in your R-code / R-session.
-
If option \code{"StartupBanner"} is not defined (default) or setting
\code{options("StartupBanner"=NULL)} or
\code{options("StartupBanner"="complete")} the complete start-up banner is
displayed.
-
For any other value of option \code{"StartupBanner"}
(i.e., not in \code{c(NULL,"off","complete")})
only the version information is displayed.
-
The same can be achieved by wrapping the \code{library} or \code{require}
call into either \code{suppressStartupMessages()} or
\code{onlytypeStartupMessages(.,atypes="version")}.
-
As for general \code{packageStartupMessage}'s, you may also suppress all
the start-up banner by wrapping the \code{library} or \code{require}
call into \code{suppressPackageStartupMessages()} from
\pkg{startupmsg}-version 0.5 on.
}
-
\section{Demos}{
Demos are available --- see \code{demo(package="distrMod")}.}
\section{Scripts}{
@@ -339,7 +282,6 @@
Example scripts are available --- see folder \file{scripts}
in the package folder to package \pkg{distrMod} in your library.
}
-
\section{Package versions}{
Note: The first two numbers of package versions do not necessarily reflect
package-individual development, but rather are chosen for the
Modified: branches/distr-2.9/pkg/distrMod/man/BetaFamily.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/BetaFamily.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/BetaFamily.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -28,7 +28,9 @@
\examples{
(B1 <- BetaFamily())
FisherInfo(B1)
+## IGNORE_RDIFF_BEGIN
checkL2deriv(B1)
+## IGNORE_RDIFF_END
}
\concept{Beta model}
\keyword{models}
Modified: branches/distr-2.9/pkg/distrMod/man/ExpScaleFamily.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/ExpScaleFamily.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/ExpScaleFamily.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -30,7 +30,9 @@
(E1 <- ExpScaleFamily())
plot(E1)
Map(L2deriv(E1)[[1]])
+## IGNORE_RDIFF_BEGIN
checkL2deriv(E1)
+## IGNORE_RDIFF_END
}
\concept{exponential scale model}
\concept{scale model}
Modified: branches/distr-2.9/pkg/distrMod/man/GammaFamily.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/GammaFamily.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/GammaFamily.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -31,7 +31,9 @@
\examples{
(G1 <- GammaFamily())
FisherInfo(G1)
+## IGNORE_RDIFF_BEGIN
checkL2deriv(G1)
+## IGNORE_RDIFF_END
}
\concept{Gamma model}
\keyword{models}
Modified: branches/distr-2.9/pkg/distrMod/man/InfoNorm.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/InfoNorm.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/InfoNorm.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -19,12 +19,12 @@
%\note{}
\seealso{\code{\link{InfoNorm-class}}}
\examples{
- ## IGNORE_RDIFF_BEGIN
+## IGNORE_RDIFF_BEGIN
InfoNorm()
## The function is currently defined as
function(){ new("InfoNorm") }
- ## IGNORE_RDIFF_END
+## IGNORE_RDIFF_END
}
\concept{asymptotic covariance}
\concept{risk}
Modified: branches/distr-2.9/pkg/distrMod/man/MCEstimator.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/MCEstimator.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/MCEstimator.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -97,8 +97,11 @@
## Total variation minimum distance estimator
## Note: you can also use function MDEstimator!
## discretize Gamma distribution
-MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
+
+## IGNORE_RDIFF_BEGIN
+MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
crit.name = "Total variation distance")
+## IGNORE_RDIFF_END
## or smooth empirical distribution (takes some time!)
#MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
Modified: branches/distr-2.9/pkg/distrMod/man/NormType-class.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/NormType-class.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/NormType-class.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -40,11 +40,11 @@
%\note{}
\seealso{\code{\link{BiasType-class}}}
\examples{
- ## IGNORE_RDIFF_BEGIN
+## IGNORE_RDIFF_BEGIN
EuclNorm <- NormType("EuclideanNorm",EuclideanNorm)
fct(EuclNorm)
name(EuclNorm)
- ## IGNORE_RDIFF_END
+## IGNORE_RDIFF_END
}
\concept{asymptotic risk}
\concept{bias}
Modified: branches/distr-2.9/pkg/distrMod/man/NormType.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/NormType.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/NormType.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -22,9 +22,9 @@
%\note{}
\seealso{\code{\link{NormType-class}}}
\examples{
- ## IGNORE_RDIFF_BEGIN
+## IGNORE_RDIFF_BEGIN
NormType()
- ## IGNORE_RDIFF_END
+## IGNORE_RDIFF_END
}
\concept{asymptotic covariance}
\concept{risk}
Modified: branches/distr-2.9/pkg/distrMod/man/QFNorm.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/QFNorm.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/QFNorm.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -24,12 +24,12 @@
%\note{}
\seealso{\code{\link{QFNorm-class}}}
\examples{
+## IGNORE_RDIFF_BEGIN
QFNorm()
## The function is currently defined as
- ## IGNORE_RDIFF_BEGIN
function(){ new("QFNorm") }
- ## IGNORE_RDIFF_END
+## IGNORE_RDIFF_END
}
\concept{asymptotic covariance}
\concept{risk}
Modified: branches/distr-2.9/pkg/distrMod/man/SelfNorm.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/SelfNorm.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/SelfNorm.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -19,12 +19,12 @@
%\note{}
\seealso{\code{\link{SelfNorm-class}}}
\examples{
- ## IGNORE_RDIFF_BEGIN
+## IGNORE_RDIFF_BEGIN
SelfNorm()
## The function is currently defined as
function(){ new("SelfNorm") }
- ## IGNORE_RDIFF_END
+## IGNORE_RDIFF_END
}
\concept{asymptotic covariance}
\concept{risk}
Modified: branches/distr-2.9/pkg/distrMod/man/returnlevelplot.Rd
===================================================================
--- branches/distr-2.9/pkg/distrMod/man/returnlevelplot.Rd 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/man/returnlevelplot.Rd 2019-05-01 12:21:36 UTC (rev 1346)
@@ -255,9 +255,13 @@
returnlevelplot(r(Norm(15,sqrt(30)))(40), Chisq(df=15))
### more could be seen after installing RobExtremes and ismev
#
- ## IGNORE_RDIFF_BEGIN
- if(require(RobExtremes) && require(ismev)){
+## IGNORE_RDIFF_BEGIN
+\donttest{ ## at R CMD check --as-cran, it does not find package cluster
+ ## when trying to attach package rrcov
+ ## so remove this from testing
+if(require(RobExtremes) && require(ismev)){
+
data(portpirie)
gevfit <- gev.fit(portpirie[,2]) ## taken from example from ismev::gev.fit
GEVF <- GEVFamily(scale=gevfit$mle[2],shape=gevfit$mle[3],loc=gevfit$mle[1])
@@ -271,7 +275,9 @@
gpdfit <- gpd.fit(rain,10) ## taken from example from ismev::gpd.fit
GPDF <- GParetoFamily(scale=gpdfit$mle[1],shape=gpdfit$mle[2],loc=10)
returnlevelplot(rain, GPDF, MaxOrPOT="POT", xlim=c(1e-1,1e3))
-}}
- ## IGNORE_RDIFF_END
+}
+}
+## IGNORE_RDIFF_END
+}
\keyword{hplot}
\keyword{distribution}
Modified: branches/distr-2.9/pkg/distrMod/tests/Examples/distrMod-Ex.Rout.save
===================================================================
--- branches/distr-2.9/pkg/distrMod/tests/Examples/distrMod-Ex.Rout.save 2019-05-01 12:18:38 UTC (rev 1345)
+++ branches/distr-2.9/pkg/distrMod/tests/Examples/distrMod-Ex.Rout.save 2019-05-01 12:21:36 UTC (rev 1346)
@@ -1,7 +1,7 @@
R Under development (unstable) (2019-03-28 r76288) -- "Unsuffered Consequences"
Copyright (C) 2019 The R Foundation for Statistical Computing
-Platform: x86_64-w64-mingw32/x64 (64-bit)
+Platform: i386-w64-mingw32/i386 (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
@@ -92,7 +92,7 @@
Loading required package: MASS
Loading required package: stats4
:distrMod> Object Oriented Implementation of Probability Models
-:distrMod> (version 2.8.1)
+:distrMod> (version 2.8.2)
:distrMod>
:distrMod> Some functions from pkg's 'base' and 'stats' are
:distrMod> intentionally masked ---see distrModMASK().
@@ -167,8 +167,9 @@
shape1 shape2
shape1 1.0000000 -0.6449341
shape2 -0.6449341 1.0000000
+> ## IGNORE_RDIFF_BEGIN
> checkL2deriv(B1)
-precision of centering: 3.96327e-05 3.963591e-05
+precision of centering: 3.963281e-05 3.963334e-05
precision of Fisher information:
shape1 shape2
shape1 -1.8511e-05 1.6483e-06
@@ -179,6 +180,7 @@
[2,] -2.56e-04 -1.85e-03
condition of Fisher information:
[1] 5.2777
+> ## IGNORE_RDIFF_END
>
>
>
@@ -419,7 +421,7 @@
dimnames = list(nms, nms0))
list(fval = fval0, mat = mat0)
}
-<bytecode: 0x000000000d34f360>
+<bytecode: 0x0785b5b8>
Trafo / derivative matrix at which estimate was produced:
scale shape
shape 0.000 1
@@ -620,10 +622,11 @@
1)/c(scale = 1)
return(y)
}
-<environment: 0x000000000e638530>
+<environment: 0x0cd99b98>
+> ## IGNORE_RDIFF_BEGIN
> checkL2deriv(E1)
-precision of centering: -2.042661e-06
+precision of centering: -2.04266e-06
precision of Fisher information:
scale
scale -3.5986e-05
@@ -632,6 +635,7 @@
[1,] -3.6e-03
condition of Fisher information:
[1] 1
+> ## IGNORE_RDIFF_END
>
>
>
@@ -749,8 +753,9 @@
scale shape
scale 1 1.000000
shape 1 1.644934
+> ## IGNORE_RDIFF_BEGIN
> checkL2deriv(G1)
-precision of centering: -2.042661e-06 1.791171e-06
+precision of centering: -2.04266e-06 1.791171e-06
precision of Fisher information:
scale shape
scale -3.5986e-05 -9.5036e-06
@@ -761,6 +766,7 @@
[2,] -9.5e-04 -2.4e-03
condition of Fisher information:
[1] 10.603
+> ## IGNORE_RDIFF_END
>
>
>
@@ -777,6 +783,7 @@
>
> ### ** Examples
>
+> ## IGNORE_RDIFF_BEGIN
> InfoNorm()
An object of class "InfoNorm"
Slot "QuadForm":
@@ -789,8 +796,8 @@
Slot "fct":
function(x) QuadFormNorm(x, A = A)
-<bytecode: 0x000000000d5f1b98>
-<environment: 0x000000000d5f5498>
+<bytecode: 0x0835e4e0>
+<environment: 0x0835d0f0>
>
> ## The function is currently defined as
@@ -799,6 +806,7 @@
{
new("InfoNorm")
}
+> ## IGNORE_RDIFF_END
>
>
>
@@ -1072,7 +1080,7 @@
1)/c(meanlog = 1)
return(y)
}
-<environment: 0x00000000065c9f70>
+<environment: 0x0d93e088>
> checkL2deriv(L1)
precision of centering: -0.003003394
@@ -1307,7 +1315,9 @@
> ## Total variation minimum distance estimator
> ## Note: you can also use function MDEstimator!
> ## discretize Gamma distribution
-> MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
+>
+> ## IGNORE_RDIFF_BEGIN
+> MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
+ crit.name = "Total variation distance")
Evaluations of Minimum Total variation distance estimate:
---------------------------------------------------------
@@ -1318,10 +1328,11 @@
samplesize: 50
estimate:
scale shape
-0.2829687 5.0197306
+0.2829625 5.0197311
Criterion:
Total variation distance
- 0.4866141
+ 0.4819868
+> ## IGNORE_RDIFF_END
>
> ## or smooth empirical distribution (takes some time!)
> #MCEstimator(x = x, ParamFamily = G, criterion = TotalVarDist,
@@ -2103,6 +2114,7 @@
>
> ### ** Examples
>
+> ## IGNORE_RDIFF_BEGIN
> EuclNorm <- NormType("EuclideanNorm",EuclideanNorm)
> fct(EuclNorm)
function(x){
@@ -2111,10 +2123,11 @@
else
return(sqrt(colSums(x^2)))
}
-<bytecode: 0x000000000d7b0298>
+<bytecode: 0x08375ea0>
<environment: namespace:distrMod>
> name(EuclNorm)
[1] "EuclideanNorm"
+> ## IGNORE_RDIFF_END
>
>
>
@@ -2131,6 +2144,7 @@
>
> ### ** Examples
>
+> ## IGNORE_RDIFF_BEGIN
> NormType()
An object of class "NormType"
Slot "name":
@@ -2143,9 +2157,10 @@
else
return(sqrt(colSums(x^2)))
}
-<bytecode: 0x000000000d7b0298>
+<bytecode: 0x08375ea0>
<environment: namespace:distrMod>
+> ## IGNORE_RDIFF_END
>
>
>
@@ -2586,6 +2601,7 @@
>
> ### ** Examples
>
+> ## IGNORE_RDIFF_BEGIN
> QFNorm()
An object of class "QFNorm"
Slot "QuadForm":
@@ -2598,8 +2614,8 @@
Slot "fct":
function(x) QuadFormNorm(x, A= A0)
-<bytecode: 0x0000000014108388>
-<environment: 0x0000000014109fd0>
+<bytecode: 0x0adda7f0>
+<environment: 0x0af18d28>
>
> ## The function is currently defined as
@@ -2608,6 +2624,7 @@
{
new("QFNorm")
}
+> ## IGNORE_RDIFF_END
>
>
>
@@ -2624,6 +2641,7 @@
>
> ### ** Examples
>
+> ## IGNORE_RDIFF_BEGIN
> SelfNorm()
An object of class "SelfNorm"
Slot "QuadForm":
@@ -2636,8 +2654,8 @@
Slot "fct":
function(x) QuadFormNorm(x, A = A)
-<bytecode: 0x00000000138dd9d0>
-<environment: 0x00000000138e33d8>
+<bytecode: 0x0b06e640>
+<environment: 0x0c6661f8>
>
> ## The function is currently defined as
@@ -2646,6 +2664,7 @@
{
new("SelfNorm")
}
+> ## IGNORE_RDIFF_END
>
>
>
@@ -2663,7 +2682,7 @@
> ### ** Examples
>
> addAlphTrsp2col(rgb(1,0.3,0.03), 25)
-[1] "#FF4D0819"
+[1] "#FF4C0819"
> addAlphTrsp2col("darkblue", 25)
[1] "#00008B19"
> addAlphTrsp2col("#AAAAAAAA",25)
@@ -3655,7 +3674,7 @@
dimnames(mat) <- list(nfval, c("mean", "sd"))
return(list(fval = fval, mat = mat))
}
-<bytecode: 0x0000000014f4f340>
+<bytecode: 0x0b7f0dc8>
> print(param(NS), show.details = "minimal")
An object of class "ParamWithScaleFamParameter"
name: location and scale
@@ -3704,7 +3723,7 @@
dimnames(mat) <- list(nfval, c("mean", "sd"))
return(list(fval = fval, mat = mat))
}
-<bytecode: 0x0000000014f4f340>
+<bytecode: 0x0b7f0dc8>
Trafo / derivative matrix:
mean sd
mu/sig 0.3668695 -0.3024814
@@ -3747,7 +3766,7 @@
dimnames(mat) <- list(nfval, c("mean", "sd"))
return(list(fval = fval, mat = mat))
}
-<bytecode: 0x0000000014f4f340>
+<bytecode: 0x0b7f0dc8>
Trafo / derivative matrix:
mean sd
mu/sig 0.3669 -0.3025
@@ -3827,132 +3846,13 @@
> returnlevelplot(r(Norm(15,sqrt(30)))(40), Chisq(df=15))
> ### more could be seen after installing RobExtremes and ismev
> #
-> if(require(RobExtremes) && require(ismev)){
-+
-+ data(portpirie)
-+ gevfit <- gev.fit(portpirie[,2]) ## taken from example from ismev::gev.fit
-+ GEVF <- GEVFamily(scale=gevfit$mle[2],shape=gevfit$mle[3],loc=gevfit$mle[1])
-+ erg <- returnlevelplot(portpirie[,2], GEVF)
-+ print(names(erg))
-+ print(names(erg$plotArgs))
-+ print(names(erg$IdLineArgs))
-+ returnlevelplot(portpirie[,2], GEVF, datax=TRUE)
-+
-+ data(rain)
-+ gpdfit <- gpd.fit(rain,10) ## taken from example from ismev::gpd.fit
-+ GPDF <- GParetoFamily(scale=gpdfit$mle[1],shape=gpdfit$mle[2],loc=10)
-+ returnlevelplot(rain, GPDF, MaxOrPOT="POT", xlim=c(1e-1,1e3))
-+ }
-Loading required package: RobExtremes
-Loading required package: ROptEst
-Loading required package: RobAStBase
-Loading required package: rrcov
-Loading required package: robustbase
-Scalable Robust Estimators with High Breakdown Point (version 1.4-7)
-
-:RobAStBase> Robust Asymptotic Statistics (version 1.2.0)
-:RobAStBase>
-:RobAStBase> Some functions from pkg's 'stats' and 'graphics'
-:RobAStBase> are intentionally masked ---see RobAStBaseMASK().
-:RobAStBase>
-:RobAStBase> Note that global options are controlled by
-:RobAStBase> RobAStBaseoptions() ---c.f. ?"RobAStBaseoptions".
-:RobAStBase>
-:RobAStBase> For more information see ?"RobAStBase",
-:RobAStBase> NEWS("RobAStBase"), as well as
-:RobAStBase> http://robast.r-forge.r-project.org/
-
-
-Attaching package: 'RobAStBase'
-
-The following object is masked from 'package:graphics':
-
- clip
-
-Loading required package: evd
-:RobExtremes> Optimally Robust Estimation for Extreme Value
-:RobExtremes> Distributions (version 1.2.0)
-:RobExtremes>
-:RobExtremes>
-:RobExtremes> For more information see ?"RobExtremes",
-:RobExtremes> NEWS("RobExtremes"), as well as
-:RobExtremes> http://robast.r-forge.r-project.org/
-
-
-Attaching package: 'RobExtremes'
-
-The following objects are masked from 'package:robustbase':
-
- Qn, Sn
-
-Loading required package: ismev
-Loading required package: mgcv
-Loading required package: nlme
-
-Attaching package: 'nlme'
-
-The following object is masked from 'package:rrcov':
-
- getData
-
-The following object is masked from 'package:RandVar':
-
- Dim
-
-This is mgcv 1.8-28. For overview type 'help("mgcv-package")'.
-
-Attaching package: 'ismev'
-
-The following objects are masked from 'package:RobExtremes':
-
- gev.diag, gev.prof, gev.profxi, gpd.diag, gpd.prof, gpd.profxi
-
-$conv
-[1] 0
-
-$nllh
-[1] -4.339058
-
-$mle
-[1] 3.87474692 0.19804120 -0.05008773
-
-$se
-[1] 0.02793211 0.02024610 0.09825633
-
- [1] "call" "dots" "args" "plotArgs" "pointArgs"
- [6] "usr" "IdLineArgs" "ret.x" "ret.y" "qqb.crit"
-[11] "qqb.err"
- [1] "x" "y" "log" "type" "xlab" "ylab" "main" "pch" "cex" "col"
-[1] "" "" "col" "lty" "lwd"
-$threshold
-[1] 10
-
-$nexc
-[1] 2003
-
-$conv
-[1] 0
-
-$nllh
-[1] 6123.465
-
-$mle
-[1] 7.43768624 0.05045225
-
-$rate
-[1] 0.1142547
-
-$se
-[1] 0.23606472 0.02256649
-
>
+> ## IGNORE_RDIFF_BEGIN
+> ## IGNORE_RDIFF_END
>
+>
+>
> cleanEx()
-
-detaching 'package:ismev', 'package:mgcv', 'package:nlme',
- 'package:RobExtremes', 'package:evd', 'package:ROptEst',
- 'package:RobAStBase', 'package:rrcov', 'package:robustbase'
-
> nameEx("symmetricBias-class")
> ### * symmetricBias-class
>
@@ -4144,7 +4044,7 @@
(0.2113366) (1.9183645)
> ## confidence interval
> confint(res)
-A[n] asymptotic (LAN-based) confidence interval:
+A[n] asymptotic (CLT-based) confidence interval:
2.5 % 97.5 %
mu/sig 0.4102811 1.238705
sig^2 3.6698685 11.189719
@@ -4189,7 +4089,7 @@
> res <- trafoEst(mtrafo, res0)
> ## confidence interval
> confint(res)
-A[n] asymptotic (LAN-based) confidence interval:
+A[n] asymptotic (CLT-based) confidence interval:
2.5 % 97.5 %
mu/sig 0.4102811 1.238705
sig^2 3.6698685 11.189719
@@ -4232,7 +4132,7 @@
> cleanEx()
> options(digits = 7L)
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
-Time elapsed: 29.69 0.99 31.06 NA NA
+Time elapsed: 29.3 0.58 30 NA NA
> grDevices::dev.off()
null device
1
More information about the Distr-commits
mailing list