R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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 an HTML browser interface to help. Type 'q()' to quit R. > x <- c(840,880,930,920,940,880,980,860,900,930,870,1000,870,860,930,980,1010,860,1140,880,800,900,900,1000,890,890,870,1000,1050,790,1160,830,730,950,980,910,840,860,880,1030,1060,770,1140,890,740,860,1050,840,810,830,920,1070,1040,740,1250,850,790,810,1080,760,840,820,900,1010,1080,780,1150,820,790,820,1130,800,890,810,950,1090,1090,850,1200,790,800,850,1230,800,930,700,1030,1040,1000,830,1190,720,810,870,1190,800,970,690,1010,1030,950,830,1150,750,840,880,1210,830) > par3 = '0.01' > par2 = '0.99' > par1 = '0.O1' > ylab = 'Aantal bezoekers Zimmertoren' > xlab = 'Percentielen' > main = 'Harrell-Davis Quantiles' > par3 <- '0.01' > par2 <- '0.99' > par1 <- '0.O1' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa, (2012), Harrell-Davis Quantile Estimator (v1.0.13) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_harrell_davis.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > par1 <- as(par1,'numeric') Warning message: In asMethod(object) : NAs introduced by coercion > par2 <- as(par2,'numeric') > par3 <- as(par3,'numeric') > library(Hmisc) Loading required package: survival Loading required package: splines Attaching package: 'Hmisc' The following object(s) are masked from 'package:survival': untangle.specials The following object(s) are masked from 'package:base': format.pval, round.POSIXt, trunc.POSIXt, units > myseq <- seq(par1, par2, par3) Error in seq.default(par1, par2, par3) : invalid (to - from)/by in seq(.) Calls: seq -> seq.default Execution halted