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 computationMon, 19 Oct 2009 11:58:27 -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/19/t1255975222wfsh8goi8m4ku03.htm/, Retrieved Mon, 29 Apr 2024 23:58:22 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=48022, Retrieved Mon, 29 Apr 2024 23:58:22 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact110
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]
- RMPD        [Harrell-Davis Quantiles] [WS3_Part1_harrell] [2009-10-19 17:58:27] [7a6d96edf94be87996de99db5f42363b] [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'RServer@AstonUniversity' @ vre.aston.ac.uk

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=48022&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=48022&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=48022&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'RServer@AstonUniversity' @ vre.aston.ac.uk







Harrell-Davis Quantiles
quantilesvaluestandard error
0.0146.25653886969270.304261396219172
0.0649.45162873967290.809645216072289
0.1153.97636206084141.28104601141963
0.1658.12917600079951.64481531240375
0.2163.06055442557282.17040906810869
0.2667.13744085426661.20690048515007
0.3169.36014632588040.941217429797304
0.3671.5259774781741.18838981582641
0.4173.9865902057651.15601832858473
0.4675.79889112239440.695862879271246
0.5177.17046227438650.793998788921974
0.5678.94056598613330.96986831333195
0.6180.84347648983750.962579419834337
0.6683.00817417055861.13524719499450
0.7185.26405741809461.01081051214317
0.7687.33161689679780.8364665900989
0.8189.73246492604131.15360250519506
0.8693.5571853692251.76646523747657
0.91100.8890584384462.61624682092327
0.96111.7805646714223.01596430595947

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 46.2565388696927 & 0.304261396219172 \tabularnewline
0.06 & 49.4516287396729 & 0.809645216072289 \tabularnewline
0.11 & 53.9763620608414 & 1.28104601141963 \tabularnewline
0.16 & 58.1291760007995 & 1.64481531240375 \tabularnewline
0.21 & 63.0605544255728 & 2.17040906810869 \tabularnewline
0.26 & 67.1374408542666 & 1.20690048515007 \tabularnewline
0.31 & 69.3601463258804 & 0.941217429797304 \tabularnewline
0.36 & 71.525977478174 & 1.18838981582641 \tabularnewline
0.41 & 73.986590205765 & 1.15601832858473 \tabularnewline
0.46 & 75.7988911223944 & 0.695862879271246 \tabularnewline
0.51 & 77.1704622743865 & 0.793998788921974 \tabularnewline
0.56 & 78.9405659861333 & 0.96986831333195 \tabularnewline
0.61 & 80.8434764898375 & 0.962579419834337 \tabularnewline
0.66 & 83.0081741705586 & 1.13524719499450 \tabularnewline
0.71 & 85.2640574180946 & 1.01081051214317 \tabularnewline
0.76 & 87.3316168967978 & 0.8364665900989 \tabularnewline
0.81 & 89.7324649260413 & 1.15360250519506 \tabularnewline
0.86 & 93.557185369225 & 1.76646523747657 \tabularnewline
0.91 & 100.889058438446 & 2.61624682092327 \tabularnewline
0.96 & 111.780564671422 & 3.01596430595947 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=48022&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.01[/C][C]46.2565388696927[/C][C]0.304261396219172[/C][/ROW]
[ROW][C]0.06[/C][C]49.4516287396729[/C][C]0.809645216072289[/C][/ROW]
[ROW][C]0.11[/C][C]53.9763620608414[/C][C]1.28104601141963[/C][/ROW]
[ROW][C]0.16[/C][C]58.1291760007995[/C][C]1.64481531240375[/C][/ROW]
[ROW][C]0.21[/C][C]63.0605544255728[/C][C]2.17040906810869[/C][/ROW]
[ROW][C]0.26[/C][C]67.1374408542666[/C][C]1.20690048515007[/C][/ROW]
[ROW][C]0.31[/C][C]69.3601463258804[/C][C]0.941217429797304[/C][/ROW]
[ROW][C]0.36[/C][C]71.525977478174[/C][C]1.18838981582641[/C][/ROW]
[ROW][C]0.41[/C][C]73.986590205765[/C][C]1.15601832858473[/C][/ROW]
[ROW][C]0.46[/C][C]75.7988911223944[/C][C]0.695862879271246[/C][/ROW]
[ROW][C]0.51[/C][C]77.1704622743865[/C][C]0.793998788921974[/C][/ROW]
[ROW][C]0.56[/C][C]78.9405659861333[/C][C]0.96986831333195[/C][/ROW]
[ROW][C]0.61[/C][C]80.8434764898375[/C][C]0.962579419834337[/C][/ROW]
[ROW][C]0.66[/C][C]83.0081741705586[/C][C]1.13524719499450[/C][/ROW]
[ROW][C]0.71[/C][C]85.2640574180946[/C][C]1.01081051214317[/C][/ROW]
[ROW][C]0.76[/C][C]87.3316168967978[/C][C]0.8364665900989[/C][/ROW]
[ROW][C]0.81[/C][C]89.7324649260413[/C][C]1.15360250519506[/C][/ROW]
[ROW][C]0.86[/C][C]93.557185369225[/C][C]1.76646523747657[/C][/ROW]
[ROW][C]0.91[/C][C]100.889058438446[/C][C]2.61624682092327[/C][/ROW]
[ROW][C]0.96[/C][C]111.780564671422[/C][C]3.01596430595947[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=48022&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=48022&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.0146.25653886969270.304261396219172
0.0649.45162873967290.809645216072289
0.1153.97636206084141.28104601141963
0.1658.12917600079951.64481531240375
0.2163.06055442557282.17040906810869
0.2667.13744085426661.20690048515007
0.3169.36014632588040.941217429797304
0.3671.5259774781741.18838981582641
0.4173.9865902057651.15601832858473
0.4675.79889112239440.695862879271246
0.5177.17046227438650.793998788921974
0.5678.94056598613330.96986831333195
0.6180.84347648983750.962579419834337
0.6683.00817417055861.13524719499450
0.7185.26405741809461.01081051214317
0.7687.33161689679780.8364665900989
0.8189.73246492604131.15360250519506
0.8693.5571853692251.76646523747657
0.91100.8890584384462.61624682092327
0.96111.7805646714223.01596430595947



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