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 computationSun, 07 Mar 2010 13:41:29 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Mar/07/t1267994577o496sa3zfswhzos.htm/, Retrieved Thu, 28 Mar 2024 16:29:05 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=74099, Retrieved Thu, 28 Mar 2024 16:29:05 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDGP1W42
Estimated Impact147
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [H-D percentielen-...] [2010-03-07 20:41:29] [17264c076bcd9ffe53c701d21c08cfe9] [Current]
Feedback Forum

Post a new message
Dataseries X:
2.06
2.06
2.05
2.04
2.04
2.03
2.03
2.03
2.02
2.02
2.02
2.02
2.01
2.01
2.01
2
2
1.99
1.99
1.99
1.98
1.97
1.97
1.96
1.96
1.96
1.96
1.95
1.94
1.94
1.93
1.93
1.93
1.92
1.92
1.9
1.9
1.9
1.9
1.89
1.88
1.88
1.87
1.86
1.86
1.85
1.83
1.82
1.8
1.8
1.79
1.78
1.78
1.78
1.77
1.77
1.76
1.75
1.75
1.75
1.75
1.75
1.74
1.74
1.74
1.74
1.73
1.73
1.73
1.73
1.73
1.72
1.71
1.71
1.71
1.7
1.7
1.7
1.7
1.7
1.69
1.68
1.67
1.65
1.61




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

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







Harrell-Davis Quantiles
quantilesvaluestandard error
0.011.624776950906150.0328569099019673
0.021.644652997371490.0252948429594947
0.031.661530272367150.0199950111139237
0.041.673825858654640.0163750426741892
0.051.682603495843880.0134658704119771
0.061.688976939855690.0108576043420088
0.071.693673333727260.00863981820099488
0.081.697213223620290.00713592575119727
0.091.700037831116370.00659832185114523

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 1.62477695090615 & 0.0328569099019673 \tabularnewline
0.02 & 1.64465299737149 & 0.0252948429594947 \tabularnewline
0.03 & 1.66153027236715 & 0.0199950111139237 \tabularnewline
0.04 & 1.67382585865464 & 0.0163750426741892 \tabularnewline
0.05 & 1.68260349584388 & 0.0134658704119771 \tabularnewline
0.06 & 1.68897693985569 & 0.0108576043420088 \tabularnewline
0.07 & 1.69367333372726 & 0.00863981820099488 \tabularnewline
0.08 & 1.69721322362029 & 0.00713592575119727 \tabularnewline
0.09 & 1.70003783111637 & 0.00659832185114523 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=74099&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.01[/C][C]1.62477695090615[/C][C]0.0328569099019673[/C][/ROW]
[ROW][C]0.02[/C][C]1.64465299737149[/C][C]0.0252948429594947[/C][/ROW]
[ROW][C]0.03[/C][C]1.66153027236715[/C][C]0.0199950111139237[/C][/ROW]
[ROW][C]0.04[/C][C]1.67382585865464[/C][C]0.0163750426741892[/C][/ROW]
[ROW][C]0.05[/C][C]1.68260349584388[/C][C]0.0134658704119771[/C][/ROW]
[ROW][C]0.06[/C][C]1.68897693985569[/C][C]0.0108576043420088[/C][/ROW]
[ROW][C]0.07[/C][C]1.69367333372726[/C][C]0.00863981820099488[/C][/ROW]
[ROW][C]0.08[/C][C]1.69721322362029[/C][C]0.00713592575119727[/C][/ROW]
[ROW][C]0.09[/C][C]1.70003783111637[/C][C]0.00659832185114523[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=74099&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=74099&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.011.624776950906150.0328569099019673
0.021.644652997371490.0252948429594947
0.031.661530272367150.0199950111139237
0.041.673825858654640.0163750426741892
0.051.682603495843880.0134658704119771
0.061.688976939855690.0108576043420088
0.071.693673333727260.00863981820099488
0.081.697213223620290.00713592575119727
0.091.700037831116370.00659832185114523



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