Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_harrell_davies.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationMon, 22 Oct 2007 07:54:44 -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/2007/Oct/22/l29j0v30ep5f8661193064657.htm/, Retrieved Mon, 06 May 2024 10:24:04 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=1517, Retrieved Mon, 06 May 2024 10:24:04 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsTaak2b.5G19
Estimated Impact203
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [95% Confidence in...] [2007-10-22 14:54:44] [6b5c00822e2ce0f7cf73539c28d95782] [Current]
Feedback Forum

Post a new message
Dataseries X:
98,91
91,91
92,54
92,56
106,67
85,11
103,07
98,90
91,30
79,27
97,82
98,75
90,94
88,36
82,75
82,52
102,33
81,54
102,62
91,47
80,73
73,91
98,97
97,22
85,39
85,69
77,57
78,95
93,14
73,97
88,52
84,24
81,85
72,76
92,32
92,27
88,38
81,83
77,10
82,59
92,45
75,52
87,47
81,60
79,98
71,26
83,30
91,03
86,91
77,92
75,03
75,86
92,09
74,40
81,63
76,96
78,40
68,33
85,05
90,07
86,84




Summary of compuational 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 compuational 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=1517&T=0

[TABLE]
[ROW][C]Summary of compuational 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=1517&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1517&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 compuational 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.0169.03716033049572.57421141308358
0.0673.3107648178421.13721456461338
0.1175.13408445078981.12620725670397
0.1676.68135234463551.29239105902577
0.2178.17236711415031.39786055262194
0.2679.62933401441411.46919688433279
0.3180.94267702754231.27062376297728
0.3682.03425122425621.15481178659679
0.4183.1335184045831.391697870781
0.4684.4313105043571.65057646107131
0.5185.88361953746651.76811912874389
0.5687.42235780610971.84389496214624
0.6189.0018476635591.81749445139327
0.6690.44933046164871.48951425425951
0.7191.56248518028171.01176342940769
0.7692.54556945111.06814266940975
0.8194.2411019571392.10219325152500
0.8696.90750236017762.3513268396712
0.9199.5377963104591.71567013900682
0.96102.7950789375441.61293967388695

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 69.0371603304957 & 2.57421141308358 \tabularnewline
0.06 & 73.310764817842 & 1.13721456461338 \tabularnewline
0.11 & 75.1340844507898 & 1.12620725670397 \tabularnewline
0.16 & 76.6813523446355 & 1.29239105902577 \tabularnewline
0.21 & 78.1723671141503 & 1.39786055262194 \tabularnewline
0.26 & 79.6293340144141 & 1.46919688433279 \tabularnewline
0.31 & 80.9426770275423 & 1.27062376297728 \tabularnewline
0.36 & 82.0342512242562 & 1.15481178659679 \tabularnewline
0.41 & 83.133518404583 & 1.391697870781 \tabularnewline
0.46 & 84.431310504357 & 1.65057646107131 \tabularnewline
0.51 & 85.8836195374665 & 1.76811912874389 \tabularnewline
0.56 & 87.4223578061097 & 1.84389496214624 \tabularnewline
0.61 & 89.001847663559 & 1.81749445139327 \tabularnewline
0.66 & 90.4493304616487 & 1.48951425425951 \tabularnewline
0.71 & 91.5624851802817 & 1.01176342940769 \tabularnewline
0.76 & 92.5455694511 & 1.06814266940975 \tabularnewline
0.81 & 94.241101957139 & 2.10219325152500 \tabularnewline
0.86 & 96.9075023601776 & 2.3513268396712 \tabularnewline
0.91 & 99.537796310459 & 1.71567013900682 \tabularnewline
0.96 & 102.795078937544 & 1.61293967388695 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=1517&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]69.0371603304957[/C][C]2.57421141308358[/C][/ROW]
[ROW][C]0.06[/C][C]73.310764817842[/C][C]1.13721456461338[/C][/ROW]
[ROW][C]0.11[/C][C]75.1340844507898[/C][C]1.12620725670397[/C][/ROW]
[ROW][C]0.16[/C][C]76.6813523446355[/C][C]1.29239105902577[/C][/ROW]
[ROW][C]0.21[/C][C]78.1723671141503[/C][C]1.39786055262194[/C][/ROW]
[ROW][C]0.26[/C][C]79.6293340144141[/C][C]1.46919688433279[/C][/ROW]
[ROW][C]0.31[/C][C]80.9426770275423[/C][C]1.27062376297728[/C][/ROW]
[ROW][C]0.36[/C][C]82.0342512242562[/C][C]1.15481178659679[/C][/ROW]
[ROW][C]0.41[/C][C]83.133518404583[/C][C]1.391697870781[/C][/ROW]
[ROW][C]0.46[/C][C]84.431310504357[/C][C]1.65057646107131[/C][/ROW]
[ROW][C]0.51[/C][C]85.8836195374665[/C][C]1.76811912874389[/C][/ROW]
[ROW][C]0.56[/C][C]87.4223578061097[/C][C]1.84389496214624[/C][/ROW]
[ROW][C]0.61[/C][C]89.001847663559[/C][C]1.81749445139327[/C][/ROW]
[ROW][C]0.66[/C][C]90.4493304616487[/C][C]1.48951425425951[/C][/ROW]
[ROW][C]0.71[/C][C]91.5624851802817[/C][C]1.01176342940769[/C][/ROW]
[ROW][C]0.76[/C][C]92.5455694511[/C][C]1.06814266940975[/C][/ROW]
[ROW][C]0.81[/C][C]94.241101957139[/C][C]2.10219325152500[/C][/ROW]
[ROW][C]0.86[/C][C]96.9075023601776[/C][C]2.3513268396712[/C][/ROW]
[ROW][C]0.91[/C][C]99.537796310459[/C][C]1.71567013900682[/C][/ROW]
[ROW][C]0.96[/C][C]102.795078937544[/C][C]1.61293967388695[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=1517&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1517&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.0169.03716033049572.57421141308358
0.0673.3107648178421.13721456461338
0.1175.13408445078981.12620725670397
0.1676.68135234463551.29239105902577
0.2178.17236711415031.39786055262194
0.2679.62933401441411.46919688433279
0.3180.94267702754231.27062376297728
0.3682.03425122425621.15481178659679
0.4183.1335184045831.391697870781
0.4684.4313105043571.65057646107131
0.5185.88361953746651.76811912874389
0.5687.42235780610971.84389496214624
0.6189.0018476635591.81749445139327
0.6690.44933046164871.48951425425951
0.7191.56248518028171.01176342940769
0.7692.54556945111.06814266940975
0.8194.2411019571392.10219325152500
0.8696.90750236017762.3513268396712
0.9199.5377963104591.71567013900682
0.96102.7950789375441.61293967388695



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