FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_harrell_davis.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationFri, 04 Dec 2009 11:29:06 -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/2009/Dec/04/t12599516348krd1sl6tte37tg.htm/, Retrieved Sun, 28 Apr 2024 08:47:06 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=64010, Retrieved Sun, 28 Apr 2024 08:47:06 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact111

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-12 13:32:37] [76963dc1903f0f612b6153510a3818cf]
- R  D  [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-17 12:14:40] [76963dc1903f0f612b6153510a3818cf]
-         [Univariate Explorative Data Analysis] [Run Sequence Plot...] [2008-12-22 18:19:51] [1ce0d16c8f4225c977b42c8fa93bc163]
- RMP       [Univariate Data Series] [Identifying Integ...] [2009-11-22 12:08:06] [b98453cac15ba1066b407e146608df68]
- RMPD        [Percentiles] [] [2009-12-03 10:35:08] [3b0db66ac8145b1be856a517e2900332]
-    D          [Percentiles] [] [2009-12-03 17:22:23] [ee35698a38947a6c6c039b1e3deafc05]
- RMPD              [Harrell-Davis Quantiles] [] [2009-12-04 18:29:06] [c5f9f441970441f2f938cd843072158d] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
5.22817665838978e-05
-0.000463820140926281
-4.40505915974873e-05
3.66701376492632e-05
8.27894181489517e-05
0.000120006006809757
0.000171017242898982
-0.000279502333089345
8.75237612852431e-06
-3.91501689793012e-05
-0.000194144893208234
-0.000313101733963577
-5.08960481814131e-05
-0.000233376280487795
-0.000191085353174356
5.24200071684507e-05
-8.82943807301144e-05
3.61935156370373e-05
-8.69486245395421e-07
-0.000275113058485018
-0.000320465698642699
-0.000253480138722013
7.7227498765784e-05
0.000267886776359312
0.000106485389860988
-0.000165627713471103
-7.2254564411878e-05
-0.000135378139756055
-0.000127539859690442
-0.000213768154000569
-0.00018342039892232
-0.000138894755312727
6.85874943347128e-05
-0.000156651547253009
-1.07738256694651e-05
0.000243449321514413
-4.74576674314878e-05
-8.55531074017408e-05
0.000172352041347206
0.000139023530626794
6.7862815596465e-05
0.000195657674159593
-0.000347965684078771
0.000294749351197024
-8.77554455976112e-05
-0.000337265407768624
0.000379267264700164
-0.000141585653849405

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64010&T=0

