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, 16 Feb 2014 11:07:15 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Feb/16/t1392566979ytnymym6bnk9h7k.htm/, Retrieved Thu, 16 May 2024 02:23:24 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=233759, Retrieved Thu, 16 May 2024 02:23:24 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact105
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [OP4-oef2- Benjami...] [2014-02-16 16:07:15] [bde58939c1a4c9e8add61ba7c57eaa76] [Current]
Feedback Forum

Post a new message
Dataseries X:
530.3
527.76
521.41
1601.93
1577.49
1551.43
1551.43
1516.88
1485.95
1438.22
1385.06
1329.49
1329.49
1276.16
1242.34
1181.59
1160.21
1135.18
1135.18
1084.96
1077.35
1061.13
1029.98
1013.08
1013.08
996.04
975.02
951.89
944.4
932.47
932.47
920.44
900.18
886.9
869.74
859.03
859.03
844.99
834.82
825.62
816.92
813.21
813.21
811.03
804.16
788.62
778.76
765.91
765.91
753.85
742.22
732.11
729.94
731.22
731.22
729.11
726.94
720.52
709.36
703.21




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net

\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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=233759&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=233759&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=233759&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'Herman Ole Andreas Wold' @ wold.wessa.net







Harrell-Davis Quantiles
quantilesvaluestandard error
0.05606.53167528531981.4370811889513
0.1706.88118916446626.4201707472485
0.15731.1159474381529.87423400916843
0.2745.72254895014518.1052736847222
0.25767.13252429994525.0275370005656
0.3791.13328195944625.3115801952841
0.35813.53442237558824.2972875102148
0.4836.33362751798328.5713039993765
0.45863.59999928859135.8579768540672
0.5896.14022703816641.3332837545999
0.55932.43461927169944.5903620044717
0.6972.4961776675949.1234011520149
0.651018.4316208834855.5991491274296
0.71073.007846471564.0489467091888
0.751140.7777144289476.4467273880547
0.81228.0516139622290.9419440848481
0.851335.1004438892195.3130047555547
0.91451.705501626782.4004495815868
0.951544.7660936805438.240685259396

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.05 & 606.531675285319 & 81.4370811889513 \tabularnewline
0.1 & 706.881189164466 & 26.4201707472485 \tabularnewline
0.15 & 731.115947438152 & 9.87423400916843 \tabularnewline
0.2 & 745.722548950145 & 18.1052736847222 \tabularnewline
0.25 & 767.132524299945 & 25.0275370005656 \tabularnewline
0.3 & 791.133281959446 & 25.3115801952841 \tabularnewline
0.35 & 813.534422375588 & 24.2972875102148 \tabularnewline
0.4 & 836.333627517983 & 28.5713039993765 \tabularnewline
0.45 & 863.599999288591 & 35.8579768540672 \tabularnewline
0.5 & 896.140227038166 & 41.3332837545999 \tabularnewline
0.55 & 932.434619271699 & 44.5903620044717 \tabularnewline
0.6 & 972.49617766759 & 49.1234011520149 \tabularnewline
0.65 & 1018.43162088348 & 55.5991491274296 \tabularnewline
0.7 & 1073.0078464715 & 64.0489467091888 \tabularnewline
0.75 & 1140.77771442894 & 76.4467273880547 \tabularnewline
0.8 & 1228.05161396222 & 90.9419440848481 \tabularnewline
0.85 & 1335.10044388921 & 95.3130047555547 \tabularnewline
0.9 & 1451.7055016267 & 82.4004495815868 \tabularnewline
0.95 & 1544.76609368054 & 38.240685259396 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=233759&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.05[/C][C]606.531675285319[/C][C]81.4370811889513[/C][/ROW]
[ROW][C]0.1[/C][C]706.881189164466[/C][C]26.4201707472485[/C][/ROW]
[ROW][C]0.15[/C][C]731.115947438152[/C][C]9.87423400916843[/C][/ROW]
[ROW][C]0.2[/C][C]745.722548950145[/C][C]18.1052736847222[/C][/ROW]
[ROW][C]0.25[/C][C]767.132524299945[/C][C]25.0275370005656[/C][/ROW]
[ROW][C]0.3[/C][C]791.133281959446[/C][C]25.3115801952841[/C][/ROW]
[ROW][C]0.35[/C][C]813.534422375588[/C][C]24.2972875102148[/C][/ROW]
[ROW][C]0.4[/C][C]836.333627517983[/C][C]28.5713039993765[/C][/ROW]
[ROW][C]0.45[/C][C]863.599999288591[/C][C]35.8579768540672[/C][/ROW]
[ROW][C]0.5[/C][C]896.140227038166[/C][C]41.3332837545999[/C][/ROW]
[ROW][C]0.55[/C][C]932.434619271699[/C][C]44.5903620044717[/C][/ROW]
[ROW][C]0.6[/C][C]972.49617766759[/C][C]49.1234011520149[/C][/ROW]
[ROW][C]0.65[/C][C]1018.43162088348[/C][C]55.5991491274296[/C][/ROW]
[ROW][C]0.7[/C][C]1073.0078464715[/C][C]64.0489467091888[/C][/ROW]
[ROW][C]0.75[/C][C]1140.77771442894[/C][C]76.4467273880547[/C][/ROW]
[ROW][C]0.8[/C][C]1228.05161396222[/C][C]90.9419440848481[/C][/ROW]
[ROW][C]0.85[/C][C]1335.10044388921[/C][C]95.3130047555547[/C][/ROW]
[ROW][C]0.9[/C][C]1451.7055016267[/C][C]82.4004495815868[/C][/ROW]
[ROW][C]0.95[/C][C]1544.76609368054[/C][C]38.240685259396[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=233759&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=233759&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.05606.53167528531981.4370811889513
0.1706.88118916446626.4201707472485
0.15731.1159474381529.87423400916843
0.2745.72254895014518.1052736847222
0.25767.13252429994525.0275370005656
0.3791.13328195944625.3115801952841
0.35813.53442237558824.2972875102148
0.4836.33362751798328.5713039993765
0.45863.59999928859135.8579768540672
0.5896.14022703816641.3332837545999
0.55932.43461927169944.5903620044717
0.6972.4961776675949.1234011520149
0.651018.4316208834855.5991491274296
0.71073.007846471564.0489467091888
0.751140.7777144289476.4467273880547
0.81228.0516139622290.9419440848481
0.851335.1004438892195.3130047555547
0.91451.705501626782.4004495815868
0.951544.7660936805438.240685259396



Parameters (Session):
par1 = 0.05 ; par2 = 0.95 ; par3 = 0.05 ;
Parameters (R input):
par1 = 0.05 ; par2 = 0.95 ; par3 = 0.05 ;
R code (references can be found in the software module):
par3 <- '0.1'
par2 <- '0.95'
par1 <- '0.05'
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')