Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_harrell_davies.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationWed, 21 Oct 2009 06:31:26 -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/Oct/21/t125612850052jitys8479i99e.htm/, Retrieved Thu, 02 May 2024 09:34:10 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=49344, Retrieved Thu, 02 May 2024 09:34:10 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact153
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Bivariate Data Series] [Bivariate dataset] [2008-01-05 23:51:08] [74be16979710d4c4e7c6647856088456]
F RMPD  [Univariate Explorative Data Analysis] [Colombia Coffee] [2008-01-07 14:21:11] [74be16979710d4c4e7c6647856088456]
F RMPD    [Univariate Data Series] [] [2009-10-14 08:30:28] [74be16979710d4c4e7c6647856088456]
- RMP       [Percentiles] [DSHW-WS3-1.3] [2009-10-20 16:49:39] [f15cfb7053d35072d573abca87df96a0]
- RMPD          [Harrell-Davis Quantiles] [workshop3.2.3] [2009-10-21 12:31:26] [ff423994c38282a6d306f7d0147a5924] [Current]
Feedback Forum

Post a new message
Dataseries X:
87.28
87.28
87.09
86.92
87.59
90.72
90.69
90.3
89.55
88.94
88.41
87.82
87.07
86.82
86.4
86.02
85.66
85.32
85
84.67
83.94
82.83
81.95
81.19
80.48
78.86
69.47
68.77
70.06
73.95
75.8
77.79
81.57
83.07
84.34
85.1
85.25
84.26
83.63
86.44
85.3
84.1
83.36
82.48
81.58
80.47
79.34
82.13
81.69
80.7
79.88
79.16
78.38
77.42
76.47
75.46
74.48
78.27
80.7
79.91
78.75
77.78
81.14
81.08
80.03
78.91
78.01
76.9
75.97
81.93
80.27
78.67
77.42
76.16
74.7
76.39
76.04
74.65
73.29
71.79
74.39
74.91
74.54
73.08
72.75
71.32
70.38
70.35
70.01
69.36
67.77
69.26
69.8
68.38
67.62
68.39
66.95
65.21
66.64
63.45
60.66
62.34
60.32
58.64
60.46
58.59
61.87
61.85
67.44
77.06
91.74
93.15
94.15
93.11
91.51
89.96
88.16
86.98
88.03
86.24
84.65
83.23
81.7
80.25
78.8
77.51
76.2
75.04
74
75.49
77.14
76.15
76.27
78.19
76.49
77.31
76.65
74.99
73.51
72.07
70.59
71.96
76.29
74.86
74.93
71.9
71.01
77.47
75.78
76.6
76.07
74.57
73.02
72.65
73.16
71.53
69.78
67.98
69.96
72.16
70.47
68.86
67.37
65.87
72.16
71.34
69.93
68.44
67.16
66.01
67.25
70.91
69.75
68.59
67.48
66.31
64.81
66.58
65.97
64.7
64.7
60.94
59.08
58.42
57.77
57.11
53.31
49.96
49.4
48.84
48.3
47.74
47.24
46.76
46.29
48.9
49.23
48.53
48.03
54.34
53.79
53.24
52.96
52.17
51.7
58.55
78.2
77.03
76.19
77.15
75.87
95.47
109.67
112.28
112.01
107.93
105.96
105.06
102.98
102.2
105.23
101.85
99.89
96.23
94.76
91.51
91.63
91.54
85.23
87.83
87.38
84.44
85.19
84.03
86.73
102.52
104.45
106.98
107.02
99.26
94.45
113.44
157.33
147.38
171.89
171.95
132.71
126.02
121.18
115.45
110.48
117.85
117.63
124.65
109.59
111.27
99.78
98.21
99.2
97.97
89.55
87.91
93.34
94.42
93.2
90.29
91.46
89.98
88.35
88.41
82.44
79.89
75.69
75.66
84.5
96.73
87.48
82.39
83.48
79.31
78.16
72.77
72.45
68.46
67.62
68.76
70.07
68.55
65.3
58.96
59.17
62.37
66.28
55.62
55.23
55.85
56.75
50.89
53.88
52.95
55.08
53.61
58.78
61.85
55.91
53.32
46.41
44.57
50
50
53.36
46.23
50.45
49.07
45.85
48.45
49.96
46.53
50.51
47.58
48.05
46.84
47.67
49.16
55.54
55.82
58.22
56.19
57.77
63.19
54.76
55.74
62.54
61.39
69.6
79.23
80
93.68
107.63
100.18
97.3
90.45
80.64
80.58
75.82
85.59
89.35
89.42
104.73
95.32
89.27
90.44
86.97
79.98
81.22
87.35
83.64
82.22
94.4
102.18




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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 & 2 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=49344&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]2 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=49344&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49344&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 time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Harrell-Davis Quantiles
quantilesvaluestandard error
0.02547.28658948631820.546341372358286
0.07550.64997651131711.27159636157667
0.12555.18175703225091.33526114349448
0.17559.493371569561.81885285863324
0.22564.56504417906651.99132137071581
0.27567.88742216051651.04016061046058
0.32569.9645984602860.976870163090224
0.37572.26654956563011.23936597244925
0.42574.63523720930491.01611315746269
0.47576.19528615312850.656027125845838
0.52577.65965146664010.869484516384336
0.57579.51438371173220.963638056287923
0.62581.44321988815121.01971894998476
0.67583.70759552908551.10249722194178
0.72585.91747146600720.985244729277902
0.77587.94074904658920.907693297263047
0.82590.6616359930371.26106525568625
0.87595.274484827822.10713833741742
0.925103.7062931697872.52592332917429
0.975119.5039686831685.54742147865435

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.025 & 47.2865894863182 & 0.546341372358286 \tabularnewline
0.075 & 50.6499765113171 & 1.27159636157667 \tabularnewline
0.125 & 55.1817570322509 & 1.33526114349448 \tabularnewline
0.175 & 59.49337156956 & 1.81885285863324 \tabularnewline
0.225 & 64.5650441790665 & 1.99132137071581 \tabularnewline
0.275 & 67.8874221605165 & 1.04016061046058 \tabularnewline
0.325 & 69.964598460286 & 0.976870163090224 \tabularnewline
0.375 & 72.2665495656301 & 1.23936597244925 \tabularnewline
0.425 & 74.6352372093049 & 1.01611315746269 \tabularnewline
0.475 & 76.1952861531285 & 0.656027125845838 \tabularnewline
0.525 & 77.6596514666401 & 0.869484516384336 \tabularnewline
0.575 & 79.5143837117322 & 0.963638056287923 \tabularnewline
0.625 & 81.4432198881512 & 1.01971894998476 \tabularnewline
0.675 & 83.7075955290855 & 1.10249722194178 \tabularnewline
0.725 & 85.9174714660072 & 0.985244729277902 \tabularnewline
0.775 & 87.9407490465892 & 0.907693297263047 \tabularnewline
0.825 & 90.661635993037 & 1.26106525568625 \tabularnewline
0.875 & 95.27448482782 & 2.10713833741742 \tabularnewline
0.925 & 103.706293169787 & 2.52592332917429 \tabularnewline
0.975 & 119.503968683168 & 5.54742147865435 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=49344&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.025[/C][C]47.2865894863182[/C][C]0.546341372358286[/C][/ROW]
[ROW][C]0.075[/C][C]50.6499765113171[/C][C]1.27159636157667[/C][/ROW]
[ROW][C]0.125[/C][C]55.1817570322509[/C][C]1.33526114349448[/C][/ROW]
[ROW][C]0.175[/C][C]59.49337156956[/C][C]1.81885285863324[/C][/ROW]
[ROW][C]0.225[/C][C]64.5650441790665[/C][C]1.99132137071581[/C][/ROW]
[ROW][C]0.275[/C][C]67.8874221605165[/C][C]1.04016061046058[/C][/ROW]
[ROW][C]0.325[/C][C]69.964598460286[/C][C]0.976870163090224[/C][/ROW]
[ROW][C]0.375[/C][C]72.2665495656301[/C][C]1.23936597244925[/C][/ROW]
[ROW][C]0.425[/C][C]74.6352372093049[/C][C]1.01611315746269[/C][/ROW]
[ROW][C]0.475[/C][C]76.1952861531285[/C][C]0.656027125845838[/C][/ROW]
[ROW][C]0.525[/C][C]77.6596514666401[/C][C]0.869484516384336[/C][/ROW]
[ROW][C]0.575[/C][C]79.5143837117322[/C][C]0.963638056287923[/C][/ROW]
[ROW][C]0.625[/C][C]81.4432198881512[/C][C]1.01971894998476[/C][/ROW]
[ROW][C]0.675[/C][C]83.7075955290855[/C][C]1.10249722194178[/C][/ROW]
[ROW][C]0.725[/C][C]85.9174714660072[/C][C]0.985244729277902[/C][/ROW]
[ROW][C]0.775[/C][C]87.9407490465892[/C][C]0.907693297263047[/C][/ROW]
[ROW][C]0.825[/C][C]90.661635993037[/C][C]1.26106525568625[/C][/ROW]
[ROW][C]0.875[/C][C]95.27448482782[/C][C]2.10713833741742[/C][/ROW]
[ROW][C]0.925[/C][C]103.706293169787[/C][C]2.52592332917429[/C][/ROW]
[ROW][C]0.975[/C][C]119.503968683168[/C][C]5.54742147865435[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=49344&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=49344&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.02547.28658948631820.546341372358286
0.07550.64997651131711.27159636157667
0.12555.18175703225091.33526114349448
0.17559.493371569561.81885285863324
0.22564.56504417906651.99132137071581
0.27567.88742216051651.04016061046058
0.32569.9645984602860.976870163090224
0.37572.26654956563011.23936597244925
0.42574.63523720930491.01611315746269
0.47576.19528615312850.656027125845838
0.52577.65965146664010.869484516384336
0.57579.51438371173220.963638056287923
0.62581.44321988815121.01971894998476
0.67583.70759552908551.10249722194178
0.72585.91747146600720.985244729277902
0.77587.94074904658920.907693297263047
0.82590.6616359930371.26106525568625
0.87595.274484827822.10713833741742
0.925103.7062931697872.52592332917429
0.975119.5039686831685.54742147865435



Parameters (Session):
par1 = 0.025 ; par2 = 0.975 ; par3 = 0.05 ;
Parameters (R input):
par1 = 0.025 ; par2 = 0.975 ; par3 = 0.05 ;
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')