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, 02 Jan 2012 08:45:24 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Jan/02/t1325511957rrdg9usj21rj7nd.htm/, Retrieved Sat, 04 May 2024 12:44:20 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=160918, Retrieved Sat, 04 May 2024 12:44:20 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKDG201142
Estimated Impact161
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [decielen] [2012-01-02 13:45:24] [5f178b5bce8a01d64692a8a5c649399b] [Current]
Feedback Forum

Post a new message
Dataseries X:
123.46
123.24
123.86
124.28
124.78
125.19
125.46
127.60
127.80
126.63
127.06
126.77
127.05
128.23
128.60
128.97
129.34
129.71
130.08
130.45
128.82
132.19
131.56
131.93
132.30
130.67
133.05
132.42
133.79
134.16
134.53
136.90
135.27
136.64
136.01
136.38
136.75
138.12
137.50
137.87
138.24
138.61
138.98
140.35
139.72
143.09
140.46
141.83
143.20
140.57
141.95
141.32
142.69
143.06
144.43
143.80
144.17
144.54
146.91
145.28
144.65
145.02
144.40
146.77
146.14
147.51
148.88
148.25
147.62
150.99
148.36
149.73
150.10




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net

\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' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=160918&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' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=160918&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160918&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' @ jenkins.wessa.net







Harrell-Davis Quantiles
quantilesvaluestandard error
0.1126.0273359751110.979200326110047
0.2128.53087951360.995160959740071
0.3131.0789131906871.39274115227044
0.4134.0876814479661.72127008287878
0.5137.0726680283471.43251744744574
0.6139.7773577435031.52919104141471
0.7142.5654715544221.25710817993663
0.8144.6995437096140.872948077379848
0.9147.4192863154830.943293560649678

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 126.027335975111 & 0.979200326110047 \tabularnewline
0.2 & 128.5308795136 & 0.995160959740071 \tabularnewline
0.3 & 131.078913190687 & 1.39274115227044 \tabularnewline
0.4 & 134.087681447966 & 1.72127008287878 \tabularnewline
0.5 & 137.072668028347 & 1.43251744744574 \tabularnewline
0.6 & 139.777357743503 & 1.52919104141471 \tabularnewline
0.7 & 142.565471554422 & 1.25710817993663 \tabularnewline
0.8 & 144.699543709614 & 0.872948077379848 \tabularnewline
0.9 & 147.419286315483 & 0.943293560649678 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=160918&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]126.027335975111[/C][C]0.979200326110047[/C][/ROW]
[ROW][C]0.2[/C][C]128.5308795136[/C][C]0.995160959740071[/C][/ROW]
[ROW][C]0.3[/C][C]131.078913190687[/C][C]1.39274115227044[/C][/ROW]
[ROW][C]0.4[/C][C]134.087681447966[/C][C]1.72127008287878[/C][/ROW]
[ROW][C]0.5[/C][C]137.072668028347[/C][C]1.43251744744574[/C][/ROW]
[ROW][C]0.6[/C][C]139.777357743503[/C][C]1.52919104141471[/C][/ROW]
[ROW][C]0.7[/C][C]142.565471554422[/C][C]1.25710817993663[/C][/ROW]
[ROW][C]0.8[/C][C]144.699543709614[/C][C]0.872948077379848[/C][/ROW]
[ROW][C]0.9[/C][C]147.419286315483[/C][C]0.943293560649678[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=160918&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=160918&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.1126.0273359751110.979200326110047
0.2128.53087951360.995160959740071
0.3131.0789131906871.39274115227044
0.4134.0876814479661.72127008287878
0.5137.0726680283471.43251744744574
0.6139.7773577435031.52919104141471
0.7142.5654715544221.25710817993663
0.8144.6995437096140.872948077379848
0.9147.4192863154830.943293560649678



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