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 computationMon, 09 Mar 2009 14:02:00 -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/09/t1236628939ky4crn3hjw0pntw.htm/, Retrieved Sun, 28 Apr 2024 19:44:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=38829, Retrieved Sun, 28 Apr 2024 19:44:16 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact149
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2009-03-09 20:02:00] [738a25b0d97c8f3fa6714f905e8e3fd3] [Current]
Feedback Forum

Post a new message
Dataseries X:
27,00
26,88
27,38
26,82
27,00
26,15
25,85
26,17
25,66
25,72
25,42
25,10
26,20
26,39
26,27
26,63
21,10
20,30
20,61
21,05
20,45
20,91
21,22
20,85
21,90
22,71
22,40
22,81
23,96
23,37
23,55
23,01
22,63
22,63
22,00
22,15
22,00
22,00
21,84
22,10
22,37
21,83
21,77
21,89
20,76
20,21
20,19
20,01
19,16
18,50
17,41
18,14
18,60
18,32
18,40
18,16
17,29
16,65
16,36
16,32
17,37
17,30
18,10
19,00
18,38
18,41
18,10
17,87
18,70
18,81
18,88
19,44
18,60
18,80
18,62
18,24
17,84
17,85
17,67
17,99
18,15
18,39
18,07
18,39




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=38829&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=38829&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=38829&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'Gwilym Jenkins' @ 72.249.127.135







Harrell-Davis Quantiles
quantilesvaluestandard error
0.117.6944311506250.236336567726454
0.218.18525088592390.129013037608781
0.318.52643494631790.197841755099074
0.419.30140680378500.623387076734050
0.520.61326585579850.611524379631726
0.621.66313904137290.426817978263806
0.722.43023220640120.423028328627039
0.824.23464456761241.11666117375964
0.926.23023696973610.331729661375538

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 17.694431150625 & 0.236336567726454 \tabularnewline
0.2 & 18.1852508859239 & 0.129013037608781 \tabularnewline
0.3 & 18.5264349463179 & 0.197841755099074 \tabularnewline
0.4 & 19.3014068037850 & 0.623387076734050 \tabularnewline
0.5 & 20.6132658557985 & 0.611524379631726 \tabularnewline
0.6 & 21.6631390413729 & 0.426817978263806 \tabularnewline
0.7 & 22.4302322064012 & 0.423028328627039 \tabularnewline
0.8 & 24.2346445676124 & 1.11666117375964 \tabularnewline
0.9 & 26.2302369697361 & 0.331729661375538 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=38829&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]17.694431150625[/C][C]0.236336567726454[/C][/ROW]
[ROW][C]0.2[/C][C]18.1852508859239[/C][C]0.129013037608781[/C][/ROW]
[ROW][C]0.3[/C][C]18.5264349463179[/C][C]0.197841755099074[/C][/ROW]
[ROW][C]0.4[/C][C]19.3014068037850[/C][C]0.623387076734050[/C][/ROW]
[ROW][C]0.5[/C][C]20.6132658557985[/C][C]0.611524379631726[/C][/ROW]
[ROW][C]0.6[/C][C]21.6631390413729[/C][C]0.426817978263806[/C][/ROW]
[ROW][C]0.7[/C][C]22.4302322064012[/C][C]0.423028328627039[/C][/ROW]
[ROW][C]0.8[/C][C]24.2346445676124[/C][C]1.11666117375964[/C][/ROW]
[ROW][C]0.9[/C][C]26.2302369697361[/C][C]0.331729661375538[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=38829&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=38829&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.117.6944311506250.236336567726454
0.218.18525088592390.129013037608781
0.318.52643494631790.197841755099074
0.419.30140680378500.623387076734050
0.520.61326585579850.611524379631726
0.621.66313904137290.426817978263806
0.722.43023220640120.423028328627039
0.824.23464456761241.11666117375964
0.926.23023696973610.331729661375538



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