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 computationThu, 12 Mar 2009 13:49:45 -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/Mar/12/t1236887477ey76pcrbq1e8b94.htm/, Retrieved Thu, 02 May 2024 06:20:11 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=39012, Retrieved Thu, 02 May 2024 06:20:11 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact116
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [aantal werkzoeken...] [2009-03-12 19:49:45] [d19dedf3d9cbe742e7fe81d57ca6286d] [Current]
Feedback Forum

Post a new message
Dataseries X:
76284
74250
71362
72008
73454
76207
75091
67610
66466
69848
70483
73302
74988
74291
72932
75576
77304
81329
79357
72187
75728
77828
78665
81640
83614
82719
83324
86625
91790
92134
90951
84957
84831
84816
85220
86629
87184
86858
86285
89664
92194
92729
91378
84915
85159
86587
88521
89488
89587
89955
89528
92159
94118
93571
92603
87175
87032
88219
90390
90269
90398
89168
87888
91301
92472
91995
90580
82880
82387
83903
83823
84757
84886




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=39012&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]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=39012&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39012&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'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Harrell-Davis Quantiles
quantilesvaluestandard error
0.0166962.14044364181341.70403517596
0.0267787.7547261711594.10990939348
0.0368709.02216712841649.67729699965
0.0469567.60896152221533.16028794887
0.0570303.29456525631372.68120589042
0.0670916.36410609351242.85933594668
0.0771431.09826815871155.90306751584
0.0871874.63919958841102.11199762333
0.0972268.79905451131071.70081836059

\begin{tabular}{lllllllll}
\hline
Harrell-Davis Quantiles \tabularnewline
quantiles & value & standard error \tabularnewline
0.01 & 66962.1404436418 & 1341.70403517596 \tabularnewline
0.02 & 67787.754726171 & 1594.10990939348 \tabularnewline
0.03 & 68709.0221671284 & 1649.67729699965 \tabularnewline
0.04 & 69567.6089615222 & 1533.16028794887 \tabularnewline
0.05 & 70303.2945652563 & 1372.68120589042 \tabularnewline
0.06 & 70916.3641060935 & 1242.85933594668 \tabularnewline
0.07 & 71431.0982681587 & 1155.90306751584 \tabularnewline
0.08 & 71874.6391995884 & 1102.11199762333 \tabularnewline
0.09 & 72268.7990545113 & 1071.70081836059 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=39012&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]66962.1404436418[/C][C]1341.70403517596[/C][/ROW]
[ROW][C]0.02[/C][C]67787.754726171[/C][C]1594.10990939348[/C][/ROW]
[ROW][C]0.03[/C][C]68709.0221671284[/C][C]1649.67729699965[/C][/ROW]
[ROW][C]0.04[/C][C]69567.6089615222[/C][C]1533.16028794887[/C][/ROW]
[ROW][C]0.05[/C][C]70303.2945652563[/C][C]1372.68120589042[/C][/ROW]
[ROW][C]0.06[/C][C]70916.3641060935[/C][C]1242.85933594668[/C][/ROW]
[ROW][C]0.07[/C][C]71431.0982681587[/C][C]1155.90306751584[/C][/ROW]
[ROW][C]0.08[/C][C]71874.6391995884[/C][C]1102.11199762333[/C][/ROW]
[ROW][C]0.09[/C][C]72268.7990545113[/C][C]1071.70081836059[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=39012&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=39012&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.0166962.14044364181341.70403517596
0.0267787.7547261711594.10990939348
0.0368709.02216712841649.67729699965
0.0469567.60896152221533.16028794887
0.0570303.29456525631372.68120589042
0.0670916.36410609351242.85933594668
0.0771431.09826815871155.90306751584
0.0871874.63919958841102.11199762333
0.0972268.79905451131071.70081836059



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