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 computationFri, 27 Mar 2009 02:13:42 -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/27/t1238141684i5i26p3xo94axen.htm/, Retrieved Tue, 30 Apr 2024 12:26:56 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=39212, Retrieved Tue, 30 Apr 2024 12:26:56 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact188
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [deciel us fedfunds] [2009-03-27 08:13:42] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
1.73
1.75
1.75
1.75
1.73
1.74
1.75
1.75
1.34
1.24
1.24
1.26
1.25
1.26
1.26
1.22
1.01
1.03
1.01
1.01
1
0.98
1
1.01
1
1
1
1.03
1.26
1.43
1.61
1.76
1.93
2.16
2.28
2.5
2.63
2.79
3
3.04
3.26
3.5
3.62
3.78
4
4.16
4.29
4.49
4.59
4.79
4.94
4.99
5.24
5.25
5.25
5.25
5.25
5.24
5.25
5.26
5.26
5.25
5.25
5.25
5.26
5.02
4.94
4.76
4.49
4.24
3.94
2.98
2.61
2.28
1.98
2
2.01
2
1.81
0.97
0.39
0.16
0.15
0.22




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39212&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.10.9797725740432980.0396157863476157
0.21.123702666416400.104296007116207
0.31.421349880659190.185431825582093
0.41.761316923841420.126126386128966
0.52.136140446616430.292752836299755
0.62.939614708313350.500665954318674
0.73.999437078756290.467275337326379
0.84.826467301505630.273215878849616
0.95.232737309673020.0391279726884938

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 0.979772574043298 & 0.0396157863476157 \tabularnewline
0.2 & 1.12370266641640 & 0.104296007116207 \tabularnewline
0.3 & 1.42134988065919 & 0.185431825582093 \tabularnewline
0.4 & 1.76131692384142 & 0.126126386128966 \tabularnewline
0.5 & 2.13614044661643 & 0.292752836299755 \tabularnewline
0.6 & 2.93961470831335 & 0.500665954318674 \tabularnewline
0.7 & 3.99943707875629 & 0.467275337326379 \tabularnewline
0.8 & 4.82646730150563 & 0.273215878849616 \tabularnewline
0.9 & 5.23273730967302 & 0.0391279726884938 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=39212&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]0.979772574043298[/C][C]0.0396157863476157[/C][/ROW]
[ROW][C]0.2[/C][C]1.12370266641640[/C][C]0.104296007116207[/C][/ROW]
[ROW][C]0.3[/C][C]1.42134988065919[/C][C]0.185431825582093[/C][/ROW]
[ROW][C]0.4[/C][C]1.76131692384142[/C][C]0.126126386128966[/C][/ROW]
[ROW][C]0.5[/C][C]2.13614044661643[/C][C]0.292752836299755[/C][/ROW]
[ROW][C]0.6[/C][C]2.93961470831335[/C][C]0.500665954318674[/C][/ROW]
[ROW][C]0.7[/C][C]3.99943707875629[/C][C]0.467275337326379[/C][/ROW]
[ROW][C]0.8[/C][C]4.82646730150563[/C][C]0.273215878849616[/C][/ROW]
[ROW][C]0.9[/C][C]5.23273730967302[/C][C]0.0391279726884938[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=39212&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39212&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.10.9797725740432980.0396157863476157
0.21.123702666416400.104296007116207
0.31.421349880659190.185431825582093
0.41.761316923841420.126126386128966
0.52.136140446616430.292752836299755
0.62.939614708313350.500665954318674
0.73.999437078756290.467275337326379
0.84.826467301505630.273215878849616
0.95.232737309673020.0391279726884938



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