* using log directory 'd:/Rcompile/CRANpkg/local/2.7/ROptEst.Rcheck' * using R version 2.7.2 (2008-08-25) * using session charset: ISO8859-1 * checking for file 'ROptEst/DESCRIPTION' ... OK * this is package 'ROptEst' version '0.6.2' * checking package name space information ... OK * checking package dependencies ... OK * checking if this is a source package ... OK * checking whether package 'ROptEst' can be installed ... OK * checking package directory ... OK * checking for portable file names ... OK * checking DESCRIPTION meta-information ... OK * checking top-level files ... OK * checking index information ... OK * checking package subdirectories ... OK * checking R files for non-ASCII characters ... OK * checking R files for syntax errors ... OK * checking whether the package can be loaded ... OK * checking whether the package can be loaded with stated dependencies ... OK * checking whether the name space can be loaded with stated dependencies ... OK * checking for unstated dependencies in R code ... OK * checking S3 generic/method consistency ... OK * checking replacement functions ... OK * checking foreign function calls ... OK * checking R code for possible problems ... OK * checking Rd files ... OK * checking Rd cross-references ... OK * checking for missing documentation entries ... OK * checking for code/documentation mismatches ... OK * checking Rd \usage sections ... OK * creating ROptEst-Ex.R ... OK * checking examples ... ERROR Running examples in 'ROptEst-Ex.R' failed. The error most likely occurred in: > ### * 0ROptEst-package > > flush(stderr()); flush(stdout()) > > ### Name: ROptEst-package > ### Title: Optimally robust estimation > ### Aliases: ROptEst-package ROptEst > ### Keywords: package > > ### ** Examples > > library(ROptEst) > > ## 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 from package distrMod > MLest <- MLEstimator(x, PoisFamily()) > MLest Evaluations of Maximum likelihood estimate: ------------------------------------------- An object of class "Estimate" generated by call MLEstimator(x = x, ParamFamily = PoisFamily()) samplesize: 2608 estimate: lambda 3.87154908 (0.03852908) asymptotic (co)variance (multiplied with samplesize): [1] 3.871549 Criterium: negative log-likelihood 5352.105 > ## confidence interval based on CLT > confint(MLest) A[n] asymptotic (CLT-based) confidence interval: 2.5 % 97.5 % [1,] 3.796033 3.947065 Type of estimator: Maximum likelihood estimate samplesize: 2608 Call by which estimate was produced: MLEstimator(x = x, ParamFamily = PoisFamily()) > > ## compute optimally (w.r.t to MSE) robust estimator (unknown contamination) > robest <- roptest(x, PoisFamily(), eps.upper = 0.1, steps = 3) > estimate(robest) lambda 3.908322 > ## check influence curve > checkIC(pIC(robest)) precision of centering: 4.790607e-17 precision of Fisher consistency: lambda lambda -1.968968e-06 maximum deviation 1.968968e-06 > ## plot influence curve > plot(pIC(robest)) > ## confidence interval based on LAN - neglecting bias > confint(robest) A[n] asymptotic (CLT-based) confidence interval: 2.5 % 97.5 % lambda 3.826170 3.990474 Type of estimator: 3-step estimate samplesize: 2608 Call by which estimate was produced: roptest(x = x, L2Fam = PoisFamily(), eps.upper = 0.1, steps = 3) > ## confidence interval based on LAN - including bias > confint(robest, method = symmetricBias()) Error in object$coefficients : $ operator not defined for this S4 class Calls: confint ... -> confint.default -> coef -> coef.default Execution halted