Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_harrell_davis.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationWed, 07 Jul 2010 15:55:54 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Jul/07/t1278518227wgydy1irem8qdkv.htm/, Retrieved Mon, 06 May 2024 07:58:28 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=77944, Retrieved Mon, 06 May 2024 07:58:28 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact206
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Kelly Janbroers -...] [2010-07-07 15:55:54] [413e0fefcf22560c5655fbc122c1a3c2] [Current]
Feedback Forum

Post a new message
Dataseries X:
33
32
31
29
49
48
33
23
24
24
25
27
24
21
15
21
49
48
35
36
51
50
61
63
61
62
58
65
93
94
86
88
102
107
121
127
125
128
117
127
160
162
153
160
177
178
196
212
212
211
204
216
248
250
240
249
275
277
286
302
290
290
277
285
311
300
291
299
332
337
343
360
353
351
341
348
381
358
353
358
399
409
407
419
418
421
414
424
463
437
430
436
474
489
482
492
502
500
493
504
538
516
502
501
541
571
559
569
576
573
562
570
597
573
562
556
600
630
624
634




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77944&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77944&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77944&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Harrell-Davis Quantiles
quantilesvaluestandard error
0.133.05523073426575.71548504746708
0.262.904521650360314.4733373233442
0.3121.82577660402727.2400030978170
0.4201.68746221480936.6331220045253
0.5279.62291407837127.303169331959
0.6341.46292915671626.3721785167921
0.7410.57971322089128.5042369249107
0.8487.08491125876523.7315989667620
0.9558.63907970206214.7886758892593

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 33.0552307342657 & 5.71548504746708 \tabularnewline
0.2 & 62.9045216503603 & 14.4733373233442 \tabularnewline
0.3 & 121.825776604027 & 27.2400030978170 \tabularnewline
0.4 & 201.687462214809 & 36.6331220045253 \tabularnewline
0.5 & 279.622914078371 & 27.303169331959 \tabularnewline
0.6 & 341.462929156716 & 26.3721785167921 \tabularnewline
0.7 & 410.579713220891 & 28.5042369249107 \tabularnewline
0.8 & 487.084911258765 & 23.7315989667620 \tabularnewline
0.9 & 558.639079702062 & 14.7886758892593 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=77944&T=1

[TABLE]
[ROW][C]Harrell-Davis Quantiles[/C][/ROW]
[ROW][C]quantiles[/C][C]value[/C][C]standard error[/C][/ROW]
[ROW][C]0.1[/C][C]33.0552307342657[/C][C]5.71548504746708[/C][/ROW]
[ROW][C]0.2[/C][C]62.9045216503603[/C][C]14.4733373233442[/C][/ROW]
[ROW][C]0.3[/C][C]121.825776604027[/C][C]27.2400030978170[/C][/ROW]
[ROW][C]0.4[/C][C]201.687462214809[/C][C]36.6331220045253[/C][/ROW]
[ROW][C]0.5[/C][C]279.622914078371[/C][C]27.303169331959[/C][/ROW]
[ROW][C]0.6[/C][C]341.462929156716[/C][C]26.3721785167921[/C][/ROW]
[ROW][C]0.7[/C][C]410.579713220891[/C][C]28.5042369249107[/C][/ROW]
[ROW][C]0.8[/C][C]487.084911258765[/C][C]23.7315989667620[/C][/ROW]
[ROW][C]0.9[/C][C]558.639079702062[/C][C]14.7886758892593[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=77944&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77944&T=1

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Harrell-Davis Quantiles
quantilesvaluestandard error
0.133.05523073426575.71548504746708
0.262.904521650360314.4733373233442
0.3121.82577660402727.2400030978170
0.4201.68746221480936.6331220045253
0.5279.62291407837127.303169331959
0.6341.46292915671626.3721785167921
0.7410.57971322089128.5042369249107
0.8487.08491125876523.7315989667620
0.9558.63907970206214.7886758892593



Parameters (Session):
par1 = 0.1 ; par2 = 0.9 ; par3 = 0.1 ;
Parameters (R input):
par1 = 0.1 ; par2 = 0.9 ; par3 = 0.1 ;
R code (references can be found in the software module):
par1 <- as(par1,'numeric')
par2 <- as(par2,'numeric')
par3 <- as(par3,'numeric')
library(Hmisc)
myseq <- seq(par1, par2, par3)
hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE)
bitmap(file='test1.png')
plot(myseq,hd,col=2,main=main,xlab=xlab,ylab=ylab)
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Harrell-Davis Quantiles',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'quantiles',header=TRUE)
a<-table.element(a,'value',header=TRUE)
a<-table.element(a,'standard error',header=TRUE)
a<-table.row.end(a)
length(hd)
for (i in 1:length(hd))
{
a<-table.row.start(a)
a<-table.element(a,as(labels(hd)[i],'numeric'),header=TRUE)
a<-table.element(a,as.matrix(hd[i])[1,1])
a<-table.element(a,as.matrix(attr(hd,'se')[i])[1,1])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')