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, 26 Feb 2008 06:58:33 -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/2008/Feb/26/t1204034404t5ojn90spqwcynt.htm/, Retrieved Sun, 05 May 2024 00:30:10 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=9135, Retrieved Sun, 05 May 2024 00:30:10 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact176
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Harrel-Davis Deci...] [2008-02-26 13:58:33] [18a9c7462de321e50920adc850709d38] [Current]
Feedback Forum

Post a new message
Dataseries X:
161,79
161,79
161,85
161,77
161,86
161,89
161,89
161,89
162,18
162,43
162,58
162,57
162,57
162,57
162,44
162,79
163,15
163,23
163,23
163,23
163,38
163,71
163,73
163,73
163,73
163,73
163,93
164,27
164,57
164,73
164,73
164,76
165,75
165,86
165,99
166,13
166,13
166,13
166,15
166,45
166,48
166,51
166,51
166,51
166,58
166,82
167,35
167,5
167,5
167,6
167,72
167,29
166,98
166,98
166,98
166,98
167,63
167,83
167,85
167,87
167,87
167,96
167,7
169,25
168,79
168,77
168,77
169
168,92
169,23
169,28
169,29




Summary of compuational 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 compuational 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=9135&T=0

[TABLE]
[ROW][C]Summary of compuational 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=9135&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=9135&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 compuational 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.1162.0499190144380.226985761610883
0.2162.7707251127400.327322076889594
0.3163.5583167229530.39899395982576
0.4164.6078506077370.743484447947587
0.5165.9032799550940.546209491318028
0.6166.6052905399040.304005371031202
0.7167.2141132476520.316120923155895
0.8167.7592226293410.215526410188946
0.9168.6659768089260.374706586086423

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 162.049919014438 & 0.226985761610883 \tabularnewline
0.2 & 162.770725112740 & 0.327322076889594 \tabularnewline
0.3 & 163.558316722953 & 0.39899395982576 \tabularnewline
0.4 & 164.607850607737 & 0.743484447947587 \tabularnewline
0.5 & 165.903279955094 & 0.546209491318028 \tabularnewline
0.6 & 166.605290539904 & 0.304005371031202 \tabularnewline
0.7 & 167.214113247652 & 0.316120923155895 \tabularnewline
0.8 & 167.759222629341 & 0.215526410188946 \tabularnewline
0.9 & 168.665976808926 & 0.374706586086423 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=9135&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]162.049919014438[/C][C]0.226985761610883[/C][/ROW]
[ROW][C]0.2[/C][C]162.770725112740[/C][C]0.327322076889594[/C][/ROW]
[ROW][C]0.3[/C][C]163.558316722953[/C][C]0.39899395982576[/C][/ROW]
[ROW][C]0.4[/C][C]164.607850607737[/C][C]0.743484447947587[/C][/ROW]
[ROW][C]0.5[/C][C]165.903279955094[/C][C]0.546209491318028[/C][/ROW]
[ROW][C]0.6[/C][C]166.605290539904[/C][C]0.304005371031202[/C][/ROW]
[ROW][C]0.7[/C][C]167.214113247652[/C][C]0.316120923155895[/C][/ROW]
[ROW][C]0.8[/C][C]167.759222629341[/C][C]0.215526410188946[/C][/ROW]
[ROW][C]0.9[/C][C]168.665976808926[/C][C]0.374706586086423[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=9135&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=9135&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.1162.0499190144380.226985761610883
0.2162.7707251127400.327322076889594
0.3163.5583167229530.39899395982576
0.4164.6078506077370.743484447947587
0.5165.9032799550940.546209491318028
0.6166.6052905399040.304005371031202
0.7167.2141132476520.316120923155895
0.8167.7592226293410.215526410188946
0.9168.6659768089260.374706586086423



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