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 computationFri, 05 Jun 2009 04:15:43 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2009/Jun/05/t1244196992j3lvfd2jcw8g70n.htm/, Retrieved Thu, 09 May 2024 20:18:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=41777, Retrieved Thu, 09 May 2024 20:18:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsDuncan Huysmans opgave 4
Estimated Impact178
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Duncan Huysmans o...] [2009-06-05 10:15:43] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
49.1
52.7
55.1
60.4
57.5
55.9
62.6
61.3
67.1
68.6
74.8
76.9
85.7
86.5
90.8
89.7
96.3
107.5
109.2
100.2
116.8
120.1
123.3
130.2
131.4
125.6
124.5
134.3
135.2
151.8
146.4
140.0
127.8
148.0
165.9
165.5
179.9
190.0
189.8
190.9
203.6
183.5
169.3
144.2
141.5
154.3
169.5
194.0
203.2
192.9
209.4
227.2
263.7
297.8
337.1
361.3
355.2
312.6
309.9
323.7
324.1
355.3
383.4
395.1
412.8
407.0
438.0
446.1
452.5
447.3
475.9
487.7




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41777&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.164.22188658050516.49082106121734
0.289.481419905376612.0478837016291
0.3116.55477259935111.5820719510960
0.4135.5500685930229.35434976553517
0.5156.67280654948413.9728402362689
0.6184.43322810199415.8132953987837
0.7236.47182485341842.7887481275058
0.8328.57377264632035.0982573073144
0.9409.8662500853328.0957683128366

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.1 & 64.2218865805051 & 6.49082106121734 \tabularnewline
0.2 & 89.4814199053766 & 12.0478837016291 \tabularnewline
0.3 & 116.554772599351 & 11.5820719510960 \tabularnewline
0.4 & 135.550068593022 & 9.35434976553517 \tabularnewline
0.5 & 156.672806549484 & 13.9728402362689 \tabularnewline
0.6 & 184.433228101994 & 15.8132953987837 \tabularnewline
0.7 & 236.471824853418 & 42.7887481275058 \tabularnewline
0.8 & 328.573772646320 & 35.0982573073144 \tabularnewline
0.9 & 409.86625008533 & 28.0957683128366 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41777&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]64.2218865805051[/C][C]6.49082106121734[/C][/ROW]
[ROW][C]0.2[/C][C]89.4814199053766[/C][C]12.0478837016291[/C][/ROW]
[ROW][C]0.3[/C][C]116.554772599351[/C][C]11.5820719510960[/C][/ROW]
[ROW][C]0.4[/C][C]135.550068593022[/C][C]9.35434976553517[/C][/ROW]
[ROW][C]0.5[/C][C]156.672806549484[/C][C]13.9728402362689[/C][/ROW]
[ROW][C]0.6[/C][C]184.433228101994[/C][C]15.8132953987837[/C][/ROW]
[ROW][C]0.7[/C][C]236.471824853418[/C][C]42.7887481275058[/C][/ROW]
[ROW][C]0.8[/C][C]328.573772646320[/C][C]35.0982573073144[/C][/ROW]
[ROW][C]0.9[/C][C]409.86625008533[/C][C]28.0957683128366[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41777&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41777&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.164.22188658050516.49082106121734
0.289.481419905376612.0478837016291
0.3116.55477259935111.5820719510960
0.4135.5500685930229.35434976553517
0.5156.67280654948413.9728402362689
0.6184.43322810199415.8132953987837
0.7236.47182485341842.7887481275058
0.8328.57377264632035.0982573073144
0.9409.8662500853328.0957683128366



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