R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for a HTML browser interface to help. Type 'q()' to quit R. > ### *
> ### > attach(NULL, name = "CheckExEnv") > assign(".CheckExEnv", as.environment(2), pos = length(search())) # base > ## add some hooks to label plot pages for base and grid graphics > setHook("plot.new", ".newplot.hook") > setHook("persp", ".newplot.hook") > setHook("grid.newpage", ".gridplot.hook") > > assign("cleanEx", + function(env = .GlobalEnv) { + rm(list = ls(envir = env, all.names = TRUE), envir = env) + RNGkind("default", "default") + set.seed(1) + options(warn = 1) + delayedAssign("T", stop("T used instead of TRUE"), + assign.env = .CheckExEnv) + delayedAssign("F", stop("F used instead of FALSE"), + assign.env = .CheckExEnv) + sch <- search() + newitems <- sch[! sch %in% .oldSearch] + for(item in rev(newitems)) + eval(substitute(detach(item), list(item=item))) + missitems <- .oldSearch[! .oldSearch %in% sch] + if(length(missitems)) + warning("items ", paste(missitems, collapse=", "), + " have been removed from the search path") + }, + env = .CheckExEnv) > assign("..nameEx", "__{must remake R-ex/*.R}__", env = .CheckExEnv) # for now > assign("ptime", proc.time(), env = .CheckExEnv) > grDevices::postscript("reldist-Examples.ps") > assign("par.postscript", graphics::par(no.readonly = TRUE), env = .CheckExEnv) > options(contrasts = c(unordered = "contr.treatment", ordered = "contr.poly")) > options(warn = 1) > library('reldist') Permanent wage growth in the NLS 1966 cohort Version NLS created on Permanent wage growth in the NLS 1966 cohort. copyright (c) 2003, Mark S. Handcock, University of Washington Martina Morris, University of Washington Type help(package="reldist") to get started. > > assign(".oldSearch", search(), env = .CheckExEnv) > assign(".oldNS", loadedNamespaces(), env = .CheckExEnv) > cleanEx(); ..nameEx <- "reldist" > > ### * reldist > > flush(stderr()); flush(stdout()) > > ### Name: reldist > ### Title: Inference for Relative Distributions > ### Aliases: reldist > ### Keywords: distribution > > ### ** Examples > > # > # First load the data. > # > > data(nls, package="reldist") > > # > # A simple example comparing permanent wages of the original to the > # recent cohort in the NLS. See H&M (1999) for details. > > reldist(y=recent$chpermwage,yo=original$chpermwage) This is mgcv 1.3-1 Smoothing using 0.01023483 > > # > # A more sophisticated version of the same. > # > > reldist(y=recent$chpermwage, yo=original$chpermwage, + yowgt=original$wgt, ywgt=recent$wgt, + bar=TRUE, + smooth=0.1, + yolabs=seq(-1, 3, by=0.5), + ylim=c(0, 3.0),cex=0.8, + ylab="Relative Density", + xlab="Proportion of the Original Cohort") Warning: non-integer x = 139.264435 Warning: non-integer x = 56.214203 Warning: non-integer x = 81.266848 Warning: non-integer x = 77.164064 Warning: non-integer x = 55.902336 Warning: non-integer x = 28.847352 Warning: non-integer x = 38.083914 Warning: non-integer x = 22.911842 Warning: non-integer x = 42.292327 Warning: non-integer x = 24.269062 Warning: non-integer x = 16.275185 Warning: non-integer x = 46.339240 Warning: non-integer x = 47.848528 Warning: non-integer x = 16.946477 Warning: non-integer x = 24.192113 Warning: non-integer x = 16.810936 Warning: non-integer x = 40.597703 Warning: non-integer x = 26.497973 Warning: non-integer x = 23.444595 Warning: non-integer x = 25.115800 Warning: non-integer x = 18.404057 Warning: non-integer x = 17.011200 Warning: non-integer x = 14.894040 Warning: non-integer x = 22.457135 Warning: non-integer x = 17.570731 Warning: non-integer x = 24.041263 Warning: non-integer x = 14.225373 Warning: non-integer x = 23.650721 Warning: non-integer x = 19.998861 Warning: non-integer x = 23.464152 Warning: non-integer x = 14.764665 Warning: non-integer x = 15.624735 Warning: non-integer x = 14.298946 Warning: non-integer x = 12.287932 Warning: non-integer x = 16.842208 Warning: non-integer x = 13.222588 Warning: non-integer x = 17.069721 Warning: non-integer x = 9.033178 Warning: non-integer x = 14.017716 Warning: non-integer x = 9.621783 Warning: non-integer x = 5.036696 Warning: non-integer x = 15.292992 Warning: non-integer x = 10.531835 Warning: non-integer x = 18.824495 Warning: non-integer x = 14.763840 Warning: non-integer x = 8.097593 Warning: non-integer x = 17.602013 Warning: non-integer x = 16.571880 Warning: non-integer x = 23.417000 Warning: non-integer x = 16.562075 Warning: non-integer x = 6.372568 Warning: non-integer x = 15.212301 Warning: non-integer x = 14.763477 Warning: non-integer x = 13.739640 Warning: non-integer x = 13.353577 Warning: non-integer x = 14.512362 Warning: non-integer x = 13.745090 Warning: non-integer x = 20.787475 Warning: non-integer x = 24.700145 Warning: non-integer x = 9.199532 Warning: non-integer x = 18.785873 Warning: non-integer x = 7.992932 Warning: non-integer x = 19.914973 Warning: non-integer x = 17.174332 Warning: non-integer x = 18.188600 Warning: non-integer x = 15.274076 Warning: non-integer x = 12.993273 Warning: non-integer x = 10.096926 Warning: non-integer x = 8.220636 Warning: non-integer x = 14.074624 Warning: non-integer x = 14.670440 Warning: non-integer x = 8.331124 Warning: non-integer x = 14.252971 Warning: non-integer x = 3.349117 Warning: non-integer x = 15.440684 Warning: non-integer x = 13.948787 Warning: non-integer x = 13.688317 Warning: non-integer x = 15.592511 Warning: non-integer x = 14.010352 Warning: non-integer x = 22.016996 Warning: non-integer x = 16.542906 Warning: non-integer x = 12.984732 Warning: non-integer x = 24.991252 Warning: non-integer x = 13.127210 Warning: non-integer x = 11.072164 Warning: non-integer x = 8.472304 Warning: non-integer x = 23.301249 Warning: non-integer x = 24.387749 Warning: non-integer x = 16.024427 Warning: non-integer x = 17.320723 Warning: non-integer x = 8.299960 Warning: non-integer x = 8.835011 Warning: non-integer x = 10.132579 Warning: non-integer x = 23.963467 Warning: non-integer x = 11.535496 Warning: non-integer x = 14.371816 Warning: non-integer x = 14.540651 Warning: non-integer x = 34.351498 Warning: non-integer x = 33.841608 Warning: non-integer x = 21.013132 Smoothing using 0.1 > > # > # A CDF version. > # > > reldist(y=recent$chpermwage, yo=original$chpermwage, + yowgt=original$wgt, ywgt=recent$wgt, + cdfplot=TRUE, + smooth=0.4, + yolabs=seq(-1,3,by=0.5), + ylabs=seq(-1,3,by=0.5), + cex=0.8, + ylab="proportion of the recent cohort", + xlab="proportion of the original cohort") Warning: non-integer x = 139.264435 Warning: non-integer x = 56.214203 Warning: non-integer x = 81.266848 Warning: non-integer x = 77.164064 Warning: non-integer x = 55.902336 Warning: non-integer x = 28.847352 Warning: non-integer x = 38.083914 Warning: non-integer x = 22.911842 Warning: non-integer x = 42.292327 Warning: non-integer x = 24.269062 Warning: non-integer x = 16.275185 Warning: non-integer x = 46.339240 Warning: non-integer x = 47.848528 Warning: non-integer x = 16.946477 Warning: non-integer x = 24.192113 Warning: non-integer x = 16.810936 Warning: non-integer x = 40.597703 Warning: non-integer x = 26.497973 Warning: non-integer x = 23.444595 Warning: non-integer x = 25.115800 Warning: non-integer x = 18.404057 Warning: non-integer x = 17.011200 Warning: non-integer x = 14.894040 Warning: non-integer x = 22.457135 Warning: non-integer x = 17.570731 Warning: non-integer x = 24.041263 Warning: non-integer x = 14.225373 Warning: non-integer x = 23.650721 Warning: non-integer x = 19.998861 Warning: non-integer x = 23.464152 Warning: non-integer x = 14.764665 Warning: non-integer x = 15.624735 Warning: non-integer x = 14.298946 Warning: non-integer x = 12.287932 Warning: non-integer x = 16.842208 Warning: non-integer x = 13.222588 Warning: non-integer x = 17.069721 Warning: non-integer x = 9.033178 Warning: non-integer x = 14.017716 Warning: non-integer x = 9.621783 Warning: non-integer x = 5.036696 Warning: non-integer x = 15.292992 Warning: non-integer x = 10.531835 Warning: non-integer x = 18.824495 Warning: non-integer x = 14.763840 Warning: non-integer x = 8.097593 Warning: non-integer x = 17.602013 Warning: non-integer x = 16.571880 Warning: non-integer x = 23.417000 Warning: non-integer x = 16.562075 Warning: non-integer x = 6.372568 Warning: non-integer x = 15.212301 Warning: non-integer x = 14.763477 Warning: non-integer x = 13.739640 Warning: non-integer x = 13.353577 Warning: non-integer x = 14.512362 Warning: non-integer x = 13.745090 Warning: non-integer x = 20.787475 Warning: non-integer x = 24.700145 Warning: non-integer x = 9.199532 Warning: non-integer x = 18.785873 Warning: non-integer x = 7.992932 Warning: non-integer x = 19.914973 Warning: non-integer x = 17.174332 Warning: non-integer x = 18.188600 Warning: non-integer x = 15.274076 Warning: non-integer x = 12.993273 Warning: non-integer x = 10.096926 Warning: non-integer x = 8.220636 Warning: non-integer x = 14.074624 Warning: non-integer x = 14.670440 Warning: non-integer x = 8.331124 Warning: non-integer x = 14.252971 Warning: non-integer x = 3.349117 Warning: non-integer x = 15.440684 Warning: non-integer x = 13.948787 Warning: non-integer x = 13.688317 Warning: non-integer x = 15.592511 Warning: non-integer x = 14.010352 Warning: non-integer x = 22.016996 Warning: non-integer x = 16.542906 Warning: non-integer x = 12.984732 Warning: non-integer x = 24.991252 Warning: non-integer x = 13.127210 Warning: non-integer x = 11.072164 Warning: non-integer x = 8.472304 Warning: non-integer x = 23.301249 Warning: non-integer x = 24.387749 Warning: non-integer x = 16.024427 Warning: non-integer x = 17.320723 Warning: non-integer x = 8.299960 Warning: non-integer x = 8.835011 Warning: non-integer x = 10.132579 Warning: non-integer x = 23.963467 Warning: non-integer x = 11.535496 Warning: non-integer x = 14.371816 Warning: non-integer x = 14.540651 Warning: non-integer x = 34.351498 Warning: non-integer x = 33.841608 Warning: non-integer x = 21.013132 Smoothing using 0.4 > > > > cleanEx(); ..nameEx <- "resplot" > > ### * resplot > > flush(stderr()); flush(stdout()) > > ### Name: resplot > ### Title: Relative distribution plot to a Standard Normal > ### Aliases: resplot > ### Keywords: hplot distribution > > ### ** Examples > > y <- rnorm(2000) > resplot(y) This is mgcv 1.3-1 Smoothing using 1.281549 Smoothing using 1.281549 > data(precip) > resplot(precip, ylab = "Precipitation [in/yr] for 70 US cities") Smoothing using 0.008829239 Smoothing using 0.008829239 > > > > ### *