[Robast-commits] r184 - branches/robast-0.6/pkg branches/robast-0.6/pkg/ROptEst.Rcheck branches/robast-0.7/pkg/ROptEst/man branches/robast-0.7/pkg/RobAStBase/R branches/robast-0.7/pkg/RobAStBase/man pkg/ROptEst/chm pkg/ROptEst/man pkg/RobAStBase/R pkg/RobAStBase/chm pkg/RobAStBase/man pkg/RobLox/chm
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Wed Oct 29 19:26:15 CET 2008
Author: ruckdeschel
Date: 2008-10-29 19:26:14 +0100 (Wed, 29 Oct 2008)
New Revision: 184
Added:
branches/robast-0.6/pkg/ROptEst.Rcheck/
branches/robast-0.6/pkg/ROptEst.Rcheck/ROptEst-manual.log
pkg/ROptEst/chm/getModifyIC.html
pkg/ROptEst/chm/roptest.html
Modified:
branches/robast-0.7/pkg/ROptEst/man/roptest.Rd
branches/robast-0.7/pkg/RobAStBase/R/AllPlot.R
branches/robast-0.7/pkg/RobAStBase/R/comparePlot.R
branches/robast-0.7/pkg/RobAStBase/R/infoPlot.R
branches/robast-0.7/pkg/RobAStBase/man/IC-class.Rd
branches/robast-0.7/pkg/RobAStBase/man/comparePlot.Rd
branches/robast-0.7/pkg/RobAStBase/man/infoPlot.Rd
pkg/ROptEst/chm/ROptEst.chm
pkg/ROptEst/man/roptest.Rd
pkg/RobAStBase/R/AllPlot.R
pkg/RobAStBase/R/comparePlot.R
pkg/RobAStBase/R/infoPlot.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/chm/comparePlot.html
pkg/RobAStBase/chm/infoPlot.html
pkg/RobAStBase/man/IC-class.Rd
pkg/RobAStBase/man/comparePlot.Rd
pkg/RobAStBase/man/infoPlot.Rd
pkg/RobLox/chm/RobLox.chm
Log:
some new plotting features (similar to plot() in distr ...)
Added: branches/robast-0.6/pkg/ROptEst.Rcheck/ROptEst-manual.log
===================================================================
--- branches/robast-0.6/pkg/ROptEst.Rcheck/ROptEst-manual.log (rev 0)
+++ branches/robast-0.6/pkg/ROptEst.Rcheck/ROptEst-manual.log 2008-10-29 18:26:14 UTC (rev 184)
@@ -0,0 +1,1142 @@
+This is pdfTeX, Version 3.141592-1.40.4 (MiKTeX 2.7) (preloaded format=pdflatex 2008.6.30) 27 AUG 2008 21:59
+entering extended mode
+**\nonstopmode{}\input{ROptEst-manual.tex}
+(ROptEst-manual.tex ("C:\Program Files\MiKTeX 2.7\tex\latex\base\article.cls"
+Document Class: article 2005/09/16 v1.4f Standard LaTeX document class
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\size10.clo"
+File: size10.clo 2005/09/16 v1.4f Standard LaTeX file (size option)
+)
+\c at part=\count79
+\c at section=\count80
+\c at subsection=\count81
+\c at subsubsection=\count82
+\c at paragraph=\count83
+\c at subparagraph=\count84
+\c at figure=\count85
+\c at table=\count86
+\abovecaptionskip=\skip41
+\belowcaptionskip=\skip42
+\bibindent=\dimen102
+)
+(C:\PROGRA~1\R\R-28~1.0DE\share\texmf\Rd.sty
+Package: Rd
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\ifthen.sty"
+Package: ifthen 2001/05/26 v1.1c Standard LaTeX ifthen package (DPC)
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\tools\longtable.sty"
+Package: longtable 2004/02/01 v4.11 Multi-page Table package (DPC)
+\LTleft=\skip43
+\LTright=\skip44
+\LTpre=\skip45
+\LTpost=\skip46
+\LTchunksize=\count87
+\LTcapwidth=\dimen103
+\LT at head=\box26
+\LT at firsthead=\box27
+\LT at foot=\box28
+\LT at lastfoot=\box29
+\LT at cols=\count88
+\LT at rows=\count89
+\c at LT@tables=\count90
+\c at LT@chunks=\count91
+\LT at p@ftn=\toks14
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\tools\bm.sty"
+Package: bm 2004/02/26 v1.1c Bold Symbol Support (DPC/FMi)
+\symboldoperators=\mathgroup4
+\symboldletters=\mathgroup5
+\symboldsymbols=\mathgroup6
+LaTeX Font Info: Redeclaring math alphabet \mathbf on input line 137.
+LaTeX Info: Redefining \bm on input line 203.
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\alltt.sty"
+Package: alltt 1997/06/16 v2.0g defines alltt environment
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\tools\verbatim.sty"
+Package: verbatim 2003/08/22 v1.5q LaTeX2e package for verbatim enhancements
+\every at verbatim=\toks15
+\verbatim at line=\toks16
+\verbatim at in@stream=\read1
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\ltxmisc\url.sty"
+\Urlmuskip=\muskip10
+Package: url 2006/04/12 ver 3.3 Verb mode for urls, etc.
+)
+(C:\PROGRA~1\R\R-28~1.0DE\share\texmf\upquote.sty
+Package: upquote 2003/08/11 v1.1 Covington's upright-quote modification to verb
+atim and verb
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\textcomp.sty"
+Package: textcomp 2005/09/27 v1.99g Standard LaTeX package
+Package textcomp Info: Sub-encoding information:
+(textcomp) 5 = only ISO-Adobe without \textcurrency
+(textcomp) 4 = 5 + \texteuro
+(textcomp) 3 = 4 + \textohm
+(textcomp) 2 = 3 + \textestimated + \textcurrency
+(textcomp) 1 = TS1 - \textcircled - \t
+(textcomp) 0 = TS1 (full)
+(textcomp) Font families with sub-encoding setting implement
+(textcomp) only a restricted character set as indicated.
+(textcomp) Family '?' is the default used for unknown fonts.
+(textcomp) See the documentation for details.
+Package textcomp Info: Setting ? sub-encoding to TS1/1 on input line 71.
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\ts1enc.def"
+File: ts1enc.def 2001/06/05 v3.0e (jk/car/fm) Standard LaTeX file
+)
+LaTeX Info: Redefining \oldstylenums on input line 266.
+Package textcomp Info: Setting cmr sub-encoding to TS1/0 on input line 281.
+Package textcomp Info: Setting cmss sub-encoding to TS1/0 on input line 282.
+Package textcomp Info: Setting cmtt sub-encoding to TS1/0 on input line 283.
+Package textcomp Info: Setting cmvtt sub-encoding to TS1/0 on input line 284.
+Package textcomp Info: Setting cmbr sub-encoding to TS1/0 on input line 285.
+Package textcomp Info: Setting cmtl sub-encoding to TS1/0 on input line 286.
+Package textcomp Info: Setting ccr sub-encoding to TS1/0 on input line 287.
+Package textcomp Info: Setting ptm sub-encoding to TS1/4 on input line 288.
+Package textcomp Info: Setting pcr sub-encoding to TS1/4 on input line 289.
+Package textcomp Info: Setting phv sub-encoding to TS1/4 on input line 290.
+Package textcomp Info: Setting ppl sub-encoding to TS1/3 on input line 291.
+Package textcomp Info: Setting pag sub-encoding to TS1/4 on input line 292.
+Package textcomp Info: Setting pbk sub-encoding to TS1/4 on input line 293.
+Package textcomp Info: Setting pnc sub-encoding to TS1/4 on input line 294.
+Package textcomp Info: Setting pzc sub-encoding to TS1/4 on input line 295.
+Package textcomp Info: Setting bch sub-encoding to TS1/4 on input line 296.
+Package textcomp Info: Setting put sub-encoding to TS1/5 on input line 297.
+Package textcomp Info: Setting uag sub-encoding to TS1/5 on input line 298.
+Package textcomp Info: Setting ugq sub-encoding to TS1/5 on input line 299.
+Package textcomp Info: Setting ul8 sub-encoding to TS1/4 on input line 300.
+Package textcomp Info: Setting ul9 sub-encoding to TS1/4 on input line 301.
+Package textcomp Info: Setting augie sub-encoding to TS1/5 on input line 302.
+Package textcomp Info: Setting dayrom sub-encoding to TS1/3 on input line 303.
+Package textcomp Info: Setting dayroms sub-encoding to TS1/3 on input line 304.
+
+Package textcomp Info: Setting pxr sub-encoding to TS1/0 on input line 305.
+Package textcomp Info: Setting pxss sub-encoding to TS1/0 on input line 306.
+Package textcomp Info: Setting pxtt sub-encoding to TS1/0 on input line 307.
+Package textcomp Info: Setting txr sub-encoding to TS1/0 on input line 308.
+Package textcomp Info: Setting txss sub-encoding to TS1/0 on input line 309.
+Package textcomp Info: Setting txtt sub-encoding to TS1/0 on input line 310.
+Package textcomp Info: Setting futs sub-encoding to TS1/4 on input line 311.
+Package textcomp Info: Setting futx sub-encoding to TS1/4 on input line 312.
+Package textcomp Info: Setting futj sub-encoding to TS1/4 on input line 313.
+Package textcomp Info: Setting hlh sub-encoding to TS1/3 on input line 314.
+Package textcomp Info: Setting hls sub-encoding to TS1/3 on input line 315.
+Package textcomp Info: Setting hlst sub-encoding to TS1/3 on input line 316.
+Package textcomp Info: Setting hlct sub-encoding to TS1/5 on input line 317.
+Package textcomp Info: Setting hlx sub-encoding to TS1/5 on input line 318.
+Package textcomp Info: Setting hlce sub-encoding to TS1/5 on input line 319.
+Package textcomp Info: Setting hlcn sub-encoding to TS1/5 on input line 320.
+Package textcomp Info: Setting hlcw sub-encoding to TS1/5 on input line 321.
+Package textcomp Info: Setting hlcf sub-encoding to TS1/5 on input line 322.
+Package textcomp Info: Setting pplx sub-encoding to TS1/3 on input line 323.
+Package textcomp Info: Setting pplj sub-encoding to TS1/3 on input line 324.
+Package textcomp Info: Setting ptmx sub-encoding to TS1/4 on input line 325.
+Package textcomp Info: Setting ptmj sub-encoding to TS1/4 on input line 326.
+))
+\ldescriptionwidth=\skip47
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\fontenc.sty"
+Package: fontenc 2005/09/27 v1.99g Standard LaTeX package
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\t1enc.def"
+File: t1enc.def 2005/09/27 v1.99g Standard LaTeX file
+LaTeX Font Info: Redeclaring font encoding T1 on input line 43.
+))
+("C:\Program Files\MiKTeX 2.7\tex\latex\ae\ae.sty"
+Package: ae 2001/02/12 1.3 Almost European Computer Modern
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\fontenc.sty"
+Package: fontenc 2005/09/27 v1.99g Standard LaTeX package
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\t1enc.def"
+File: t1enc.def 2005/09/27 v1.99g Standard LaTeX file
+LaTeX Font Info: Redeclaring font encoding T1 on input line 43.
+)
+LaTeX Font Info: Try loading font information for T1+aer on input line 100.
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\ae\t1aer.fd"
+File: t1aer.fd 1997/11/16 Font definitions for T1/aer.
+)))
+("C:\Program Files\MiKTeX 2.7\tex\latex\ae\t1aett.fd"
+File: t1aett.fd 1997/11/16 Font definitions for T1/aett.
+) NOT loading times
+NOT loading lmodern ("C:\Program Files\MiKTeX 2.7\tex\latex\graphics\color.sty"
+Package: color 2005/11/14 v1.0j Standard LaTeX Color (DPC)
+("C:\Program Files\MiKTeX 2.7\tex\latex\00miktex\color.cfg"
+File: color.cfg 2007/01/18 v1.5 color configuration of teTeX/TeXLive
+)
+Package color Info: Driver file: pdftex.def on input line 130.
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\graphics\pdftex.def"
+File: pdftex.def 2007/06/12 v0.04h Graphics/color for pdfTeX
+\Gread at gobject=\count92
+))
+("C:\Program Files\MiKTeX 2.7\tex\latex\hyperref\hyperref.sty"
+Package: hyperref 2008/01/09 v6.77i Hypertext links for LaTeX
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\graphics\keyval.sty"
+Package: keyval 1999/03/16 v1.13 key=value parser (DPC)
+\KV at toks@=\toks17
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\oberdiek\hycolor.sty"
+Package: hycolor 2007/04/11 v1.1 Code for color options of hyperref/bookmark (H
+O)
+)
+\@linkdim=\dimen104
+\Hy at linkcounter=\count93
+\Hy at pagecounter=\count94
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\hyperref\pd1enc.def"
+File: pd1enc.def 2008/01/09 v6.77i Hyperref: PDFDocEncoding definition (HO)
+)
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\etexcmds.sty"
+Package: etexcmds 2007/12/12 v1.2 Prefix for e-TeX command names (HO)
+
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\infwarerr.sty"
+Package: infwarerr 2007/09/09 v1.2 Providing info/warning/message (HO)
+)
+Package etexcmds Info: Could not find \expanded.
+(etexcmds) That can mean that you are not using pdfTeX 1.50 or
+(etexcmds) that some package has redefined \expanded.
+(etexcmds) In the latter case, load this package earlier.
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\00miktex\hyperref.cfg"
+File: hyperref.cfg 2002/06/06 v1.2 hyperref configuration of TeXLive
+)
+("C:\Program Files\MiKTeX 2.7\tex\latex\oberdiek\kvoptions.sty"
+Package: kvoptions 2007/10/18 v3.0 Keyval support for LaTeX options (HO)
+)
+Package hyperref Info: Hyper figures OFF on input line 2722.
+Package hyperref Info: Link nesting OFF on input line 2727.
+Package hyperref Info: Hyper index ON on input line 2730.
+Package hyperref Info: Plain pages OFF on input line 2737.
+Package hyperref Info: Backreferencing OFF on input line 2742.
+
+Implicit mode ON; LaTeX internals redefined
+Package hyperref Info: Bookmarks ON on input line 2898.
+LaTeX Info: Redefining \href on input line 3049.
+LaTeX Info: Redefining \url on input line 3068.
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\bitset.sty"
+Package: bitset 2007/09/28 v1.0 Data type bit set (HO)
+
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\intcalc.sty"
+Package: intcalc 2007/09/27 v1.1 Expandable integer calculations (HO)
+)
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\bigintcalc.sty"
+Package: bigintcalc 2007/11/11 v1.1 Expandable big integer calculations (HO)
+
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\pdftexcmds.sty"
+Package: pdftexcmds 2007/12/12 v0.3 LuaTeX support for pdfTeX utility functions
+ (HO)
+Package pdftexcmds Info: LuaTeX not detected on input line 139.
+)))
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\kvsetkeys.sty"
+Package: kvsetkeys 2007/09/29 v1.3 Key value parser with default handler suppor
+t (HO)
+)
+\Fld at menulength=\count95
+\Field at Width=\dimen105
+\Fld at charsize=\dimen106
+\Field at toks=\toks18
+Package hyperref Info: Hyper figures OFF on input line 3949.
+Package hyperref Info: Link nesting OFF on input line 3954.
+Package hyperref Info: Hyper index ON on input line 3957.
+Package hyperref Info: backreferencing OFF on input line 3964.
+Package hyperref Info: Link coloring OFF on input line 3969.
+
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\atbegshi.sty"
+Package: atbegshi 2007/09/09 v1.6 At begin shipout hook (HO)
+
+("C:\Program Files\MiKTeX 2.7\tex\generic\oberdiek\ifpdf.sty"
+Package: ifpdf 2007/12/12 v1.6 Provides the ifpdf switch (HO)
+Package ifpdf Info: pdfTeX in pdf mode detected.
+))
+\Hy at abspage=\count96
+\c at Item=\count97
+\c at Hfootnote=\count98
+)
+*hyperref using default driver hpdftex*
+("C:\Program Files\MiKTeX 2.7\tex\latex\hyperref\hpdftex.def"
+File: hpdftex.def 2008/01/09 v6.77i Hyperref driver for pdfTeX
+\Fld at listcount=\count99
+)
+
+Package hyperref Warning: Option `hyperindex' has already been used,
+(hyperref) setting the option has no effect on input line 315.
+
+Package hyperref Info: Option `colorlinks' set `true' on input line 315.
+
+Package hyperref Warning: Option `pagebackref' has already been used,
+(hyperref) setting the option has no effect on input line 315.
+
+Package hyperref Info: Option `linktocpage' set `true' on input line 315.
+Package hyperref Info: Option `plainpages' set `false' on input line 315.
+)
+No file ROptEst-manual.aux.
+LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for TS1/cmr/m/n on input line 3.
+LaTeX Font Info: Try loading font information for TS1+cmr on input line 3.
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\ts1cmr.fd"
+File: ts1cmr.fd 1999/05/25 v2.5h Standard LaTeX font definitions
+)
+LaTeX Font Info: ... okay on input line 3.
+LaTeX Font Info: Checking defaults for PD1/pdf/m/n on input line 3.
+LaTeX Font Info: ... okay on input line 3.
+
+("C:\Program Files\MiKTeX 2.7\tex\context\base\supp-pdf.tex"
+[Loading MPS to PDF converter (version 2006.09.02).]
+\scratchcounter=\count100
+\scratchdimen=\dimen107
+\scratchbox=\box30
+\nofMPsegments=\count101
+\nofMParguments=\count102
+\MPscratchCnt=\count103
+\MPscratchDim=\dimen108
+\MPnumerator=\count104
+\everyMPtoPDFconversion=\toks19
+)
+Package hyperref Info: Link coloring ON on input line 3.
+ ("C:\Program Files\MiKTeX 2.7\tex\latex\hyperref\nameref.sty"
+Package: nameref 2007/05/29 v2.31 Cross-referencing by name of section
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\oberdiek\refcount.sty"
+Package: refcount 2006/02/20 v3.0 Data extraction from references (HO)
+)
+\c at section@level=\count105
+)
+LaTeX Info: Redefining \ref on input line 3.
+LaTeX Info: Redefining \pageref on input line 3.
+\@outlinefile=\write3
+\AtBeginShipoutBox=\box31
+LaTeX Font Info: External font `cmex10' loaded for size
+(Font) <7> on input line 4.
+LaTeX Font Info: External font `cmex10' loaded for size
+(Font) <5> on input line 4.
+LaTeX Font Info: Try loading font information for TS1+aett on input line 30.
+
+
+(C:\PROGRA~1\R\R-28~1.0DE\share\texmf\ts1aett.fd
+File: ts1aett.fd
+)
+LaTeX Font Info: Try loading font information for TS1+cmtt on input line 30.
+
+
+("C:\Program Files\MiKTeX 2.7\tex\latex\base\ts1cmtt.fd"
+File: ts1cmtt.fd 1999/05/25 v2.5h Standard LaTeX font definitions
+)
+LaTeX Font Info: Font shape `TS1/aett/m/n' in size <10> not available
+(Font) Font shape `TS1/cmtt/m/n' tried instead on input line 30.
+
+Overfull \hbox (39.70007pt too wide) in paragraph at lines 33--33
+ [] \T1/aett/m/n/10 neighbor, biastype, cl
+ip, cent, stand, trafo)[]
+ []
+
+
+Overfull \hbox (60.66003pt too wide) in paragraph at lines 38--38
+ [] \T1/aett/m/n/10 biastype
+, clip, cent, stand, trafo)[]
+ []
+
+
+Overfull \hbox (65.90002pt too wide) in paragraph at lines 43--43
+ [] \T1/aett/m/n/1
+0 neighbor, biastype, trafo)[]
+ []
+
+
+Overfull \hbox (71.14001pt too wide) in paragraph at lines 48--48
+ [] \T1/aett/m/n/10 L2deri
+v, neighbor, biastype, trafo)[]
+ []
+
+
+Overfull \hbox (81.62pt too wide) in paragraph at lines 53--53
+ [] \T1/aett/m/n/10 L2de
+riv, neighbor, biastype, trafo)[]
+ []
+
+
+Overfull \hbox (86.85999pt too wide) in paragraph at lines 58--58
+ [] \T1/aett/m
+/n/10 neighbor, biastype, trafo)[]
+ []
+
+
+Overfull \hbox (65.90002pt too wide) in paragraph at lines 63--63
+ [] \T1/aett/m/n/10 biastype, Distr, DistrSymm, L2
+derivSymm, L2derivDistrSymm,[]
+ []
+
+[1
+
+{C:/Documents and Settings/All Users/Application Data/MiKTeX/2.7/pdftex/config/
+pdftex.map}]
+Overfull \hbox (65.90002pt too wide) in paragraph at lines 69--69
+ [] \T1/aett/m/n/10 neighbor,
+biastype, clip, cent, stand)[]
+ []
+
+
+Overfull \hbox (76.38pt too wide) in paragraph at lines 74--74
+ [] \T1/aett/m/n/10 neighbor
+, biastype, clip, cent, stand)[]
+ []
+
+
+Overfull \hbox (92.09998pt too wide) in paragraph at lines 80--80
+ [] \T1/aett/m/n/10 V.comp = matrix(TRUE, ncol = nro
+w(stand), nrow = nrow(stand)), w)[]
+ []
+
+
+Overfull \hbox (86.85999pt too wide) in paragraph at lines 85--85
+ [] \T1/aett/m/n/10 neighb
+or, biastype, clip, cent, stand)[]
+ []
+
+
+Overfull \hbox (92.09998pt too wide) in paragraph at lines 90--90
+ [] \T1/aett/m/n/10 biastype, Di
+str, clip, cent, stand, normtype)[]
+ []
+
+
+Overfull \hbox (102.57996pt too wide) in paragraph at lines 95--95
+ [] \T1/aett/m/n/10 neighbor,
+biastype, clip, cent, stand, trafo)[]
+ []
+
+
+Overfull \hbox (113.05994pt too wide) in paragraph at lines 100--100
+ [] \T1/aett/m/n/10 neighbor
+, biastype, clip, cent, stand, trafo)[]
+ []
+
+[2]
+Overfull \hbox (127.80415pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (104.54443pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (149.03398pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (187.29309pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (202.26332pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (168.54395pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (126.21434pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (147.33415pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (166.84412pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (124.51451pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (173.88379pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+[3]
+Overfull \hbox (136.81439pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (200.72327pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (181.39322pt too wide) in paragraph at lines 172--172
+[]
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 186--188
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[4]
+Overfull \hbox (2.1563pt too wide) in paragraph at lines 230--230
+[]\T1/aer/m/n/10 reads off the as. bias from
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 246--248
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[5]
+Overfull \hbox (21.62582pt too wide) in paragraph at lines 308--308
+[]\T1/aer/m/n/10 computes
+ []
+
+[6]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 317--319
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+
+Overfull \hbox (0.6396pt too wide) in paragraph at lines 320--322
+[]\T1/aer/m/n/10 Ruckdeschel, P. and Kohl, M. (2005) Com-pu-ta-tion of the Fi-n
+ite Sam-ple Risk of M-estimators
+ []
+
+[7]
+Overfull \hbox (71.85513pt too wide) in paragraph at lines 367--367
+[]
+ []
+
+
+Overfull \hbox (91.3651pt too wide) in paragraph at lines 367--367
+[]
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 376--378
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[8]
+Overfull \hbox (35.87558pt too wide) in paragraph at lines 430--430
+[]\T1/aer/m/n/10 computes
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 441--443
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[9] [10]
+Overfull \hbox (130.78477pt too wide) in paragraph at lines 504--504
+[]
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 519--521
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[11] [12]
+Overfull \hbox (83.38472pt too wide) in paragraph at lines 603--603
+[]
+ []
+
+
+Overfull \hbox (102.89468pt too wide) in paragraph at lines 603--603
+[]
+ []
+
+
+Overfull \hbox (60.56508pt too wide) in paragraph at lines 603--603
+[]
+ []
+
+
+Overfull \hbox (100.82428pt too wide) in paragraph at lines 603--603
+[]
+ []
+
+
+Overfull \hbox (115.79451pt too wide) in paragraph at lines 603--603
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 606--608
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 615--617
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[13]
+Overfull \hbox (86.85999pt too wide) in paragraph at lines 642--642
+ [] \T1/aett/m/n/10 risk, neighb
+or, biastype, cent, symm, trafo)[]
+ []
+
+
+Overfull \hbox (107.81995pt too wide) in paragraph at lines 647--647
+ [] \T1/aett/m/n/10 risk, ne
+ighbor, biastype, cent, symm, trafo)[]
+ []
+
+
+Overfull \hbox (86.85999pt too wide) in paragraph at lines 652--652
+ [] \T1/aett/m/n/10 neighbor, bias
+type, Distr, stand, cent, trafo)[]
+ []
+
+
+Overfull \hbox (128.7799pt too wide) in paragraph at lines 658--658
+ [] \T1/aett/m/n/10 risk
+, neighbor, biastype, cent, symm, trafo)[]
+ []
+
+
+Overfull \hbox (107.81995pt too wide) in paragraph at lines 663--663
+ [] \T1/aett/m/n/1
+0 risk, neighbor, cent, symm, trafo)[]
+ []
+
+[14]
+Overfull \hbox (139.15411pt too wide) in paragraph at lines 702--702
+[]
+ []
+
+
+Overfull \hbox (158.66408pt too wide) in paragraph at lines 702--702
+[]
+ []
+
+
+Overfull \hbox (115.8944pt too wide) in paragraph at lines 702--702
+[]
+ []
+
+
+Overfull \hbox (187.49344pt too wide) in paragraph at lines 702--702
+[]
+ []
+
+
+Overfull \hbox (154.48407pt too wide) in paragraph at lines 702--702
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 705--707
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 716--718
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[15] [16]
+Overfull \hbox (170.01355pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (196.1235pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (147.19391pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (204.10312pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (187.45311pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (202.42334pt too wide) in paragraph at lines 815--815
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 818--820
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 832--834
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[17]
+Overfull \hbox (149.73987pt too wide) in paragraph at lines 884--884
+ [] \T1/aett/m/n/10 L2derivDistrSymm, z.start, A.start, Finfo, traf
+o, maxiter, tol, warn, verbose = FALSE, ...)[]
+ []
+
+
+Overfull \hbox (254.53967pt too wide) in paragraph at lines 894--894
+ [] \T1/aett/m/n/10 L2derivDistrSymm, Finfo, trafo, onesetLM = FALS
+E, z.start, A.start, upper, maxiter, tol, warn, verbose = FALSE)[]
+ []
+
+[18]
+Overfull \hbox (254.53967pt too wide) in paragraph at lines 904--904
+ [] \T1/aett/m/n/10 L2derivDistrSymm, Finfo, trafo, onesetLM = FALS
+E, z.start, A.start, upper, maxiter, tol, warn, verbose = FALSE)[]
+ []
+
+
+Overfull \hbox (42.74551pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+[19]
+Overfull \hbox (62.25548pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (19.92587pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (58.6951pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (21.6257pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (76.51492pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (39.44553pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (68.64505pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (31.57565pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+
+Overfull \hbox (96.13463pt too wide) in paragraph at lines 989--989
+[]
+ []
+
+[20]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1005--1007
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[21]
+Overfull \hbox (83.38472pt too wide) in paragraph at lines 1090--1090
+[]
+ []
+
+
+Overfull \hbox (102.89468pt too wide) in paragraph at lines 1090--1090
+[]
+ []
+
+
+Overfull \hbox (60.56508pt too wide) in paragraph at lines 1090--1090
+[]
+ []
+
+
+Overfull \hbox (100.82428pt too wide) in paragraph at lines 1090--1090
+[]
+ []
+
+
+Overfull \hbox (115.79451pt too wide) in paragraph at lines 1090--1090
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 1093--1095
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1102--1104
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[22] [23]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1178--1180
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[24] [25] [26]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1307--1309
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+
+Overfull \hbox (39.92572pt too wide) in paragraph at lines 1364--1364
+[]\T1/aer/m/n/10 asymptotic
+ []
+
+[27]
+Overfull \hbox (34.30647pt too wide) in paragraph at lines 1364--1364
+[]
+ []
+
+
+Overfull \hbox (47.85645pt too wide) in paragraph at lines 1364--1364
+[]
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1380--1382
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+
+Overfull \hbox (0.6396pt too wide) in paragraph at lines 1383--1385
+[]\T1/aer/m/n/10 Ruckdeschel, P. and Kohl, M. (2005) Com-pu-ta-tion of the Fi-n
+ite Sam-ple Risk of M-estimators
+ []
+
+[28]
+Overfull \hbox (29.90625pt too wide) in paragraph at lines 1435--1435
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 1438--1440
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+[29]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1454--1456
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+
+Overfull \hbox (131.27475pt too wide) in paragraph at lines 1518--1518
+[]
+ []
+
+
+Overfull \hbox (150.78471pt too wide) in paragraph at lines 1518--1518
+[]
+ []
+
+
+Overfull \hbox (148.71431pt too wide) in paragraph at lines 1518--1518
+[]
+ []
+
+[30]
+Overfull \hbox (163.68454pt too wide) in paragraph at lines 1518--1518
+[]
+ []
+
+
+Overfull \hbox (183.27342pt too wide) in paragraph at lines 1518--1518
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 1521--1523
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1525--1527
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[31] [32]
+Overfull \hbox (83.38472pt too wide) in paragraph at lines 1628--1628
+[]
+ []
+
+
+Overfull \hbox (115.79451pt too wide) in paragraph at lines 1628--1628
+[]
+ []
+
+
+Overfull \hbox (102.89468pt too wide) in paragraph at lines 1628--1628
+[]
+ []
+
+
+Overfull \hbox (60.56508pt too wide) in paragraph at lines 1628--1628
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 1631--1633
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1642--1644
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[33]
+Overfull \hbox (23.9801pt too wide) in paragraph at lines 1665--1665
+ [] \T1/aett/m/n/10 tol = .Machine$double.e
+ps^0.4, warn = TRUE,[]
+ []
+
+
+Overfull \hbox (50.18005pt too wide) in paragraph at lines 1670--1670
+ [] \T1/aett/m/n/10 tol = .Machine$dou
+ble.eps^0.4, warn = TRUE)[]
+ []
+
+
+Overfull \hbox (55.42004pt too wide) in paragraph at lines 1674--1674
+ [] \T1/aett/m/n/10 maxiter = 50, tol
+= .Machine$double.eps^0.4,[]
+ []
+
+
+Overfull \hbox (29.2201pt too wide) in paragraph at lines 1675--1675
+ [] \T1/aett/m/n/10 warn = TRUE, Algo
+= "A", cont = "left",[]
+ []
+
+[34]
+Overfull \hbox (11.55566pt too wide) in paragraph at lines 1719--1719
+\T1/aer/m/n/10 ence curve for ro-bust mod-els with fixed neigh-bor-hoods and fi
+nite-sample under-/overshoot
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1732--1734
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[35]
+Overfull \hbox (81.62pt too wide) in paragraph at lines 1769--1769
+ [] \T1/aett/m/n/10 maxiter = 50, tol = .Machine$double.eps^0
+.4, warn = TRUE, noLow = FALSE)[]
+ []
+
+
+Overfull \hbox (39.70007pt too wide) in paragraph at lines 1773--1773
+ [] \T1/aett/m/n/10 tol = .Machine$double.eps^0.4, warn = TRUE, Alg
+o = "A", cont = "left")[]
+ []
+
+[36]
+Overfull \hbox (5.41653pt too wide) in paragraph at lines 1811--1811
+[]\T1/aer/m/n/10 asymptotic co-vari-ance of L2 dif-fer-en-tiable
+ []
+
+
+Overfull \hbox (3.31567pt too wide) in paragraph at lines 1811--1811
+[]\T1/aer/m/n/10 finite-sample under-/overshoot risk
+ []
+
+
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1824--1826
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+
+Overfull \hbox (3.02014pt too wide) in paragraph at lines 1850--1850
+ [] \T1/aett/m/n/10 L2Fam, neighbor, risk, loRad, upRad, z.start = NULL,
+ A.start = NULL,[]
+ []
+
+[37]
+Overfull \hbox (33.09619pt too wide) in paragraph at lines 1884--1884
+[]
+ []
+
+
+Overfull \hbox (58.06248pt too wide) in paragraph at lines 1887--1889
+[]\T1/aer/m/n/10 Matthias Kohl $\OMS/cmsy/m/n/10 h$\T1/aer/m/n/10 Matthias.Kohl
+ at stamats.de$\OMS/cmsy/m/n/10 i$\T1/aer/m/n/10 , Pe-ter Ruck-de-schel $\OMS/cmsy
+/m/n/10 h$\T1/aer/m/n/10 Peter.Ruckdeschel at itwm.fraunhofer.de$\OMS/cmsy/m/n/10
+i$
+ []
+
+[38]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 1900--1902
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[39] [40]
+Overfull \hbox (4.02951pt too wide) in paragraph at lines 2007--2009
+[]\T1/aer/m/n/10 Kohl, M. (2005) \T1/aer/m/it/10 Nu-mer-i-cal Con-tri-bu-tions
+to the Asymp-totic The-ory of Ro-bust-ness\T1/aer/m/n/10 . Bayreuth:
+ []
+
+[41] [42]
+Overfull \hbox (19.3521pt too wide) in paragraph at lines 2124--2124
+[]\T1/aett/m/n/10 signature(normtype = "SelfNorm")\T1/aer/m/n/10 : udates the n
+orm in the self-standardized
+ []
+
+[43] (ROptEst-manual.aux) )
+Here is how much of TeX's memory you used:
+ 5264 strings out of 95327
+ 78035 string characters out of 1183974
+ 130192 words of memory out of 1500000
+ 7978 multiletter control sequences out of 110000
+ 17143 words of font info for 40 fonts, out of 1200000 for 2000
+ 14 hyphenation exceptions out of 8191
+ 46i,6n,27p,264b,544s stack positions out of 5000i,500n,10000p,200000b,5000s
+pdfTeX warning (dest): name{Rfn.NormType.Rdash.class
+} has been referenced but does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.UncondNeighborhood.Rdash.class} has been refere
+nced but does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.roblox} has been referenced but does not exist,
+ replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.kStepEstimator} has been referenced but does no
+t exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.RiskType.Rdash.class} has been referenced but d
+oes not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.InfluenceCurve.Rdash.class} has been referenced
+ but does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.Neighborhood.Rdash.class} has been referenced b
+ut does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.L2ParamFamily.Rdash.class} has been referenced
+but does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.IC.Rdash.class} has been referenced but does no
+t exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.asUnOvShoot.Rdash.class} has been referenced bu
+t does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.asMSE.Rdash.class} has been referenced but does
+ not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.asGRisk.Rdash.class} has been referenced but do
+es not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.FixRobModel.Rdash.class} has been referenced bu
+t does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.TotalVarIC.Rdash.class} has been referenced but
+ does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.ContIC.Rdash.class} has been referenced but doe
+s not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.fiRisk.Rdash.class} has been referenced but doe
+s not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.InfRobModel.Rdash.class} has been referenced bu
+t does not exist, replaced by a fixed one
+
+pdfTeX warning (dest): name{Rfn.asRisk.Rdash.class} has been referenced but doe
+s not exist, replaced by a fixed one
+
+ <C:\Documents and Settings\All Users\Application Data\MiKTeX\2.7\Fonts\Pk\ljfo
+ur\jknappen\ec\dpi600\tctt1000.pk><C:/Program Files/MiKTeX 2.7/fonts/type1/blue
+sky/cm/cmbx10.pfb><C:/Program Files/MiKTeX 2.7/fonts/type1/bluesky/cm/cmmi10.pf
+b><C:/Program Files/MiKTeX 2.7/fonts/type1/bluesky/cm/cmr10.pfb><C:/Program Fil
+es/MiKTeX 2.7/fonts/type1/bluesky/cm/cmsy10.pfb><C:/Program Files/MiKTeX 2.7/fo
+nts/type1/bluesky/cm/cmti10.pfb><C:/Program Files/MiKTeX 2.7/fonts/type1/bluesk
+y/cm/cmtt10.pfb><C:/Program Files/MiKTeX 2.7/fonts/type1/bluesky/cm/cmtt9.pfb>
+Output written on ROptEst-manual.pdf (43 pages, 232772 bytes).
+PDF statistics:
+ 738 PDF objects out of 1000 (max. 8388607)
+ 423 named destinations out of 1000 (max. 131072)
+ 17 words of extra memory for PDF output out of 10000 (max. 10000000)
+
Modified: branches/robast-0.7/pkg/ROptEst/man/roptest.Rd
===================================================================
--- branches/robast-0.7/pkg/ROptEst/man/roptest.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/ROptEst/man/roptest.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -117,7 +117,7 @@
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)
+roptest(x, BinomFamily(size = 25), eps.lower = 0, eps.upper = 0.2, steps = 3)
#############################
@@ -132,7 +132,7 @@
MLEstimator(x, PoisFamily())
## compute optimally robust estimator (unknown contamination)
-roptest(x, PoisFamily(), eps.upper = 0.05, steps = 3)
+roptest(x, PoisFamily(), eps.upper = 0.1, steps = 3)
#############################
## 3. Normal (Gaussian) location and scale
@@ -150,7 +150,7 @@
## compute optimally robust estimator (unknown contamination)
## takes some time
-roptest(x, NormLocationScaleFamily(), eps.upper = 0.1, steps = 3)
+roptest(x, NormLocationScaleFamily(), eps.upper = 0.2, steps = 3)
}
\concept{k-step construction}
\concept{optimally robust estimation}
Modified: branches/robast-0.7/pkg/RobAStBase/R/AllPlot.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/AllPlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/R/AllPlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,10 +1,22 @@
-setMethod("plot", "IC",
- function(x,y=NULL,...){
+setMethod("plot", signature(x = "IC", y = "missing"),
+ function(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){
+ xc <- match.call(call = sys.call(sys.parent(1)))$x
dots <- match.call(call = sys.call(sys.parent(1)),
expand.dots = FALSE)$"..."
+ if(!is.logical(inner)){
+ if(!is.list(inner)||length(inner) != 4)
+ stop("Argument 'inner' must either be 'logical' or a 'list' vector of length 4")
+ innerD <- inner[1:3]
+ innerL <- inner[4]
+ }else{innerD <- innerL <- inner}
+
L2Fam <- eval(x at CallL2Fam)
e1 <- L2Fam at distribution
if(!is(e1, "UnivariateDistribution")) stop("not yet implemented")
@@ -32,23 +44,97 @@
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
dims <- nrow(L2Fam at param@trafo)
IC1 <- as(diag(dims) %*% x at Curve, "EuclRandVariable")
+ mainL <- FALSE
+ subL <- FALSE
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c("%C", "%D", "%A"),
+ c(as.character(class(x)[1]),
+ as.character(date()),
+ as.character(deparse(xc))))
+
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- gettextf("Plot for IC %%A") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+
+ if(is.logical(innerL)){
+ innerT <- paste(gettextf("Component "), 1:dims,
+ gettextf(" of (partial) IC\nfor"),
+ name(L2Fam)[1],
+ gettextf("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3), collapse = ", "),")")
+ if(!is.null(L2Fam at param@nuisance))
+ innerT <- paste(innerT,
+ gettextf("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
+ ")",
+ sep="" )
+ if(!is.null(L2Fam at param@fixed))
+ innerT <- paste(innerT,
+ gettextf("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3), collapse = ", "),
+ ")",
+ sep="" )
+ }else{
+ innerT <- rep(sapply(inner, .mpresubs), length.out=dims)
+ }
+
+
w0 <- options("warn")
options(warn = -1)
opar <- par()
+ if (!withSweave)
+ devNew()
nrows <- trunc(sqrt(dims))
ncols <- ceiling(dims/nrows)
- par(mfrow = c(nrows, ncols))
+
+ if(mfColRow)
+ parArgs <- list(mfrow = c(nrows, ncols))
- if(is.null(dots[["cex.main"]])) dots["cex.main"] <- 0.8
+ omar <- par("mar")
+ parArgs <- c(parArgs,list(mar = c(bmar,omar[2],tmar,omar[4])))
+ do.call(par,args=parArgs)
+
+ dotsT <- dots
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+
+
for(i in 1:dims){
do.call(plot, args=c(list(x.vec, sapply(x.vec, IC1 at Map[[i]]),
type = plty, lty = lty,
@@ -59,21 +145,21 @@
do.call(lines,args=c(list(x.vec1, sapply(x.vec1, IC1 at Map[[i]]),
lty = "dotted"), dots))
}
- if(is.null(L2Fam at param@nuisance))
- do.call(title,args=c(list(paste("Component", i,
- "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")")), dots))
- else
- do.call(title,args=c(list(paste("Component", i,
- "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
- ")")), dots))
+ do.call(title,args=c(list(main = innerT[i]), dotsT, line = lineT,
+ cex.main = cex.inner, col.main = col.inner))
}
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
par(opar)
options(w0)
invisible()
Modified: branches/robast-0.7/pkg/RobAStBase/R/comparePlot.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/comparePlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/R/comparePlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,5 +1,10 @@
setMethod("comparePlot", signature("IC","IC"),
- function(obj1,obj2, obj3 = NULL, obj4 = NULL, ...){
+ function(obj1,obj2, obj3 = NULL, obj4 = NULL,
+ ..., 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){
xc1 <- as.character(deparse(match.call(call = sys.call(sys.parent(1)))$obj1))
xc2 <- as.character(deparse(match.call(call = sys.call(sys.parent(1)))$obj2))
@@ -14,15 +19,18 @@
ncomp <- 2+ !is.null(obj3) + !is.null(obj4)
+ if(is.null(dots[["col"]])) dots$"col" <- 1:ncomp
+ if(is.null(dots[["lwd"]])) dots$"lwd" <- 1
+
+ col <- dots[["col"]]
+ lwd <- dots[["lwd"]]
+
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
- if(is.null(dots[["col"]])) dots$"col" <- 1:ncomp
- if(is.null(dots[["cex.main"]])) dots$"cex.main" <- 0.8
- if(is.null(dots[["lwd"]])) dots$"lwd" <- 1
+
+ dotsP <- dotsL <- dotsT <- dots
L2Fam <- eval(obj1 at CallL2Fam)
@@ -58,6 +66,7 @@
IC1 <- as(diag(dims) %*% obj1 at Curve, "EuclRandVariable")
IC2 <- as(diag(dims) %*% obj2 at Curve, "EuclRandVariable")
+
obj <- obj3
if(is(obj, "IC"))
{
@@ -74,6 +83,83 @@
IC4 <- as(diag(dims) %*% obj4 at Curve, "EuclRandVariable")
}
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c(paste("%C",1:ncomp,sep=""),
+ "%D",
+ paste("%A",1:ncomp,sep="")),
+ c(as.character(class(obj1)[1]),
+ as.character(class(obj2)[1]),
+ if(is.null(obj3))NULL else as.character(class(obj3)[1]),
+ if(is.null(obj4))NULL else as.character(class(obj4)[1]),
+ as.character(date()),
+ xc))
+
+ mainL <- FALSE
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- paste(gettextf("Plot for ICs"),
+ paste("%A", 1:ncomp, sep="", collapse=", "),
+ sep=" ") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ subL <- FALSE
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+ innerParam <- paste(gettext("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3),
+ collapse = ", "),
+ ")", sep = "")
+ if(!is.null(L2Fam at param@nuisance))
+ innerParam <- paste(innerParam,
+ gettext("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3),
+ collapse = ", "),
+ ")", sep ="")
+ if(!is.null(L2Fam at param@fixed))
+ innerParam <- paste(innerParam,
+ gettext("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3),
+ collapse = ", "),
+ ")", sep ="")
+
+ if(!is.logical(inner)){
+ if(!is.character(inner))
+ stop("Argument 'inner' must either be 'logical' or a character vector")
+ innerT <- rep(inner,length.out=dims)
+ innerL <- TRUE
+ }else{if(any(is.na(inner))||any(!inner)) {
+ innerT <- ""; innerL <- FALSE
+ }else{innerL <- TRUE
+ innerT <- paste(paste(gettext("Component "), 1:dims,
+ gettext(" of (partial) IC\nfor "),
+ name(L2Fam)[1], sep =""), innerParam)
+ }
+ }
+
+
w0 <- options("warn")
options(warn = -1)
opar <- par()
@@ -84,6 +170,11 @@
if(is(e1, "DiscreteDistribution"))
x.vec1 <- seq(from = min(x.vec), to = max(x.vec), length = 1000)
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+
for(i in 1:dims){
matp <- cbind(sapply(x.vec, IC1 at Map[[i]]),sapply(x.vec, IC2 at Map[[i]]))
if(is(obj3, "IC"))
@@ -93,31 +184,39 @@
do.call(matplot, args=c(list( x= x.vec, y=matp,
type = plty, lty = lty,
- xlab = "x", ylab = "(partial) IC"), dots))
+ xlab = "x", ylab = "(partial) IC"), dotsP))
+
if(is(e1, "DiscreteDistribution")){
matp1 <- cbind(sapply(x.vec1, IC1 at Map[[i]]),sapply(x.vec1, IC2 at Map[[i]]))
if(is(obj3, "IC"))
matp1 <- cbind(matp1,sapply(x.vec1, IC3 at Map[[i]]))
if(is(obj4, "IC"))
matp1 <- cbind(matp1,sapply(x.vec1, IC4 at Map[[i]]))
- do.call(matlines, c(list(x.vec1, matp1, lty = "dotted"),dots))
+ do.call(matlines, c(list(x.vec1, matp1, lty = "dotted"),dotsL))
}
- if(is.null(L2Fam at param@nuisance))
- do.call(title, c(list(paste("Component", i, "of (partial) ICs\nfor", name(L2Fam)[1],
- "\nwith main parameter (", paste(round(L2Fam at param@main, 3), collapse = ", "), ")")),
- dots))
- else
- do.call(title, c(list(paste("Component", i, "of (partial) ICs\nfor", name(L2Fam)[1],
- "\nwith main parameter (", paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (", paste(round(L2Fam at param@nuisance, 3), collapse = ", "), ")")),
- dots))
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[i]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
}
legend("bottomright",
legend = xc, col = eval(dots[["col"]]),
cex=0.75, lwd=eval(dots[["lwd"]])*1.5)
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
+
par(opar)
options(w0)
invisible()
Modified: branches/robast-0.7/pkg/RobAStBase/R/infoPlot.R
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/R/infoPlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/R/infoPlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,26 +1,30 @@
setMethod("infoPlot", "IC",
- function(object, ...){
+ function(object, ..., withSweave = getdistrOption("withSweave"),
+ colI = grey(0.5), lwdI = 0.7*par("lwd"),
+ 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){
+
+ objectc <- match.call(call = sys.call(sys.parent(1)))$object
dots <- match.call(call = sys.call(sys.parent(1)),
expand.dots = FALSE)$"..."
+
L2Fam <- eval(object at CallL2Fam)
-
+
+ if(!hasArg(col)) col <- par("col") else col <- dots$col
+ if(!hasArg(lwd)) lwd <- par("lwd") else lwd <- dots$lwd
+
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
- if(!is.null(dots[["xlim"]])) dots["xlim"] <- NULL
- if(!is.null(dots[["ylim"]])) dots["ylim"] <- NULL
- if(is.null(dots[["colA"]])) dots$"colA" <- grey(0.5)
- if(is.null(dots[["colB"]])) dots$"colB" <- par("col")
- if(is.null(dots[["lwdA"]])) dots$"lwdA" <- par("lwd")
- if(is.null(dots[["lwdB"]])) dots$"lwdB" <- 2
- if(is.null(dots[["cex.main"]])) dots$"cex.main" <- 0.8
+ dotsP <- dotsL <- dotsT <- dots
-
+ print(dots)
+ print(lwd)
e1 <- L2Fam at distribution
if(!is(e1, "UnivariateDistribution") | is(e1, "CondDistribution"))
stop("not yet implemented")
@@ -42,7 +46,81 @@
trafo <- L2Fam at param@trafo
dims <- nrow(trafo)
+
+
+ mainL <- FALSE
+ subL <- FALSE
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c("%C", "%D", "%A"),
+ c(as.character(class(object)[1]),
+ as.character(date()),
+ as.character(deparse(objectc))))
+
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- gettextf("Plot for IC %%A") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+ innerParam <- paste(gettext("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3),
+ collapse = ", "),
+ ")", sep = "")
+ if(!is.null(L2Fam at param@nuisance))
+ innerParam <- paste(innerParam,
+ gettext("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3),
+ collapse = ", "),
+ ")", sep ="")
+ if(!is.null(L2Fam at param@fixed))
+ innerParam <- paste(innerParam,
+ gettext("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3),
+ collapse = ", "),
+ ")", sep ="")
+
+ if(!is.logical(inner)){
+ if(!is.character(inner))
+ stop("Argument 'inner' must either be 'logical' or a character vector")
+ innerT <- rep(inner,length.out=1+dims)
+ innerL <- TRUE
+ }else{if(any(is.na(inner))||any(!inner)) {
+ innerT <- rep("",1+dims); innerL <- FALSE
+ }else{innerL <- TRUE
+ innerT <- paste(c( paste(gettext("Absolute information of (partial) IC for"),
+ name(L2Fam)[1], sep =""),
+ paste(gettext("Relative information of \ncomponent "),
+ 1:dims,
+ gettext("of (partial) IC\nfor "),
+ name(L2Fam)[1], sep ="")), innerParam)
+ }
+ }
+
QFc <- diag(dims)
if(is(object,"ContIC") & dims>1 )
{if (is(normtype(object),"QFNorm")) QFc <- QuadForm(normtype(object))
@@ -65,80 +143,92 @@
absInfo <- t(IC1) %*% QF %*% IC1
absInfo <- sapply(x.vec, absInfo at Map[[1]])
- dots["col"] <- dots[["colA"]]
- dots["lwd"] <- dots[["lwdA"]]
+
+ w0 <- options("warn")
+ options(warn = -1)
+ opar <- par()
+ if (!withSweave)
+ devNew()
+
+ omar <- par("mar")
+ parArgs <- list(mar = c(bmar,omar[2],tmar,omar[4]))
+ do.call(par,args=parArgs)
+
+
+
+ dotsP["col"] <- NULL
+ dotsP["lwd"] <- NULL
+ if(!hasArg(ylim)) dots["ylim"] <- c(0, 2*max(absInfo, na.rm = TRUE))
+
do.call(plot, args=c(list(x.vec, absInfoClass, type = plty,
- lty = "dashed",
- ylim = c(0, 2*max(absInfo, na.rm = TRUE)), xlab = "x",
- ylab = "absolute information"), dots))
- dots["col"] <- dots[["colB"]]
- dots["lwd"] <- dots[["lwdB"]]
+ lty = "dashed", col = colI, lwd = lwdI,
+ xlab = "x",
+ ylab = "absolute information"), dotsP))
do.call(lines, args=c(list(x.vec, absInfo, type = plty, lty = lty),
- dots))
+ dotsL))
legend(max(x.vec), 0, xjust = 1, yjust = 0,
- legend = c("class. opt. IC"), lty = "dashed",
- col = c(dots[["colA"]]), cex=0.75)
+ legend = c("class. opt. IC", objectc),
+ lty = c(lty,"dashed"), col = c(colI, col),
+ lwd=c(lwdI, lwd), cex = 0.75)
- if(is.null(L2Fam at param@nuisance))
- do.call(title, args=c(list(paste("Absolute information of (partial) IC for",
- name(L2Fam)[1], "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "), ")")),
- dots))
- else
- do.call(title, args=c(list(paste("Absolute information of (partial) IC for",
- name(L2Fam)[1], "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
- ")")), dots))
-
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[1]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
+
if(dims > 1){
+ dotsP["ylim"] <- NULL
+ dotsL["ylim"] <- NULL
+ dotsT["ylim"] <- NULL
nrows <- trunc(sqrt(dims))
ncols <- ceiling(dims/nrows)
- w0 <- options("warn")
- options(warn = -1)
- opar <- par()
- devNew()
- par(mfrow = c(nrows, ncols))
+ if (!withSweave)
+ devNew()
+ if(mfColRow)
+ parArgs <- c(parArgs,list(mfrow = c(nrows, ncols)))
+
+ do.call(par,args=parArgs)
+
IC1.i.5 <- QF.5%*%IC1
classIC.i.5 <- QFc.5%*%classIC
for(i in 1:dims){
y.vec <- sapply(x.vec, IC1.i.5 at Map[[i]])^2/absInfo
- dots["col"] <- dots[["colB"]]
- dots["lwd"] <- dots[["lwdA"]]
do.call(plot, args=c(list(x.vec, y.vec, type = plty,
lty = lty, xlab = "x",
ylab = "relative information",
- ylim = c(0, 1.1)), dots))
+ ylim = c(0, 1.1),
+ col = colI, lwd = lwdI), dotsP))
yc.vec <- sapply(x.vec, classIC.i.5 at Map[[i]])^2/absInfoClass
- dots["col"] <- dots[["colA"]]
- dots["lwd"] <- dots[["lwdB"]]
do.call(lines, args=c(list(x.vec, yc.vec, type = plty,
- lty = "dashed"),dots))
- legend(max(x.vec), 1.1, xjust = 1, cex = 0.6,
- legend = c("class. opt. IC"), lty = "dashed",
- col = c(dots[["colA"]]))
- if(is.null(L2Fam at param@nuisance))
- do.call(title, args=c(list(paste("Relative information of\ncomponent",
- i, "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3),
- collapse = ", "), ")")),
- dots))
- else
- do.call(title, args=c(list(paste("Relative information of\ncomponent",
- i, "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3),
- collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3),
- collapse = ", "), ")")),
- dots))
+ lty = "dashed"), dotsL))
+ legend(max(x.vec), 1.1, xjust = 1,
+ legend = c("class. opt. IC", objectc), lty = c(lty,"dashed"),
+ col = c(colI, col), lwd=c(lwdI, lwd),
+ cex = 0.6)
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[1+i]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
}
}
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
+
par(opar)
options(w0)
}
})
+
\ No newline at end of file
Modified: branches/robast-0.7/pkg/RobAStBase/man/IC-class.Rd
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/man/IC-class.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/man/IC-class.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -12,7 +12,6 @@
\alias{evalIC,IC,numeric-method}
\alias{evalIC,IC,matrix-method}
\alias{infoPlot,IC-method}
-\alias{plot,IC,ANY-method}
\alias{show,IC-method}
\title{Influence curve}
Modified: branches/robast-0.7/pkg/RobAStBase/man/comparePlot.Rd
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/man/comparePlot.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/man/comparePlot.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -10,15 +10,61 @@
}
\usage{
comparePlot(obj1, obj2, ... )
-\S4method{comparePlot}{IC,IC}(obj1, obj2, obj3 = NULL, obj4 = NULL, ... )
+\S4method{comparePlot}{IC,IC}(obj1, obj2, obj3 = NULL, obj4 = NULL,
+ ..., 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{obj1}{ object of class \code{"InfluenceCurve"} }
\item{obj2}{ object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
\item{obj3}{ optional: object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
\item{obj4}{ optional: object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
+ \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 arguments to be passed to \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{"\%C1"},\code{"\%C2"},[\code{"\%C3"},][\code{"\%C4"}]}{class of argument
+ \code{obj<i>}, i=1,..4}
+\item{\code{"\%A1"},\code{"\%A2"},[\code{"\%A3"},][\code{"\%A4"}]}{deparsed argument
+ \code{obj<i>}, i=1,..4}
+\item{\code{"\%D"}}{time/date-string when the plot was generated}
+}
+}
+
%\value{}
\references{
Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}.
Modified: branches/robast-0.7/pkg/RobAStBase/man/infoPlot.Rd
===================================================================
--- branches/robast-0.7/pkg/RobAStBase/man/infoPlot.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ branches/robast-0.7/pkg/RobAStBase/man/infoPlot.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,15 +1,41 @@
\name{infoPlot}
\alias{infoPlot}
+\alias{infoPlot-methods}
+\alias{infoPlot,IC-method}
\title{Plot absolute and relative information}
\description{
Plot absolute and relative information of influence curves.
}
\usage{
-infoPlot(object, ...)
+infoPlot(object, ...)
+\S4method{infoPlot}{IC}(object, ..., withSweave = getdistrOption("withSweave"),
+ colI = grey(0.5), lwdI = 0.7*par("lwd"),
+ 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{object}{ object of class \code{"InfluenceCurve"} }
+ \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{colI}{color of the classically optimal IC}
+ \item{lwdI}{linewidth of the classically optimal IC}
+ \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{
@@ -17,8 +43,29 @@
of an IC. The relative information is defined as the
absolute information of one component with respect to the
absolute information of the whole IC; confer Section 8.1
- of Kohl (2005).
+ of Kohl (2005). \cr
+
+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}
}
+}
%\value{}
\references{
Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}.
Modified: pkg/ROptEst/chm/ROptEst.chm
===================================================================
(Binary files differ)
Added: pkg/ROptEst/chm/getModifyIC.html
===================================================================
--- pkg/ROptEst/chm/getModifyIC.html (rev 0)
+++ pkg/ROptEst/chm/getModifyIC.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -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: pkg/ROptEst/chm/roptest.html
===================================================================
--- pkg/ROptEst/chm/roptest.html (rev 0)
+++ pkg/ROptEst/chm/roptest.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -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.2, 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.1, 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.2, 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>
Modified: pkg/ROptEst/man/roptest.Rd
===================================================================
--- pkg/ROptEst/man/roptest.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/ROptEst/man/roptest.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -117,7 +117,7 @@
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)
+roptest(x, BinomFamily(size = 25), eps.lower = 0, eps.upper = 0.2, steps = 3)
#############################
@@ -132,7 +132,7 @@
MLEstimator(x, PoisFamily())
## compute optimally robust estimator (unknown contamination)
-roptest(x, PoisFamily(), eps.upper = 0.05, steps = 3)
+roptest(x, PoisFamily(), eps.upper = 0.1, steps = 3)
#############################
## 3. Normal (Gaussian) location and scale
@@ -150,7 +150,7 @@
## compute optimally robust estimator (unknown contamination)
## takes some time
-roptest(x, NormLocationScaleFamily(), eps.upper = 0.1, steps = 3)
+roptest(x, NormLocationScaleFamily(), eps.upper = 0.2, steps = 3)
}
\concept{k-step construction}
\concept{optimally robust estimation}
Modified: pkg/RobAStBase/R/AllPlot.R
===================================================================
--- pkg/RobAStBase/R/AllPlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/R/AllPlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,10 +1,22 @@
-setMethod("plot", "IC",
- function(x,y=NULL,...){
+setMethod("plot", signature(x = "IC", y = "missing"),
+ function(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){
+ xc <- match.call(call = sys.call(sys.parent(1)))$x
dots <- match.call(call = sys.call(sys.parent(1)),
expand.dots = FALSE)$"..."
+ if(!is.logical(inner)){
+ if(!is.list(inner)||length(inner) != 4)
+ stop("Argument 'inner' must either be 'logical' or a 'list' vector of length 4")
+ innerD <- inner[1:3]
+ innerL <- inner[4]
+ }else{innerD <- innerL <- inner}
+
L2Fam <- eval(x at CallL2Fam)
e1 <- L2Fam at distribution
if(!is(e1, "UnivariateDistribution")) stop("not yet implemented")
@@ -32,23 +44,97 @@
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
dims <- nrow(L2Fam at param@trafo)
IC1 <- as(diag(dims) %*% x at Curve, "EuclRandVariable")
+ mainL <- FALSE
+ subL <- FALSE
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c("%C", "%D", "%A"),
+ c(as.character(class(x)[1]),
+ as.character(date()),
+ as.character(deparse(xc))))
+
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- gettextf("Plot for IC %%A") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+
+ if(is.logical(innerL)){
+ innerT <- paste(gettextf("Component "), 1:dims,
+ gettextf(" of (partial) IC\nfor"),
+ name(L2Fam)[1],
+ gettextf("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3), collapse = ", "),")")
+ if(!is.null(L2Fam at param@nuisance))
+ innerT <- paste(innerT,
+ gettextf("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
+ ")",
+ sep="" )
+ if(!is.null(L2Fam at param@fixed))
+ innerT <- paste(innerT,
+ gettextf("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3), collapse = ", "),
+ ")",
+ sep="" )
+ }else{
+ innerT <- rep(sapply(inner, .mpresubs), length.out=dims)
+ }
+
+
w0 <- options("warn")
options(warn = -1)
opar <- par()
+ if (!withSweave)
+ devNew()
nrows <- trunc(sqrt(dims))
ncols <- ceiling(dims/nrows)
- par(mfrow = c(nrows, ncols))
+
+ if(mfColRow)
+ parArgs <- list(mfrow = c(nrows, ncols))
- if(is.null(dots[["cex.main"]])) dots["cex.main"] <- 0.8
+ omar <- par("mar")
+ parArgs <- c(parArgs,list(mar = c(bmar,omar[2],tmar,omar[4])))
+ do.call(par,args=parArgs)
+
+ dotsT <- dots
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+
+
for(i in 1:dims){
do.call(plot, args=c(list(x.vec, sapply(x.vec, IC1 at Map[[i]]),
type = plty, lty = lty,
@@ -59,21 +145,21 @@
do.call(lines,args=c(list(x.vec1, sapply(x.vec1, IC1 at Map[[i]]),
lty = "dotted"), dots))
}
- if(is.null(L2Fam at param@nuisance))
- do.call(title,args=c(list(paste("Component", i,
- "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")")), dots))
- else
- do.call(title,args=c(list(paste("Component", i,
- "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
- ")")), dots))
+ do.call(title,args=c(list(main = innerT[i]), dotsT, line = lineT,
+ cex.main = cex.inner, col.main = col.inner))
}
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
par(opar)
options(w0)
invisible()
Modified: pkg/RobAStBase/R/comparePlot.R
===================================================================
--- pkg/RobAStBase/R/comparePlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/R/comparePlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,5 +1,10 @@
setMethod("comparePlot", signature("IC","IC"),
- function(obj1,obj2, obj3 = NULL, obj4 = NULL, ...){
+ function(obj1,obj2, obj3 = NULL, obj4 = NULL,
+ ..., 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){
xc1 <- as.character(deparse(match.call(call = sys.call(sys.parent(1)))$obj1))
xc2 <- as.character(deparse(match.call(call = sys.call(sys.parent(1)))$obj2))
@@ -14,15 +19,18 @@
ncomp <- 2+ !is.null(obj3) + !is.null(obj4)
+ if(is.null(dots[["col"]])) dots$"col" <- 1:ncomp
+ if(is.null(dots[["lwd"]])) dots$"lwd" <- 1
+
+ col <- dots[["col"]]
+ lwd <- dots[["lwd"]]
+
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
- if(is.null(dots[["col"]])) dots$"col" <- 1:ncomp
- if(is.null(dots[["cex.main"]])) dots$"cex.main" <- 0.8
- if(is.null(dots[["lwd"]])) dots$"lwd" <- 1
+
+ dotsP <- dotsL <- dotsT <- dots
L2Fam <- eval(obj1 at CallL2Fam)
@@ -58,6 +66,7 @@
IC1 <- as(diag(dims) %*% obj1 at Curve, "EuclRandVariable")
IC2 <- as(diag(dims) %*% obj2 at Curve, "EuclRandVariable")
+
obj <- obj3
if(is(obj, "IC"))
{
@@ -74,6 +83,83 @@
IC4 <- as(diag(dims) %*% obj4 at Curve, "EuclRandVariable")
}
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c(paste("%C",1:ncomp,sep=""),
+ "%D",
+ paste("%A",1:ncomp,sep="")),
+ c(as.character(class(obj1)[1]),
+ as.character(class(obj2)[1]),
+ if(is.null(obj3))NULL else as.character(class(obj3)[1]),
+ if(is.null(obj4))NULL else as.character(class(obj4)[1]),
+ as.character(date()),
+ xc))
+
+ mainL <- FALSE
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- paste(gettextf("Plot for ICs"),
+ paste("%A", 1:ncomp, sep="", collapse=", "),
+ sep=" ") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ subL <- FALSE
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+ innerParam <- paste(gettext("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3),
+ collapse = ", "),
+ ")", sep = "")
+ if(!is.null(L2Fam at param@nuisance))
+ innerParam <- paste(innerParam,
+ gettext("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3),
+ collapse = ", "),
+ ")", sep ="")
+ if(!is.null(L2Fam at param@fixed))
+ innerParam <- paste(innerParam,
+ gettext("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3),
+ collapse = ", "),
+ ")", sep ="")
+
+ if(!is.logical(inner)){
+ if(!is.character(inner))
+ stop("Argument 'inner' must either be 'logical' or a character vector")
+ innerT <- rep(inner,length.out=dims)
+ innerL <- TRUE
+ }else{if(any(is.na(inner))||any(!inner)) {
+ innerT <- ""; innerL <- FALSE
+ }else{innerL <- TRUE
+ innerT <- paste(paste(gettext("Component "), 1:dims,
+ gettext(" of (partial) IC\nfor "),
+ name(L2Fam)[1], sep =""), innerParam)
+ }
+ }
+
+
w0 <- options("warn")
options(warn = -1)
opar <- par()
@@ -84,6 +170,11 @@
if(is(e1, "DiscreteDistribution"))
x.vec1 <- seq(from = min(x.vec), to = max(x.vec), length = 1000)
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+
for(i in 1:dims){
matp <- cbind(sapply(x.vec, IC1 at Map[[i]]),sapply(x.vec, IC2 at Map[[i]]))
if(is(obj3, "IC"))
@@ -93,31 +184,39 @@
do.call(matplot, args=c(list( x= x.vec, y=matp,
type = plty, lty = lty,
- xlab = "x", ylab = "(partial) IC"), dots))
+ xlab = "x", ylab = "(partial) IC"), dotsP))
+
if(is(e1, "DiscreteDistribution")){
matp1 <- cbind(sapply(x.vec1, IC1 at Map[[i]]),sapply(x.vec1, IC2 at Map[[i]]))
if(is(obj3, "IC"))
matp1 <- cbind(matp1,sapply(x.vec1, IC3 at Map[[i]]))
if(is(obj4, "IC"))
matp1 <- cbind(matp1,sapply(x.vec1, IC4 at Map[[i]]))
- do.call(matlines, c(list(x.vec1, matp1, lty = "dotted"),dots))
+ do.call(matlines, c(list(x.vec1, matp1, lty = "dotted"),dotsL))
}
- if(is.null(L2Fam at param@nuisance))
- do.call(title, c(list(paste("Component", i, "of (partial) ICs\nfor", name(L2Fam)[1],
- "\nwith main parameter (", paste(round(L2Fam at param@main, 3), collapse = ", "), ")")),
- dots))
- else
- do.call(title, c(list(paste("Component", i, "of (partial) ICs\nfor", name(L2Fam)[1],
- "\nwith main parameter (", paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (", paste(round(L2Fam at param@nuisance, 3), collapse = ", "), ")")),
- dots))
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[i]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
}
legend("bottomright",
legend = xc, col = eval(dots[["col"]]),
cex=0.75, lwd=eval(dots[["lwd"]])*1.5)
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
+
par(opar)
options(w0)
invisible()
Modified: pkg/RobAStBase/R/infoPlot.R
===================================================================
--- pkg/RobAStBase/R/infoPlot.R 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/R/infoPlot.R 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,26 +1,30 @@
setMethod("infoPlot", "IC",
- function(object, ...){
+ function(object, ..., withSweave = getdistrOption("withSweave"),
+ colI = grey(0.5), lwdI = 0.7*par("lwd"),
+ 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){
+
+ objectc <- match.call(call = sys.call(sys.parent(1)))$object
dots <- match.call(call = sys.call(sys.parent(1)),
expand.dots = FALSE)$"..."
+
L2Fam <- eval(object at CallL2Fam)
-
+
+ if(!hasArg(col)) col <- par("col") else col <- dots$col
+ if(!hasArg(lwd)) lwd <- par("lwd") else lwd <- dots$lwd
+
if(!is.null(dots[["lty"]])) dots["lty"] <- NULL
if(!is.null(dots[["type"]])) dots["type"] <- NULL
- if(!is.null(dots[["main"]])) dots["main"] <- NULL
- if(!is.null(dots[["sub"]])) dots["sub"] <- NULL
if(!is.null(dots[["xlab"]])) dots["xlab"] <- NULL
if(!is.null(dots[["ylab"]])) dots["ylab"] <- NULL
- if(!is.null(dots[["xlim"]])) dots["xlim"] <- NULL
- if(!is.null(dots[["ylim"]])) dots["ylim"] <- NULL
- if(is.null(dots[["colA"]])) dots$"colA" <- grey(0.5)
- if(is.null(dots[["colB"]])) dots$"colB" <- par("col")
- if(is.null(dots[["lwdA"]])) dots$"lwdA" <- par("lwd")
- if(is.null(dots[["lwdB"]])) dots$"lwdB" <- 2
- if(is.null(dots[["cex.main"]])) dots$"cex.main" <- 0.8
+ dotsP <- dotsL <- dotsT <- dots
-
+ print(dots)
+ print(lwd)
e1 <- L2Fam at distribution
if(!is(e1, "UnivariateDistribution") | is(e1, "CondDistribution"))
stop("not yet implemented")
@@ -42,7 +46,81 @@
trafo <- L2Fam at param@trafo
dims <- nrow(trafo)
+
+
+ mainL <- FALSE
+ subL <- FALSE
+ lineT <- NA
+
+ .mpresubs <- function(inx)
+ distr:::.presubs(inx, c("%C", "%D", "%A"),
+ c(as.character(class(object)[1]),
+ as.character(date()),
+ as.character(deparse(objectc))))
+
+ if (hasArg(main)){
+ mainL <- TRUE
+ if (is.logical(main)){
+ if (!main) mainL <- FALSE
+ else
+ main <- gettextf("Plot for IC %%A") ###
+ ### double %% as % is special for gettextf
+ }
+ main <- .mpresubs(main)
+ if (mainL) {
+ if(missing(tmar))
+ tmar <- 5
+ if(missing(cex.inner))
+ cex.inner <- .65
+ lineT <- 0.6
+ }
+ }
+ if (hasArg(sub)){
+ subL <- TRUE
+ if (is.logical(sub)){
+ if (!sub) subL <- FALSE
+ else sub <- gettextf("generated %%D")
+ ### double %% as % is special for gettextf
+ }
+ sub <- .mpresubs(sub)
+ if (subL)
+ if (missing(bmar)) bmar <- 6
+ }
+ innerParam <- paste(gettext("\nwith main parameter ("),
+ paste(round(L2Fam at param@main, 3),
+ collapse = ", "),
+ ")", sep = "")
+ if(!is.null(L2Fam at param@nuisance))
+ innerParam <- paste(innerParam,
+ gettext("\nand nuisance parameter ("),
+ paste(round(L2Fam at param@nuisance, 3),
+ collapse = ", "),
+ ")", sep ="")
+ if(!is.null(L2Fam at param@fixed))
+ innerParam <- paste(innerParam,
+ gettext("\nand fixed known parameter ("),
+ paste(round(L2Fam at param@fixed, 3),
+ collapse = ", "),
+ ")", sep ="")
+
+ if(!is.logical(inner)){
+ if(!is.character(inner))
+ stop("Argument 'inner' must either be 'logical' or a character vector")
+ innerT <- rep(inner,length.out=1+dims)
+ innerL <- TRUE
+ }else{if(any(is.na(inner))||any(!inner)) {
+ innerT <- rep("",1+dims); innerL <- FALSE
+ }else{innerL <- TRUE
+ innerT <- paste(c( paste(gettext("Absolute information of (partial) IC for"),
+ name(L2Fam)[1], sep =""),
+ paste(gettext("Relative information of \ncomponent "),
+ 1:dims,
+ gettext("of (partial) IC\nfor "),
+ name(L2Fam)[1], sep ="")), innerParam)
+ }
+ }
+
QFc <- diag(dims)
if(is(object,"ContIC") & dims>1 )
{if (is(normtype(object),"QFNorm")) QFc <- QuadForm(normtype(object))
@@ -65,80 +143,92 @@
absInfo <- t(IC1) %*% QF %*% IC1
absInfo <- sapply(x.vec, absInfo at Map[[1]])
- dots["col"] <- dots[["colA"]]
- dots["lwd"] <- dots[["lwdA"]]
+
+ w0 <- options("warn")
+ options(warn = -1)
+ opar <- par()
+ if (!withSweave)
+ devNew()
+
+ omar <- par("mar")
+ parArgs <- list(mar = c(bmar,omar[2],tmar,omar[4]))
+ do.call(par,args=parArgs)
+
+
+
+ dotsP["col"] <- NULL
+ dotsP["lwd"] <- NULL
+ if(!hasArg(ylim)) dots["ylim"] <- c(0, 2*max(absInfo, na.rm = TRUE))
+
do.call(plot, args=c(list(x.vec, absInfoClass, type = plty,
- lty = "dashed",
- ylim = c(0, 2*max(absInfo, na.rm = TRUE)), xlab = "x",
- ylab = "absolute information"), dots))
- dots["col"] <- dots[["colB"]]
- dots["lwd"] <- dots[["lwdB"]]
+ lty = "dashed", col = colI, lwd = lwdI,
+ xlab = "x",
+ ylab = "absolute information"), dotsP))
do.call(lines, args=c(list(x.vec, absInfo, type = plty, lty = lty),
- dots))
+ dotsL))
legend(max(x.vec), 0, xjust = 1, yjust = 0,
- legend = c("class. opt. IC"), lty = "dashed",
- col = c(dots[["colA"]]), cex=0.75)
+ legend = c("class. opt. IC", objectc),
+ lty = c(lty,"dashed"), col = c(colI, col),
+ lwd=c(lwdI, lwd), cex = 0.75)
- if(is.null(L2Fam at param@nuisance))
- do.call(title, args=c(list(paste("Absolute information of (partial) IC for",
- name(L2Fam)[1], "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "), ")")),
- dots))
- else
- do.call(title, args=c(list(paste("Absolute information of (partial) IC for",
- name(L2Fam)[1], "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3), collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3), collapse = ", "),
- ")")), dots))
-
+ dotsT["main"] <- NULL
+ dotsT["cex.main"] <- NULL
+ dotsT["col.main"] <- NULL
+ dotsT["line"] <- NULL
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[1]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
+
if(dims > 1){
+ dotsP["ylim"] <- NULL
+ dotsL["ylim"] <- NULL
+ dotsT["ylim"] <- NULL
nrows <- trunc(sqrt(dims))
ncols <- ceiling(dims/nrows)
- w0 <- options("warn")
- options(warn = -1)
- opar <- par()
- devNew()
- par(mfrow = c(nrows, ncols))
+ if (!withSweave)
+ devNew()
+ if(mfColRow)
+ parArgs <- c(parArgs,list(mfrow = c(nrows, ncols)))
+
+ do.call(par,args=parArgs)
+
IC1.i.5 <- QF.5%*%IC1
classIC.i.5 <- QFc.5%*%classIC
for(i in 1:dims){
y.vec <- sapply(x.vec, IC1.i.5 at Map[[i]])^2/absInfo
- dots["col"] <- dots[["colB"]]
- dots["lwd"] <- dots[["lwdA"]]
do.call(plot, args=c(list(x.vec, y.vec, type = plty,
lty = lty, xlab = "x",
ylab = "relative information",
- ylim = c(0, 1.1)), dots))
+ ylim = c(0, 1.1),
+ col = colI, lwd = lwdI), dotsP))
yc.vec <- sapply(x.vec, classIC.i.5 at Map[[i]])^2/absInfoClass
- dots["col"] <- dots[["colA"]]
- dots["lwd"] <- dots[["lwdB"]]
do.call(lines, args=c(list(x.vec, yc.vec, type = plty,
- lty = "dashed"),dots))
- legend(max(x.vec), 1.1, xjust = 1, cex = 0.6,
- legend = c("class. opt. IC"), lty = "dashed",
- col = c(dots[["colA"]]))
- if(is.null(L2Fam at param@nuisance))
- do.call(title, args=c(list(paste("Relative information of\ncomponent",
- i, "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3),
- collapse = ", "), ")")),
- dots))
- else
- do.call(title, args=c(list(paste("Relative information of\ncomponent",
- i, "of (partial) IC\nfor", name(L2Fam)[1],
- "\nwith main parameter (",
- paste(round(L2Fam at param@main, 3),
- collapse = ", "),
- ")\nand nuisance parameter (",
- paste(round(L2Fam at param@nuisance, 3),
- collapse = ", "), ")")),
- dots))
+ lty = "dashed"), dotsL))
+ legend(max(x.vec), 1.1, xjust = 1,
+ legend = c("class. opt. IC", objectc), lty = c(lty,"dashed"),
+ col = c(colI, col), lwd=c(lwdI, lwd),
+ cex = 0.6)
+ if(innerL)
+ do.call(title, args=c(list(main = innerT[1+i]), dotsT,
+ line = lineT, cex.main = cex.inner, col.main = col.inner))
}
}
+ if(!hasArg(cex.main)) cex.main <- par("cex.main") else cex.main <- dots$"cex.main"
+ if(!hasArg(col.main)) col.main <- par("col.main") else col.main <- dots$"col.main"
+ if (mainL)
+ mtext(text = main, side = 3, cex = cex.main, adj = .5,
+ outer = TRUE, padj = 1.4, col = col.main)
+
+ if(!hasArg(cex.sub)) cex.sub <- par("cex.sub") else cex.sub <- dots$"cex.sub"
+ if(!hasArg(col.sub)) col.sub <- par("col.sub") else col.sub <- dots$"col.sub"
+ if (subL)
+ mtext(text = sub, side = 1, cex = cex.sub, adj = .5,
+ outer = TRUE, line = -1.6, col = col.sub)
+
+
par(opar)
options(w0)
}
})
+
\ No newline at end of file
Modified: pkg/RobAStBase/chm/00Index.html
===================================================================
--- pkg/RobAStBase/chm/00Index.html 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/00Index.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -10,7 +10,7 @@
<param name="keyword" value=".. contents">
</object>
-<h2>Help pages for package ‘RobAStBase’ version 0.1.0</h2>
+<h2>Help pages for package ‘RobAStBase’ version 0.1.3</h2>
<p align="center">
<a href="#A">A</a>
@@ -252,8 +252,10 @@
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="infoPlot.html">infoPlot</a></td>
<td>Plot absolute and relative information</td></tr>
-<tr><td width="25%"><a href="IC-class.html">infoPlot,IC-method</a></td>
-<td>Influence curve</td></tr>
+<tr><td width="25%"><a href="infoPlot.html">infoPlot,IC-method</a></td>
+<td>Plot absolute and relative information</td></tr>
+<tr><td width="25%"><a href="infoPlot.html">infoPlot-methods</a></td>
+<td>Plot absolute and relative information</td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">Infos</a></td>
<td>Influence curve</td></tr>
<tr><td width="25%"><a href="InfluenceCurve-class.html">Infos,InfluenceCurve-method</a></td>
@@ -416,8 +418,12 @@
<td>ALEstimate-class.</td></tr>
<tr><td width="25%"><a href="ALEstimate-class.html">pIC,ALEstimate-method</a></td>
<td>ALEstimate-class.</td></tr>
-<tr><td width="25%"><a href="IC-class.html">plot,IC,ANY-method</a></td>
-<td>Influence curve</td></tr>
+<tr><td width="25%"><a href="plot-methods.html">plot</a></td>
+<td>Methods for Function plot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="plot-methods.html">plot,IC,missing-method</a></td>
+<td>Methods for Function plot in Package ‘RobAStBase’ </td></tr>
+<tr><td width="25%"><a href="plot-methods.html">plot-methods</a></td>
+<td>Methods for Function plot in Package ‘RobAStBase’ </td></tr>
</table>
<h2><a name="R">-- R --</a></h2>
Modified: pkg/RobAStBase/chm/IC-class.html
===================================================================
--- pkg/RobAStBase/chm/IC-class.html 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/IC-class.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -16,8 +16,6 @@
<param name="keyword" value="R: checkIC,IC,L2ParamFamily-method">
<param name="keyword" value="R: evalIC,IC,numeric-method">
<param name="keyword" value="R: evalIC,IC,matrix-method">
-<param name="keyword" value="R: infoPlot,IC-method">
-<param name="keyword" value="R: plot,IC,ANY-method">
<param name="keyword" value="R: show,IC-method">
<param name="keyword" value=" Influence curve">
</object>
@@ -160,6 +158,6 @@
</script>
-<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.3 <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-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/RobAStBase.hhp 2008-10-29 18:26:14 UTC (rev 184)
@@ -56,3 +56,4 @@
makeIC-methods.html
oneStepEstimator.html
optIC.html
+plot-methods.html
Modified: pkg/RobAStBase/chm/RobAStBase.toc
===================================================================
--- pkg/RobAStBase/chm/RobAStBase.toc 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/RobAStBase.toc 2008-10-29 18:26:14 UTC (rev 184)
@@ -395,9 +395,13 @@
</OBJECT>
<LI> <OBJECT type="text/sitemap">
<param name="Name" value="infoPlot,IC-method">
-<param name="Local" value="IC-class.html">
+<param name="Local" value="infoPlot.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="infoPlot-methods">
+<param name="Local" value="infoPlot.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Infos">
<param name="Local" value="InfluenceCurve-class.html">
</OBJECT>
@@ -682,10 +686,18 @@
<param name="Local" value="ALEstimate-class.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
-<param name="Name" value="plot,IC,ANY-method">
-<param name="Local" value="IC-class.html">
+<param name="Name" value="plot">
+<param name="Local" value="plot-methods.html">
</OBJECT>
<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="plot,IC,missing-method">
+<param name="Local" value="plot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
+<param name="Name" value="plot-methods">
+<param name="Local" value="plot-methods.html">
+</OBJECT>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="radius">
<param name="Local" value="Neighborhood-class.html">
</OBJECT>
@@ -983,6 +995,10 @@
<param name="Local" value="MEstimate-class.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>
+<LI> <OBJECT type="text/sitemap">
<param name="Name" value="Neighborhood">
<param name="Local" value="Neighborhood-class.html">
</OBJECT>
Modified: pkg/RobAStBase/chm/comparePlot.html
===================================================================
--- pkg/RobAStBase/chm/comparePlot.html 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/comparePlot.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -27,7 +27,12 @@
<pre>
comparePlot(obj1, obj2, ... )
## S4 method for signature 'IC, IC':
-comparePlot(obj1, obj2, obj3 = NULL, obj4 = NULL, ... )
+comparePlot(obj1, obj2, obj3 = NULL, obj4 = NULL,
+ ..., 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>
@@ -46,11 +51,72 @@
<tr valign="top"><td><code>obj4</code></td>
<td>
optional: object of class <code>"InfluenceCurve"</code> to be compared with <code>obj1</code></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>
further arguments to be passed to <code>plot</code></td></tr>
</table>
+<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>"%C1"</code>,<code>"%C2"</code>,[<code>"%C3"</code>,][<code>"%C4"</code>]</dt><dd>class of argument
+<code>obj<i></code>, i=1,..4</dd>
+<dt><code>"%A1"</code>,<code>"%A2"</code>,[<code>"%A3"</code>,][<code>"%A4"</code>]</dt><dd>deparsed argument
+<code>obj<i></code>, i=1,..4</dd>
+<dt><code>"%D"</code></dt><dd>time/date-string when the plot was generated</dd>
+</ul>
+
<h3>Author(s)</h3>
<p>
@@ -99,6 +165,6 @@
</script>
-<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.3 <a href="00Index.html">Index]</a></div>
</body></html>
Modified: pkg/RobAStBase/chm/infoPlot.html
===================================================================
--- pkg/RobAStBase/chm/infoPlot.html 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/chm/infoPlot.html 2008-10-29 18:26:14 UTC (rev 184)
@@ -6,6 +6,8 @@
<table width="100%"><tr><td>infoPlot(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: infoPlot">
+<param name="keyword" value="R: infoPlot-methods">
+<param name="keyword" value="R: infoPlot,IC-method">
<param name="keyword" value=" Plot absolute and relative information">
</object>
@@ -23,7 +25,14 @@
<h3>Usage</h3>
<pre>
-infoPlot(object, ...)
+infoPlot(object, ...)
+## S4 method for signature 'IC':
+infoPlot(object, ..., withSweave = getdistrOption("withSweave"),
+ colI = grey(0.5), lwdI = 0.7*par("lwd"),
+ 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>
@@ -33,6 +42,46 @@
<tr valign="top"><td><code>object</code></td>
<td>
object of class <code>"InfluenceCurve"</code> </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>colI</code></td>
+<td>
+color of the classically optimal IC</td></tr>
+<tr valign="top"><td><code>lwdI</code></td>
+<td>
+linewidth of the classically optimal IC</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>
@@ -47,10 +96,31 @@
of an IC. The relative information is defined as the
absolute information of one component with respect to the
absolute information of the whole IC; confer Section 8.1
-of Kohl (2005).
+of Kohl (2005). <br>
</p>
+<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>Author(s)</h3>
<p>
@@ -92,6 +162,6 @@
</script>
-<hr><div align="center">[Package <em>RobAStBase</em> version 0.1.0 <a href="00Index.html">Index]</a></div>
+<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-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/man/IC-class.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -11,8 +11,6 @@
\alias{checkIC,IC,L2ParamFamily-method}
\alias{evalIC,IC,numeric-method}
\alias{evalIC,IC,matrix-method}
-\alias{infoPlot,IC-method}
-\alias{plot,IC,ANY-method}
\alias{show,IC-method}
\title{Influence curve}
Modified: pkg/RobAStBase/man/comparePlot.Rd
===================================================================
--- pkg/RobAStBase/man/comparePlot.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/man/comparePlot.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -10,15 +10,61 @@
}
\usage{
comparePlot(obj1, obj2, ... )
-\S4method{comparePlot}{IC,IC}(obj1, obj2, obj3 = NULL, obj4 = NULL, ... )
+\S4method{comparePlot}{IC,IC}(obj1, obj2, obj3 = NULL, obj4 = NULL,
+ ..., 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{obj1}{ object of class \code{"InfluenceCurve"} }
\item{obj2}{ object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
\item{obj3}{ optional: object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
\item{obj4}{ optional: object of class \code{"InfluenceCurve"} to be compared with \code{obj1}}
+ \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 arguments to be passed to \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{"\%C1"},\code{"\%C2"},[\code{"\%C3"},][\code{"\%C4"}]}{class of argument
+ \code{obj<i>}, i=1,..4}
+\item{\code{"\%A1"},\code{"\%A2"},[\code{"\%A3"},][\code{"\%A4"}]}{deparsed argument
+ \code{obj<i>}, i=1,..4}
+\item{\code{"\%D"}}{time/date-string when the plot was generated}
+}
+}
+
%\value{}
\references{
Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}.
Modified: pkg/RobAStBase/man/infoPlot.Rd
===================================================================
--- pkg/RobAStBase/man/infoPlot.Rd 2008-10-24 09:33:45 UTC (rev 183)
+++ pkg/RobAStBase/man/infoPlot.Rd 2008-10-29 18:26:14 UTC (rev 184)
@@ -1,15 +1,41 @@
\name{infoPlot}
\alias{infoPlot}
+\alias{infoPlot-methods}
+\alias{infoPlot,IC-method}
\title{Plot absolute and relative information}
\description{
Plot absolute and relative information of influence curves.
}
\usage{
-infoPlot(object, ...)
+infoPlot(object, ...)
+\S4method{infoPlot}{IC}(object, ..., withSweave = getdistrOption("withSweave"),
+ colI = grey(0.5), lwdI = 0.7*par("lwd"),
+ 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{object}{ object of class \code{"InfluenceCurve"} }
+ \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{colI}{color of the classically optimal IC}
+ \item{lwdI}{linewidth of the classically optimal IC}
+ \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{
@@ -17,8 +43,29 @@
of an IC. The relative information is defined as the
absolute information of one component with respect to the
absolute information of the whole IC; confer Section 8.1
- of Kohl (2005).
+ of Kohl (2005). \cr
+
+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}
}
+}
%\value{}
\references{
Kohl, M. (2005) \emph{Numerical Contributions to the Asymptotic Theory of Robustness}.
Modified: pkg/RobLox/chm/RobLox.chm
===================================================================
(Binary files differ)
More information about the Robast-commits
mailing list