FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_harrell_davis.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationFri, 18 May 2012 11:09:42 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/May/18/t1337353800i6x3ine8afbrsf1.htm/, Retrieved Fri, 03 May 2024 15:21:20 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=166685, Retrieved Fri, 03 May 2024 15:21:20 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG20114(1)
Estimated Impact134

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2012-05-18 15:09:42] [606e5654d317e57bd58e5a48c9e4e9a9] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
20
25
15
15
25
25
25
21
30
25
20
40
13
30
25
20
25
20
25
20
20
15
15
12
20
5
20
15
25
22
20
22
25
20
20
35
30
25
20
20
20
25
25
15
20
35
25
25
30
23
10
22
25
25
22
30
20
25
25
22
25
25
25
22
25
12
18
20
20
22
30
25
22
20
50
30
25
20
30
22
25
30
22
25
22
22
25
25
25
20
22
15
20
30
20
25
30
35
22
12
30
15
10
30
9
25
20
20
35
25
35
30
12
25
15
25
25
20
20
6
15
40
20
40
25
25
20
15
15
22
24
22
20
25
25
25
35
40
20
22
22
20
25
25
18
25
20
25
30
20
22
35
22
25
25
25
25
22
23
35
15
25
18
22
25
25
28
30
20
25
25
30
22
30
10
10
25
20
22
25
25
15
22
25
25
28
22
30
25
20
25
25
20
30
20
30
50
19
20
28
20
25
35
25
25
15
16
20
20
25
30
20
25
25
25
20
20
25
25
30
22
20
25
25
18
18
20
25
25
30
25
20
25
20
20
20
22
18
22
20
15
25
25
20
25
15
22
25
25
15
12
25
30
22
15
22
25
12
18
30
25
25
40
24
25
15
25
20
25
25
25
20
30
20
25
30
22
25
25
25
50
19
50
25
35
20
20
20
20
20
25
25
25
20
20
20
20
25
18
25
22
22
30
30
8
20
25
30
50
22
20
10
25
25
25
25
18
25
20
25
30
18
20
25
22
22
20
20
25
20
20
20
20
25
20
10
20
25
30
25
50
30
30
50
15
25
25
22
20
22
30
25
18
22
22
30
40
25
20
10
20
9
15
20
15
20
30
12
15
12
20
15
12
25
20
25
25
25
30
20
25
15
15
22
10
15
10
20
25
20
20
38
20
20
20
40
25
25
30
25
10
20
25
12
15
25
20
22
22
20
25
25
25
15
40
20
20
16
25
15
20
25
20
30
50
20
25
20
30
30
25
25
12
25
25
25
20
20
20
15
20
25
15
25
50
30
20
20
25
12
15
20
20
35
22
15
18
30
22
12
12
20
20
15
25
15
20
20
25
18
30
20
25
25
25
20
20
25
20
22
15
15
22
20
10
25
20
20
15
12
20
5
20
15
15
25
25
25
15
25
22
25
20
18
22
25
35
25
25
25
35
30
22
30
50
15
25
24
20
25
25
25
12
15
22
25
25
25
25
15
20
20
15
35
30
20
22
65
20
25
22
20
25
25
20
25
15
20
12
15
10
25
15
30
35
25
25
25
25
25
40
40
25
25
20
25
25
22
25
30
25
25
30
25
25
30
25
25
20
22
22
20
25
22
25
22
40
25
25
25
22
20
35
20
35
25
22
25
25
25
25
25
40
25
30
25
20
25
25
30
22
22
20
15
15
25
25
20
20
15
25
15
20
22
25
15
15
18
5
15
25
18
40
25
25
20
30
20
25
25
25
22
22
25
25
30
25
25
25
25
20
20
25
25
25
25
20
30
25
22
30
20
20
30
25
25
30
20
25
25
24
25
30
18
15
22
22
25
22
22
25
15
20
22
18
35
20
20
20
25
25
30
15
25
22
26
25
20
25
25
25
22
25
25
20
22
30
15
30
25
20
25
25
35
22
20
25
20
20
18
20
22
25
10
20
25
20
20
30
25
20
15
20
25
10
20
25
22
22
25
25
15
25
20
10
25
16
25
35
25
15
25
25
30
25
10
22
20
25
20
20
25
22
18
30
19
25
20
25
20
25
20
22
12
30
12
22
25
25
25
25
30
30
10
22
22
25
20
22
20
25
20
15
25
20
25
20
30
15
40
25
20
22
22
30
20
40
20
25
20
25
20
50
50
25
25
40
30
22
30
20
25
25
30
25
25
20
18
18
28
25
22
15
40
40
12
12
18
12
25
26
18
25
22
15
25
15
15
15
25
15
12
22
20
20
25
20
12
9
15
12
15
25
20
20
15
15
30
21
25
22
22
50
15
25
15
25
22
18
50
20
50
20
20
30
25
20
22
25
50
40
25
25
25
25
30
40
25
30
20

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net

\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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=166685&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=166685&T=0

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

As an alternative you can also use a QR Code:  
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'Herman Ole Andreas Wold' @ wold.wessa.net






