R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing 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 an HTML browser interface to help. Type 'q()' to quit R. > x <- c(10926126 + ,204514194 + ,20655116 + ,7365385 + ,3918657 + ,38494822 + ,90838469 + ,8652220 + ,47482511 + ,5438191 + ,5986958 + ,59447386 + ,37113055 + ,23669000 + ,5572851 + ,7784176 + ,5270786 + ,43583746 + ,19903814 + ,7799064 + ,35989363 + ,29552113 + ,21506498 + ,7859684 + ,13006528 + ,89426565 + ,40100910 + ,38650820 + ,17818636 + ,51985948 + ,24571878 + ,29266469 + ,45758250 + ,22712439 + ,86459749 + ,49703065 + ,65187751 + ,39091325 + ,34955804 + ,33411147 + ,15575837 + ,36116313 + ,27405375 + ,55382004 + ,7478767 + ,72881801 + ,37479848 + ,27030251 + ,109851111 + ,28356456 + ,48560722 + ,23550558 + ,52319234 + ,105318425 + ,17345510 + ,71519979 + ,59398481 + ,36287185 + ,44559709 + ,114256076 + ,23174724 + ,38795625 + ,94992548 + ,48389051 + ,46785833 + ,56063795 + ,27622103 + ,11848002 + ,55125816 + ,42484133 + ,49639784 + ,30006747 + ,19776143 + ,41176416 + ,34933018 + ,39696259 + ,38610419 + ,44404281 + ,44628459 + ,20428394 + ,33494104 + ,15246622 + ,10078039 + ,36293305 + ,11594321 + ,34818494 + ,36718379 + ,19690322 + ,44432615 + ,36047082 + ,40379563 + ,23458900 + ,24477773 + ,21472076 + ,27623825 + ,42834376 + ,21692408 + ,2490957 + ,5804447 + ,4149313 + ,22772412 + ,8743945 + ,6818743 + ,30796414 + ,17697315 + ,62514975 + ,17618542 + ,37171201 + ,29747133 + ,39734547 + ,18636718 + ,30751381 + ,58198919 + ,35172258 + ,30568494 + ,32732200 + ,33931175 + ,38495426 + ,39639276 + ,18386312 + ,49064295 + ,34517203 + ,54053990 + ,56632146 + ,46131933 + ,31138610 + ,37573144 + ,32922160 + ,71766977 + ,21574297 + ,46733768 + ,36375026 + ,50499581 + ,44157910 + ,26780541 + ,71675809 + ,17411180 + ,138402202 + ,89520475 + ,53764838 + ,79645460 + ,45232765 + ,60335913 + ,48862499 + ,56443228 + ,149401818 + ,55893936 + ,39185258 + ,54999796 + ,62938759 + ,130610415 + ,3578335 + ,119330610 + ,43612479 + ,67573324 + ,94235304 + ,52529046 + ,38148241 + ,41788652 + ,21748127 + ,16815679 + ,29140007 + ,80020386 + ,10008346 + ,13150180 + ,33070836 + ,44343752 + ,38761463 + ,63242438 + ,10776650 + ,13079087 + ,42254065 + ,45334695 + ,75653353 + ,37723883 + ,35383003 + ,42169175 + ,37836247 + ,16157265 + ,25483137 + ,24758200 + ,25931415 + ,35341058 + ,59793935 + ,46240034 + ,25935511 + ,29348057 + ,16621841 + ,34421545 + ,30201523 + ,68504357 + ,21760376 + ,68420735 + ,34576094 + ,80009790 + ,23896896 + ,43959624 + ,31004430 + ,96339703 + ,78115427 + ,13773684 + ,79714473 + ,30300608 + ,47400000 + ,29685182 + ,72574056 + ,26667187 + ,56654962 + ,34032684 + ,37819725 + ,25253403 + ,55674674 + ,40214078 + ,32066001 + ,26916456 + ,31925699 + ,36677072 + ,22865662 + ,33057141 + ,87224667 + ,30046064 + ,33705054 + ,94643577 + ,74500253 + ,13705457 + ,81920371 + ,75534971 + ,17984250 + ,26564680 + ,75039636 + ,31345912 + ,156181130 + ,37274013 + ,60344095 + ,19480849 + ,15213215 + ,36646854 + ,37800225 + ,35192619 + ,51894705 + ,43116866 + ,14343902 + ,21908691 + ,34532932 + ,38781483 + ,32523543 + ,46755676 + ,68059999 + ,34481666 + ,23934252 + ,56464201 + ,30381482 + ,41957576 + ,11196762 + ,26125759 + ,45626240 + ,33816982 + ,150181287 + ,33424206 + ,67497380 + ,54983068 + ,83797137 + ,44798290 + ,39918703 + ,22068324 + ,78296142 + ,41187726 + ,87854574 + ,56492262 + ,36294788 + ,102238796 + ,42832378 + ,22490585 + ,53384345 + ,22461067 + ,24979034 + ,106897458 + ,51351879 + ,39339354 + ,32067212 + ,32690247 + ,53227096 + ,59342371 + ,64960155 + ,54764596 + ,102374548 + ,76913992 + ,42777761 + ,40847167 + ,86951363 + ,43659732 + ,76698629 + ,78268805 + ,91261974 + ,65606487 + ,28139063 + ,75895280 + ,51122832 + ,53066011 + ,1354614 + ,30021100 + ,33899609 + ,58995770 + ,57876331 + ,87561067 + ,26724300 + ,76278583 + ,50626962) > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > 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 > m <- mean(x) > e <- median(x) > postscript(file="/var/www/html/rcomp/tmp/14emm1290779069.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > op <- par(mfrow=c(2,1)) > mydensity1 <- density(x,kernel='gaussian',na.rm=TRUE) > plot(mydensity1,main='Density Plot - Gaussian Kernel',xlab='Median (0 -> full line) | Mean (0 -> dashed line)',ylab='density') > abline(v=e,lty=1) > abline(v=m,lty=5) > grid() > myseq <- seq(0.01, 0.99, 0.01) > hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE) > plot(myseq,hd,col=2,main='Harrell-Davis Quantiles',xlab='quantiles',ylab='Median (0 -> full) | Mean (0 -> dashed)') > abline(h=m,lty=5) > abline(h=e,lty=1) > grid() > par(op) > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Median versus Mean',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'mean',header=TRUE) > a<-table.element(a,mean(x)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'median',header=TRUE) > a<-table.element(a,median(x)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/27wls1290779069.tab") > > try(system("convert tmp/14emm1290779069.ps tmp/14emm1290779069.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.951 0.256 2.009