Free Statistics

of Irreproducible Research!

Author's title

Q9: Compute a confidence interval for the random component of the model in ...

Author*Unverified author*
R Software Modulerwasp_harrell_davies.wasp
Title produced by softwareHarrell-Davis Quantiles
Date of computationSun, 04 Nov 2007 11:13:54 -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/Nov/04/mgrmhl019h42x5x1194200205.htm/, Retrieved Sun, 05 May 2024 15:08:07 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=334, Retrieved Sun, 05 May 2024 15:08:07 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsconfidence interval, harrell-davis quantiles, Q9
Estimated Impact195
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Q9: Compute a con...] [2007-11-04 18:13:54] [bd7b8d7754bcf95ad80b21f541dc6b78] [Current]
Feedback Forum

Post a new message
Dataseries X:
-5,20
-5,04
-4,79
-4,76
-4,61
-4,88
-4,57
-4,60
-4,33
-4,44
-4,52
-4,47
-4,04
-3,37
-3,09
-3,34
-3,61
-3,33
-3,14
-2,87
-2,59
-2,91
-2,73
-2,79
-2,50
-2,13
-1,99
-1,46
-1,07
-1,12
-0,78
-0,66
-0,59
-0,10
-0,22
-0,55
-0,27
0,41
1,06
1,29
1,44
1,75
2,37
2,47
2,59
2,40
2,37
2,35
2,38
2,83
2,79
3,30
3,69
3,67
4,03
4,16
3,86
3,65
3,91
4,05
4,10
4,67
4,68
5,16
5,03
5,02
5,47
5,34
5,44




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=334&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-5.2NA
0.025-5.023012203262250.152318025716172
0.05-4.832296869990890.136367583898141
0.075-4.695457005370250.129908713489573
0.1-4.578777256181720.155614363289524
0.125-4.443945140303920.244579987389161
0.15-4.257148246880860.373297369127306
0.175-4.01233526131570.47314210800336
0.2-3.737860499918480.497128157416909
0.225-3.471901530904820.460808783514607
0.25-3.233623987044050.418484939799904
0.275-3.015678729913060.415682769462612
0.3-2.795021362208050.468638350718464
0.325-2.547214448449870.566829635375835
0.35-2.257462593201900.67882909400745
0.375-1.926943851641100.76507445825357
0.4-1.571392313430030.800135451188219
0.425-1.210700246388850.791011510055775
0.45-0.8552970858566010.775237809930013
0.475-0.4991548846147590.79529792830234
0.5-0.1251661703479710.863680643069782
0.5250.2804208689147760.94951900563004
0.550.7145001550906131.00019570895379
0.5751.154015732210010.97527094060258
0.61.565871918577620.870300807804545
0.6251.924146577886010.7204413731757
0.652.224204414294950.583692229395327
0.6752.484506100316320.510408644055965
0.72.733874841899780.507184376096636
0.7252.991927253037790.528665023106946
0.753.256790076535290.520275493865977
0.7753.510119392844520.461966451949402
0.83.735862278055090.379021840749663
0.8253.938044655535590.321409852453113
0.854.143235215869520.326012323503082
0.8754.381263791744260.370280770393109
0.94.654372210838190.37993743445695
0.9254.93090586291730.317281950180866
0.955.180439812152340.229962245314860
0.9755.378991897988630.130896578643583
15.47NA

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0 & -5.2 & NA \tabularnewline
0.025 & -5.02301220326225 & 0.152318025716172 \tabularnewline
0.05 & -4.83229686999089 & 0.136367583898141 \tabularnewline
0.075 & -4.69545700537025 & 0.129908713489573 \tabularnewline
0.1 & -4.57877725618172 & 0.155614363289524 \tabularnewline
0.125 & -4.44394514030392 & 0.244579987389161 \tabularnewline
0.15 & -4.25714824688086 & 0.373297369127306 \tabularnewline
0.175 & -4.0123352613157 & 0.47314210800336 \tabularnewline
0.2 & -3.73786049991848 & 0.497128157416909 \tabularnewline
0.225 & -3.47190153090482 & 0.460808783514607 \tabularnewline
0.25 & -3.23362398704405 & 0.418484939799904 \tabularnewline
0.275 & -3.01567872991306 & 0.415682769462612 \tabularnewline
0.3 & -2.79502136220805 & 0.468638350718464 \tabularnewline
0.325 & -2.54721444844987 & 0.566829635375835 \tabularnewline
0.35 & -2.25746259320190 & 0.67882909400745 \tabularnewline
0.375 & -1.92694385164110 & 0.76507445825357 \tabularnewline
0.4 & -1.57139231343003 & 0.800135451188219 \tabularnewline
0.425 & -1.21070024638885 & 0.791011510055775 \tabularnewline
0.45 & -0.855297085856601 & 0.775237809930013 \tabularnewline
0.475 & -0.499154884614759 & 0.79529792830234 \tabularnewline
0.5 & -0.125166170347971 & 0.863680643069782 \tabularnewline
0.525 & 0.280420868914776 & 0.94951900563004 \tabularnewline
0.55 & 0.714500155090613 & 1.00019570895379 \tabularnewline
0.575 & 1.15401573221001 & 0.97527094060258 \tabularnewline
0.6 & 1.56587191857762 & 0.870300807804545 \tabularnewline
0.625 & 1.92414657788601 & 0.7204413731757 \tabularnewline
0.65 & 2.22420441429495 & 0.583692229395327 \tabularnewline
0.675 & 2.48450610031632 & 0.510408644055965 \tabularnewline
0.7 & 2.73387484189978 & 0.507184376096636 \tabularnewline
0.725 & 2.99192725303779 & 0.528665023106946 \tabularnewline
0.75 & 3.25679007653529 & 0.520275493865977 \tabularnewline
0.775 & 3.51011939284452 & 0.461966451949402 \tabularnewline
0.8 & 3.73586227805509 & 0.379021840749663 \tabularnewline
0.825 & 3.93804465553559 & 0.321409852453113 \tabularnewline
0.85 & 4.14323521586952 & 0.326012323503082 \tabularnewline
0.875 & 4.38126379174426 & 0.370280770393109 \tabularnewline
0.9 & 4.65437221083819 & 0.37993743445695 \tabularnewline
0.925 & 4.9309058629173 & 0.317281950180866 \tabularnewline
0.95 & 5.18043981215234 & 0.229962245314860 \tabularnewline
0.975 & 5.37899189798863 & 0.130896578643583 \tabularnewline
1 & 5.47 & NA \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=334&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[/C][C]-5.2[/C][C]NA[/C][/ROW]
[ROW][C]0.025[/C][C]-5.02301220326225[/C][C]0.152318025716172[/C][/ROW]
[ROW][C]0.05[/C][C]-4.83229686999089[/C][C]0.136367583898141[/C][/ROW]
[ROW][C]0.075[/C][C]-4.69545700537025[/C][C]0.129908713489573[/C][/ROW]
[ROW][C]0.1[/C][C]-4.57877725618172[/C][C]0.155614363289524[/C][/ROW]
[ROW][C]0.125[/C][C]-4.44394514030392[/C][C]0.244579987389161[/C][/ROW]
[ROW][C]0.15[/C][C]-4.25714824688086[/C][C]0.373297369127306[/C][/ROW]
[ROW][C]0.175[/C][C]-4.0123352613157[/C][C]0.47314210800336[/C][/ROW]
[ROW][C]0.2[/C][C]-3.73786049991848[/C][C]0.497128157416909[/C][/ROW]
[ROW][C]0.225[/C][C]-3.47190153090482[/C][C]0.460808783514607[/C][/ROW]
[ROW][C]0.25[/C][C]-3.23362398704405[/C][C]0.418484939799904[/C][/ROW]
[ROW][C]0.275[/C][C]-3.01567872991306[/C][C]0.415682769462612[/C][/ROW]
[ROW][C]0.3[/C][C]-2.79502136220805[/C][C]0.468638350718464[/C][/ROW]
[ROW][C]0.325[/C][C]-2.54721444844987[/C][C]0.566829635375835[/C][/ROW]
[ROW][C]0.35[/C][C]-2.25746259320190[/C][C]0.67882909400745[/C][/ROW]
[ROW][C]0.375[/C][C]-1.92694385164110[/C][C]0.76507445825357[/C][/ROW]
[ROW][C]0.4[/C][C]-1.57139231343003[/C][C]0.800135451188219[/C][/ROW]
[ROW][C]0.425[/C][C]-1.21070024638885[/C][C]0.791011510055775[/C][/ROW]
[ROW][C]0.45[/C][C]-0.855297085856601[/C][C]0.775237809930013[/C][/ROW]
[ROW][C]0.475[/C][C]-0.499154884614759[/C][C]0.79529792830234[/C][/ROW]
[ROW][C]0.5[/C][C]-0.125166170347971[/C][C]0.863680643069782[/C][/ROW]
[ROW][C]0.525[/C][C]0.280420868914776[/C][C]0.94951900563004[/C][/ROW]
[ROW][C]0.55[/C][C]0.714500155090613[/C][C]1.00019570895379[/C][/ROW]
[ROW][C]0.575[/C][C]1.15401573221001[/C][C]0.97527094060258[/C][/ROW]
[ROW][C]0.6[/C][C]1.56587191857762[/C][C]0.870300807804545[/C][/ROW]
[ROW][C]0.625[/C][C]1.92414657788601[/C][C]0.7204413731757[/C][/ROW]
[ROW][C]0.65[/C][C]2.22420441429495[/C][C]0.583692229395327[/C][/ROW]
[ROW][C]0.675[/C][C]2.48450610031632[/C][C]0.510408644055965[/C][/ROW]
[ROW][C]0.7[/C][C]2.73387484189978[/C][C]0.507184376096636[/C][/ROW]
[ROW][C]0.725[/C][C]2.99192725303779[/C][C]0.528665023106946[/C][/ROW]
[ROW][C]0.75[/C][C]3.25679007653529[/C][C]0.520275493865977[/C][/ROW]
[ROW][C]0.775[/C][C]3.51011939284452[/C][C]0.461966451949402[/C][/ROW]
[ROW][C]0.8[/C][C]3.73586227805509[/C][C]0.379021840749663[/C][/ROW]
[ROW][C]0.825[/C][C]3.93804465553559[/C][C]0.321409852453113[/C][/ROW]
[ROW][C]0.85[/C][C]4.14323521586952[/C][C]0.326012323503082[/C][/ROW]
[ROW][C]0.875[/C][C]4.38126379174426[/C][C]0.370280770393109[/C][/ROW]
[ROW][C]0.9[/C][C]4.65437221083819[/C][C]0.37993743445695[/C][/ROW]
[ROW][C]0.925[/C][C]4.9309058629173[/C][C]0.317281950180866[/C][/ROW]
[ROW][C]0.95[/C][C]5.18043981215234[/C][C]0.229962245314860[/C][/ROW]
[ROW][C]0.975[/C][C]5.37899189798863[/C][C]0.130896578643583[/C][/ROW]
[ROW][C]1[/C][C]5.47[/C][C]NA[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=334&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=334&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-5.2NA
0.025-5.023012203262250.152318025716172
0.05-4.832296869990890.136367583898141
0.075-4.695457005370250.129908713489573
0.1-4.578777256181720.155614363289524
0.125-4.443945140303920.244579987389161
0.15-4.257148246880860.373297369127306
0.175-4.01233526131570.47314210800336
0.2-3.737860499918480.497128157416909
0.225-3.471901530904820.460808783514607
0.25-3.233623987044050.418484939799904
0.275-3.015678729913060.415682769462612
0.3-2.795021362208050.468638350718464
0.325-2.547214448449870.566829635375835
0.35-2.257462593201900.67882909400745
0.375-1.926943851641100.76507445825357
0.4-1.571392313430030.800135451188219
0.425-1.210700246388850.791011510055775
0.45-0.8552970858566010.775237809930013
0.475-0.4991548846147590.79529792830234
0.5-0.1251661703479710.863680643069782
0.5250.2804208689147760.94951900563004
0.550.7145001550906131.00019570895379
0.5751.154015732210010.97527094060258
0.61.565871918577620.870300807804545
0.6251.924146577886010.7204413731757
0.652.224204414294950.583692229395327
0.6752.484506100316320.510408644055965
0.72.733874841899780.507184376096636
0.7252.991927253037790.528665023106946
0.753.256790076535290.520275493865977
0.7753.510119392844520.461966451949402
0.83.735862278055090.379021840749663
0.8253.938044655535590.321409852453113
0.854.143235215869520.326012323503082
0.8754.381263791744260.370280770393109
0.94.654372210838190.37993743445695
0.9254.93090586291730.317281950180866
0.955.180439812152340.229962245314860
0.9755.378991897988630.130896578643583
15.47NA



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