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, 19 May 2012 05:39:10 -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/19/t1337420390gxeb74gv7m54cy5.htm/, Retrieved Sun, 05 May 2024 16:58:57 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=166716, Retrieved Sun, 05 May 2024 16:58:57 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG201142
Estimated Impact139
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2012-05-19 09:39:10] [21ca38b5c207d9fbe3078e977c625188] [Current]
Feedback Forum

Post a new message
Dataseries X:
31.1
31.8
32.5
34.4
35.5
35.5
36.6
37.1
37.9
38.1
39
41.5
41.8
41.9
44.6
46.1
46.4
47.2
47.7
49.2
49.3
49.3
49.5
50.1
51.9
52.6
53.2
53.5
53.7
53.7
53.9
54.1
54.8
55.4
55.9
56.8
58.4
59.3
60.3
60.5
60.8
61
61.1
61.3
61.4
61.5
63.9
63.9
64
64.1
64.5
64.5
65.9
66.8
68.7
69.2
69.6
70.2
70.6
70.7
70.7
71
72.1
73.7
77.4
79.7
91.6
93.6
94.3
97.3
101.7
103
103.1
104.6
107.2
107.7
108.3
108.8
113.1
113.8
113.8
116.5
116.9
117.6




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=166716&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=166716&T=0

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







Harrell-Davis Quantiles
quantilesvaluestandard error
0.137.97623905153592.1319887776153
0.246.0764821088492.89600703306002
0.351.69515231667442.12851266227144
0.456.13037268202772.50477637296161
0.560.91651695172822.04767038440182
0.665.16887846716412.50179002729203
0.771.28056160328423.81699410664255
0.889.471716924459710.4621034709853
0.9107.6290215694493.60519960497615

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 37.9762390515359 & 2.1319887776153 \tabularnewline
0.2 & 46.076482108849 & 2.89600703306002 \tabularnewline
0.3 & 51.6951523166744 & 2.12851266227144 \tabularnewline
0.4 & 56.1303726820277 & 2.50477637296161 \tabularnewline
0.5 & 60.9165169517282 & 2.04767038440182 \tabularnewline
0.6 & 65.1688784671641 & 2.50179002729203 \tabularnewline
0.7 & 71.2805616032842 & 3.81699410664255 \tabularnewline
0.8 & 89.4717169244597 & 10.4621034709853 \tabularnewline
0.9 & 107.629021569449 & 3.60519960497615 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=166716&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]37.9762390515359[/C][C]2.1319887776153[/C][/ROW]
[ROW][C]0.2[/C][C]46.076482108849[/C][C]2.89600703306002[/C][/ROW]
[ROW][C]0.3[/C][C]51.6951523166744[/C][C]2.12851266227144[/C][/ROW]
[ROW][C]0.4[/C][C]56.1303726820277[/C][C]2.50477637296161[/C][/ROW]
[ROW][C]0.5[/C][C]60.9165169517282[/C][C]2.04767038440182[/C][/ROW]
[ROW][C]0.6[/C][C]65.1688784671641[/C][C]2.50179002729203[/C][/ROW]
[ROW][C]0.7[/C][C]71.2805616032842[/C][C]3.81699410664255[/C][/ROW]
[ROW][C]0.8[/C][C]89.4717169244597[/C][C]10.4621034709853[/C][/ROW]
[ROW][C]0.9[/C][C]107.629021569449[/C][C]3.60519960497615[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=166716&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=166716&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.137.97623905153592.1319887776153
0.246.0764821088492.89600703306002
0.351.69515231667442.12851266227144
0.456.13037268202772.50477637296161
0.560.91651695172822.04767038440182
0.665.16887846716412.50179002729203
0.771.28056160328423.81699410664255
0.889.471716924459710.4621034709853
0.9107.6290215694493.60519960497615



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):
par3 <- '0.1'
par2 <- '0.9'
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')