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 09:48:33 -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/gkd1nu3wjyarxxa1193071490.htm/, Retrieved Mon, 06 May 2024 02:12:18 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=1468, Retrieved Mon, 06 May 2024 02:12:18 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords80%, Confidence Interval, Q6, Total Production
Estimated Impact191
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [80% confidence in...] [2007-10-22 16:48:33] [9fe578921d87f9af8e79a90d6142ba02] [Current]
Feedback Forum

Post a new message
Dataseries X:
110,40
96,40
101,90
106,20
81,00
94,70
101,00
109,40
102,30
90,70
96,20
96,10
106,00
103,10
102,00
104,70
86,00
92,10
106,90
112,60
101,70
92,00
97,40
97,00
105,40
102,70
98,10
104,50
87,40
89,90
109,80
111,70
98,60
96,90
95,10
97,00
112,70
102,90
97,40
111,40
87,40
96,80
114,10
110,30
103,90
101,60
94,60
95,90
104,70
102,80
98,10
113,90
80,90
95,70
113,20
105,90
108,80
102,30
99,00
100,70
115,50




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1468&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.189.99141564358462.37418160585810
0.295.00132410074661.27198204922879
0.396.8134051839340.725863145593432
0.498.655450860861.4799064365811
0.5101.2318360290211.30308488044260
0.6103.0240763527630.994073412489012
0.7105.1479979277841.35610783414605
0.8108.4636416800831.83726345527739
0.9111.9836800237171.10775542432169

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 89.9914156435846 & 2.37418160585810 \tabularnewline
0.2 & 95.0013241007466 & 1.27198204922879 \tabularnewline
0.3 & 96.813405183934 & 0.725863145593432 \tabularnewline
0.4 & 98.65545086086 & 1.4799064365811 \tabularnewline
0.5 & 101.231836029021 & 1.30308488044260 \tabularnewline
0.6 & 103.024076352763 & 0.994073412489012 \tabularnewline
0.7 & 105.147997927784 & 1.35610783414605 \tabularnewline
0.8 & 108.463641680083 & 1.83726345527739 \tabularnewline
0.9 & 111.983680023717 & 1.10775542432169 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=1468&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]89.9914156435846[/C][C]2.37418160585810[/C][/ROW]
[ROW][C]0.2[/C][C]95.0013241007466[/C][C]1.27198204922879[/C][/ROW]
[ROW][C]0.3[/C][C]96.813405183934[/C][C]0.725863145593432[/C][/ROW]
[ROW][C]0.4[/C][C]98.65545086086[/C][C]1.4799064365811[/C][/ROW]
[ROW][C]0.5[/C][C]101.231836029021[/C][C]1.30308488044260[/C][/ROW]
[ROW][C]0.6[/C][C]103.024076352763[/C][C]0.994073412489012[/C][/ROW]
[ROW][C]0.7[/C][C]105.147997927784[/C][C]1.35610783414605[/C][/ROW]
[ROW][C]0.8[/C][C]108.463641680083[/C][C]1.83726345527739[/C][/ROW]
[ROW][C]0.9[/C][C]111.983680023717[/C][C]1.10775542432169[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=1468&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1468&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.189.99141564358462.37418160585810
0.295.00132410074661.27198204922879
0.396.8134051839340.725863145593432
0.498.655450860861.4799064365811
0.5101.2318360290211.30308488044260
0.6103.0240763527630.994073412489012
0.7105.1479979277841.35610783414605
0.8108.4636416800831.83726345527739
0.9111.9836800237171.10775542432169



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