[Robast-commits] r210 - branches/robast-0.6/pkg/ROptEst/chm branches/robast-0.6/pkg/RobAStBase/chm branches/robast-0.7/pkg/RobAStBase/man pkg/ROptEst/chm pkg/ROptRegTS pkg/RobAStBase pkg/RobAStBase/R pkg/RobAStBase/chm pkg/RobAStBase/man pkg/RobLox pkg/RobRex
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Wed Nov 12 22:33:46 CET 2008
Author: ruckdeschel
Date: 2008-11-12 22:33:45 +0100 (Wed, 12 Nov 2008)
New Revision: 210
Added:
branches/robast-0.6/pkg/ROptEst/chm/getModifyIC.html
branches/robast-0.6/pkg/ROptEst/chm/roptest.html
branches/robast-0.6/pkg/RobAStBase/chm/ALEstimate-class.html
branches/robast-0.6/pkg/RobAStBase/chm/MEstimate-class.html
branches/robast-0.6/pkg/RobAStBase/chm/RobAStBaseOptions.html
branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimate-class.html
branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimator.html
branches/robast-0.7/pkg/RobAStBase/man/plot-methods.Rd
pkg/RobAStBase/R/cutoff-class.R
pkg/RobAStBase/R/ddPlot.R
pkg/RobAStBase/R/ddPlot_utils.R
pkg/RobAStBase/R/outlyingPlot.R
pkg/RobAStBase/chm/0RobAStBase-package.html
pkg/RobAStBase/chm/ALEstimate-class.html
pkg/RobAStBase/chm/MEstimate-class.html
pkg/RobAStBase/chm/RobAStBaseOptions.html
pkg/RobAStBase/chm/cutoff-class.html
pkg/RobAStBase/chm/cutoff.html
pkg/RobAStBase/chm/ddPlot-methods.html
pkg/RobAStBase/chm/internals_ddPlot.html
pkg/RobAStBase/chm/kStepEstimate-class.html
pkg/RobAStBase/chm/kStepEstimator.html
pkg/RobAStBase/chm/outlyingPlotIC.html
pkg/RobAStBase/chm/plot-methods.html
pkg/RobAStBase/man/cutoff-class.Rd
pkg/RobAStBase/man/cutoff.Rd
pkg/RobAStBase/man/ddPlot-methods.Rd
pkg/RobAStBase/man/internals_ddPlot.Rd
pkg/RobAStBase/man/outlyingPlotIC.Rd
pkg/RobAStBase/man/plot-methods.Rd
Modified:
pkg/ROptEst/chm/00Index.html
pkg/ROptRegTS/DESCRIPTION
pkg/RobAStBase/NAMESPACE
pkg/RobAStBase/R/AllClass.R
pkg/RobAStBase/R/AllGeneric.R
pkg/RobAStBase/R/bALEstimate.R
pkg/RobAStBase/chm/00Index.html
pkg/RobAStBase/chm/IC-class.html
pkg/RobAStBase/chm/RobAStBase.chm
pkg/RobAStBase/chm/RobAStBase.hhp
pkg/RobAStBase/chm/RobAStBase.toc
pkg/RobAStBase/man/IC-class.Rd
pkg/RobLox/DESCRIPTION
pkg/RobRex/DESCRIPTION
Log:
corrected a minor bug in confint (for asym. Bias)
included new plotting functionality : ddPlot and outlyingPlot
Added: branches/robast-0.6/pkg/ROptEst/chm/getModifyIC.html
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/getModifyIC.html (rev 0)
+++ branches/robast-0.6/pkg/ROptEst/chm/getModifyIC.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,118 @@
+<html><head><title>Generic Function for the Computation of Functions for Slot modifyIC</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>getModifyIC(ROptEst)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: getModifyIC">
+<param name="keyword" value="R: getModifyIC-methods">
+<param name="keyword" value="R: getModifyIC,L2ParamFamily,Neighborhood,asRisk-method">
+<param name="keyword" value="R: getModifyIC,L2LocationFamily,UncondNeighborhood,asGRisk-method">
+<param name="keyword" value="R: getModifyIC,L2LocationFamily,UncondNeighborhood,fiUnOvShoot-method">
+<param name="keyword" value="R: getModifyIC,L2ScaleFamily,ContNeighborhood,asGRisk-method">
+<param name="keyword" value="R: getModifyIC,L2ScaleFamily,TotalVarNeighborhood,asGRisk-method">
+<param name="keyword" value="R: getModifyIC,L2LocationScaleFamily,ContNeighborhood,asGRisk-method">
+<param name="keyword" value=" Generic Function for the Computation of Functions for Slot modifyIC">
+</object>
+
+
+<h2>Generic Function for the Computation of Functions for Slot modifyIC</h2>
+
+
+<h3>Description</h3>
+
+<p>
+This function is used by internal computations and is rarely called directly.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2ParamFamily, Neighborhood,
+## asRisk':
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2LocationFamily,
+## UncondNeighborhood, asGRisk':
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2LocationFamily,
+## UncondNeighborhood, fiUnOvShoot':
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2ScaleFamily,
+## ContNeighborhood, asGRisk':
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2ScaleFamily,
+## TotalVarNeighborhood, asGRisk':
+getModifyIC(L2FamIC, neighbor, risk)
+## S4 method for signature 'L2LocationScaleFamily,
+## ContNeighborhood, asGRisk':
+getModifyIC(L2FamIC, neighbor, risk)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>L2FamIC</code></td>
+<td>
+object of class <code>L2ParamFamily</code>. </td></tr>
+<tr valign="top"><td><code>neighbor</code></td>
+<td>
+object of class <code>"Neighborhood"</code>. </td></tr>
+<tr valign="top"><td><code>risk</code></td>
+<td>
+object of class <code>"RiskType"</code> </td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+This function is used for internal computations.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+Function for slot <code>modifyIC</code> of <code>IC</code>s</p>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>References</h3>
+
+<p>
+Rieder, H. (1994) <EM>Robust Asymptotic Statistics</EM>. New York: Springer.
+</p>
+<p>
+Kohl, M. (2005) <EM>Numerical Contributions to the Asymptotic Theory of Robustness</EM>.
+Bayreuth: Dissertation.
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="optIC.html">optIC</a></code>, <code><a onclick="findlink('RobAStBase', 'IC-class.html')" style="text-decoration: underline; color: blue; cursor: hand">IC-class</a></code>
+</p>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>ROptEst</em> version 0.6.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/ROptEst/chm/roptest.html
===================================================================
--- branches/robast-0.6/pkg/ROptEst/chm/roptest.html (rev 0)
+++ branches/robast-0.6/pkg/ROptEst/chm/roptest.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,254 @@
+<html><head><title>Optimally robust estimation</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>roptest(ROptEst)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: roptest">
+<param name="keyword" value=" Optimally robust estimation">
+</object>
+
+
+<h2>Optimally robust estimation</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Function to compute optimally robust estimates for L2-differentiable
+parametric families via k-step construction.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+roptest(x, L2Fam, eps, eps.lower, eps.upper, initial.est,
+ neighbor = ContNeighborhood(), risk = asMSE(), steps = 1,
+ distance = CvMDist, startPar = NULL, verbose = FALSE,
+ useLast = getRobAStBaseOption("kStepUseLast"), ...)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>x</code></td>
+<td>
+sample </td></tr>
+<tr valign="top"><td><code>L2Fam</code></td>
+<td>
+object of class <code>"L2ParamFamily"</code> </td></tr>
+<tr valign="top"><td><code>eps</code></td>
+<td>
+positive real (0 < <code>eps</code> <= 0.5): amount of gross errors.
+See details below. </td></tr>
+<tr valign="top"><td><code>eps.lower</code></td>
+<td>
+positive real (0 <= <code>eps.lower</code> <= <code>eps.upper</code>):
+lower bound for the amount of gross errors. See details below. </td></tr>
+<tr valign="top"><td><code>eps.upper</code></td>
+<td>
+positive real (<code>eps.lower</code> <= <code>eps.upper</code> <= 0.5):
+upper bound for the amount of gross errors. See details below. </td></tr>
+<tr valign="top"><td><code>initial.est</code></td>
+<td>
+initial estimate for unknown parameter. If missing
+minimum distance estimator is computed. </td></tr>
+<tr valign="top"><td><code>neighbor</code></td>
+<td>
+object of class <code>"UncondNeighborhood"</code> </td></tr>
+<tr valign="top"><td><code>risk</code></td>
+<td>
+object of class <code>"RiskType"</code> </td></tr>
+<tr valign="top"><td><code>steps</code></td>
+<td>
+positive integer: number of steps used for k-steps construction </td></tr>
+<tr valign="top"><td><code>distance</code></td>
+<td>
+distance function </td></tr>
+<tr valign="top"><td><code>startPar</code></td>
+<td>
+initial information used by <code>optimize</code> resp. <code>optim</code>;
+i.e; if (total) parameter is of length 1, <code>startPar</code> is
+a search interval, else it is an initial parameter value; if <code>NULL</code>
+slot <code>startPar</code> of <code>ParamFamily</code> is used to produce it;
+in the multivariate case, <code>startPar</code> may also be of class <code>Estimate</code>,
+in which case slot <code>untransformed.estimate</code> is used.</td></tr>
+<tr valign="top"><td><code>verbose</code></td>
+<td>
+logical: if <code>TRUE</code>, some messages are printed </td></tr>
+<tr valign="top"><td><code>useLast</code></td>
+<td>
+which parameter estimate (initial estimate or
+k-step estimate) shall be used to fill the slots <code>pIC</code>,
+<code>asvar</code> and <code>asbias</code> of the return value. </td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+further arguments </td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+Computes the optimally robust estimator for a given L2 differentiable
+parametric family. The computation uses a k-step construction with an
+appropriate initial estimate; cf. also <code><a onclick="findlink('RobAStBase', 'kStepEstimator.html')" style="text-decoration: underline; color: blue; cursor: hand">kStepEstimator</a></code>.
+Valid candidates are e.g. Kolmogorov(-Smirnov) or von Mises minimum
+distance estimators (default); cf. Rieder (1994) and Kohl (2005).
+</p>
+<p>
+If the amount of gross errors (contamination) is known, it can be
+specified by <code>eps</code>. The radius of the corresponding infinitesimal
+contamination neighborhood is obtained by multiplying <code>eps</code>
+by the square root of the sample size.
+</p>
+<p>
+If the amount of gross errors (contamination) is unknown, try to find a
+rough estimate for the amount of gross errors, such that it lies
+between <code>eps.lower</code> and <code>eps.upper</code>.
+</p>
+<p>
+In case <code>eps.lower</code> is specified and <code>eps.upper</code> is missing,
+<code>eps.upper</code> is set to 0.5. In case <code>eps.upper</code> is specified and
+<code>eps.lower</code> is missing, <code>eps.lower</code> is set to 0.
+</p>
+<p>
+If neither <code>eps</code> nor <code>eps.lower</code> and/or <code>eps.upper</code> is
+specified, <code>eps.lower</code> and <code>eps.upper</code> are set to 0 and 0.5,
+respectively.
+</p>
+<p>
+If <code>eps</code> is missing, the radius-minimax estimator in sense of
+Rieder et al. (2001, 2008), respectively Section 2.2 of Kohl (2005) is returned.
+</p>
+<p>
+The default value of argument <code>useLast</code> is set by the
+global option <code>kStepUseLast</code> which by default is set to
+<code>FALSE</code>. In case of general models <code>useLast</code>
+remains unchanged during the computations. However, if
+slot <code>CallL2Fam</code> of <code>IC</code> generates an object of
+class <code>"L2GroupParamFamily"</code> the value of <code>useLast</code>
+is changed to <code>TRUE</code>.
+Explicitly setting <code>useLast</code> to <code>TRUE</code> should
+be done with care as in this situation the influence curve
+is re-computed using the value of the one-step estimate
+which may take quite a long time depending on the model.
+</p>
+<p>
+If <code>useLast</code> is set to <code>TRUE</code> the computation of <code>asvar</code>,
+<code>asbias</code> and <code>IC</code> is based on the k-step estimate.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+Object of class <code>"kStepEstimate"</code>.</p>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>References</h3>
+
+<p>
+Kohl, M. (2005) <EM>Numerical Contributions to the Asymptotic Theory of Robustness</EM>.
+Bayreuth: Dissertation.
+</p>
+<p>
+Rieder, H. (1994) <EM>Robust Asymptotic Statistics</EM>. New York: Springer.
+</p>
+<p>
+Rieder, H., Kohl, M. and Ruckdeschel, P. (2008) The Costs of not Knowing
+the Radius. Statistical Methods and Applications <EM>17</EM>(1) 13-40.
+</p>
+<p>
+Rieder, H., Kohl, M. and Ruckdeschel, P. (2001) The Costs of not Knowing
+the Radius. Appeared as discussion paper Nr. 81.
+SFB 373 (Quantification and Simulation of Economic Processes),
+Humboldt University, Berlin; also available under
+<a href="www.uni-bayreuth.de/departments/math/org/mathe7/RIEDER/pubs/RR.pdf">www.uni-bayreuth.de/departments/math/org/mathe7/RIEDER/pubs/RR.pdf</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('RobLox', 'roblox.html')" style="text-decoration: underline; color: blue; cursor: hand">roblox</a></code>,
+<code><a onclick="findlink('distrMod', 'L2ParamFamily-class.html')" style="text-decoration: underline; color: blue; cursor: hand">L2ParamFamily-class</a></code>
+<code><a onclick="findlink('RobAStBase', 'UncondNeighborhood-class.html')" style="text-decoration: underline; color: blue; cursor: hand">UncondNeighborhood-class</a></code>,
+<code><a onclick="findlink('distrMod', 'RiskType-class.html')" style="text-decoration: underline; color: blue; cursor: hand">RiskType-class</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+#############################
+## 1. Binomial data
+#############################
+## generate a sample of contaminated data
+ind <- rbinom(100, size=1, prob=0.05)
+x <- rbinom(100, size=25, prob=(1-ind)*0.25 + ind*0.9)
+
+## ML-estimate
+MLEstimator(x, BinomFamily(size = 25))
+
+## compute optimally robust estimator (known contamination)
+roptest(x, BinomFamily(size = 25), eps = 0.05, steps = 3)
+
+## compute optimally robust estimator (unknown contamination)
+roptest(x, BinomFamily(size = 25), eps.lower = 0, eps.upper = 0.1, steps = 3)
+
+#############################
+## 2. Poisson data
+#############################
+## Example: Rutherford-Geiger (1910); cf. Feller~(1968), Section VI.7 (a)
+x <- c(rep(0, 57), rep(1, 203), rep(2, 383), rep(3, 525), rep(4, 532),
+ rep(5, 408), rep(6, 273), rep(7, 139), rep(8, 45), rep(9, 27),
+ rep(10, 10), rep(11, 4), rep(12, 0), rep(13, 1), rep(14, 1))
+
+## ML-estimate
+MLEstimator(x, PoisFamily())
+
+## compute optimally robust estimator (unknown contamination)
+roptest(x, PoisFamily(), eps.upper = 0.05, steps = 3)
+
+#############################
+## 3. Normal (Gaussian) location and scale
+#############################
+## Generate a contaminated sample
+ind <- rbinom(100, size=1, prob=0.05)
+x <- rnorm(100, mean=0, sd=(1-ind) + ind*9)
+
+## ML-estimate
+MLEstimator(x, NormLocationScaleFamily())
+
+## compute optimally robust estimator (known contamination)
+## takes some time
+roptest(x, NormLocationScaleFamily(), eps = 0.05, steps = 3)
+
+## compute optimally robust estimator (unknown contamination)
+## takes some time
+roptest(x, NormLocationScaleFamily(), eps.upper = 0.1, steps = 3)
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>ROptEst</em> version 0.6.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/RobAStBase/chm/ALEstimate-class.html
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/chm/ALEstimate-class.html (rev 0)
+++ branches/robast-0.6/pkg/RobAStBase/chm/ALEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,108 @@
+<html><head><title>ALEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>ALEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: ALEstimate-class">
+<param name="keyword" value="R: pIC">
+<param name="keyword" value="R: pIC,ALEstimate-method">
+<param name="keyword" value="R: asbias">
+<param name="keyword" value="R: asbias,ALEstimate-method">
+<param name="keyword" value="R: show,ALEstimate-method">
+<param name="keyword" value=" ALEstimate-class.">
+</object>
+
+
+<h2>ALEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("ALEstimate", ...)</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"Estimate"</code>, directly.
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>pIC</dt><dd><code>signature(object = "ALEstimate")</code>:
+accessor function for slot <code>pIC</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "ALEstimate")</code> </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('distrMod', 'Estimate-class.html')" style="text-decoration: underline; color: blue; cursor: hand">Estimate-class</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+## prototype
+new("ALEstimate")
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/RobAStBase/chm/MEstimate-class.html
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/chm/MEstimate-class.html (rev 0)
+++ branches/robast-0.6/pkg/RobAStBase/chm/MEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,102 @@
+<html><head><title>MEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>MEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: MEstimate-class">
+<param name="keyword" value="R: Mroot">
+<param name="keyword" value="R: Mroot,MEstimate-method">
+<param name="keyword" value="R: show,MEstimate-method">
+<param name="keyword" value=" MEstimate-class.">
+</object>
+
+
+<h2>MEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("MEstimate", ...)</code>.
+More frequently they are created via the generating function
+<code>locMEstimator</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>Mroot</code>:</dt><dd>Object of class <code>"numeric"</code>: value of
+the M equation at the estimate. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"ALEstimate"</code>, directly.<br>
+Class <code>"Estimate"</code>, by class <code>"ALEstimate"</code>.
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>Mroot</dt><dd><code>signature(object = "MEstimate")</code>:
+accessor function for slot <code>Mroot</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "MEstimate")</code> </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="ALEstimate-class.html">ALEstimate-class</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+## prototype
+new("MEstimate")
+</pre>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/RobAStBase/chm/RobAStBaseOptions.html
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/chm/RobAStBaseOptions.html (rev 0)
+++ branches/robast-0.6/pkg/RobAStBase/chm/RobAStBaseOptions.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,102 @@
+<html><head><title>Function to change the global variables of the package ‘RobAStBase’</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>RobAStBaseOptions(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: RobAStBaseOptions">
+<param name="keyword" value="R: getRobAStBaseOption">
+<param name="keyword" value="R: kStepUseLast">
+<param name="keyword" value=" Function to change the global variables of the package ‘RobAStBase’">
+</object>
+
+
+<h2>Function to change the global variables of the package ‘RobAStBase’</h2>
+
+
+<h3>Description</h3>
+
+<p>
+With <code>RobAStBaseOptions</code> you can inspect and change
+the global variables of the package <span class="pkg">RobAStBase</span>.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+RobAStBaseOptions(...)
+getRobAStBaseOption(x)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>...</code></td>
+<td>
+any options can be defined, using name = value or by passing a list of such tagged values. </td></tr>
+<tr valign="top"><td><code>x</code></td>
+<td>
+a character string holding an option name.</td></tr>
+</table>
+
+<h3>Value</h3>
+
+<p>
+<code>RobAStBaseOptions()</code> returns a list of the global variables.newline
+<code>RobAStBaseOptions(x)</code> returns the global variable <VAR>x</VAR>.newline
+<code>getRobAStBaseOption(x)</code> returns the global variable <VAR>x</VAR>.newline
+<code>RobAStBaseOptions(x=y)</code> sets the value of the global variable <VAR>x</VAR> to <VAR>y</VAR>.</p>
+
+<h3>Global Options</h3>
+
+<dl>
+<dt>kStepUseLast:</dt><dd>The default value of argument <code>kStepUseLast</code> is
+<code>FALSE</code>. Explicitly setting <code>kStepUseLast</code> to <code>TRUE</code> should
+be done with care as in this situation the influence curve in case of
+<code>oneStepEstimator</code> and <code>kStepEstimator</code> is re-computed using
+the value of the one- resp. k-step estimate which may take quite a long
+time depending on the model. </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('base', 'options.html')" style="text-decoration: underline; color: blue; cursor: hand">options</a></code>, <code><a onclick="findlink('base', 'getOption.html')" style="text-decoration: underline; color: blue; cursor: hand">getOption</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+RobAStBaseOptions()
+RobAStBaseOptions("kStepUseLast")
+RobAStBaseOptions("kStepUseLast" = TRUE)
+# or
+RobAStBaseOptions(kStepUseLast = 1e-6)
+getRobAStBaseOption("kStepUseLast")
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimate-class.html
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimate-class.html (rev 0)
+++ branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,94 @@
+<html><head><title>kStepEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>kStepEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: kStepEstimate-class">
+<param name="keyword" value="R: steps">
+<param name="keyword" value="R: steps,kStepEstimate-method">
+<param name="keyword" value="R: show,kStepEstimate-method">
+<param name="keyword" value=" kStepEstimate-class.">
+</object>
+
+
+<h2>kStepEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("kStepEstimate", ...)</code>.
+More frequently they are created via the generating function
+<code>kStepEstimator</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>steps</code>:</dt><dd>Object of class <code>"integer"</code>: number
+of steps. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"ALEstimate"</code>, directly.<br>
+Class <code>"Estimate"</code>, by class <code>"ALEstimate"</code>
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>steps</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>steps</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "kStepEstimate")</code> </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="ALEstimate-class.html">ALEstimate-class</a></code>
+</p>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimator.html
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimator.html (rev 0)
+++ branches/robast-0.6/pkg/RobAStBase/chm/kStepEstimator.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,144 @@
+<html><head><title>Generic function for the computation of k-step estimates</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>kStepEstimator(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: kStepEstimator">
+<param name="keyword" value="R: kStepEstimator-methods">
+<param name="keyword" value="R: kStepEstimator,numeric,IC,numeric-method">
+<param name="keyword" value="R: kStepEstimator,matrix,IC,numeric-method">
+<param name="keyword" value="R: kStepEstimator,numeric,IC,Estimate-method">
+<param name="keyword" value="R: kStepEstimator,matrix,IC,Estimate-method">
+<param name="keyword" value=" Generic function for the computation of k-step estimates">
+</object>
+
+
+<h2>Generic function for the computation of k-step estimates</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Generic function for the computation of k-step estimates.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+kStepEstimator(x, IC, start, ...)
+
+## S4 method for signature 'numeric, IC, numeric':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'matrix, IC, numeric':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'numeric, IC, Estimate':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'matrix, IC, Estimate':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>x</code></td>
+<td>
+sample </td></tr>
+<tr valign="top"><td><code>IC</code></td>
+<td>
+object of class <code>"IC"</code> </td></tr>
+<tr valign="top"><td><code>start</code></td>
+<td>
+initial estimate </td></tr>
+<tr valign="top"><td><code>steps</code></td>
+<td>
+integer: number of steps </td></tr>
+<tr valign="top"><td><code>useLast</code></td>
+<td>
+which parameter estimate (initial estimate or
+k-step estimate) shall be used to fill the slots <code>pIC</code>,
+<code>asvar</code> and <code>asbias</code> of the return value. </td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+additional parameters </td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+Given an initial estimation <code>start</code>, a sample <code>x</code>
+and an influence curve <code>IC</code> the corresponding k-step
+estimator is computed.
+</p>
+<p>
+The default value of argument <code>useLast</code> is set by the
+global option <code>kStepUseLast</code> which by default is set to
+<code>FALSE</code>. In case of general models <code>useLast</code>
+remains unchanged during the computations. However, if
+slot <code>CallL2Fam</code> of <code>IC</code> generates an object of
+class <code>"L2GroupParamFamily"</code> the value of <code>useLast</code>
+is changed to <code>TRUE</code>.
+Explicitly setting <code>useLast</code> to <code>TRUE</code> should
+be done with care as in this situation the influence curve
+is re-computed using the value of the one-step estimate
+which may take quite a long time depending on the model.
+</p>
+<p>
+If <code>useLast</code> is set to <code>TRUE</code> and slot <code>modifyIC</code>
+of <code>IC</code> is filled with some function (which can be
+used to re-compute the IC for a different parameter), the
+computation of <code>asvar</code>, <code>asbias</code> and <code>IC</code> is
+based on the k-step estimate.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+Object of class <code>"kStepEstimate"</code>.</p>
+
+<h3>Methods</h3>
+
+<dl>
+<dt>x = "numeric", IC = "IC", start = "numeric"</dt><dd>univariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "numeric"</dt><dd>multivariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "Estimate"</dt><dd>multivariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "Estimate"</dt><dd>multivariate samples. </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>References</h3>
+
+<p>
+Rieder, H. (1994) <EM>Robust Asymptotic Statistics</EM>. New York: Springer.
+</p>
+<p>
+Kohl, M. (2005) <EM>Numerical Contributions to the Asymptotic Theory of Robustness</EM>.
+Bayreuth: Dissertation.
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="IC-class.html">IC-class</a></code>, <code><a href="kStepEstimate-class.html">kStepEstimate-class</a></code>
+</p>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: branches/robast-0.7/pkg/RobAStBase/man/plot-methods.Rd
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/man/plot-methods.Rd (rev 0)
+++ branches/robast-0.7/pkg/RobAStBase/man/plot-methods.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,66 @@
+\name{plot-methods}
+\docType{methods}
+\alias{plot}
+\alias{plot-methods}
+\alias{plot,IC,missing-method}
+\title{ Methods for Function plot in Package `RobAStBase' }
+\description{plot-methods}
+\usage{plot(x, y, ...)
+\S4method{plot}{IC,missing}(x, ..., withSweave = getdistrOption("withSweave"),
+ main = FALSE, inner = TRUE, sub = FALSE,
+ col.inner = par("col.main"), cex.inner = 0.8,
+ bmar = par("mar")[1], tmar = par("mar")[3],
+ mfColRow = TRUE)
+}
+\arguments{
+ \item{x}{object of class \code{"IC"}: IC to be plotted }
+ \item{y}{ missing }
+ \item{withSweave}{logical: if \code{TRUE} (for working with \command{Sweave})
+ no extra device is opened}
+ \item{main}{logical: is a main title to be used? or \cr
+ just as argument \code{main} in \code{\link{plot.default}}.}
+ \item{inner}{logical: panels have their own titles? or \cr
+ character vector of / cast to length number of comparands:
+ \code{main} in \code{\link{plot.default}})}
+ \item{sub}{logical: is a sub-title to be used? or \cr
+ just as argument \code{sub} in \code{\link{plot.default}}.}
+ \item{tmar}{top margin -- useful for non-standard main title sizes}
+ \item{bmar}{bottom margin -- useful for non-standard sub title sizes}
+ \item{cex.inner}{magnification to be used for inner titles relative
+ to the current setting of \code{cex}; as in
+ \code{\link[stats]{par}}}
+ \item{col.inner}{character or integer code; color for the inner title}
+ \item{mfColRow}{shall default partition in panels be used --- defaults to \code{TRUE}}
+ \item{\dots} {further parameters for \code{plot}}
+}
+\details{
+Any parameters of \code{plot.default} may be passed on to this particular
+\code{plot} method.
+
+For main-, inner, and subtitles given as arguments \code{main},
+\code{inner}, and \code{sub}, top and bottom margins are enlarged to 5 resp.
+6 by default but may also be specified by \code{tmar} / \code{bmar} arguments.
+If \code{main} / \code{inner} / \code{sub} are
+logical then if the respective argument is \code{FALSE} nothing is done/plotted,
+but if it is \code{TRUE}, we use a default main title taking up the calling
+arguments in case of \code{main}, default inner titles taking up the
+class and (named) parameter slots of arguments in case of \code{inner},
+and a "generated on <data>"-tag in case of \code{sub}.
+Of course, if \code{main} / \code{inner} / \code{sub} are \code{character}, this
+is used for the title; in case of \code{inner} it is then checked whether it
+has correct length. In all title arguments, the following patterns are substituted:
+\itemize{
+\item{\code{"\%C"}}{class of argument \code{object}}
+\item{\code{"\%A"}}{deparsed argument \code{object}}
+\item{\code{"\%D"}}{time/date-string when the plot was generated}
+}
+}
+\examples{
+IC1 <- new("IC")
+plot(IC1)
+plot(IC1, main = TRUE, panel.first= grid(),
+ col = "blue", cex.main = 2, cex.inner = 1)
+}
+\keyword{methods}
+\keyword{distribution}
+
Modified: pkg/ROptEst/chm/00Index.html
===================================================================
--- pkg/ROptEst/chm/00Index.html 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/ROptEst/chm/00Index.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -10,7 +10,7 @@
<param name="keyword" value=".. contents">
</object>
-<h2>Help pages for package ‘ROptEst’ version 0.6.0</h2>
+<h2>Help pages for package ‘ROptEst’ version 0.6.1</h2>
<p align="center">
<a href="#G">G</a>
Modified: pkg/ROptRegTS/DESCRIPTION
===================================================================
--- pkg/ROptRegTS/DESCRIPTION 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/ROptRegTS/DESCRIPTION 2008-11-12 21:33:45 UTC (rev 210)
@@ -2,12 +2,13 @@
Version: 0.6.0
Date: 2008-07-21
Title: Optimally robust estimation for regression-type models
-Description: Optimally robust estimation for regression-type
- models using S4 classes and methods
-Depends: R (>= 2.4.0), methods, distr(>= 1.9), distrEx(>= 1.9), RandVar(>= 0.6), ROptEst(>= 0.5.0)
+Description: Optimally robust estimation for regression-type models
+ using S4 classes and methods
+Depends: R (>= 2.4.0), methods, distr(>= 1.9), distrEx(>= 1.9),
+ RandVar(>= 0.6), ROptEst(>= 0.5.0)
Author: Matthias Kohl <Matthias.Kohl at stamats.de>, Peter Ruckdeschel
Maintainer: Matthias Kohl <Matthias.Kohl at stamats.de>
LazyLoad: yes
-License: GPL version 2 or later
+License: LGPL-3
URL: http://robast.r-forge.r-project.org/
-Packaged: Mon Aug 6 08:09:59 2007; kohl
+Packaged: Mon Aug 6 08:09:59 2007; kohl
Modified: pkg/RobAStBase/NAMESPACE
===================================================================
--- pkg/RobAStBase/NAMESPACE 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/NAMESPACE 2008-11-12 21:33:45 UTC (rev 210)
@@ -17,6 +17,7 @@
exportClasses("RobAStControl", "RobWeight", "BoundedWeight",
"BdStWeight", "HampelWeight")
exportClasses("ALEstimate", "kStepEstimate", "MEstimate")
+exportClasses("cutoff")
exportMethods("show",
"plot")
exportMethods("type", "radius")
@@ -54,8 +55,12 @@
exportMethods("pIC", "asbias",
"steps",
"Mroot")
+exportMethods("ddPlot")
+exportMethods("cutoff.quantile", "cutoff.quantile<-")
export("ContNeighborhood", "TotalVarNeighborhood")
export("FixRobModel", "InfRobModel")
export("InfluenceCurve", "IC", "ContIC", "TotalVarIC")
export(".eq", ".getDistr")
export("RobAStBaseOptions", "getRobAStBaseOption")
+export("cutoff","cutoff.chisq","cutoff.sememp")
+export("outlyingPlotIC")
Modified: pkg/RobAStBase/R/AllClass.R
===================================================================
--- pkg/RobAStBase/R/AllClass.R 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/R/AllClass.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -255,3 +255,12 @@
untransformed.estimate = NULL,
untransformed.asvar = NULL),
contains = "ALEstimate")
+#################################################
+## "cutoff" class
+#################################################
+setClass("cutoff", representation = representation(name = "character",
+ fct = "function",
+ cutoff.quantile = "numeric"),
+ prototype = prototype(name = "empirical",
+ fct = function(data) quantile(data),
+ cutoff.quantile = 0.95))
Modified: pkg/RobAStBase/R/AllGeneric.R
===================================================================
--- pkg/RobAStBase/R/AllGeneric.R 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/R/AllGeneric.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -179,3 +179,14 @@
if(!isGeneric("modifyIC")){
setGeneric("modifyIC", function(object) standardGeneric("modifyIC"))
}
+if(!isGeneric("cutoff.quantile")){
+ setGeneric("cutoff.quantile", function(object) standardGeneric("cutoff.quantile"))
+}
+if(!isGeneric("cutoff.quantile<-")){
+ setGeneric("cutoff.quantile<-", function(object,value)
+ standardGeneric("cutoff.quantile<-"))
+}
+if(!isGeneric("ddPlot")){
+ setGeneric("ddPlot", function(data, dist.x, dist.y, cutoff.x, cutoff.y,...)
+ standardGeneric("ddPlot"))
+}
Modified: pkg/RobAStBase/R/bALEstimate.R
===================================================================
--- pkg/RobAStBase/R/bALEstimate.R 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/R/bALEstimate.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -65,7 +65,10 @@
ci[] <- main(object) + sd0 %o% fac
### end of borrowed code
- new("Confint", type = gettext("asymptotic (LAN-based)"),
+ new("Confint", type = c(
+ gettext("asymptotic (LAN-based), uniform (bias-aware)\n"),
+ gettextf("for %s", name(method))
+ ),
samplesize.estimate = object at samplesize,
call.estimate = object at estimate.call,
name.estimate = object at name,
@@ -105,7 +108,10 @@
ci[] <- main(object) + sd0 %o% fac
### end of borrowed code
- new("Confint", type = gettext("asymptotic (LAN-based)"),
+ new("Confint", type = c(
+ gettext("asymptotic (LAN-based), uniform (bias-aware)\n"),
+ gettextf("for %s", name(method))
+ ),
samplesize.estimate = object at samplesize,
call.estimate = object at estimate.call,
name.estimate = object at name,
@@ -135,13 +141,18 @@
a <- (1 - level)/2
a <- c(a, 1 - a)
pct <- stats:::format.perc(a, 3)
- fac <- qnorm(a, mean = c(-object at asbias, object at asbias)*method at nu)
+ fac <- qnorm(a, mean = c(-object at asbias, object at asbias)/method at nu)
ci <- array(NA, dim = c(length(object at estimate), 2),
dimnames = list(names(object at estimate), pct))
ci[] <- main(object) + sd0 %o% fac
### end of borrowed code
- new("Confint", type = gettext("asymptotic (LAN-based)"),
+ nuround <- round(nu,3)
+ new("Confint", type = c(
+ gettext("asymptotic (LAN-based), uniform (bias-aware)\n"),
+ gettextf("for %s with nu =(%f,%f)",
+ name(method), nuround[1], nuround[2])
+ ),
samplesize.estimate = object at samplesize,
call.estimate = object at estimate.call,
name.estimate = object at name,
Added: pkg/RobAStBase/R/cutoff-class.R
===================================================================
--- pkg/RobAStBase/R/cutoff-class.R (rev 0)
+++ pkg/RobAStBase/R/cutoff-class.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,52 @@
+
+
+setMethod("cutoff.quantile", "cutoff", function(object)object at cutoff.quantile)
+setReplaceMethod("cutoff.quantile", "cutoff", function(object, value){
+ object at cutoff.quantile<-value; object})
+
+setMethod("fct", "cutoff", function(object)object at fct)
+setMethod("name", "cutoff", function(object)object at name)
+
+
+
+cutoff <- function(name = "empirical",
+ body.fct0,
+ cutoff.quantile = 0.95,
+ norm = NormType(), QF, nsim = 100000){
+ mc <- match.call()
+ if (missing(body.fct0)) body.fct0 <- quote(quantile(slot(norm,"fct")(data), cutoff.quantile))
+ arglist <- list()
+ arglist$norm <- if (!is.null(mc$norm)) mc$norm else NormType()
+ QFsub <- if(is.null(mc$QF))
+ quote({QF <- if(is(norm,"QFNorm"))
+ QuadForm(norm) else diag(nrow(data))})
+ else quote({})
+ arglist$cutoff.quantile <- if (!is.null(mc$cutoff.quantile))
+ mc$cutoff.quantile else 0.95
+ arglist$nsim <- if (!is.null(mc$nsim)) mc$nsim else 100000
+ arglist$QFsub0 <- substitute(QFsub)
+ arglist$body.fct1 <- body.fct0
+ fct0 <- function(data){}
+ body(fct0) <- #eval(
+ substitute({
+ QFsub0
+ body.fct1
+ }, arglist)
+ #)
+ new("cutoff", fct = fct0, name = name, cutoff.quantile = cutoff.quantile)
+}
+
+cutoff.sememp <- function(){cutoff(name = "semi-empirical",
+ body.fct0 = substitute({n.05 <- chol(QF)
+# print(QF)
+ N0 <- matrix(rnorm(nsim*nrow(QF)),ncol=ncol(QF))
+ N0 <- N0 %*% n.05
+ quantile((rowSums(N0^2))^.5,cutoff.quantile)
+ }),
+ cutoff.quantile = 0.95)}
+
+cutoff.chisq <- function(){cutoff(name = "chisq",
+ body.fct0 = substitute({dim = nrow(data)
+ qchisq(df=dim,cutoff.quantile)^.5
+ }),
+ cutoff.quantile = 0.95)}
Added: pkg/RobAStBase/R/ddPlot.R
===================================================================
--- pkg/RobAStBase/R/ddPlot.R (rev 0)
+++ pkg/RobAStBase/R/ddPlot.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,44 @@
+setMethod("ddPlot", signature = signature(data = "matrix"),
+ function(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff){
+ mc <- match.call(expand.dots = FALSE, call = sys.call(sys.parent(1)))
+ dots <- mc$"..."
+
+do.call(.ddPlot.MatNtNtCoCo, args = c(list(data=data), dots, list(dist.x = mc$dist.x,
+ dist.y = mc$dist.y, cutoff.x = mc$cutoff.x, cutoff.y = mc$cutoff.y,
+ cutoff.quantile.x = mc$cutoff.quantile.x, cutoff.quantile.y = mc$cutoff.quantile.y,
+ transform.x = mc$transform.x, transform.y = mc$transform.y,
+ id.n = mc$id.n, lab.pts = mc$lab.pts, adj = mc$adj, cex.idn = mc$cex.idn,
+ col.idn = mc$col.idn, lty.cutoff = mc$lty.cutoff,
+ lwd.cutoff = mc$lwd.cutoff, col.cutoff = mc$col.cutoff)))
+})
+
+setMethod("ddPlot", signature = signature(data = "data.frame"),
+ function(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff){
+
+ mc <- match.call(call = sys.call(sys.parent(1)))
+ mc$data <- t(as.matrix(data))
+ do.call(ddPlot,args = as.list(mc[-1]))
+ })
+
+setMethod("ddPlot", signature = signature(data = "numeric"),
+ function(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff){
+
+ mc <- match.call(call = sys.call(sys.parent(1)))
+ mc$data <- matrix(data,nrow=1)
+ do.call(ddPlot,args = as.list(mc[-1]))
+ })
Added: pkg/RobAStBase/R/ddPlot_utils.R
===================================================================
--- pkg/RobAStBase/R/ddPlot_utils.R (rev 0)
+++ pkg/RobAStBase/R/ddPlot_utils.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,153 @@
+.ddPlot.MatNtNtCoCo <- function(data, ..., dist.x = NormType(), dist.y = NormType(),
+ cutoff.x = cutoff(norm = dist.x, cutoff.quantile = cutoff.quantile.x),
+ cutoff.y = cutoff(norm = dist.y, cutoff.quantile = cutoff.quantile.y),
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj =0, cex.idn = 1,
+ col.idn = par("col"), lty.cutoff,
+ lwd.cutoff, col.cutoff = "red"){
+
+ dots <- match.call(expand.dots = FALSE)$"..."
+
+ id.n1 <- 1:ncol(data)
+
+ if(missing(id.n) || is.null(id.n))
+ id.n <- id.n1
+
+
+ if(missing(lab.pts)|| is.null(lab.pts)){
+ lab.pts <- if(!is.null(colnames(data))) colnames(data) else 1:ncol(data)
+ }
+
+ data <- data[,id.n, drop = FALSE]
+ lab.pts <- lab.pts[id.n]
+ id.n1 <- id.n1[id.n]
+
+ id.n0 <- 1:length(id.n)
+
+ data.x <- if(missing(transform.x)||is.null(transform.x)) data else
+ transform.x(data)
+ #if(!is(IC.x,"IC")) stop("Argument 'IC.x' of 'ddPlot' must be an 'IC'")
+ # apply(data,2,function(xx) evalIC(IC.x,xx))}
+ data.y <- if(missing(transform.y)||is.null(transform.y)) data else
+ transform.y(data)
+ #if(!is(IC.y,"IC")) stop("Argument 'IC.y' of 'ddPlot' must be an 'IC'")
+ # apply(data,2,function(xx) evalIC(IC.y,xx))}
+
+
+ if(is.null(dist.x)) dist.x <- NormType()
+ if(is.null(dist.y)) dist.y <- NormType()
+ if(is.null(dots$xlab)) dots$xlab <- name(dist.x)
+ if(is.null(dots$ylab)) dots$ylab <- name(dist.y)
+
+ if(is.null(cutoff.quantile.x))
+ cutoff.quantile.x <- 0.95
+
+ if(is.null(cutoff.quantile.y))
+ cutoff.quantile.y <- cutoff.quantile.x
+
+ if(is.null(cutoff.x))
+ cutoff.x <- cutoff(norm = dist.x, cutoff.quantile = cutoff.quantile.x)
+ else {assign("norm", dist.x, environment(fct(cutoff.x)))
+ assign("cutoff.quantile", cutoff.quantile.x, environment(fct(cutoff.x)))}
+
+ if(is.null(cutoff.y))
+ cutoff.y <- cutoff(norm = dist.y, cutoff.quantile = cutoff.quantile.y)
+ else {assign("norm", dist.y, environment(fct(cutoff.y)))
+ assign("cutoff.quantile", cutoff.quantile.y, environment(fct(cutoff.y)))}
+
+ if(!is(dist.x, "NormType")) stop("Argument 'dist.x' of 'ddPlot' must be of class 'NormType'")
+ if(!is(dist.y, "NormType")) stop("Argument 'dist.y' of 'ddPlot' must be of class 'NormType'")
+ if(!is(cutoff.x, "cutoff")) stop("Argument 'cutoff.x' of 'ddPlot' must be of class 'cutoff'")
+ if(!is(cutoff.y, "cutoff")) stop("Argument 'cutoff.y' of 'ddPlot' must be of class 'cutoff'")
+
+ ndata.x <- fct(dist.x)(data.x)
+ ndata.y <- fct(dist.y)(data.y)
+
+
+ if(is.null(adj)) adj <- 0
+ if(is.null(cex.idn)) cex.idn <- 1
+ if(is.null(col.idn)) col.idn <- par("col")
+ if(is.null(col.cutoff)) col.cutoff <- "red"
+
+ if(is.null(dots$lwd)) dots$lwd <- par("lwd")
+ if(is.null(dots$lty)) dots$lty <- par("lty")
+
+
+ pdots <- dots
+ pdots$type <- NULL
+ pdots$x <- NULL
+ pdots$y <- NULL
+ pdots$offset <- NULL
+ pdots$pos <- NULL
+ pdots$log <- NULL
+ pdots$untf <- NULL
+
+ abdots <- pdots
+ abdots$col <- col.cutoff
+ if(!missing(lwd.cutoff)) abdots$lwd <- lwd.cutoff
+ if(!missing(lty.cutoff)) abdots$lty <- lty.cutoff
+ abdots$pos <- NULL
+ abdots$untf <- dots$untf
+ abdots$adj <- NULL
+
+ adots <- pdots
+ adots$col <- pdots$col.axis
+ adots$lty <- pdots$lty.axis
+ adots$adj <- par("adj")
+
+ tdots <- pdots
+ tdots$cex <- cex.idn
+ tdots$col <- col.idn
+ tdots$offset <- dots$offset
+ tdots$pos <- dots$pos
+ tdots$adj <- adj
+
+ pdots$axes <- FALSE
+ pdots$log <- dots$log
+ pdots$adj <- par("adj")
+
+ ####
+
+ co.x <- fct(cutoff.x)(data.x)
+ co.y <- fct(cutoff.y)(data.y)
+# print(quantile(ndata.x))
+# print(co.x)
+# print(fct(cutoff.x))
+# print(dist.x)
+# print(get("norm", environment(fct(cutoff.x))))
+# print(quantile(ndata.y))
+# print(co.y)
+# print(fct(cutoff.y))
+# print(get("norm", environment(fct(cutoff.y))))
+##
+
+ if(!is.null(dots$xlim))
+ if(is.logical(dots$xlim))
+ if(dots$xlim) pdots$xlim <- c(min(ndata.x)/1.03,max(ndata.x,co.x)*1.03)
+ if(!is.null(dots$ylim))
+ if(is.logical(dots$ylim))
+ if(dots$ylim) pdots$ylim <- c(min(ndata.y)/1.03,max(ndata.y,co.y)*1.03)
+
+
+ id.x <- id.n0[ndata.x >= co.x*.999]
+ id.y <- id.n0[ndata.y >= co.y*.999]
+ id.xy <- intersect(id.x,id.y)
+
+
+ id0.xy <- id.n1[id.xy]
+ id0.x <- id.n1[id.x]
+ id0.y <- id.n1[id.y]
+
+ do.call(plot, args = c(list(x = ndata.x,ndata.y, type = "p"), pdots))
+ do.call(box,args=c(adots))
+ do.call(abline, args = c(list(h=co.y), abdots))
+ do.call(abline, args = c(list(v=co.x), abdots))
+ if(length(id.xy))
+ do.call(text, args = c(list(ndata.x[id.xy], ndata.y[id.xy],
+ labels=lab.pts[id.xy]), tdots))
+ return(list(id.x=id0.x, id.y= id0.y, id.xy = id0.xy,
+ qtx = quantile(ndata.x), qty = quantile(ndata.y),
+ cutoff.x.v = co.x, cutoff.y.v = co.y
+ ))
+}
Added: pkg/RobAStBase/R/outlyingPlot.R
===================================================================
--- pkg/RobAStBase/R/outlyingPlot.R (rev 0)
+++ pkg/RobAStBase/R/outlyingPlot.R 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,48 @@
+outlyingPlotIC <- function(data, IC.x, IC.y = IC.x, dist.x = NormType(),
+ dist.y, cutoff.y = cutoff.chisq(), cutoff.x = cutoff.sememp(), ...,
+ cutoff.quantile.x = 0.95,
+ cutoff.quantile.y = cutoff.quantile.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff,
+ main = gettext("Outlyingness by means of a distance-distance plot")
+ ){
+ mc <- as.list(match.call(expand.dots = FALSE))[-1]
+ dots <- mc$"..."
+ if(is.null(dots$xlim)) dots$xlim <- TRUE
+ if(is.null(dots$ylim)) dots$ylim <- TRUE
+ if(is.null(mc$cutoff.quantile.x)) mc$cutoff.quantile.x <- 0.95
+ if(is.null(mc$cutoff.quantile.y)) mc$cutoff.quantile.y <- cutoff.quantile.x
+ if(is.null(mc$cutoff.x)) mc$cutoff.x <- cutoff.sememp()
+ if(is.null(mc$cutoff.y)) mc$cutoff.y <- cutoff.chisq()
+ if(missing(IC.x)) stop("Argument 'IC.x' must be given as argument to 'outlyingPlot'")
+ if(missing(data)) stop("Argument 'data' must be given as argument to 'outlyingPlot'")
+
+ if(missing(dist.y)){
+ if("asCov" %in% names(Risks(IC.y)))
+ if(is.matrix(Risks(IC.y)$asCov) || length(Risks(IC.y)$asCov) == 1)
+ asVar <- Risks(IC.y)$asCov
+ else
+ asVar <- Risks(IC.y)$asCov$value
+ else
+ asVar <- getRiskIC(IC.y, risk = asCov())$asCov$value
+
+ asVar <- PosSemDefSymmMatrix(solve(asVar))
+ mc$dist.y <- QFNorm(name = gettext("Mahalonobis-Norm"), QuadForm = asVar)
+ }
+ if(missing(dist.x))
+ mc$dist.x <- NormType()
+
+ tf.x <- function(x) apply(x,2,function(xx) evalIC(IC.x,xx))
+ tf.y <- function(x) apply(x,2,function(xx) evalIC(IC.y,xx))
+
+
+ do.call(ddPlot,args=c(list(data=data), dots, list(dist.x = mc$dist.x,
+ dist.y = mc$dist.y, cutoff.x = mc$cutoff.x, cutoff.y = mc$cutoff.y,
+ cutoff.quantile.x = mc$cutoff.quantile.x, cutoff.quantile.y = mc$cutoff.quantile.y,
+ transform.x = tf.x, transform.y = tf.y,
+ id.n = mc$id.n, lab.pts = mc$lab.pts, adj = mc$adj, cex.idn = mc$cex.idn,
+ col.idn = mc$col.idn, lty.cutoff = mc$lty.cutoff,
+ lwd.cutoff = mc$lwd.cutoff, col.cutoff = mc$col.cutoff, main = main)))
+
+ }
+
Modified: pkg/RobAStBase/chm/00Index.html
===================================================================
--- pkg/RobAStBase/chm/00Index.html 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/chm/00Index.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -10,9 +10,10 @@
<param name="keyword" value=".. contents">
</object>
-<h2>Help pages for package ‘RobAStBase’ version 0.1.3</h2>
+<h2>Help pages for package ‘RobAStBase’ version 0.1.2</h2>
<p align="center">
+<a href="# "> </a>
<a href="#A">A</a>
<a href="#B">B</a>
<a href="#C">C</a>
@@ -37,6 +38,13 @@
<table width="100%">
</table>
+<h2><a name=" ">-- --</a></h2>
+
+<table width="100%">
+<tr><td width="25%"><a href="0RobAStBase-package.html">RobAStBase-package</a></td>
+<td>Robust Asymptotic Statistics </td></tr>
+</table>
+
<h2><a name="A">-- A --</a></h2>
<table width="100%">
@@ -126,6 +134,14 @@
<td>Compare - Plots</td></tr>
<tr><td width="25%"><a href="comparePlot.html">comparePlot-methods</a></td>
<td>Compare - Plots</td></tr>
+<tr><td width="25%"><a href="ALEstimate-class.html">confint,ALEstimate,asymmetricBias-method</a></td>
+<td>ALEstimate-class.</td></tr>
+<tr><td width="25%"><a href="ALEstimate-class.html">confint,ALEstimate,missing-method</a></td>
+<td>ALEstimate-class.</td></tr>
+<tr><td width="25%"><a href="ALEstimate-class.html">confint,ALEstimate,onesidedBias-method</a></td>
+<td>ALEstimate-class.</td></tr>
+<tr><td width="25%"><a href="ALEstimate-class.html">confint,ALEstimate,symmetricBias-method</a></td>
+<td>ALEstimate-class.</td></tr>
<tr><td width="25%"><a href="ContIC.html">ContIC</a></td>
<td>Generating function for ContIC-class</td></tr>
<tr><td width="25%"><a href="ContIC-class.html">ContIC-class</a></td>
@@ -138,11 +154,35 @@
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">Curve,InfluenceCurve-method</a></td>
<td>Influence curve</td></tr>
+<tr><td width="25%"><a href="cutoff.html">cutoff</a></td>
+<td>Generating function(s) for class 'cutoff'</td></tr>
+<tr><td width="25%"><a href="cutoff-class.html">cutoff-class</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
+<tr><td width="25%"><a href="cutoff.html">cutoff.chisq</a></td>
+<td>Generating function(s) for class 'cutoff'</td></tr>
+<tr><td width="25%"><a href="cutoff-class.html">cutoff.quantile</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
+<tr><td width="25%"><a href="cutoff-class.html">cutoff.quantile,cutoff-method</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
+<tr><td width="25%"><a href="cutoff-class.html">cutoff.quantile<-,cutoff-method</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
+<tr><td width="25%"><a href="cutoff.html">cutoff.sememp</a></td>
+<td>Generating function(s) for class 'cutoff'</td></tr>
</table>
<h2><a name="D">-- D --</a></h2>
<table width="100%">
+<tr><td width="25%"><a href="ddPlot-methods.html">ddPlot</a></td>
+<td>Methods for Function ddPlot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="ddPlot-methods.html">ddPlot,data.frame-method</a></td>
+<td>Methods for Function ddPlot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="ddPlot-methods.html">ddPlot,matrix-method</a></td>
+<td>Methods for Function ddPlot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="ddPlot-methods.html">ddPlot,numeric-method</a></td>
+<td>Methods for Function ddPlot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="ddPlot-methods.html">ddPlot-methods</a></td>
+<td>Methods for Function ddPlot in Package ‘RobAStBase’ </td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">Domain,InfluenceCurve-method</a></td>
<td>Influence curve</td></tr>
</table>
@@ -161,6 +201,8 @@
<h2><a name="F">-- F --</a></h2>
<table width="100%">
+<tr><td width="25%"><a href="cutoff-class.html">fct,cutoff-method</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
<tr><td width="25%"><a href="FixRobModel.html">FixRobModel</a></td>
<td>Generating function for FixRobModel-class</td></tr>
<tr><td width="25%"><a href="FixRobModel-class.html">FixRobModel-class</a></td>
@@ -348,6 +390,8 @@
<h2><a name="N">-- N --</a></h2>
<table width="100%">
+<tr><td width="25%"><a href="cutoff-class.html">name,cutoff-method</a></td>
+<td>Cutoff class for distance-distance plots</td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">name,InfluenceCurve-method</a></td>
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="RobAStControl-class.html">name,RobAStControl-method</a></td>
@@ -409,6 +453,8 @@
<td>Generic function for the computation of optimally robust ICs</td></tr>
<tr><td width="25%"><a href="optIC.html">optIC-methods</a></td>
<td>Generic function for the computation of optimally robust ICs</td></tr>
+<tr><td width="25%"><a href="outlyingPlotIC.html">outlyingPlotIC</a></td>
+<td>Function outlyingPlotIC in Package ‘RobAStBase’ </td></tr>
</table>
<h2><a name="P">-- P --</a></h2>
@@ -441,6 +487,8 @@
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">Risks<-,InfluenceCurve-method</a></td>
<td>Influence curve</td></tr>
+<tr><td width="25%"><a href="0RobAStBase-package.html">RobAStBase</a></td>
+<td>Robust Asymptotic Statistics </td></tr>
<tr><td width="25%"><a href="RobAStBaseOptions.html">RobAStBaseOptions</a></td>
<td>Function to change the global variables of the package ‘RobAStBase’ </td></tr>
<tr><td width="25%"><a href="RobAStControl-class.html">RobAStControl-class</a></td>
Added: pkg/RobAStBase/chm/0RobAStBase-package.html
===================================================================
--- pkg/RobAStBase/chm/0RobAStBase-package.html (rev 0)
+++ pkg/RobAStBase/chm/0RobAStBase-package.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,109 @@
+<html><head><title>Robust Asymptotic Statistics</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>RobAStBase-package(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: RobAStBase-package">
+<param name="keyword" value="R: RobAStBase">
+<param name="keyword" value=" Robust Asymptotic Statistics">
+</object>
+
+
+<h2>Robust Asymptotic Statistics</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Base S4-classes and functions for robust asymptotic statistics.
+</p>
+
+
+<h3>Details</h3>
+
+<p>
+<table summary="Rd table">
+<tr>
+ <td align="left">Package: </td> <td align="left"> RobAStBase</td>
+</tr>
+<tr>
+ <td align="left"> Type: </td> <td align="left"> Package</td>
+</tr>
+<tr>
+ <td align="left"> Version: </td> <td align="left"> 0.1.2</td>
+</tr>
+<tr>
+ <td align="left"> Date: </td> <td align="left"> 2008-11-05</td>
+</tr>
+<tr>
+ <td align="left"> Depends: </td> <td align="left"> R(>= 2.7.0), methods, distr(>= 2.0), distrEx(>= 2.0), distrMod(>= 2.0), RandVar(>= 0.6.3)</td>
+</tr>
+<tr>
+ <td align="left"> LazyLoad: </td> <td align="left"> yes</td>
+</tr>
+<tr>
+ <td align="left"> License: </td> <td align="left"> LGPL-3</td>
+</tr>
+<tr>
+ <td align="left"> URL: </td> <td align="left"> http://robast.r-forge.r-project.org/</td>
+</tr>
+</table>
+
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>,<br>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a><br>
+</p>
+<p>
+Maintainer: Matthias Kohl <a href="mailto:matthias.kohl at stamats.de">matthias.kohl at stamats.de</a>
+</p>
+
+
+<h3>References</h3>
+
+<p>
+M. Kohl (2005). Numerical Contributions to the Asymptotic Theory of Robustness.
+Dissertation. University of Bayreuth.
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('distr', 'distr-package.html')" style="text-decoration: underline; color: blue; cursor: hand">distr-package</a></code>, <code><a onclick="findlink('distrEx', 'distrEx-package.html')" style="text-decoration: underline; color: blue; cursor: hand">distrEx-package</a></code>,
+<code><a onclick="findlink('distrMod', 'distrMod-package.html')" style="text-decoration: underline; color: blue; cursor: hand">distrMod-package</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+library(RobAStBase)
+
+## some L2 differentiable parametric family from package distrMod, e.g.
+B <- BinomFamily(size = 25, prob = 0.25)
+
+## classical optimal IC
+IC0 <- optIC(model = B, risk = asCov())
+plot(IC0) # plot IC
+checkIC(IC0, B)
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/ALEstimate-class.html
===================================================================
--- pkg/RobAStBase/chm/ALEstimate-class.html (rev 0)
+++ pkg/RobAStBase/chm/ALEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,131 @@
+<html><head><title>ALEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>ALEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: ALEstimate-class">
+<param name="keyword" value="R: pIC">
+<param name="keyword" value="R: pIC,ALEstimate-method">
+<param name="keyword" value="R: asbias">
+<param name="keyword" value="R: asbias,ALEstimate-method">
+<param name="keyword" value="R: show,ALEstimate-method">
+<param name="keyword" value="R: confint,ALEstimate,missing-method">
+<param name="keyword" value="R: confint,ALEstimate,symmetricBias-method">
+<param name="keyword" value="R: confint,ALEstimate,onesidedBias-method">
+<param name="keyword" value="R: confint,ALEstimate,asymmetricBias-method">
+<param name="keyword" value=" ALEstimate-class.">
+</object>
+
+
+<h2>ALEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("ALEstimate", ...)</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"Estimate"</code>, directly.
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>pIC</dt><dd><code>signature(object = "ALEstimate")</code>:
+accessor function for slot <code>pIC</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "ALEstimate")</code> </dd>
+
+
+<dt>confint</dt><dd><code>signature(object = "ALEstimate", method = "missing")</code>:
+compute asymptotic (LAN-based) confidence interval neglecting any bias. </dd>
+
+
+<dt>confint</dt><dd><code>signature(object = "ALEstimate", method = "symmetricBias")</code>:
+compute asymptotic (LAN-based) confidence interval incorporating bias
+symmetrically. </dd>
+
+
+<dt>confint</dt><dd><code>signature(object = "ALEstimate", method = "onesidedBias")</code>:
+compute asymptotic (LAN-based) confidence interval incorporating bias
+one-sided; i.e., positive or negative, respectively. </dd>
+
+
+<dt>confint</dt><dd><code>signature(object = "ALEstimate", method = "asymmetricBias")</code>:
+compute asymptotic (LAN-based) confidence interval incorporating bias
+asymmetrically. </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('distrMod', 'Estimate-class.html')" style="text-decoration: underline; color: blue; cursor: hand">Estimate-class</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+## prototype
+new("ALEstimate")
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Modified: pkg/RobAStBase/chm/IC-class.html
===================================================================
--- pkg/RobAStBase/chm/IC-class.html 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/chm/IC-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -105,7 +105,7 @@
Plot absolute and relative information of <code>IC</code>. </dd>
-<dt>plot</dt><dd><code>signature(x = "IC")</code></dd>
+<dt>plot</dt><dd><code>signature(x = "IC", y = "missing")</code></dd>
<dt>show</dt><dd><code>signature(object = "IC")</code></dd>
@@ -158,6 +158,6 @@
</script>
-<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.3 <a href="00Index.html">Index]</a></div>
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
</body></html>
Added: pkg/RobAStBase/chm/MEstimate-class.html
===================================================================
--- pkg/RobAStBase/chm/MEstimate-class.html (rev 0)
+++ pkg/RobAStBase/chm/MEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,102 @@
+<html><head><title>MEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>MEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: MEstimate-class">
+<param name="keyword" value="R: Mroot">
+<param name="keyword" value="R: Mroot,MEstimate-method">
+<param name="keyword" value="R: show,MEstimate-method">
+<param name="keyword" value=" MEstimate-class.">
+</object>
+
+
+<h2>MEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("MEstimate", ...)</code>.
+More frequently they are created via the generating function
+<code>locMEstimator</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>Mroot</code>:</dt><dd>Object of class <code>"numeric"</code>: value of
+the M equation at the estimate. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"ALEstimate"</code>, directly.<br>
+Class <code>"Estimate"</code>, by class <code>"ALEstimate"</code>.
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>Mroot</dt><dd><code>signature(object = "MEstimate")</code>:
+accessor function for slot <code>Mroot</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "MEstimate")</code> </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="ALEstimate-class.html">ALEstimate-class</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+## prototype
+new("MEstimate")
+</pre>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Modified: pkg/RobAStBase/chm/RobAStBase.chm
===================================================================
(Binary files differ)
Modified: pkg/RobAStBase/chm/RobAStBase.hhp
===================================================================
--- pkg/RobAStBase/chm/RobAStBase.hhp 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/chm/RobAStBase.hhp 2008-11-12 21:33:45 UTC (rev 210)
@@ -12,6 +12,7 @@
[FILES]
00Index.html
+0RobAStBase-package.html
ALEstimate-class.html
BdStWeight-class.html
BoundedWeight-class.html
@@ -42,6 +43,9 @@
UncondNeighborhood-class.html
checkIC.html
comparePlot.html
+cutoff-class.html
+cutoff.html
+ddPlot-methods.html
evalIC.html
generateIC.html
generateICfct.html
@@ -50,10 +54,12 @@
getweight.html
infoPlot.html
internals.html
+internals_ddPlot.html
kStepEstimate-class.html
kStepEstimator.html
locMEstimator.html
makeIC-methods.html
oneStepEstimator.html
optIC.html
+outlyingPlotIC.html
plot-methods.html
Modified: pkg/RobAStBase/chm/RobAStBase.toc
===================================================================
--- pkg/RobAStBase/chm/RobAStBase.toc 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/chm/RobAStBase.toc 2008-11-12 21:33:45 UTC (rev 210)
@@ -10,6 +10,14 @@
</OBJECT>
<UL>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="RobAStBase-package">
+<param name="Local" value="0RobAStBase-package.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value=".ddPlot.MatNtNtCoCo">
+<param name="Local" value="internals_ddPlot.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value=".eq">
<param name="Local" value="internals.html">
</OBJECT>
@@ -202,6 +210,22 @@
<param name="Local" value="comparePlot.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="confint,ALEstimate,asymmetricBias-method">
+<param name="Local" value="ALEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="confint,ALEstimate,missing-method">
+<param name="Local" value="ALEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="confint,ALEstimate,onesidedBias-method">
+<param name="Local" value="ALEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="confint,ALEstimate,symmetricBias-method">
+<param name="Local" value="ALEstimate-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="ContIC">
<param name="Local" value="ContIC.html">
</OBJECT>
@@ -226,6 +250,58 @@
<param name="Local" value="InfluenceCurve-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff">
+<param name="Local" value="cutoff.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff-class">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.chisq">
+<param name="Local" value="cutoff.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.quantile">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.quantile,cutoff-method">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.quantile<-">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.quantile<-,cutoff-method">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="cutoff.sememp">
+<param name="Local" value="cutoff.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ddPlot">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ddPlot,data.frame-method">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ddPlot,matrix-method">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ddPlot,numeric-method">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="ddPlot-methods">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Domain,InfluenceCurve-method">
<param name="Local" value="InfluenceCurve-class.html">
</OBJECT>
@@ -242,6 +318,10 @@
<param name="Local" value="IC-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="fct,cutoff-method">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="FixRobModel">
<param name="Local" value="FixRobModel.html">
</OBJECT>
@@ -430,6 +510,10 @@
<param name="Local" value="internals.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="internals_for_RobAStBase_ddPlot">
+<param name="Local" value="internals_ddPlot.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="kStepEstimate-class">
<param name="Local" value="kStepEstimate-class.html">
</OBJECT>
@@ -558,6 +642,10 @@
<param name="Local" value="MEstimate-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="name,cutoff-method">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="name,InfluenceCurve-method">
<param name="Local" value="InfluenceCurve-class.html">
</OBJECT>
@@ -678,6 +766,10 @@
<param name="Local" value="optIC.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="outlyingPlotIC">
+<param name="Local" value="outlyingPlotIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="pIC">
<param name="Local" value="ALEstimate-class.html">
</OBJECT>
@@ -726,6 +818,10 @@
<param name="Local" value="InfluenceCurve-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="RobAStBase">
+<param name="Local" value="0RobAStBase-package.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="RobAStBaseOptions">
<param name="Local" value="RobAStBaseOptions.html">
</OBJECT>
@@ -887,6 +983,14 @@
<param name="Local" value="RobAStControl-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Cutoff class for distance-distance plots">
+<param name="Local" value="cutoff-class.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Function outlyingPlotIC in Package `RobAStBase' ">
+<param name="Local" value="outlyingPlotIC.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Function to change the global variables of the package `RobAStBase' ">
<param name="Local" value="RobAStBaseOptions.html">
</OBJECT>
@@ -923,6 +1027,10 @@
<param name="Local" value="TotalVarNeighborhood.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Generating function(s) for class 'cutoff'">
+<param name="Local" value="cutoff.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Generating weights">
<param name="Local" value="getweight.html">
</OBJECT>
@@ -987,6 +1095,10 @@
<param name="Local" value="internals.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Internal / Helper functions of package RobAStBase for ddPlot">
+<param name="Local" value="internals_ddPlot.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="kStepEstimate-class.">
<param name="Local" value="kStepEstimate-class.html">
</OBJECT>
@@ -995,6 +1107,10 @@
<param name="Local" value="MEstimate-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Methods for Function ddPlot in Package `RobAStBase' ">
+<param name="Local" value="ddPlot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Methods for Function plot in Package `RobAStBase' ">
<param name="Local" value="plot-methods.html">
</OBJECT>
@@ -1007,6 +1123,10 @@
<param name="Local" value="infoPlot.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="Robust Asymptotic Statistics ">
+<param name="Local" value="0RobAStBase-package.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Robust model">
<param name="Local" value="RobModel-class.html">
</OBJECT>
Added: pkg/RobAStBase/chm/RobAStBaseOptions.html
===================================================================
--- pkg/RobAStBase/chm/RobAStBaseOptions.html (rev 0)
+++ pkg/RobAStBase/chm/RobAStBaseOptions.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,102 @@
+<html><head><title>Function to change the global variables of the package ‘RobAStBase’</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>RobAStBaseOptions(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: RobAStBaseOptions">
+<param name="keyword" value="R: getRobAStBaseOption">
+<param name="keyword" value="R: kStepUseLast">
+<param name="keyword" value=" Function to change the global variables of the package ‘RobAStBase’">
+</object>
+
+
+<h2>Function to change the global variables of the package ‘RobAStBase’</h2>
+
+
+<h3>Description</h3>
+
+<p>
+With <code>RobAStBaseOptions</code> you can inspect and change
+the global variables of the package <span class="pkg">RobAStBase</span>.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+RobAStBaseOptions(...)
+getRobAStBaseOption(x)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>...</code></td>
+<td>
+any options can be defined, using name = value or by passing a list of such tagged values. </td></tr>
+<tr valign="top"><td><code>x</code></td>
+<td>
+a character string holding an option name.</td></tr>
+</table>
+
+<h3>Value</h3>
+
+<p>
+<code>RobAStBaseOptions()</code> returns a list of the global variables.newline
+<code>RobAStBaseOptions(x)</code> returns the global variable <VAR>x</VAR>.newline
+<code>getRobAStBaseOption(x)</code> returns the global variable <VAR>x</VAR>.newline
+<code>RobAStBaseOptions(x=y)</code> sets the value of the global variable <VAR>x</VAR> to <VAR>y</VAR>.</p>
+
+<h3>Global Options</h3>
+
+<dl>
+<dt>kStepUseLast:</dt><dd>The default value of argument <code>kStepUseLast</code> is
+<code>FALSE</code>. Explicitly setting <code>kStepUseLast</code> to <code>TRUE</code> should
+be done with care as in this situation the influence curve in case of
+<code>oneStepEstimator</code> and <code>kStepEstimator</code> is re-computed using
+the value of the one- resp. k-step estimate which may take quite a long
+time depending on the model. </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('base', 'options.html')" style="text-decoration: underline; color: blue; cursor: hand">options</a></code>, <code><a onclick="findlink('base', 'getOption.html')" style="text-decoration: underline; color: blue; cursor: hand">getOption</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+RobAStBaseOptions()
+RobAStBaseOptions("kStepUseLast")
+RobAStBaseOptions("kStepUseLast" = TRUE)
+# or
+RobAStBaseOptions(kStepUseLast = 1e-6)
+getRobAStBaseOption("kStepUseLast")
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/cutoff-class.html
===================================================================
--- pkg/RobAStBase/chm/cutoff-class.html (rev 0)
+++ pkg/RobAStBase/chm/cutoff-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,110 @@
+<html><head><title>Cutoff class for distance-distance plots</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>cutoff-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: cutoff-class">
+<param name="keyword" value="R: cutoff.quantile<-,cutoff-method">
+<param name="keyword" value="R: cutoff.quantile<-">
+<param name="keyword" value="R: cutoff.quantile,cutoff-method">
+<param name="keyword" value="R: cutoff.quantile">
+<param name="keyword" value="R: name,cutoff-method">
+<param name="keyword" value="R: fct,cutoff-method">
+<param name="keyword" value=" Cutoff class for distance-distance plots">
+</object>
+
+
+<h2>Cutoff class for distance-distance plots</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of methods to determine cutoff point for distance-distance plots;
+used to derive other cutoff methods later by method dispatch.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects could in principle be created by calls of the form <code>new("cutoff", ...)</code>.
+More frequently they are created via the generating function
+<code><a href="cutoff.html">cutoff</a></code>, respectively via the helper functions
+<code>cutoff.sememp</code> and <code>cutoff.chisq</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>object of class <code>"character"</code>; defaults to <code>"empirical"</code>
+in prototype; </dd>
+
+
+<dt><code>fct</code>:</dt><dd>an object of of class <code>"function"</code>;
+for this class layer, this function
+must only have one argument <code>data</code> (which may
+but need not be used to determine the cutoff point empirically);
+in derived classes this restriction could be dropped, if corresponding
+special methods for <code>ddPlot</code> are derived. Defaults to
+<code>function(data) quantile(data)</code>.</dd>
+
+
+<dt><code>cutoff.quantile</code>:</dt><dd>Object of class <code>"numeric"</code>:
+a quantile (empirical or theoretical) to plot the cutoff line; defaults
+to <code>0.95</code> in prototype;</dd>
+
+<p>
+</dl>
+
+<h3>Methods</h3>
+
+<dl>
+</p>
+
+<dt>cutoff.quantile</dt><dd><code>signature(object = "cutoff")</code>:
+accessor function for slot <code>cutoff.quantile</code>. </dd>
+
+
+<dt>cutoff.quantile<-</dt><dd><code>signature(object = "cutoff")</code>:
+replacement function for slot <code>cutoff.quantile</code>. </dd>
+
+
+<dt>fct</dt><dd><code>signature(object = "cutoff")</code>:
+accessor function for slot <code>fct</code>. </dd>
+
+
+<dt>name</dt><dd><code>signature(object = "cutoff")</code>:
+accessor function for slot <code>name</code>. </dd>
+
+<p>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="ddPlot-methods.html">ddPlot</a></code>, <code><a href="outlyingPlotIC.html">outlyingPlotIC</a></code> <code><a href="cutoff.html">cutoff</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+cutoff()
+</pre>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/cutoff.html
===================================================================
--- pkg/RobAStBase/chm/cutoff.html (rev 0)
+++ pkg/RobAStBase/chm/cutoff.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,122 @@
+<html><head><title>Generating function(s) for class 'cutoff'</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>cutoff(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: cutoff">
+<param name="keyword" value="R: cutoff.sememp">
+<param name="keyword" value="R: cutoff.chisq">
+<param name="keyword" value=" Generating function(s) for class 'cutoff'">
+</object>
+
+
+<h2>Generating function(s) for class 'cutoff'</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Generating function(s) for class <code>cutoff</code>.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+cutoff(name = "empirical", body.fct0,
+ cutoff.quantile = 0.95,
+ norm = NormType(), QF, nsim = 100000)
+cutoff.sememp()
+cutoff.chisq()
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>name</code></td>
+<td>
+argument for name slot of <code>cutoff</code> object</td></tr>
+<tr valign="top"><td><code>body.fct0</code></td>
+<td>
+a call generated by code wrapped to <code>substitute</code> resp. <code>quote</code>;
+the body of the <code>fct</code> slot of the <code>cutoff</code> object</td></tr>
+<tr valign="top"><td><code>cutoff.quantile</code></td>
+<td>
+numeric; the corresponding slot value for the <code>cutoff</code> object</td></tr>
+<tr valign="top"><td><code>norm</code></td>
+<td>
+an object of class <code>NormType</code> – the norm/distance by which to produce
+the cutoff - value.</td></tr>
+<tr valign="top"><td><code>nsim</code></td>
+<td>
+integer: the sample size used for determining the quantiles
+of <i>x'Qx^.5</i> for <i>x</i> multivariate standard
+normal and <i>Q</i> a corresponding quadratic form</td></tr>
+<tr valign="top"><td><code>QF</code></td>
+<td>
+a quadratic (positive semidefinite, symmetric) matrix used
+as quadratic form </td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+<code>cutoff</code> generates a valid object of class <code>"cutoff"</code>.
+As function slot <code>fct</code> may only have a formal argument <code>data</code>,
+the other arguments to determine the cutoff value, i.e. <code>norm</code>, <code>QF</code>, <code>nsim</code>,
+<code>cutoff.quantile</code>, <code>nsim</code> have to enter the scope of this
+function by lexical scoping; now <code>cutoff.quantile</code>, <code>norm</code>, <code>QF</code>
+are to be taken from the calling environment (not from the defining one),
+so we have delay evaluation of the function body, which is why we assume
+it to be given wrapped into <code>substitute</code> resp. <code>quote</code>.
+<code>body.fct0</code> is by default (i.e. if argument <code>body.fct0</code> is missing) set
+to <code>quote(quantile(slot(norm,"fct")(data), cutoff.quantile))</code>, internally,
+i.e.; to an empirical quantile of the corresponding norms.<br>
+</p>
+<p>
+<code>cutoff.sememp()</code> is a helper function generating the theoretical (asymptotic)
+quantile of (the square root of) a corresponding quadratic form, assuming multivariate
+normality; to determine this quantile <code>nsim</code> simulations are used.<br>
+</p>
+<p>
+<code>cutoff.chisq()</code> is a helper function generating the theoretical (asymptotic)
+quantile of (the square root of) a (self-standardized) quadratic form, assuming multivariate
+normality; i.e.; a corresponding quantile of a Chi-Square distribution.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+Object of class <code>"cutoff"</code>.</p>
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="cutoff-class.html">cutoff-class</a></code>, <code><a href="ddPlot-methods.html">ddPlot</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+cutoff()
+cutoff.sememp()
+cutoff.chisq()
+</pre>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/ddPlot-methods.html
===================================================================
--- pkg/RobAStBase/chm/ddPlot-methods.html (rev 0)
+++ pkg/RobAStBase/chm/ddPlot-methods.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,192 @@
+<html><head><title>Methods for Function ddPlot in Package ‘RobAStBase’</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>ddPlot-methods(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: ddPlot">
+<param name="keyword" value="R: ddPlot-methods">
+<param name="keyword" value="R: ddPlot,matrix-method">
+<param name="keyword" value="R: ddPlot,numeric-method">
+<param name="keyword" value="R: ddPlot,data.frame-method">
+<param name="keyword" value=" Methods for Function ddPlot in Package ‘RobAStBase’">
+</object>
+
+
+<h2>Methods for Function ddPlot in Package ‘RobAStBase’</h2>
+
+
+<h3>Description</h3>
+
+<p>
+ddPlot-methods
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+ddPlot(data, dist.x, dist.y, cutoff.x, cutoff.y, ...)
+## S4 method for signature 'matrix':
+ddPlot(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+## S4 method for signature 'numeric':
+ddPlot(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+## S4 method for signature 'data.frame':
+ddPlot(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>data</code></td>
+<td>
+data coercable to <code>matrix</code>; the data at which to produce the <code>ddPlot</code>.</td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+further arguments to be passed to <code>plot.default</code>, <code>text</code>, and <code>abline</code></td></tr>
+<tr valign="top"><td><code>dist.x</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>dist.y</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.x</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>x</code> axis
+(the vertical line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.y</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>y</code> axis
+(the horizontal line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.x</code></td>
+<td>
+numeric; the cutoff quantile for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.y</code></td>
+<td>
+numeric; the cutoff quantile for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>transform.x</code></td>
+<td>
+function; a transformation to be performed before determining the
+distances of the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>transform.y</code></td>
+<td>
+function; a transformation to be performed before determining the
+distances of the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>id.n</code></td>
+<td>
+a set of indices (or a corresponding logical vector); to select a subset
+of the data in argument <code>data</code>.</td></tr>
+<tr valign="top"><td><code>lab.pts</code></td>
+<td>
+a vector of labels for the (unsubsetted) <code>data</code>.</td></tr>
+<tr valign="top"><td><code>adj</code></td>
+<td>
+the corresponding argument for <code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for
+labelling the outliers.</td></tr>
+<tr valign="top"><td><code>cex.idn</code></td>
+<td>
+the corresponding <code>cex</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>col.idn</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>lty.cutoff</code></td>
+<td>
+the corresponding <code>lty</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>lwd.cutoff</code></td>
+<td>
+the corresponding <code>lwd</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>col.cutoff</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+The <code>matrix</code>-method calls <code>.ddPlot.MatNtNtCoCo</code>,
+the <code>numeric</code>- and <code>data.frame</code>-methods coerce argument <code>data</code>
+to <code>matrix</code> — the <code>numeric</code>-method by a call to <code>matrix(data, nrow=1)</code>,
+in the <code>data.frame</code>-methods by a call to <code>t(as.matrix(data))</code>.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+a list with items
+</p>
+<table summary="R argblock">
+<tr valign="top"><td><code>id.x</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.y</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.xy</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff and the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>qtx</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>qty</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>y</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.x.v</code></td>
+<td>
+the cutoff value in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.y.v</code></td>
+<td>
+the cutoff value in <code>y</code> direction</td></tr>
+</table>
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+MX <- matrix(rnorm(1500),nrow=6)
+QM <- matrix(rnorm(36),nrow=6); QM <- QM %*% t(QM)
+ddPlot(data=MX, dist.y=QFNorm(QuadF=PosSemDefSymmMatrix(QM)))
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/internals_ddPlot.html
===================================================================
--- pkg/RobAStBase/chm/internals_ddPlot.html (rev 0)
+++ pkg/RobAStBase/chm/internals_ddPlot.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,193 @@
+<html><head><title>Internal / Helper functions of package RobAStBase for ddPlot</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>internals_for_RobAStBase_ddPlot(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: internals_for_RobAStBase_ddPlot">
+<param name="keyword" value="R: .ddPlot.MatNtNtCoCo">
+<param name="keyword" value=" Internal / Helper functions of package RobAStBase for ddPlot">
+</object>
+
+
+<h2>Internal / Helper functions of package RobAStBase for ddPlot</h2>
+
+
+<h3>Description</h3>
+
+<p>
+This function is an internally used helper function for <code><a href="ddPlot-methods.html">ddPlot</a></code> in package <span class="pkg">RobAStBase</span>.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+.ddPlot.MatNtNtCoCo(data, ..., dist.x = NormType(), dist.y = NormType(),
+ cutoff.x = cutoff(norm = dist.x, cutoff.quantile = cutoff.quantile.x),
+ cutoff.y = cutoff(norm = dist.y, cutoff.quantile = cutoff.quantile.y),
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj =0, cex.idn = 1,
+ col.idn = par("col"), lty.cutoff,
+ lwd.cutoff, col.cutoff = "red")</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>data</code></td>
+<td>
+data in <code>matrix</code> form (columns are observations; rows are variable
+dimensions) at which to produce the <code>ddPlot</code>.</td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+further arguments to be passed to <code>plot.default</code>, <code>text</code>, and <code>abline</code></td></tr>
+<tr valign="top"><td><code>dist.x</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>dist.y</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.x</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>x</code> axis
+(the vertical line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.y</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>y</code> axis
+(the horizontal line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.x</code></td>
+<td>
+numeric; the cutoff quantile for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.y</code></td>
+<td>
+numeric; the cutoff quantile for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>transform.x</code></td>
+<td>
+function; a transformation to be performed before determining the
+distances of the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>transform.y</code></td>
+<td>
+function; a transformation to be performed before determining the
+distances of the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>id.n</code></td>
+<td>
+a set of indices (or a corresponding logical vector); to select a subset
+of the data in argument <code>data</code>.</td></tr>
+<tr valign="top"><td><code>lab.pts</code></td>
+<td>
+a vector of labels for the (unsubsetted) <code>data</code>.</td></tr>
+<tr valign="top"><td><code>adj</code></td>
+<td>
+the corresponding argument for <code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for
+labelling the outliers.</td></tr>
+<tr valign="top"><td><code>cex.idn</code></td>
+<td>
+the corresponding <code>cex</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>col.idn</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>lty.cutoff</code></td>
+<td>
+the corresponding <code>lty</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>lwd.cutoff</code></td>
+<td>
+the corresponding <code>lwd</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>col.cutoff</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+performs the plotting for <code><a href="ddPlot-methods.html">ddPlot</a></code> and <code><a href="outlyingPlotIC.html">outlyingPlotIC</a></code>;
+all arguments except for <code>data</code> are optional. In case they are missing
+default values are used as usual; for those arguments without default arguments,
+we do
+<ul>
+<dt>transform.x</dt><dd>defaults to identity, internally</dd>.
+<dt>id.n</dt><dd>defaults to <code>1:ncol(data)</code>, internally</dd>.
+<dt>lab.pts</dt><dd>defaults to <code>(1:ncol(data))[id.n]</code>, internally</dd>.
+<dt>lwd.cutoff</dt><dd>defaults to argument <code>lwd</code>, if given, else to <code>par{lwd}</code>, internally</dd>.
+<dt>lty.cutoff</dt><dd>defaults to argument <code>lty</code>, if given, else to <code>par{lty}</code>, internally</dd>.
+</ul>
+
+<h3>Value</h3>
+
+<p>
+a list with items
+</p>
+<table summary="R argblock">
+<tr valign="top"><td><code>id.x</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.y</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.xy</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff and the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>qtx</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>qty</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>y</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.x.v</code></td>
+<td>
+the cutoff value in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.y.v</code></td>
+<td>
+the cutoff value in <code>y</code> direction</td></tr>
+</table>
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a onclick="findlink('graphics', 'plot.default.html')" style="text-decoration: underline; color: blue; cursor: hand">plot.default</a></code>, <code><a onclick="findlink('graphics', 'par.html')" style="text-decoration: underline; color: blue; cursor: hand">par</a></code>,
+<code><a href="ddPlot-methods.html">ddPlot</a></code>, <code><a href="outlyingPlotIC.html">outlyingPlotIC</a></code>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+MX <- matrix(rnorm(1500),nrow=6)
+QM <- matrix(rnorm(36),nrow=6); QM <- QM %*% t(QM)
+RobAStBase:::.ddPlot.MatNtNtCoCo(data=MX,
+ dist.y=QFNorm(QuadF=PosSemDefSymmMatrix(QM)),
+ xlab="Norm.x",ylab="Norm.y", cex.idn = 1.3, offset=0,
+ lwd=2, lwd.cutoff=4, lty=2, col.cutoff =2, col.idn="green",
+ col = "blue", adj=0.4, pos=4,id.n = sample(1:200,size=100),
+ lab.pts=letters,log="x", main="GA", sub="NO",cex.sub=0.2)
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/kStepEstimate-class.html
===================================================================
--- pkg/RobAStBase/chm/kStepEstimate-class.html (rev 0)
+++ pkg/RobAStBase/chm/kStepEstimate-class.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,94 @@
+<html><head><title>kStepEstimate-class.</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>kStepEstimate-class(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: kStepEstimate-class">
+<param name="keyword" value="R: steps">
+<param name="keyword" value="R: steps,kStepEstimate-method">
+<param name="keyword" value="R: show,kStepEstimate-method">
+<param name="keyword" value=" kStepEstimate-class.">
+</object>
+
+
+<h2>kStepEstimate-class.</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Class of asymptotically linear estimates.
+</p>
+
+
+<h3>Objects from the Class</h3>
+
+<p>
+Objects can be created by calls of the form <code>new("kStepEstimate", ...)</code>.
+More frequently they are created via the generating function
+<code>kStepEstimator</code>.
+</p>
+
+
+<h3>Slots</h3>
+
+<dl>
+<dt><code>name</code>:</dt><dd>Object of class <code>"character"</code>:
+name of the estimator. </dd>
+<dt><code>estimate</code>:</dt><dd>Object of class <code>"ANY"</code>:
+estimate. </dd>
+<dt><code>samplesize</code>:</dt><dd>Object of class <code>"numeric"</code>:
+sample size. </dd>
+<dt><code>asvar</code>:</dt><dd>Optional object of class <code>"matrix"</code>:
+asymptotic variance. </dd>
+<dt><code>asbias</code>:</dt><dd>Optional object of class <code>"numeric"</code>:
+asymptotic bias. </dd>
+<dt><code>pIC</code>:</dt><dd>Optional object of class <code>InfluenceCurve</code>:
+influence curve. </dd>
+<dt><code>nuis.idx</code>:</dt><dd>object of class <code>"OptionalNumeric"</code>:
+indices of <code>estimate</code> belonging to the nuisance part. </dd>
+<dt><code>steps</code>:</dt><dd>Object of class <code>"integer"</code>: number
+of steps. </dd>
+<dt><code>Infos</code>:</dt><dd>object of class <code>"matrix"</code>
+with two columns named <code>method</code> and <code>message</code>:
+additional informations. </dd>
+</dl>
+
+<h3>Extends</h3>
+
+<p>
+Class <code>"ALEstimate"</code>, directly.<br>
+Class <code>"Estimate"</code>, by class <code>"ALEstimate"</code>
+</p>
+
+
+<h3>Methods</h3>
+
+<dl>
+<dt>steps</dt><dd><code>signature(object = "kStepEstimate")</code>:
+accessor function for slot <code>steps</code>. </dd>
+
+
+<dt>show</dt><dd><code>signature(object = "kStepEstimate")</code> </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="ALEstimate-class.html">ALEstimate-class</a></code>
+</p>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/kStepEstimator.html
===================================================================
--- pkg/RobAStBase/chm/kStepEstimator.html (rev 0)
+++ pkg/RobAStBase/chm/kStepEstimator.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,144 @@
+<html><head><title>Generic function for the computation of k-step estimates</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>kStepEstimator(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: kStepEstimator">
+<param name="keyword" value="R: kStepEstimator-methods">
+<param name="keyword" value="R: kStepEstimator,numeric,IC,numeric-method">
+<param name="keyword" value="R: kStepEstimator,matrix,IC,numeric-method">
+<param name="keyword" value="R: kStepEstimator,numeric,IC,Estimate-method">
+<param name="keyword" value="R: kStepEstimator,matrix,IC,Estimate-method">
+<param name="keyword" value=" Generic function for the computation of k-step estimates">
+</object>
+
+
+<h2>Generic function for the computation of k-step estimates</h2>
+
+
+<h3>Description</h3>
+
+<p>
+Generic function for the computation of k-step estimates.
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+kStepEstimator(x, IC, start, ...)
+
+## S4 method for signature 'numeric, IC, numeric':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'matrix, IC, numeric':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'numeric, IC, Estimate':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+## S4 method for signature 'matrix, IC, Estimate':
+kStepEstimator(x, IC, start, steps = 1L,
+ useLast = getRobAStBaseOption("kStepUseLast"))
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>x</code></td>
+<td>
+sample </td></tr>
+<tr valign="top"><td><code>IC</code></td>
+<td>
+object of class <code>"IC"</code> </td></tr>
+<tr valign="top"><td><code>start</code></td>
+<td>
+initial estimate </td></tr>
+<tr valign="top"><td><code>steps</code></td>
+<td>
+integer: number of steps </td></tr>
+<tr valign="top"><td><code>useLast</code></td>
+<td>
+which parameter estimate (initial estimate or
+k-step estimate) shall be used to fill the slots <code>pIC</code>,
+<code>asvar</code> and <code>asbias</code> of the return value. </td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+additional parameters </td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+Given an initial estimation <code>start</code>, a sample <code>x</code>
+and an influence curve <code>IC</code> the corresponding k-step
+estimator is computed.
+</p>
+<p>
+The default value of argument <code>useLast</code> is set by the
+global option <code>kStepUseLast</code> which by default is set to
+<code>FALSE</code>. In case of general models <code>useLast</code>
+remains unchanged during the computations. However, if
+slot <code>CallL2Fam</code> of <code>IC</code> generates an object of
+class <code>"L2GroupParamFamily"</code> the value of <code>useLast</code>
+is changed to <code>TRUE</code>.
+Explicitly setting <code>useLast</code> to <code>TRUE</code> should
+be done with care as in this situation the influence curve
+is re-computed using the value of the one-step estimate
+which may take quite a long time depending on the model.
+</p>
+<p>
+If <code>useLast</code> is set to <code>TRUE</code> and slot <code>modifyIC</code>
+of <code>IC</code> is filled with some function (which can be
+used to re-compute the IC for a different parameter), the
+computation of <code>asvar</code>, <code>asbias</code> and <code>IC</code> is
+based on the k-step estimate.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+Object of class <code>"kStepEstimate"</code>.</p>
+
+<h3>Methods</h3>
+
+<dl>
+<dt>x = "numeric", IC = "IC", start = "numeric"</dt><dd>univariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "numeric"</dt><dd>multivariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "Estimate"</dt><dd>multivariate samples. </dd>
+<dt>x = "matrix", IC = "IC", start = "Estimate"</dt><dd>multivariate samples. </dd>
+</dl>
+
+<h3>Author(s)</h3>
+
+<p>
+Matthias Kohl <a href="mailto:Matthias.Kohl at stamats.de">Matthias.Kohl at stamats.de</a>
+</p>
+
+
+<h3>References</h3>
+
+<p>
+Rieder, H. (1994) <EM>Robust Asymptotic Statistics</EM>. New York: Springer.
+</p>
+<p>
+Kohl, M. (2005) <EM>Numerical Contributions to the Asymptotic Theory of Robustness</EM>.
+Bayreuth: Dissertation.
+</p>
+
+
+<h3>See Also</h3>
+
+<p>
+<code><a href="IC-class.html">IC-class</a></code>, <code><a href="kStepEstimate-class.html">kStepEstimate-class</a></code>
+</p>
+
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/outlyingPlotIC.html
===================================================================
--- pkg/RobAStBase/chm/outlyingPlotIC.html (rev 0)
+++ pkg/RobAStBase/chm/outlyingPlotIC.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,182 @@
+<html><head><title>Function outlyingPlotIC in Package ‘RobAStBase’</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>outlyingPlotIC(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: outlyingPlotIC">
+<param name="keyword" value=" Function outlyingPlotIC in Package ‘RobAStBase’">
+</object>
+
+
+<h2>Function outlyingPlotIC in Package ‘RobAStBase’</h2>
+
+
+<h3>Description</h3>
+
+<p>
+outlyingPlotIC produces an outlyingness plot based on distances applied
+to ICs
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>
+outlyingPlotIC(data, IC.x, IC.y = IC.x, dist.x = NormType(),
+ dist.y, cutoff.y = cutoff.chisq(), cutoff.x = cutoff.sememp(), ...,
+ cutoff.quantile.x = 0.95,
+ cutoff.quantile.y = cutoff.quantile.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff,
+ main = gettext("Outlyingness by means of a distance-distance plot")
+ )
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>data</code></td>
+<td>
+data coercable to <code>matrix</code>; the data at which to produce the <code>ddPlot</code>.</td></tr>
+<tr valign="top"><td><code>IC.x</code></td>
+<td>
+object of class <code>IC</code> the influence curve to produce
+the distances for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>IC.y</code></td>
+<td>
+object of class <code>IC</code> the influence curve to produce
+the distances for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+further arguments to be passed to <code>plot.default</code>, <code>text</code>, and <code>abline</code></td></tr>
+<tr valign="top"><td><code>dist.x</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>dist.y</code></td>
+<td>
+object of class <code>NormType</code>; the distance for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.x</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>x</code> axis
+(the vertical line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.y</code></td>
+<td>
+object of class <code>cutoff</code>; the cutoff information for the <code>y</code> axis
+(the horizontal line discriminating 'good' and 'bad' points).</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.x</code></td>
+<td>
+numeric; the cutoff quantile for the <code>x</code> axis.</td></tr>
+<tr valign="top"><td><code>cutoff.quantile.y</code></td>
+<td>
+numeric; the cutoff quantile for the <code>y</code> axis.</td></tr>
+<tr valign="top"><td><code>id.n</code></td>
+<td>
+a set of indices (or a corresponding logical vector); to select a subset
+of the data in argument <code>data</code>.</td></tr>
+<tr valign="top"><td><code>lab.pts</code></td>
+<td>
+a vector of labels for the (unsubsetted) <code>data</code>.</td></tr>
+<tr valign="top"><td><code>adj</code></td>
+<td>
+the corresponding argument for <code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for
+labelling the outliers.</td></tr>
+<tr valign="top"><td><code>cex.idn</code></td>
+<td>
+the corresponding <code>cex</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>col.idn</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'text.html')" style="text-decoration: underline; color: blue; cursor: hand">text</a></code> for labelling the outliers.</td></tr>
+<tr valign="top"><td><code>lty.cutoff</code></td>
+<td>
+the corresponding <code>lty</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>lwd.cutoff</code></td>
+<td>
+the corresponding <code>lwd</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>col.cutoff</code></td>
+<td>
+the corresponding <code>col</code> argument for
+<code><a onclick="findlink('graphics', 'abline.html')" style="text-decoration: underline; color: blue; cursor: hand">abline</a></code> for drawing the cutoff lines.</td></tr>
+<tr valign="top"><td><code>main</code></td>
+<td>
+the main title.</td></tr>
+</table>
+
+<h3>Details</h3>
+
+<p>
+calls a corresponding <code><a href="ddPlot-methods.html">ddPlot</a></code> method to produce the plot.
+</p>
+
+
+<h3>Value</h3>
+
+<p>
+a list with items
+</p>
+<table summary="R argblock">
+<tr valign="top"><td><code>id.x</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.y</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>id.xy</code></td>
+<td>
+the indices of (possibly transformed) data (within subset <code>id.n</code>) beyond the <code>x</code>-cutoff and the <code>y</code>-cutoff</td></tr>
+<tr valign="top"><td><code>qtx</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>qty</code></td>
+<td>
+the quantiles of the distances of the (possibly transformed) data in <code>y</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.x.v</code></td>
+<td>
+the cutoff value in <code>x</code> direction</td></tr>
+<tr valign="top"><td><code>cutoff.y.v</code></td>
+<td>
+the cutoff value in <code>y</code> direction</td></tr>
+</table>
+
+<h3>Author(s)</h3>
+
+<p>
+Peter Ruckdeschel <a href="mailto:Peter.Ruckdeschel at itwm.fraunhofer.de">Peter.Ruckdeschel at itwm.fraunhofer.de</a>
+</p>
+
+
+<h3>Examples</h3>
+
+<pre>
+if(require(ROptEst)){
+## generates normal location and scale family with mean = -2 and sd = 3
+N0 <- NormLocationScaleFamily()
+N0.IC0 <- optIC(model = N0, risk = asCov())
+N0.Rob1 <- InfRobModel(center = N0, neighbor = ContNeighborhood(radius = 0.5))
+N0.IC1 <- optIC(model = N0.Rob1, risk = asMSE())
+xn <- c(rnorm(100),rcauchy(20)+20)
+outlyingPlotIC(xn, IC.x=N0.IC0)
+outlyingPlotIC(xn, IC.x=N0.IC1)
+}
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.2 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Added: pkg/RobAStBase/chm/plot-methods.html
===================================================================
--- pkg/RobAStBase/chm/plot-methods.html (rev 0)
+++ pkg/RobAStBase/chm/plot-methods.html 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,134 @@
+<html><head><title>Methods for Function plot in Package ‘RobAStBase’</title>
+<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
+<link rel="stylesheet" type="text/css" href="Rchm.css">
+</head>
+<body>
+
+<table width="100%"><tr><td>plot-methods(RobAStBase)</td><td align="right">R Documentation</td></tr></table><object type="application/x-oleobject" classid="clsid:1e2a7bd0-dab9-11d0-b93a-00c04fc99f9e">
+<param name="keyword" value="R: plot">
+<param name="keyword" value="R: plot-methods">
+<param name="keyword" value="R: plot,IC,missing-method">
+<param name="keyword" value=" Methods for Function plot in Package ‘RobAStBase’">
+</object>
+
+
+<h2>Methods for Function plot in Package ‘RobAStBase’</h2>
+
+
+<h3>Description</h3>
+
+<p>
+plot-methods
+</p>
+
+
+<h3>Usage</h3>
+
+<pre>plot(x, y, ...)
+## S4 method for signature 'IC, missing':
+plot(x, ..., withSweave = getdistrOption("withSweave"),
+ main = FALSE, inner = TRUE, sub = FALSE,
+ col.inner = par("col.main"), cex.inner = 0.8,
+ bmar = par("mar")[1], tmar = par("mar")[3],
+ mfColRow = TRUE)
+</pre>
+
+
+<h3>Arguments</h3>
+
+<table summary="R argblock">
+<tr valign="top"><td><code>x</code></td>
+<td>
+object of class <code>"IC"</code>: IC to be plotted </td></tr>
+<tr valign="top"><td><code>y</code></td>
+<td>
+missing </td></tr>
+<tr valign="top"><td><code>withSweave</code></td>
+<td>
+logical: if <code>TRUE</code> (for working with <CODE>Sweave</CODE>)
+no extra device is opened</td></tr>
+<tr valign="top"><td><code>main</code></td>
+<td>
+logical: is a main title to be used? or <br>
+just as argument <code>main</code> in <code><a onclick="findlink('graphics', 'plotdefault.html')" style="text-decoration: underline; color: blue; cursor: hand">plot.default</a></code>.</td></tr>
+<tr valign="top"><td><code>inner</code></td>
+<td>
+logical: panels have their own titles? or <br>
+character vector of / cast to length number of comparands:
+<code>main</code> in <code><a onclick="findlink('graphics', 'plotdefault.html')" style="text-decoration: underline; color: blue; cursor: hand">plot.default</a></code>)</td></tr>
+<tr valign="top"><td><code>sub</code></td>
+<td>
+logical: is a sub-title to be used? or <br>
+just as argument <code>sub</code> in <code><a onclick="findlink('graphics', 'plotdefault.html')" style="text-decoration: underline; color: blue; cursor: hand">plot.default</a></code>.</td></tr>
+<tr valign="top"><td><code>tmar</code></td>
+<td>
+top margin – useful for non-standard main title sizes</td></tr>
+<tr valign="top"><td><code>bmar</code></td>
+<td>
+bottom margin – useful for non-standard sub title sizes</td></tr>
+<tr valign="top"><td><code>cex.inner</code></td>
+<td>
+magnification to be used for inner titles relative
+to the current setting of <code>cex</code>; as in
+<code><a onclick="findlink('stats', 'par.html')" style="text-decoration: underline; color: blue; cursor: hand">par</a></code></td></tr>
+<tr valign="top"><td><code>col.inner</code></td>
+<td>
+character or integer code; color for the inner title</td></tr>
+<tr valign="top"><td><code>mfColRow</code></td>
+<td>
+shall default partition in panels be used — defaults to <code>TRUE</code></td></tr>
+<tr valign="top"><td><code>...</code></td>
+<td>
+</td></tr>
+</table>
+<p>
+ {further parameters for <code>plot</code>}</p>
+
+<h3>Details</h3>
+
+<p>
+Any parameters of <code>plot.default</code> may be passed on to this particular
+<code>plot</code> method.
+</p>
+<p>
+For main-, inner, and subtitles given as arguments <code>main</code>,
+<code>inner</code>, and <code>sub</code>, top and bottom margins are enlarged to 5 resp.
+6 by default but may also be specified by <code>tmar</code> / <code>bmar</code> arguments.
+If <code>main</code> / <code>inner</code> / <code>sub</code> are
+logical then if the respective argument is <code>FALSE</code> nothing is done/plotted,
+but if it is <code>TRUE</code>, we use a default main title taking up the calling
+arguments in case of <code>main</code>, default inner titles taking up the
+class and (named) parameter slots of arguments in case of <code>inner</code>,
+and a "generated on <data>"-tag in case of <code>sub</code>.
+Of course, if <code>main</code> / <code>inner</code> / <code>sub</code> are <code>character</code>, this
+is used for the title; in case of <code>inner</code> it is then checked whether it
+has correct length. In all title arguments, the following patterns are substituted:
+<ul>
+<dt><code>"%C"</code></dt><dd>class of argument <code>object</code></dd>
+<dt><code>"%A"</code></dt><dd>deparsed argument <code>object</code></dd>
+<dt><code>"%D"</code></dt><dd>time/date-string when the plot was generated</dd>
+</ul>
+
+<h3>Examples</h3>
+
+<pre>
+IC1 <- new("IC")
+plot(IC1)
+plot(IC1, main = TRUE, panel.first= grid(),
+ col = "blue", cex.main = 2, cex.inner = 1)
+</pre>
+
+<script Language="JScript">
+function findlink(pkg, fn) {
+var Y, link;
+Y = location.href.lastIndexOf("\\") + 1;
+link = location.href.substring(0, Y);
+link = link + "../../" + pkg + "/chtml/" + pkg + ".chm::/" + fn;
+location.href = link;
+}
+</script>
+
+
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.3 <a href="00Index.html">Index]</a></div>
+
+</body></html>
Modified: pkg/RobAStBase/man/IC-class.Rd
===================================================================
--- pkg/RobAStBase/man/IC-class.Rd 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobAStBase/man/IC-class.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -11,7 +11,6 @@
\alias{checkIC,IC,L2ParamFamily-method}
\alias{evalIC,IC,numeric-method}
\alias{evalIC,IC,matrix-method}
-\alias{plot,IC,missing-method}
\alias{show,IC-method}
\title{Influence curve}
Added: pkg/RobAStBase/man/cutoff-class.Rd
===================================================================
--- pkg/RobAStBase/man/cutoff-class.Rd (rev 0)
+++ pkg/RobAStBase/man/cutoff-class.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,64 @@
+\name{cutoff-class}
+\docType{class}
+\alias{cutoff-class}
+\alias{cutoff.quantile<-,cutoff-method}
+\alias{cutoff.quantile<-}
+\alias{cutoff.quantile,cutoff-method}
+\alias{cutoff.quantile}
+\alias{name,cutoff-method}
+\alias{fct,cutoff-method}
+
+\title{Cutoff class for distance-distance plots}
+\description{Class of methods to determine cutoff point for distance-distance plots;
+ used to derive other cutoff methods later by method dispatch.
+}
+\section{Objects from the Class}{
+ Objects could in principle be created by calls of the form \code{new("cutoff", ...)}.
+ More frequently they are created via the generating function
+ \code{\link{cutoff}}, respectively via the helper functions
+ \code{cutoff.sememp} and \code{cutoff.chisq}.
+}
+\section{Slots}{
+ \describe{
+ \item{\code{name}:}{ object of class \code{"character"}; defaults to \code{"empirical"}
+ in prototype; }
+
+ \item{\code{fct}:}{ an object of of class \code{"function"};
+ for this class layer, this function
+ must only have one argument \code{data} (which may
+ but need not be used to determine the cutoff point empirically);
+ in derived classes this restriction could be dropped, if corresponding
+ special methods for \code{ddPlot} are derived. Defaults to
+ \code{function(data) quantile(data)}.}
+
+ \item{\code{cutoff.quantile}:}{ Object of class \code{"numeric"}:
+ a quantile (empirical or theoretical) to plot the cutoff line; defaults
+ to \code{0.95} in prototype;}
+
+ }
+}
+\section{Methods}{
+ \describe{
+
+ \item{cutoff.quantile}{\code{signature(object = "cutoff")}:
+ accessor function for slot \code{cutoff.quantile}. }
+
+ \item{cutoff.quantile<-}{\code{signature(object = "cutoff")}:
+ replacement function for slot \code{cutoff.quantile}. }
+
+ \item{fct}{\code{signature(object = "cutoff")}:
+ accessor function for slot \code{fct}. }
+
+ \item{name}{\code{signature(object = "cutoff")}:
+ accessor function for slot \code{name}. }
+
+
+ }
+}
+\author{Peter Ruckdeschel \email{Peter.Ruckdeschel at itwm.fraunhofer.de}}
+%\note{}
+\seealso{\code{\link{ddPlot}}, \code{\link{outlyingPlotIC}} \code{\link{cutoff}}}
+\examples{
+cutoff()
+}
+\keyword{classes}
Added: pkg/RobAStBase/man/cutoff.Rd
===================================================================
--- pkg/RobAStBase/man/cutoff.Rd (rev 0)
+++ pkg/RobAStBase/man/cutoff.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,62 @@
+\name{cutoff}
+\alias{cutoff}
+\alias{cutoff.sememp}
+\alias{cutoff.chisq}
+
+\title{Generating function(s) for class 'cutoff'}
+\description{
+ Generating function(s) for class \code{cutoff}.
+}
+\usage{
+cutoff(name = "empirical", body.fct0,
+ cutoff.quantile = 0.95,
+ norm = NormType(), QF, nsim = 100000)
+cutoff.sememp()
+cutoff.chisq()
+}
+\arguments{
+ \item{name}{argument for name slot of \code{cutoff} object}
+ \item{body.fct0}{a call generated by code wrapped to \code{substitute} resp. \code{quote};
+ the body of the \code{fct} slot of the \code{cutoff} object}
+ \item{cutoff.quantile}{numeric; the corresponding slot value for the \code{cutoff} object}
+ \item{norm}{an object of class \code{NormType} -- the norm/distance by which to produce
+ the cutoff - value.}
+ \item{nsim}{ integer: the sample size used for determining the quantiles
+ of \eqn{(x^\tau Q x)^{1/2}}{x'Qx\^.5} for \eqn{x} multivariate standard
+ normal and \eqn{Q} a corresponding quadratic form}
+ \item{QF}{ a quadratic (positive semidefinite, symmetric) matrix used
+ as quadratic form }
+}
+\details{
+\code{cutoff} generates a valid object of class \code{"cutoff"}.
+As function slot \code{fct} may only have a formal argument \code{data},
+the other arguments to determine the cutoff value, i.e. \code{norm}, \code{QF}, \code{nsim},
+\code{cutoff.quantile}, \code{nsim} have to enter the scope of this
+function by lexical scoping; now \code{cutoff.quantile}, \code{norm}, \code{QF}
+are to be taken from the calling environment (not from the defining one),
+so we have delay evaluation of the function body, which is why we assume
+it to be given wrapped into \code{substitute} resp. \code{quote}.
+\code{body.fct0} is by default (i.e. if argument \code{body.fct0} is missing) set
+to \code{quote(quantile(slot(norm,"fct")(data), cutoff.quantile))}, internally,
+i.e.; to an empirical quantile of the corresponding norms.\cr
+
+\code{cutoff.sememp()} is a helper function generating the theoretical (asymptotic)
+quantile of (the square root of) a corresponding quadratic form, assuming multivariate
+normality; to determine this quantile \code{nsim} simulations are used.\cr
+
+\code{cutoff.chisq()} is a helper function generating the theoretical (asymptotic)
+quantile of (the square root of) a (self-standardized) quadratic form, assuming multivariate
+normality; i.e.; a corresponding quantile of a Chi-Square distribution.
+}
+\value{Object of class \code{"cutoff"}.}
+\author{
+ Peter Ruckdeschel \email{Peter.Ruckdeschel at itwm.fraunhofer.de}
+ }
+%\note{}
+\seealso{\code{\link{cutoff-class}}, \code{\link{ddPlot}} }
+\examples{
+cutoff()
+cutoff.sememp()
+cutoff.chisq()
+}
+\keyword{hplot}
Added: pkg/RobAStBase/man/ddPlot-methods.Rd
===================================================================
--- pkg/RobAStBase/man/ddPlot-methods.Rd (rev 0)
+++ pkg/RobAStBase/man/ddPlot-methods.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,89 @@
+\name{ddPlot-methods}
+\docType{methods}
+\alias{ddPlot}
+\alias{ddPlot-methods}
+\alias{ddPlot,matrix-method}
+\alias{ddPlot,numeric-method}
+\alias{ddPlot,data.frame-method}
+\title{ Methods for Function ddPlot in Package `RobAStBase' }
+\description{ddPlot-methods}
+\usage{
+ddPlot(data, dist.x, dist.y, cutoff.x, cutoff.y, ...)
+\S4method{ddPlot}{matrix}(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+\S4method{ddPlot}{numeric}(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+\S4method{ddPlot}{data.frame}(data, dist.x = NormType(), dist.y = NormType(),
+ cutoff.x, cutoff.y, ...,
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff)
+}
+\arguments{
+ \item{data}{data coercable to \code{matrix}; the data at which to produce the \code{ddPlot}.}
+ \item{\dots}{further arguments to be passed to \code{plot.default}, \code{text}, and \code{abline}}
+ \item{dist.x}{object of class \code{NormType}; the distance for the \code{x} axis.}
+ \item{dist.y}{object of class \code{NormType}; the distance for the \code{y} axis.}
+ \item{cutoff.x}{object of class \code{cutoff}; the cutoff information for the \code{x} axis
+ (the vertical line discriminating 'good' and 'bad' points).}
+ \item{cutoff.y}{object of class \code{cutoff}; the cutoff information for the \code{y} axis
+ (the horizontal line discriminating 'good' and 'bad' points).}
+ \item{cutoff.quantile.x}{numeric; the cutoff quantile for the \code{x} axis.}
+ \item{cutoff.quantile.y}{numeric; the cutoff quantile for the \code{y} axis.}
+ \item{transform.x}{function; a transformation to be performed before determining the
+ distances of the \code{x} axis.}
+ \item{transform.y}{function; a transformation to be performed before determining the
+ distances of the \code{y} axis.}
+ \item{id.n}{a set of indices (or a corresponding logical vector); to select a subset
+ of the data in argument \code{data}.}
+ \item{lab.pts}{a vector of labels for the (unsubsetted) \code{data}.}
+ \item{adj}{the corresponding argument for \code{\link[graphics]{text}} for
+ labelling the outliers.}
+ \item{cex.idn}{the corresponding \code{cex} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{col.idn}{the corresponding \code{col} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{lty.cutoff}{the corresponding \code{lty} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{lwd.cutoff}{the corresponding \code{lwd} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{col.cutoff}{the corresponding \code{col} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+}
+\details{
+The \code{matrix}-method calls \code{.ddPlot.MatNtNtCoCo},
+the \code{numeric}- and \code{data.frame}-methods coerce argument \code{data}
+to \code{matrix} --- the \code{numeric}-method by a call to \code{matrix(data, nrow=1)},
+in the \code{data.frame}-methods by a call to \code{t(as.matrix(data))}.
+}
+\value{
+a list with items
+\item{id.x}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff}
+\item{id.y}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{y}-cutoff}
+\item{id.xy}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff and the \code{y}-cutoff}
+\item{qtx}{the quantiles of the distances of the (possibly transformed) data in \code{x} direction}
+\item{qty}{the quantiles of the distances of the (possibly transformed) data in \code{y} direction}
+\item{cutoff.x.v}{the cutoff value in \code{x} direction}
+\item{cutoff.y.v}{the cutoff value in \code{y} direction}
+}
+\author{
+ Peter Ruckdeschel \email{Peter.Ruckdeschel at itwm.fraunhofer.de}
+ }
+
+\examples{
+MX <- matrix(rnorm(1500),nrow=6)
+QM <- matrix(rnorm(36),nrow=6); QM <- QM \%*\% t(QM)
+ddPlot(data=MX, dist.y=QFNorm(QuadF=PosSemDefSymmMatrix(QM)))
+}
+\keyword{methods}
+\keyword{hplot}
+
Added: pkg/RobAStBase/man/internals_ddPlot.Rd
===================================================================
--- pkg/RobAStBase/man/internals_ddPlot.Rd (rev 0)
+++ pkg/RobAStBase/man/internals_ddPlot.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,98 @@
+\name{internals_for_RobAStBase_ddPlot}
+\alias{internals_for_RobAStBase_ddPlot}
+\alias{.ddPlot.MatNtNtCoCo}
+
+\title{Internal / Helper functions of package RobAStBase for ddPlot}
+
+\description{
+This function is an internally used helper function for \code{\link{ddPlot}} in package \pkg{RobAStBase}.}
+
+\usage{
+.ddPlot.MatNtNtCoCo(data, ..., dist.x = NormType(), dist.y = NormType(),
+ cutoff.x = cutoff(norm = dist.x, cutoff.quantile = cutoff.quantile.x),
+ cutoff.y = cutoff(norm = dist.y, cutoff.quantile = cutoff.quantile.y),
+ cutoff.quantile.x = 0.95, cutoff.quantile.y = cutoff.quantile.x,
+ transform.x, transform.y = transform.x,
+ id.n, lab.pts, adj =0, cex.idn = 1,
+ col.idn = par("col"), lty.cutoff,
+ lwd.cutoff, col.cutoff = "red")}
+
+\arguments{
+ \item{data}{data in \code{matrix} form (columns are observations; rows are variable
+ dimensions) at which to produce the \code{ddPlot}.}
+ \item{\dots}{further arguments to be passed to \code{plot.default}, \code{text}, and \code{abline}}
+ \item{dist.x}{object of class \code{NormType}; the distance for the \code{x} axis.}
+ \item{dist.y}{object of class \code{NormType}; the distance for the \code{y} axis.}
+ \item{cutoff.x}{object of class \code{cutoff}; the cutoff information for the \code{x} axis
+ (the vertical line discriminating 'good' and 'bad' points).}
+ \item{cutoff.y}{object of class \code{cutoff}; the cutoff information for the \code{y} axis
+ (the horizontal line discriminating 'good' and 'bad' points).}
+ \item{cutoff.quantile.x}{numeric; the cutoff quantile for the \code{x} axis.}
+ \item{cutoff.quantile.y}{numeric; the cutoff quantile for the \code{y} axis.}
+ \item{transform.x}{function; a transformation to be performed before determining the
+ distances of the \code{x} axis.}
+ \item{transform.y}{function; a transformation to be performed before determining the
+ distances of the \code{y} axis.}
+ \item{id.n}{a set of indices (or a corresponding logical vector); to select a subset
+ of the data in argument \code{data}.}
+ \item{lab.pts}{a vector of labels for the (unsubsetted) \code{data}.}
+ \item{adj}{the corresponding argument for \code{\link[graphics]{text}} for
+ labelling the outliers.}
+ \item{cex.idn}{the corresponding \code{cex} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{col.idn}{the corresponding \code{col} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{lty.cutoff}{the corresponding \code{lty} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{lwd.cutoff}{the corresponding \code{lwd} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{col.cutoff}{the corresponding \code{col} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+}
+
+\details{
+performs the plotting for \code{\link{ddPlot}} and \code{\link{outlyingPlotIC}};
+all arguments except for \code{data} are optional. In case they are missing
+default values are used as usual; for those arguments without default arguments,
+we do
+\itemize{
+\item{transform.x}{defaults to identity, internally}.
+\item{id.n}{defaults to \code{1:ncol(data)}, internally}.
+\item{lab.pts}{defaults to \code{(1:ncol(data))[id.n]}, internally}.
+\item{lwd.cutoff}{defaults to argument \code{lwd}, if given, else to \code{par{lwd}}, internally}.
+\item{lty.cutoff}{defaults to argument \code{lty}, if given, else to \code{par{lty}}, internally}.
+}
+}
+
+
+\value{
+a list with items
+\item{id.x}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff}
+\item{id.y}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{y}-cutoff}
+\item{id.xy}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff and the \code{y}-cutoff}
+\item{qtx}{the quantiles of the distances of the (possibly transformed) data in \code{x} direction}
+\item{qty}{the quantiles of the distances of the (possibly transformed) data in \code{y} direction}
+\item{cutoff.x.v}{the cutoff value in \code{x} direction}
+\item{cutoff.y.v}{the cutoff value in \code{y} direction}
+}
+
+
+\author{
+ Peter Ruckdeschel \email{Peter.Ruckdeschel at itwm.fraunhofer.de}
+ }
+\seealso{\code{\link[graphics]{plot.default}}, \code{\link[graphics]{par}},
+ \code{\link{ddPlot}}, \code{\link{outlyingPlotIC}}}
+
+\examples{
+MX <- matrix(rnorm(1500),nrow=6)
+QM <- matrix(rnorm(36),nrow=6); QM <- QM \%*\% t(QM)
+RobAStBase:::.ddPlot.MatNtNtCoCo(data=MX,
+ dist.y=QFNorm(QuadF=PosSemDefSymmMatrix(QM)),
+ xlab="Norm.x",ylab="Norm.y", cex.idn = 1.3, offset=0,
+ lwd=2, lwd.cutoff=4, lty=2, col.cutoff =2, col.idn="green",
+ col = "blue", adj=0.4, pos=4,id.n = sample(1:200,size=100),
+ lab.pts=letters,log="x", main="GA", sub="NO",cex.sub=0.2)
+}
+\keyword{internal}
+\concept{utilities}
+\keyword{hplot}
Added: pkg/RobAStBase/man/outlyingPlotIC.Rd
===================================================================
--- pkg/RobAStBase/man/outlyingPlotIC.Rd (rev 0)
+++ pkg/RobAStBase/man/outlyingPlotIC.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,79 @@
+\name{outlyingPlotIC}
+\alias{outlyingPlotIC}
+\title{Function outlyingPlotIC in Package `RobAStBase' }
+\description{outlyingPlotIC produces an outlyingness plot based on distances applied
+to ICs}
+\usage{
+outlyingPlotIC(data, IC.x, IC.y = IC.x, dist.x = NormType(),
+ dist.y, cutoff.y = cutoff.chisq(), cutoff.x = cutoff.sememp(), ...,
+ cutoff.quantile.x = 0.95,
+ cutoff.quantile.y = cutoff.quantile.x,
+ id.n, lab.pts, adj, cex.idn,
+ col.idn, lty.cutoff, lwd.cutoff, col.cutoff,
+ main = gettext("Outlyingness by means of a distance-distance plot")
+ )
+}
+\arguments{
+ \item{data}{data coercable to \code{matrix}; the data at which to produce the \code{ddPlot}.}
+ \item{IC.x}{object of class \code{IC} the influence curve to produce
+ the distances for the \code{x} axis.}
+ \item{IC.y}{object of class \code{IC} the influence curve to produce
+ the distances for the \code{y} axis.}
+ \item{\dots}{further arguments to be passed to \code{plot.default}, \code{text}, and \code{abline}}
+ \item{dist.x}{object of class \code{NormType}; the distance for the \code{x} axis.}
+ \item{dist.y}{object of class \code{NormType}; the distance for the \code{y} axis.}
+ \item{cutoff.x}{object of class \code{cutoff}; the cutoff information for the \code{x} axis
+ (the vertical line discriminating 'good' and 'bad' points).}
+ \item{cutoff.y}{object of class \code{cutoff}; the cutoff information for the \code{y} axis
+ (the horizontal line discriminating 'good' and 'bad' points).}
+ \item{cutoff.quantile.x}{numeric; the cutoff quantile for the \code{x} axis.}
+ \item{cutoff.quantile.y}{numeric; the cutoff quantile for the \code{y} axis.}
+ \item{id.n}{a set of indices (or a corresponding logical vector); to select a subset
+ of the data in argument \code{data}.}
+ \item{lab.pts}{a vector of labels for the (unsubsetted) \code{data}.}
+ \item{adj}{the corresponding argument for \code{\link[graphics]{text}} for
+ labelling the outliers.}
+ \item{cex.idn}{the corresponding \code{cex} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{col.idn}{the corresponding \code{col} argument for
+ \code{\link[graphics]{text}} for labelling the outliers.}
+ \item{lty.cutoff}{the corresponding \code{lty} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{lwd.cutoff}{the corresponding \code{lwd} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{col.cutoff}{the corresponding \code{col} argument for
+ \code{\link[graphics]{abline}} for drawing the cutoff lines.}
+ \item{main}{the main title.}
+}
+\details{
+calls a corresponding \code{\link{ddPlot}} method to produce the plot.
+}
+\author{
+ Peter Ruckdeschel \email{Peter.Ruckdeschel at itwm.fraunhofer.de}
+ }
+
+\value{
+a list with items
+\item{id.x}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff}
+\item{id.y}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{y}-cutoff}
+\item{id.xy}{the indices of (possibly transformed) data (within subset \code{id.n}) beyond the \code{x}-cutoff and the \code{y}-cutoff}
+\item{qtx}{the quantiles of the distances of the (possibly transformed) data in \code{x} direction}
+\item{qty}{the quantiles of the distances of the (possibly transformed) data in \code{y} direction}
+\item{cutoff.x.v}{the cutoff value in \code{x} direction}
+\item{cutoff.y.v}{the cutoff value in \code{y} direction}
+}
+
+\examples{
+if(require(ROptEst)){
+## generates normal location and scale family with mean = -2 and sd = 3
+N0 <- NormLocationScaleFamily()
+N0.IC0 <- optIC(model = N0, risk = asCov())
+N0.Rob1 <- InfRobModel(center = N0, neighbor = ContNeighborhood(radius = 0.5))
+N0.IC1 <- optIC(model = N0.Rob1, risk = asMSE())
+xn <- c(rnorm(100),rcauchy(20)+20)
+outlyingPlotIC(xn, IC.x=N0.IC0)
+outlyingPlotIC(xn, IC.x=N0.IC1)
+}
+}
+\keyword{hplot}
+
Added: pkg/RobAStBase/man/plot-methods.Rd
===================================================================
--- pkg/RobAStBase/man/plot-methods.Rd (rev 0)
+++ pkg/RobAStBase/man/plot-methods.Rd 2008-11-12 21:33:45 UTC (rev 210)
@@ -0,0 +1,66 @@
+\name{plot-methods}
+\docType{methods}
+\alias{plot}
+\alias{plot-methods}
+\alias{plot,IC,missing-method}
+\title{ Methods for Function plot in Package `RobAStBase' }
+\description{plot-methods}
+\usage{plot(x, y, ...)
+\S4method{plot}{IC,missing}(x, ..., withSweave = getdistrOption("withSweave"),
+ main = FALSE, inner = TRUE, sub = FALSE,
+ col.inner = par("col.main"), cex.inner = 0.8,
+ bmar = par("mar")[1], tmar = par("mar")[3],
+ mfColRow = TRUE)
+}
+\arguments{
+ \item{x}{object of class \code{"IC"}: IC to be plotted }
+ \item{y}{ missing }
+ \item{withSweave}{logical: if \code{TRUE} (for working with \command{Sweave})
+ no extra device is opened}
+ \item{main}{logical: is a main title to be used? or \cr
+ just as argument \code{main} in \code{\link{plot.default}}.}
+ \item{inner}{logical: panels have their own titles? or \cr
+ character vector of / cast to length number of comparands:
+ \code{main} in \code{\link{plot.default}})}
+ \item{sub}{logical: is a sub-title to be used? or \cr
+ just as argument \code{sub} in \code{\link{plot.default}}.}
+ \item{tmar}{top margin -- useful for non-standard main title sizes}
+ \item{bmar}{bottom margin -- useful for non-standard sub title sizes}
+ \item{cex.inner}{magnification to be used for inner titles relative
+ to the current setting of \code{cex}; as in
+ \code{\link[stats]{par}}}
+ \item{col.inner}{character or integer code; color for the inner title}
+ \item{mfColRow}{shall default partition in panels be used --- defaults to \code{TRUE}}
+ \item{\dots} {further parameters for \code{plot}}
+}
+\details{
+Any parameters of \code{plot.default} may be passed on to this particular
+\code{plot} method.
+
+For main-, inner, and subtitles given as arguments \code{main},
+\code{inner}, and \code{sub}, top and bottom margins are enlarged to 5 resp.
+6 by default but may also be specified by \code{tmar} / \code{bmar} arguments.
+If \code{main} / \code{inner} / \code{sub} are
+logical then if the respective argument is \code{FALSE} nothing is done/plotted,
+but if it is \code{TRUE}, we use a default main title taking up the calling
+arguments in case of \code{main}, default inner titles taking up the
+class and (named) parameter slots of arguments in case of \code{inner},
+and a "generated on <data>"-tag in case of \code{sub}.
+Of course, if \code{main} / \code{inner} / \code{sub} are \code{character}, this
+is used for the title; in case of \code{inner} it is then checked whether it
+has correct length. In all title arguments, the following patterns are substituted:
+\itemize{
+\item{\code{"\%C"}}{class of argument \code{object}}
+\item{\code{"\%A"}}{deparsed argument \code{object}}
+\item{\code{"\%D"}}{time/date-string when the plot was generated}
+}
+}
+\examples{
+IC1 <- new("IC")
+plot(IC1)
+plot(IC1, main = TRUE, panel.first= grid(),
+ col = "blue", cex.main = 2, cex.inner = 1)
+}
+\keyword{methods}
+\keyword{distribution}
+
Modified: pkg/RobLox/DESCRIPTION
===================================================================
--- pkg/RobLox/DESCRIPTION 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobLox/DESCRIPTION 2008-11-12 21:33:45 UTC (rev 210)
@@ -2,9 +2,8 @@
Version: 0.6.1
Date: 2008-08-04
Title: Optimally robust influence curves for location and scale
-Description: functions for the determination of optimally
- robust influence curves in case of normal
- location with unknown scale
+Description: functions for the determination of optimally robust
+ influence curves in case of normal location with unknown scale
Depends: R(>= 2.7.0), stats, distrMod(>= 2.0.1), RobAStBase(>= 0.1.1)
Suggests: Biobase
Author: Matthias Kohl <Matthias.Kohl at stamats.de>
Modified: pkg/RobRex/DESCRIPTION
===================================================================
--- pkg/RobRex/DESCRIPTION 2008-11-05 18:28:07 UTC (rev 209)
+++ pkg/RobRex/DESCRIPTION 2008-11-12 21:33:45 UTC (rev 210)
@@ -2,14 +2,14 @@
Version: 0.6.0
Date: 2008-07-21
Title: Optimally robust influence curves for regression and scale
-Description: functions for the determination of optimally
- robust influence curves in case of linear regression
- with unknown scale and standard normal distributed errors
- where the regressor is random.
+Description: functions for the determination of optimally robust
+ influence curves in case of linear regression with unknown
+ scale and standard normal distributed errors where the
+ regressor is random.
Depends: R (>= 2.4.0), ROptRegTS(>= 0.5.0)
Author: Matthias Kohl <Matthias.Kohl at stamats.de>
Maintainer: Matthias Kohl <Matthias.Kohl at stamats.de>
LazyLoad: yes
-License: GPL version 2 or later
+License: LGPL-3
URL: http://www.stamats.de/RobASt.htm
-Packaged: Mon Aug 6 08:18:03 2007; kohl
+Packaged: Mon Aug 6 08:18:03 2007; kohl
More information about the Robast-commits
mailing list