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 computationThu, 19 Aug 2010 15:39:20 +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/2010/Aug/19/t1282232525hx6o0gl4mlet5gd.htm/, Retrieved Fri, 03 May 2024 06:09:02 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=79319, Retrieved Fri, 03 May 2024 06:09:02 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsellen aerts
Estimated Impact119
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] ["Tijdreeks A stap...] [2010-08-19 15:39:20] [6e43eada780a1520be8ab5bc59456d41] [Current]
Feedback Forum

Post a new message
Dataseries X:
25
24
23
21
41
40
25
15
16
16
17
19
18
19
20
21
46
47
30
16
15
18
30
31
32
36
30
31
61
57
45
33
31
36
46
49
34
40
41
48
75
77
71
54
50
56
66
66
48
63
71
70
88
92
91
80
81
81
98
106
85
93
96
92
115
109
119
107
107
106
132
143
120
123
132
136
158
151
155
138
143
139
168
182
154
158
167
170
197
190
196
174
180
171
200
215
184
186
197
186
211
205
218
199
213
207
236
248
211
220
235
223
245
236
253
246
255
248
274
288




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=79319&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'Gwilym Jenkins' @ 72.249.127.135







Harrell-Davis Quantiles
quantilesvaluestandard error
0.121.70504036585072.92800073135328
0.233.6512566870974.13986045204726
0.348.10863456154626.57987215166813
0.469.737945509574710.3638298556173
0.594.825976469611.5225901957908
0.6126.56782009854915.2689449203561
0.7162.33024540563913.6683014077816
0.8193.9896247831489.74164791122706
0.9225.37801273890610.6383875997951

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 21.7050403658507 & 2.92800073135328 \tabularnewline
0.2 & 33.651256687097 & 4.13986045204726 \tabularnewline
0.3 & 48.1086345615462 & 6.57987215166813 \tabularnewline
0.4 & 69.7379455095747 & 10.3638298556173 \tabularnewline
0.5 & 94.8259764696 & 11.5225901957908 \tabularnewline
0.6 & 126.567820098549 & 15.2689449203561 \tabularnewline
0.7 & 162.330245405639 & 13.6683014077816 \tabularnewline
0.8 & 193.989624783148 & 9.74164791122706 \tabularnewline
0.9 & 225.378012738906 & 10.6383875997951 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=79319&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]21.7050403658507[/C][C]2.92800073135328[/C][/ROW]
[ROW][C]0.2[/C][C]33.651256687097[/C][C]4.13986045204726[/C][/ROW]
[ROW][C]0.3[/C][C]48.1086345615462[/C][C]6.57987215166813[/C][/ROW]
[ROW][C]0.4[/C][C]69.7379455095747[/C][C]10.3638298556173[/C][/ROW]
[ROW][C]0.5[/C][C]94.8259764696[/C][C]11.5225901957908[/C][/ROW]
[ROW][C]0.6[/C][C]126.567820098549[/C][C]15.2689449203561[/C][/ROW]
[ROW][C]0.7[/C][C]162.330245405639[/C][C]13.6683014077816[/C][/ROW]
[ROW][C]0.8[/C][C]193.989624783148[/C][C]9.74164791122706[/C][/ROW]
[ROW][C]0.9[/C][C]225.378012738906[/C][C]10.6383875997951[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=79319&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=79319&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.121.70504036585072.92800073135328
0.233.6512566870974.13986045204726
0.348.10863456154626.57987215166813
0.469.737945509574710.3638298556173
0.594.825976469611.5225901957908
0.6126.56782009854915.2689449203561
0.7162.33024540563913.6683014077816
0.8193.9896247831489.74164791122706
0.9225.37801273890610.6383875997951



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')