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 computationSun, 22 Mar 2009 15:23:45 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Mar/22/t1237757081o772pmhqcch4d8r.htm/, Retrieved Wed, 24 Apr 2024 23:12:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=39156, Retrieved Wed, 24 Apr 2024 23:12:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact220
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Gabriels Wim Deci...] [2009-03-22 21:23:45] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
5
5
5
6
8
9
9
9
10
10
10
10
10
10
10
10
10
10
10
10
10
10
10
10
10
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
13
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
15
16
16
16
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
18
19
19
19
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
20
21
21
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
22
23
23
24
24
24
24
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
25
26
26
28
28
28
28
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
30
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
35
38
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
40
50
50
50
50
50
50
50
50
50
50
50
50
50
50
50
50
65




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 2 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=39156&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=39156&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39156&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 time2 seconds
R Server'George Udny Yule' @ 72.249.76.132







Harrell-Davis Quantiles
quantilesvaluestandard error
0.115.00007139374610.000330841003470901
0.219.93090219759930.146281932400588
0.320.00000000243060
0.421.21087483190720.766741332207266
0.522.99177021477981.07013689280987
0.624.99999994553770
0.7250
0.825.00180702541750.00600848841895861
0.930.00014886562620.000560591507156429

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 15.0000713937461 & 0.000330841003470901 \tabularnewline
0.2 & 19.9309021975993 & 0.146281932400588 \tabularnewline
0.3 & 20.0000000024306 & 0 \tabularnewline
0.4 & 21.2108748319072 & 0.766741332207266 \tabularnewline
0.5 & 22.9917702147798 & 1.07013689280987 \tabularnewline
0.6 & 24.9999999455377 & 0 \tabularnewline
0.7 & 25 & 0 \tabularnewline
0.8 & 25.0018070254175 & 0.00600848841895861 \tabularnewline
0.9 & 30.0001488656262 & 0.000560591507156429 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=39156&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]15.0000713937461[/C][C]0.000330841003470901[/C][/ROW]
[ROW][C]0.2[/C][C]19.9309021975993[/C][C]0.146281932400588[/C][/ROW]
[ROW][C]0.3[/C][C]20.0000000024306[/C][C]0[/C][/ROW]
[ROW][C]0.4[/C][C]21.2108748319072[/C][C]0.766741332207266[/C][/ROW]
[ROW][C]0.5[/C][C]22.9917702147798[/C][C]1.07013689280987[/C][/ROW]
[ROW][C]0.6[/C][C]24.9999999455377[/C][C]0[/C][/ROW]
[ROW][C]0.7[/C][C]25[/C][C]0[/C][/ROW]
[ROW][C]0.8[/C][C]25.0018070254175[/C][C]0.00600848841895861[/C][/ROW]
[ROW][C]0.9[/C][C]30.0001488656262[/C][C]0.000560591507156429[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=39156&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39156&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.115.00007139374610.000330841003470901
0.219.93090219759930.146281932400588
0.320.00000000243060
0.421.21087483190720.766741332207266
0.522.99177021477981.07013689280987
0.624.99999994553770
0.7250
0.825.00180702541750.00600848841895861
0.930.00014886562620.000560591507156429



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')