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 computationMon, 08 Mar 2010 10:04:04 -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/2010/Mar/08/t12680679624jq2b4h9tlwaher.htm/, Retrieved Fri, 29 Mar 2024 05:30:04 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=74157, Retrieved Fri, 29 Mar 2024 05:30:04 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact161
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [] [2010-03-08 17:04:04] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
10698
11514
13253
13807
13971
15548
17351
17409
17537
17704
17970
18004
18089
18479
18730
19169
19684
19785
20366
20506
20764
21382
21705
21725
21790
22205
22228
22397
22405
22414
22679
22780
22782
23843
23970
24319
24561
24693
25316
25591
25814
26105
26482
26887
27044
27071
27192
28029
28352
28806
29096
29133
29383
29462
29743
29965
30212
30364
31514
32034
32609
34777
35022
35338
35953
36438
36845
37702
37863
39690
41086
43129




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=74157&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.116553.14030344201261.71904227439
0.218904.9333628657936.419353564739
0.321077.7807959141936.211321502927
0.422580.9310873708749.572903183123
0.524414.4514770231206.61986807374
0.626621.53178853821161.24460010690
0.728764.48045991421015.35169792374
0.831361.37668237891665.03966210523
0.936034.96789361461391.99079053831

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 16553.1403034420 & 1261.71904227439 \tabularnewline
0.2 & 18904.9333628657 & 936.419353564739 \tabularnewline
0.3 & 21077.7807959141 & 936.211321502927 \tabularnewline
0.4 & 22580.9310873708 & 749.572903183123 \tabularnewline
0.5 & 24414.451477023 & 1206.61986807374 \tabularnewline
0.6 & 26621.5317885382 & 1161.24460010690 \tabularnewline
0.7 & 28764.4804599142 & 1015.35169792374 \tabularnewline
0.8 & 31361.3766823789 & 1665.03966210523 \tabularnewline
0.9 & 36034.9678936146 & 1391.99079053831 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=74157&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]16553.1403034420[/C][C]1261.71904227439[/C][/ROW]
[ROW][C]0.2[/C][C]18904.9333628657[/C][C]936.419353564739[/C][/ROW]
[ROW][C]0.3[/C][C]21077.7807959141[/C][C]936.211321502927[/C][/ROW]
[ROW][C]0.4[/C][C]22580.9310873708[/C][C]749.572903183123[/C][/ROW]
[ROW][C]0.5[/C][C]24414.451477023[/C][C]1206.61986807374[/C][/ROW]
[ROW][C]0.6[/C][C]26621.5317885382[/C][C]1161.24460010690[/C][/ROW]
[ROW][C]0.7[/C][C]28764.4804599142[/C][C]1015.35169792374[/C][/ROW]
[ROW][C]0.8[/C][C]31361.3766823789[/C][C]1665.03966210523[/C][/ROW]
[ROW][C]0.9[/C][C]36034.9678936146[/C][C]1391.99079053831[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=74157&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=74157&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.116553.14030344201261.71904227439
0.218904.9333628657936.419353564739
0.321077.7807959141936.211321502927
0.422580.9310873708749.572903183123
0.524414.4514770231206.61986807374
0.626621.53178853821161.24460010690
0.728764.48045991421015.35169792374
0.831361.37668237891665.03966210523
0.936034.96789361461391.99079053831



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