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 computationWed, 04 Aug 2010 15:23:11 +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/04/t1280935535r901f0c79ewojid.htm/, Retrieved Fri, 03 May 2024 09:27:23 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=78342, Retrieved Fri, 03 May 2024 09:27:23 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsMathias Goossenaerts
Estimated Impact97
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Tijdreeks 1 - Sta...] [2010-08-04 15:23:11] [f7fc4e1bbbe57039ee5ebdd2c8b864c0] [Current]
Feedback Forum

Post a new message
Dataseries X:
36
35
34
32
52
51
36
26
27
27
28
30
28
29
25
28
55
53
42
32
37
41
37
38
39
32
36
39
83
83
66
53
72
77
69
72
81
71
63
66
114
116
109
97
111
120
110
106
115
110
103
112
163
166
156
140
166
176
163
162
171
167
163
168
222
216
197
178
204
220
196
195
213
218
216
225
280
272
252
230
248
259
240
237
252
250
255
255
313
291
271
247
268
283
259
259
267
270
279
269
334
326
301
276
301
313
291
287
289
298
320
312
385
380
351
322
350
363
344
345




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=78342&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.133.17335970553872.42157966224297
0.247.47263669717719.05267354907778
0.378.746084028684614.6795690611279
0.4119.67380489466619.7707340953180
0.5168.95290790107819.0800859305708
0.6215.45860430753318.8995766198814
0.7252.52837158232511.0416332030123
0.8277.9601231542429.29912263352767
0.9314.29112927489811.6732392528071

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 33.1733597055387 & 2.42157966224297 \tabularnewline
0.2 & 47.4726366971771 & 9.05267354907778 \tabularnewline
0.3 & 78.7460840286846 & 14.6795690611279 \tabularnewline
0.4 & 119.673804894666 & 19.7707340953180 \tabularnewline
0.5 & 168.952907901078 & 19.0800859305708 \tabularnewline
0.6 & 215.458604307533 & 18.8995766198814 \tabularnewline
0.7 & 252.528371582325 & 11.0416332030123 \tabularnewline
0.8 & 277.960123154242 & 9.29912263352767 \tabularnewline
0.9 & 314.291129274898 & 11.6732392528071 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=78342&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]33.1733597055387[/C][C]2.42157966224297[/C][/ROW]
[ROW][C]0.2[/C][C]47.4726366971771[/C][C]9.05267354907778[/C][/ROW]
[ROW][C]0.3[/C][C]78.7460840286846[/C][C]14.6795690611279[/C][/ROW]
[ROW][C]0.4[/C][C]119.673804894666[/C][C]19.7707340953180[/C][/ROW]
[ROW][C]0.5[/C][C]168.952907901078[/C][C]19.0800859305708[/C][/ROW]
[ROW][C]0.6[/C][C]215.458604307533[/C][C]18.8995766198814[/C][/ROW]
[ROW][C]0.7[/C][C]252.528371582325[/C][C]11.0416332030123[/C][/ROW]
[ROW][C]0.8[/C][C]277.960123154242[/C][C]9.29912263352767[/C][/ROW]
[ROW][C]0.9[/C][C]314.291129274898[/C][C]11.6732392528071[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=78342&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=78342&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.133.17335970553872.42157966224297
0.247.47263669717719.05267354907778
0.378.746084028684614.6795690611279
0.4119.67380489466619.7707340953180
0.5168.95290790107819.0800859305708
0.6215.45860430753318.8995766198814
0.7252.52837158232511.0416332030123
0.8277.9601231542429.29912263352767
0.9314.29112927489811.6732392528071



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