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, 26 May 2012 05:30:31 -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/26/t133802472672twi6y3pynu522.htm/, Retrieved Thu, 02 May 2024 21:38:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=167563, Retrieved Thu, 02 May 2024 21:38:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG201142
Estimated Impact113
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Decielen inschrij...] [2012-05-26 09:30:31] [a2709f7ab86902c8a6d752383a828f6c] [Current]
Feedback Forum

Post a new message
Dataseries X:
67,66
68
68,02
68,11
68,41
68,4
68,4
68,55
68,54
68,99
68,97
68,98
68,98
68,94
69,21
69,21
69,67
69,66
69,66
69,66
69,77
70,32
70,34
70,38
70,38
70,29
70,42
70,29
70,59
70,64
70,64
70,68
70,78
70,9
71,04
71,15
71,15
71,15
71,07
71,17
71,24
71,23
71,23
71,23
71,24
71,28
71,52
71,52
71,52
71,6
71,61
71,78
71,66
71,86
71,86
71,82
71,8
72,22
72,51
72,56




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=167563&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.168.42490038549650.203297874966722
0.269.00518226817130.262937421641831
0.369.63657161924930.375763069064088
0.470.24284833632940.277806969824116
0.570.65026221857080.23316875960649
0.671.02291823935170.182823070606844
0.771.24457096517080.101780556236766
0.871.50697939285570.147488942400825
0.971.81863610949430.112697170147362

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 68.4249003854965 & 0.203297874966722 \tabularnewline
0.2 & 69.0051822681713 & 0.262937421641831 \tabularnewline
0.3 & 69.6365716192493 & 0.375763069064088 \tabularnewline
0.4 & 70.2428483363294 & 0.277806969824116 \tabularnewline
0.5 & 70.6502622185708 & 0.23316875960649 \tabularnewline
0.6 & 71.0229182393517 & 0.182823070606844 \tabularnewline
0.7 & 71.2445709651708 & 0.101780556236766 \tabularnewline
0.8 & 71.5069793928557 & 0.147488942400825 \tabularnewline
0.9 & 71.8186361094943 & 0.112697170147362 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=167563&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]68.4249003854965[/C][C]0.203297874966722[/C][/ROW]
[ROW][C]0.2[/C][C]69.0051822681713[/C][C]0.262937421641831[/C][/ROW]
[ROW][C]0.3[/C][C]69.6365716192493[/C][C]0.375763069064088[/C][/ROW]
[ROW][C]0.4[/C][C]70.2428483363294[/C][C]0.277806969824116[/C][/ROW]
[ROW][C]0.5[/C][C]70.6502622185708[/C][C]0.23316875960649[/C][/ROW]
[ROW][C]0.6[/C][C]71.0229182393517[/C][C]0.182823070606844[/C][/ROW]
[ROW][C]0.7[/C][C]71.2445709651708[/C][C]0.101780556236766[/C][/ROW]
[ROW][C]0.8[/C][C]71.5069793928557[/C][C]0.147488942400825[/C][/ROW]
[ROW][C]0.9[/C][C]71.8186361094943[/C][C]0.112697170147362[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=167563&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=167563&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.168.42490038549650.203297874966722
0.269.00518226817130.262937421641831
0.369.63657161924930.375763069064088
0.470.24284833632940.277806969824116
0.570.65026221857080.23316875960649
0.671.02291823935170.182823070606844
0.771.24457096517080.101780556236766
0.871.50697939285570.147488942400825
0.971.81863610949430.112697170147362



Parameters (Session):
par1 = 0.01 ; par2 = 0.99 ; par3 = 0.01 ;
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')