Harrell-Davis Quantiles
quantilesvaluestandard error
0.0512.56952089122410.844603824618881
0.115.00007139374610.000330841003470901
0.1516.94784722434971.18642295313764
0.219.93090219759930.146281932400588
0.2519.99999996010870
0.320.00000000243060
0.3520.00484364569820.0152774948701164
0.421.21087483190720.766741332207266
0.4522.00004691628750.00651503412376393
0.522.99177021477981.07013689280987
0.5524.98267716247840.0485757220106025
0.624.99999994553770
0.65250
0.7250
0.7525.00000000000380
0.825.00180702541750.00600848841895861
0.8528.53236855662451.69252291226154
0.930.00014886562620.000560591507156429
0.9535.47741364078961.28406811724589

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.05 & 12.5695208912241 & 0.844603824618881 \tabularnewline
0.1 & 15.0000713937461 & 0.000330841003470901 \tabularnewline
0.15 & 16.9478472243497 & 1.18642295313764 \tabularnewline
0.2 & 19.9309021975993 & 0.146281932400588 \tabularnewline
0.25 & 19.9999999601087 & 0 \tabularnewline
0.3 & 20.0000000024306 & 0 \tabularnewline
0.35 & 20.0048436456982 & 0.0152774948701164 \tabularnewline
0.4 & 21.2108748319072 & 0.766741332207266 \tabularnewline
0.45 & 22.0000469162875 & 0.00651503412376393 \tabularnewline
0.5 & 22.9917702147798 & 1.07013689280987 \tabularnewline
0.55 & 24.9826771624784 & 0.0485757220106025 \tabularnewline
0.6 & 24.9999999455377 & 0 \tabularnewline
0.65 & 25 & 0 \tabularnewline
0.7 & 25 & 0 \tabularnewline
0.75 & 25.0000000000038 & 0 \tabularnewline
0.8 & 25.0018070254175 & 0.00600848841895861 \tabularnewline
0.85 & 28.5323685566245 & 1.69252291226154 \tabularnewline
0.9 & 30.0001488656262 & 0.000560591507156429 \tabularnewline
0.95 & 35.4774136407896 & 1.28406811724589 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=166685&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.05[/C][C]12.5695208912241[/C][C]0.844603824618881[/C][/ROW]
[ROW][C]0.1[/C][C]15.0000713937461[/C][C]0.000330841003470901[/C][/ROW]
[ROW][C]0.15[/C][C]16.9478472243497[/C][C]1.18642295313764[/C][/ROW]
[ROW][C]0.2[/C][C]19.9309021975993[/C][C]0.146281932400588[/C][/ROW]
[ROW][C]0.25[/C][C]19.9999999601087[/C][C]0[/C][/ROW]
[ROW][C]0.3[/C][C]20.0000000024306[/C][C]0[/C][/ROW]
[ROW][C]0.35[/C][C]20.0048436456982[/C][C]0.0152774948701164[/C][/ROW]
[ROW][C]0.4[/C][C]21.2108748319072[/C][C]0.766741332207266[/C][/ROW]
[ROW][C]0.45[/C][C]22.0000469162875[/C][C]0.00651503412376393[/C][/ROW]
[ROW][C]0.5[/C][C]22.9917702147798[/C][C]1.07013689280987[/C][/ROW]
[ROW][C]0.55[/C][C]24.9826771624784[/C][C]0.0485757220106025[/C][/ROW]
[ROW][C]0.6[/C][C]24.9999999455377[/C][C]0[/C][/ROW]
[ROW][C]0.65[/C][C]25[/C][C]0[/C][/ROW]
[ROW][C]0.7[/C][C]25[/C][C]0[/C][/ROW]
[ROW][C]0.75[/C][C]25.0000000000038[/C][C]0[/C][/ROW]
[ROW][C]0.8[/C][C]25.0018070254175[/C][C]0.00600848841895861[/C][/ROW]
[ROW][C]0.85[/C][C]28.5323685566245[/C][C]1.69252291226154[/C][/ROW]
[ROW][C]0.9[/C][C]30.0001488656262[/C][C]0.000560591507156429[/C][/ROW]
[ROW][C]0.95[/C][C]35.4774136407896[/C][C]1.28406811724589[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=166685&T=1

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

As an alternative you can also use a QR Code:  
QR Code


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

Harrell-Davis Quantiles
quantilesvaluestandard error
0.0512.56952089122410.844603824618881
0.115.00007139374610.000330841003470901
0.1516.94784722434971.18642295313764
0.219.93090219759930.146281932400588
0.2519.99999996010870
0.320.00000000243060
0.3520.00484364569820.0152774948701164
0.421.21087483190720.766741332207266
0.4522.00004691628750.00651503412376393
0.522.99177021477981.07013689280987
0.5524.98267716247840.0485757220106025
0.624.99999994553770
0.65250
0.7250
0.7525.00000000000380
0.825.00180702541750.00600848841895861
0.8528.53236855662451.69252291226154
0.930.00014886562620.000560591507156429
0.9535.47741364078961.28406811724589


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
par1 = 0.05 ; par2 = 0.95 ; par3 = 0.05 ;
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')