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 computationWed, 12 Oct 2011 14:37:21 -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/2011/Oct/12/t131844509970ha25aipj4sz3v.htm/, Retrieved Thu, 09 May 2024 21:19:54 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=128494, Retrieved Thu, 09 May 2024 21:19:54 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG201142
Estimated Impact80
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2011-10-12 18:37:21] [5143da885265528a100ac524801fc397] [Current]
Feedback Forum

Post a new message
Dataseries X:
31.956
29.506
34.506
27.165
26.736
23.691
18.157
17.328
18.205
20.995
17.382
9.367
31.124
26.551
30.651
25.859
25.100
25.778
20.418
18.688
20.424
24.776
19.814
12.738
31.566
30.111
30.019
31.934
25.826
26.835
20.205
17.789
20.520
22.518
15.572
11.509
25.447
24.090
27.786
26.195
20.516
22.759
19.028
16.971
20.036
22.485
18.730
14.538
27.561
25.985
34.670
32.066
27.186
29.586
21.359
21.553
19.573
24.256
22.380
16.167




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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 & 1 seconds \tabularnewline
R Server & 'George Udny Yule' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=128494&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]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=128494&T=0

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







Harrell-Davis Quantiles
quantilesvaluestandard error
0.116.15426798684291.20326716217538
0.218.4450344108250.775665243775775
0.319.96507247632720.7081293448511
0.421.36608039159920.988113272021589
0.523.34401465533691.288360741083
0.625.28185043442960.926446419781788
0.726.65330783162450.744962398065101
0.828.69347322726441.26806015510051
0.931.2421975821250.79595820710452

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 16.1542679868429 & 1.20326716217538 \tabularnewline
0.2 & 18.445034410825 & 0.775665243775775 \tabularnewline
0.3 & 19.9650724763272 & 0.7081293448511 \tabularnewline
0.4 & 21.3660803915992 & 0.988113272021589 \tabularnewline
0.5 & 23.3440146553369 & 1.288360741083 \tabularnewline
0.6 & 25.2818504344296 & 0.926446419781788 \tabularnewline
0.7 & 26.6533078316245 & 0.744962398065101 \tabularnewline
0.8 & 28.6934732272644 & 1.26806015510051 \tabularnewline
0.9 & 31.242197582125 & 0.79595820710452 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=128494&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]16.1542679868429[/C][C]1.20326716217538[/C][/ROW]
[ROW][C]0.2[/C][C]18.445034410825[/C][C]0.775665243775775[/C][/ROW]
[ROW][C]0.3[/C][C]19.9650724763272[/C][C]0.7081293448511[/C][/ROW]
[ROW][C]0.4[/C][C]21.3660803915992[/C][C]0.988113272021589[/C][/ROW]
[ROW][C]0.5[/C][C]23.3440146553369[/C][C]1.288360741083[/C][/ROW]
[ROW][C]0.6[/C][C]25.2818504344296[/C][C]0.926446419781788[/C][/ROW]
[ROW][C]0.7[/C][C]26.6533078316245[/C][C]0.744962398065101[/C][/ROW]
[ROW][C]0.8[/C][C]28.6934732272644[/C][C]1.26806015510051[/C][/ROW]
[ROW][C]0.9[/C][C]31.242197582125[/C][C]0.79595820710452[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=128494&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=128494&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.116.15426798684291.20326716217538
0.218.4450344108250.775665243775775
0.319.96507247632720.7081293448511
0.421.36608039159920.988113272021589
0.523.34401465533691.288360741083
0.625.28185043442960.926446419781788
0.726.65330783162450.744962398065101
0.828.69347322726441.26806015510051
0.931.2421975821250.79595820710452



Parameters (Session):
par1 = 0.1 ; par2 = 1 ; 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')