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 computationTue, 17 May 2011 13:27:44 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/May/17/t1305638810tzk2rzmmut609tn.htm/, Retrieved Fri, 10 May 2024 08:29:54 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=121723, Retrieved Fri, 10 May 2024 08:29:54 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDGP1W42
Estimated Impact124
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2011-05-17 13:27:44] [7ec7c2b78d31d17737f91280a6e28d4a] [Current]
Feedback Forum

Post a new message
Dataseries X:
90
51
47
59
54
79
59
80
46
62
55
77
72
72
71
50
66
78
59
52
71
98
70
84
90
98
98
78
59
0
58
55
62
80
91
86
61
49
61
56
73
85
82
32
39
30
51
48
57
59
32
56
54
74
62
78
72
48
59
61
80
69
58
63
27
23
34
45
51
51
73
37
35
66
54
30
66
61
37
55
64
53
63
70
72
52
53
50
60
73
66
78




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=121723&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=121723&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=121723&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 time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk







Harrell-Davis Quantiles
quantilesvaluestandard error
0.136.71184366028064.35481513099062
0.248.96648056160712.27588328339022
0.353.17568102970541.64850077638698
0.457.00583904693531.82675966987016
0.560.29983410898271.53363813470093
0.664.52803812970072.7389323439615
0.770.56157923869732.3320413007541
0.875.96253356570332.48374265263302
0.983.24257696748823.1385315281003

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 36.7118436602806 & 4.35481513099062 \tabularnewline
0.2 & 48.9664805616071 & 2.27588328339022 \tabularnewline
0.3 & 53.1756810297054 & 1.64850077638698 \tabularnewline
0.4 & 57.0058390469353 & 1.82675966987016 \tabularnewline
0.5 & 60.2998341089827 & 1.53363813470093 \tabularnewline
0.6 & 64.5280381297007 & 2.7389323439615 \tabularnewline
0.7 & 70.5615792386973 & 2.3320413007541 \tabularnewline
0.8 & 75.9625335657033 & 2.48374265263302 \tabularnewline
0.9 & 83.2425769674882 & 3.1385315281003 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=121723&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]36.7118436602806[/C][C]4.35481513099062[/C][/ROW]
[ROW][C]0.2[/C][C]48.9664805616071[/C][C]2.27588328339022[/C][/ROW]
[ROW][C]0.3[/C][C]53.1756810297054[/C][C]1.64850077638698[/C][/ROW]
[ROW][C]0.4[/C][C]57.0058390469353[/C][C]1.82675966987016[/C][/ROW]
[ROW][C]0.5[/C][C]60.2998341089827[/C][C]1.53363813470093[/C][/ROW]
[ROW][C]0.6[/C][C]64.5280381297007[/C][C]2.7389323439615[/C][/ROW]
[ROW][C]0.7[/C][C]70.5615792386973[/C][C]2.3320413007541[/C][/ROW]
[ROW][C]0.8[/C][C]75.9625335657033[/C][C]2.48374265263302[/C][/ROW]
[ROW][C]0.9[/C][C]83.2425769674882[/C][C]3.1385315281003[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=121723&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=121723&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.136.71184366028064.35481513099062
0.248.96648056160712.27588328339022
0.353.17568102970541.64850077638698
0.457.00583904693531.82675966987016
0.560.29983410898271.53363813470093
0.664.52803812970072.7389323439615
0.770.56157923869732.3320413007541
0.875.96253356570332.48374265263302
0.983.24257696748823.1385315281003



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