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 computationSat, 29 Dec 2012 09:31:28 -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/2012/Dec/29/t1356791636itxvce9jxhi435l.htm/, Retrieved Thu, 02 May 2024 07:54:09 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=204873, Retrieved Thu, 02 May 2024 07:54:09 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact119
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2012-12-29 14:31:28] [4fb632b17dd9b46a28f1faaf764882b5] [Current]
Feedback Forum

Post a new message
Dataseries X:
33,7
34,59
35,1
35,87
37,15
37,61
37,97
38,94
39,18
39,49
39,86
40,02
40,2
40,85
41,45
41,7
41,92
41,97
42,31
42,61
42,82
43,07
43,51
43,57
43,86
44,49
45,99
48,22
49,46
50,39
50,4
50,59
51,32
51,86
52,47
52,73
52,73
53,59
54,11
54,8
55,72
56,06
56,66
57,05
57,31
57,89
58,32
58,72
59,02
59,54
61,49
62,26
63,49
64,36
65,93
66,82
68,85
71,27
72,27
73,4
73,58
74,84
75,74
77,81
78,74
79,06
79,48
81,19
85,11
86,64
88,48
89,2




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=204873&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' @ yule.wessa.net







Harrell-Davis Quantiles
quantilesvaluestandard error
0.138.34951069580821.11431252606113
0.241.09931029312411.01292923192178
0.343.59823911135811.72758800885515
0.448.47785244521172.97398569372822
0.552.99002602092662.07471609684058
0.656.79342886738562.04470189502389
0.761.76260089163733.35454908030552
0.870.33642568716254.03726470243248
0.978.68857912161772.59775018116913

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 38.3495106958082 & 1.11431252606113 \tabularnewline
0.2 & 41.0993102931241 & 1.01292923192178 \tabularnewline
0.3 & 43.5982391113581 & 1.72758800885515 \tabularnewline
0.4 & 48.4778524452117 & 2.97398569372822 \tabularnewline
0.5 & 52.9900260209266 & 2.07471609684058 \tabularnewline
0.6 & 56.7934288673856 & 2.04470189502389 \tabularnewline
0.7 & 61.7626008916373 & 3.35454908030552 \tabularnewline
0.8 & 70.3364256871625 & 4.03726470243248 \tabularnewline
0.9 & 78.6885791216177 & 2.59775018116913 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=204873&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]38.3495106958082[/C][C]1.11431252606113[/C][/ROW]
[ROW][C]0.2[/C][C]41.0993102931241[/C][C]1.01292923192178[/C][/ROW]
[ROW][C]0.3[/C][C]43.5982391113581[/C][C]1.72758800885515[/C][/ROW]
[ROW][C]0.4[/C][C]48.4778524452117[/C][C]2.97398569372822[/C][/ROW]
[ROW][C]0.5[/C][C]52.9900260209266[/C][C]2.07471609684058[/C][/ROW]
[ROW][C]0.6[/C][C]56.7934288673856[/C][C]2.04470189502389[/C][/ROW]
[ROW][C]0.7[/C][C]61.7626008916373[/C][C]3.35454908030552[/C][/ROW]
[ROW][C]0.8[/C][C]70.3364256871625[/C][C]4.03726470243248[/C][/ROW]
[ROW][C]0.9[/C][C]78.6885791216177[/C][C]2.59775018116913[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=204873&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=204873&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.138.34951069580821.11431252606113
0.241.09931029312411.01292923192178
0.343.59823911135811.72758800885515
0.448.47785244521172.97398569372822
0.552.99002602092662.07471609684058
0.656.79342886738562.04470189502389
0.761.76260089163733.35454908030552
0.870.33642568716254.03726470243248
0.978.68857912161772.59775018116913



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