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 13:53:08 -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/t1236628432l87bilxffpxmzrm.htm/, Retrieved Tue, 07 May 2024 22:38:42 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=38826, Retrieved Tue, 07 May 2024 22:38:42 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact106
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Thomas van den Bo...] [2009-03-09 19:53:08] [50e97696ebad247f45d73cd9926afb25] [Current]
Feedback Forum

Post a new message
Dataseries X:
36.80
35.40
33.00
28.73
26.70
26.46
24.60
28.00
31.60
33.50
34.50
35.00
34.76
33.50
32.74
34.40
31.93
29.24
25.75
26.03
26.08
23.80
20.61
19.70
18.18
19.60
20.60
20.03
23.00
23.60
22.56
22.55
23.75
24.92
24.50
30.58
28.07
27.70
27.00
25.23
26.86
25.60
24.55
23.96
23.50
23.64
21.55
21.05
21.89
21.98
21.45
22.15
22.58
23.80
23.30
22.38
23.00
21.96
22.40
20.80
20.40
16.00
12.78
9.75
7.50
11.24
12.24
12.75
12.52
14.49
14.21
14.32
22.15
22.58
23.80
23.30
22.38
23.00
21.96
22.40
20.80
20.40
16.00
12.78
9.75
7.50
11.24
12.24
12.75
12.52
14.49
14.21
14.32




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 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 & 0 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=38826&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]0 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=38826&T=0

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







Harrell-Davis Quantiles
quantilesvaluestandard error
0.112.47405019129320.55325483083052
0.215.36739753474361.85308515591103
0.320.03740119023261.27102098901072
0.421.79638391697490.526780433999596
0.522.64159745541500.383458620624374
0.623.56064859257930.478671834849329
0.724.98171173695150.861980820073833
0.827.4292590756121.23781833973568
0.932.39165434294891.4842859586931

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 12.4740501912932 & 0.55325483083052 \tabularnewline
0.2 & 15.3673975347436 & 1.85308515591103 \tabularnewline
0.3 & 20.0374011902326 & 1.27102098901072 \tabularnewline
0.4 & 21.7963839169749 & 0.526780433999596 \tabularnewline
0.5 & 22.6415974554150 & 0.383458620624374 \tabularnewline
0.6 & 23.5606485925793 & 0.478671834849329 \tabularnewline
0.7 & 24.9817117369515 & 0.861980820073833 \tabularnewline
0.8 & 27.429259075612 & 1.23781833973568 \tabularnewline
0.9 & 32.3916543429489 & 1.4842859586931 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=38826&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]12.4740501912932[/C][C]0.55325483083052[/C][/ROW]
[ROW][C]0.2[/C][C]15.3673975347436[/C][C]1.85308515591103[/C][/ROW]
[ROW][C]0.3[/C][C]20.0374011902326[/C][C]1.27102098901072[/C][/ROW]
[ROW][C]0.4[/C][C]21.7963839169749[/C][C]0.526780433999596[/C][/ROW]
[ROW][C]0.5[/C][C]22.6415974554150[/C][C]0.383458620624374[/C][/ROW]
[ROW][C]0.6[/C][C]23.5606485925793[/C][C]0.478671834849329[/C][/ROW]
[ROW][C]0.7[/C][C]24.9817117369515[/C][C]0.861980820073833[/C][/ROW]
[ROW][C]0.8[/C][C]27.429259075612[/C][C]1.23781833973568[/C][/ROW]
[ROW][C]0.9[/C][C]32.3916543429489[/C][C]1.4842859586931[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=38826&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=38826&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.112.47405019129320.55325483083052
0.215.36739753474361.85308515591103
0.320.03740119023261.27102098901072
0.421.79638391697490.526780433999596
0.522.64159745541500.383458620624374
0.623.56064859257930.478671834849329
0.724.98171173695150.861980820073833
0.827.4292590756121.23781833973568
0.932.39165434294891.4842859586931



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