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 computationThu, 21 Feb 2013 11:09:50 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/Feb/21/t13614630982rcevjkwwhxw0ul.htm/, Retrieved Thu, 02 May 2024 19:46:30 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=206727, Retrieved Thu, 02 May 2024 19:46:30 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsyasmien.naciri@student.kdg.be
Estimated Impact91

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] [] [2013-02-21 16:09:50] [d06e8713ea83045a022ab0926c74dd0b] [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 time3 seconds
R Server'Gwilym Jenkins' @ jenkins.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 & 3 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=206727&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=206727&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=206727&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 time3 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net






Harrell-Davis Quantiles
quantilesvaluestandard error
0.019.481460198747260.531630278100433
0.1115.00253827858880.00915781008035258
0.2119.98981180703890.027716502034098
0.3120.00000009820240
0.4121.62036795182850.538667594485544
0.5123.6838233510181.16373117950554
0.6124.99999999849280
0.71250
0.8125.02210490027560.0598605898389705

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 9.48146019874726 & 0.531630278100433 \tabularnewline
0.11 & 15.0025382785888 & 0.00915781008035258 \tabularnewline
0.21 & 19.9898118070389 & 0.027716502034098 \tabularnewline
0.31 & 20.0000000982024 & 0 \tabularnewline
0.41 & 21.6203679518285 & 0.538667594485544 \tabularnewline
0.51 & 23.683823351018 & 1.16373117950554 \tabularnewline
0.61 & 24.9999999984928 & 0 \tabularnewline
0.71 & 25 & 0 \tabularnewline
0.81 & 25.0221049002756 & 0.0598605898389705 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=206727&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.01[/C][C]9.48146019874726[/C][C]0.531630278100433[/C][/ROW]
[ROW][C]0.11[/C][C]15.0025382785888[/C][C]0.00915781008035258[/C][/ROW]
[ROW][C]0.21[/C][C]19.9898118070389[/C][C]0.027716502034098[/C][/ROW]
[ROW][C]0.31[/C][C]20.0000000982024[/C][C]0[/C][/ROW]
[ROW][C]0.41[/C][C]21.6203679518285[/C][C]0.538667594485544[/C][/ROW]
[ROW][C]0.51[/C][C]23.683823351018[/C][C]1.16373117950554[/C][/ROW]
[ROW][C]0.61[/C][C]24.9999999984928[/C][C]0[/C][/ROW]
[ROW][C]0.71[/C][C]25[/C][C]0[/C][/ROW]
[ROW][C]0.81[/C][C]25.0221049002756[/C][C]0.0598605898389705[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=206727&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=206727&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.019.481460198747260.531630278100433
0.1115.00253827858880.00915781008035258
0.2119.98981180703890.027716502034098
0.3120.00000009820240
0.4121.62036795182850.538667594485544
0.5123.6838233510181.16373117950554
0.6124.99999999849280
0.71250
0.8125.02210490027560.0598605898389705


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 0.01 ; par2 = 0.9 ; par3 = 0.1 ;
Parameters (R input):
par1 = 0.01 ; par2 = 0.9 ; par3 = 0.1 ;
R code (references can be found in the software module):
par3 <- '0.01'
par2 <- '0.99'
par1 <- '0.01'
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')