As an alternative you can also use a QR Code:  
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.01-0.0004457039511345099.85690005025227e-05
0.06-0.0003379852987238472.84256673443792e-05
0.11-0.000291411368769073.36629513032892e-05
0.16-0.0002503624624524593.81145151488093e-05
0.21-0.0002127630554667293.59102417614131e-05
0.26-0.0001813200098318763.20549700954041e-05
0.31-0.0001542085545017223.12265821859620e-05
0.36-0.0001281805256327783.31139125901475e-05
0.41-0.0001016333706702963.48040918012818e-05
0.46-7.49527711125178e-053.4659436198585e-05
0.51-4.8295889238851e-053.5223220069855e-05
0.56-2.05585359587772e-053.73460625248096e-05
0.618.1447052094873e-063.70131163285573e-05
0.663.53108551814767e-053.21159645084506e-05
0.715.96964880508208e-052.75342585993309e-05
0.768.50211110670112e-052.9905958215845e-05
0.810.0001174986315674253.71331099220778e-05
0.860.0001602640484587084.18296650065675e-05
0.910.0002167044909387844.63333705239767e-05
0.960.0002984926902514254.82752290806702e-05

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & -0.000445703951134509 & 9.85690005025227e-05 \tabularnewline
0.06 & -0.000337985298723847 & 2.84256673443792e-05 \tabularnewline
0.11 & -0.00029141136876907 & 3.36629513032892e-05 \tabularnewline
0.16 & -0.000250362462452459 & 3.81145151488093e-05 \tabularnewline
0.21 & -0.000212763055466729 & 3.59102417614131e-05 \tabularnewline
0.26 & -0.000181320009831876 & 3.20549700954041e-05 \tabularnewline
0.31 & -0.000154208554501722 & 3.12265821859620e-05 \tabularnewline
0.36 & -0.000128180525632778 & 3.31139125901475e-05 \tabularnewline
0.41 & -0.000101633370670296 & 3.48040918012818e-05 \tabularnewline
0.46 & -7.49527711125178e-05 & 3.4659436198585e-05 \tabularnewline
0.51 & -4.8295889238851e-05 & 3.5223220069855e-05 \tabularnewline
0.56 & -2.05585359587772e-05 & 3.73460625248096e-05 \tabularnewline
0.61 & 8.1447052094873e-06 & 3.70131163285573e-05 \tabularnewline
0.66 & 3.53108551814767e-05 & 3.21159645084506e-05 \tabularnewline
0.71 & 5.96964880508208e-05 & 2.75342585993309e-05 \tabularnewline
0.76 & 8.50211110670112e-05 & 2.9905958215845e-05 \tabularnewline
0.81 & 0.000117498631567425 & 3.71331099220778e-05 \tabularnewline
0.86 & 0.000160264048458708 & 4.18296650065675e-05 \tabularnewline
0.91 & 0.000216704490938784 & 4.63333705239767e-05 \tabularnewline
0.96 & 0.000298492690251425 & 4.82752290806702e-05 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=64010&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]-0.000445703951134509[/C][C]9.85690005025227e-05[/C][/ROW]
[ROW][C]0.06[/C][C]-0.000337985298723847[/C][C]2.84256673443792e-05[/C][/ROW]
[ROW][C]0.11[/C][C]-0.00029141136876907[/C][C]3.36629513032892e-05[/C][/ROW]
[ROW][C]0.16[/C][C]-0.000250362462452459[/C][C]3.81145151488093e-05[/C][/ROW]
[ROW][C]0.21[/C][C]-0.000212763055466729[/C][C]3.59102417614131e-05[/C][/ROW]
[ROW][C]0.26[/C][C]-0.000181320009831876[/C][C]3.20549700954041e-05[/C][/ROW]
[ROW][C]0.31[/C][C]-0.000154208554501722[/C][C]3.12265821859620e-05[/C][/ROW]
[ROW][C]0.36[/C][C]-0.000128180525632778[/C][C]3.31139125901475e-05[/C][/ROW]
[ROW][C]0.41[/C][C]-0.000101633370670296[/C][C]3.48040918012818e-05[/C][/ROW]
[ROW][C]0.46[/C][C]-7.49527711125178e-05[/C][C]3.4659436198585e-05[/C][/ROW]
[ROW][C]0.51[/C][C]-4.8295889238851e-05[/C][C]3.5223220069855e-05[/C][/ROW]
[ROW][C]0.56[/C][C]-2.05585359587772e-05[/C][C]3.73460625248096e-05[/C][/ROW]
[ROW][C]0.61[/C][C]8.1447052094873e-06[/C][C]3.70131163285573e-05[/C][/ROW]
[ROW][C]0.66[/C][C]3.53108551814767e-05[/C][C]3.21159645084506e-05[/C][/ROW]
[ROW][C]0.71[/C][C]5.96964880508208e-05[/C][C]2.75342585993309e-05[/C][/ROW]
[ROW][C]0.76[/C][C]8.50211110670112e-05[/C][C]2.9905958215845e-05[/C][/ROW]
[ROW][C]0.81[/C][C]0.000117498631567425[/C][C]3.71331099220778e-05[/C][/ROW]
[ROW][C]0.86[/C][C]0.000160264048458708[/C][C]4.18296650065675e-05[/C][/ROW]
[ROW][C]0.91[/C][C]0.000216704490938784[/C][C]4.63333705239767e-05[/C][/ROW]
[ROW][C]0.96[/C][C]0.000298492690251425[/C][C]4.82752290806702e-05[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=64010&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=64010&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Harrell-Davis Quantiles
quantilesvaluestandard error
0.01-0.0004457039511345099.85690005025227e-05
0.06-0.0003379852987238472.84256673443792e-05
0.11-0.000291411368769073.36629513032892e-05
0.16-0.0002503624624524593.81145151488093e-05
0.21-0.0002127630554667293.59102417614131e-05
0.26-0.0001813200098318763.20549700954041e-05
0.31-0.0001542085545017223.12265821859620e-05
0.36-0.0001281805256327783.31139125901475e-05
0.41-0.0001016333706702963.48040918012818e-05
0.46-7.49527711125178e-053.4659436198585e-05
0.51-4.8295889238851e-053.5223220069855e-05
0.56-2.05585359587772e-053.73460625248096e-05
0.618.1447052094873e-063.70131163285573e-05
0.663.53108551814767e-053.21159645084506e-05
0.715.96964880508208e-052.75342585993309e-05
0.768.50211110670112e-052.9905958215845e-05
0.810.0001174986315674253.71331099220778e-05
0.860.0001602640484587084.18296650065675e-05
0.910.0002167044909387844.63333705239767e-05
0.960.0002984926902514254.82752290806702e-05


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 36 ; par2 = -1.8 ; par3 = 2 ; par4 = 1 ; par5 = 12 ; par6 = MA ; par7 = 0.95 ;
Parameters (R input):
par1 = 0.01 ; par2 = 0.99 ; par3 = 0.05 ;
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')