Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_harrell_davies.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationFri, 19 Oct 2007 06:57:17 -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/2007/Oct/19/cx2782qti5qtx5e1192802027.htm/, Retrieved Wed, 08 May 2024 13:20:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=1020, Retrieved Wed, 08 May 2024 13:20:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsQ5: Show that the 95% confidence interval of the prediction error of the model in Question 3 (assuming that the model would be valid) is equal to [-0.156036650732365, 0.180290511044764].
Estimated Impact216
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Q5 Workshop 2b] [2007-10-19 13:57:17] [dffb4ee60db878ec13fafd952698fec9] [Current]
Feedback Forum

Post a new message
Dataseries X:
109,20
88,60
94,30
98,30
86,40
80,60
104,10
108,20
93,40
71,90
94,10
94,90
96,40
91,10
84,40
86,40
88,00
75,10
109,70
103,00
82,10
68,00
96,40
94,30
90,00
88,00
76,10
82,50
81,40
66,50
97,20
94,10
80,70
70,50
87,80
89,50
99,60
84,20
75,10
92,00
80,80
73,10
99,80
90,00
83,10
72,40
78,80
87,30
91,00
80,10
73,60
86,40
74,50
71,20
92,40
81,50
85,30
69,90
84,20
90,70
100,30




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1020&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.0166.93552886464051.50511131042239
0.0670.38872235001741.31575761513653
0.1172.59817503509491.44826178181682
0.1674.74472999439272.05984925490576
0.2177.28390767866032.56568833999510
0.2679.69653030221722.11444854425374
0.3181.48578678140711.64332585943454
0.3682.99094157370961.69967575099472
0.4184.5093959633421.77327081904245
0.4685.99395203514081.68549250903118
0.5187.37696330829751.57220071469155
0.5688.70511327533851.55014908485959
0.6190.05489991274721.57676934616907
0.6691.47275198710961.62458617851067
0.7192.9390938270361.56589672406022
0.7694.41566821443721.49251938093099
0.8196.1407734522571.75730465039923
0.8698.42876858585792.04037373562957
0.91101.8103874647642.6930702561356
0.96107.1686316139792.78944884233111

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 66.9355288646405 & 1.50511131042239 \tabularnewline
0.06 & 70.3887223500174 & 1.31575761513653 \tabularnewline
0.11 & 72.5981750350949 & 1.44826178181682 \tabularnewline
0.16 & 74.7447299943927 & 2.05984925490576 \tabularnewline
0.21 & 77.2839076786603 & 2.56568833999510 \tabularnewline
0.26 & 79.6965303022172 & 2.11444854425374 \tabularnewline
0.31 & 81.4857867814071 & 1.64332585943454 \tabularnewline
0.36 & 82.9909415737096 & 1.69967575099472 \tabularnewline
0.41 & 84.509395963342 & 1.77327081904245 \tabularnewline
0.46 & 85.9939520351408 & 1.68549250903118 \tabularnewline
0.51 & 87.3769633082975 & 1.57220071469155 \tabularnewline
0.56 & 88.7051132753385 & 1.55014908485959 \tabularnewline
0.61 & 90.0548999127472 & 1.57676934616907 \tabularnewline
0.66 & 91.4727519871096 & 1.62458617851067 \tabularnewline
0.71 & 92.939093827036 & 1.56589672406022 \tabularnewline
0.76 & 94.4156682144372 & 1.49251938093099 \tabularnewline
0.81 & 96.140773452257 & 1.75730465039923 \tabularnewline
0.86 & 98.4287685858579 & 2.04037373562957 \tabularnewline
0.91 & 101.810387464764 & 2.6930702561356 \tabularnewline
0.96 & 107.168631613979 & 2.78944884233111 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=1020&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]66.9355288646405[/C][C]1.50511131042239[/C][/ROW]
[ROW][C]0.06[/C][C]70.3887223500174[/C][C]1.31575761513653[/C][/ROW]
[ROW][C]0.11[/C][C]72.5981750350949[/C][C]1.44826178181682[/C][/ROW]
[ROW][C]0.16[/C][C]74.7447299943927[/C][C]2.05984925490576[/C][/ROW]
[ROW][C]0.21[/C][C]77.2839076786603[/C][C]2.56568833999510[/C][/ROW]
[ROW][C]0.26[/C][C]79.6965303022172[/C][C]2.11444854425374[/C][/ROW]
[ROW][C]0.31[/C][C]81.4857867814071[/C][C]1.64332585943454[/C][/ROW]
[ROW][C]0.36[/C][C]82.9909415737096[/C][C]1.69967575099472[/C][/ROW]
[ROW][C]0.41[/C][C]84.509395963342[/C][C]1.77327081904245[/C][/ROW]
[ROW][C]0.46[/C][C]85.9939520351408[/C][C]1.68549250903118[/C][/ROW]
[ROW][C]0.51[/C][C]87.3769633082975[/C][C]1.57220071469155[/C][/ROW]
[ROW][C]0.56[/C][C]88.7051132753385[/C][C]1.55014908485959[/C][/ROW]
[ROW][C]0.61[/C][C]90.0548999127472[/C][C]1.57676934616907[/C][/ROW]
[ROW][C]0.66[/C][C]91.4727519871096[/C][C]1.62458617851067[/C][/ROW]
[ROW][C]0.71[/C][C]92.939093827036[/C][C]1.56589672406022[/C][/ROW]
[ROW][C]0.76[/C][C]94.4156682144372[/C][C]1.49251938093099[/C][/ROW]
[ROW][C]0.81[/C][C]96.140773452257[/C][C]1.75730465039923[/C][/ROW]
[ROW][C]0.86[/C][C]98.4287685858579[/C][C]2.04037373562957[/C][/ROW]
[ROW][C]0.91[/C][C]101.810387464764[/C][C]2.6930702561356[/C][/ROW]
[ROW][C]0.96[/C][C]107.168631613979[/C][C]2.78944884233111[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=1020&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=1020&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.0166.93552886464051.50511131042239
0.0670.38872235001741.31575761513653
0.1172.59817503509491.44826178181682
0.1674.74472999439272.05984925490576
0.2177.28390767866032.56568833999510
0.2679.69653030221722.11444854425374
0.3181.48578678140711.64332585943454
0.3682.99094157370961.69967575099472
0.4184.5093959633421.77327081904245
0.4685.99395203514081.68549250903118
0.5187.37696330829751.57220071469155
0.5688.70511327533851.55014908485959
0.6190.05489991274721.57676934616907
0.6691.47275198710961.62458617851067
0.7192.9390938270361.56589672406022
0.7694.41566821443721.49251938093099
0.8196.1407734522571.75730465039923
0.8698.42876858585792.04037373562957
0.91101.8103874647642.6930702561356
0.96107.1686316139792.78944884233111



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