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 computationWed, 22 Feb 2012 17:02:34 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Feb/22/t1329948175t8b6smbhnccg21f.htm/, Retrieved Mon, 06 May 2024 16:31:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=163235, Retrieved Mon, 06 May 2024 16:31:26 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact73
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Harrell-Davis Quantiles] [Decielen gemiddel...] [2012-02-22 22:02:34] [3b1daca0bbde728cac84422b4a85488f] [Current]
Feedback Forum

Post a new message
Dataseries X:
62.11
62.15
62.2
62.22
62.02
62.02
62.02
62.07
62.31
62.71
62.77
62.82
62.82
62.82
62.55
62.6
62.47
62.47
62.47
62.72
63.13
64.09
64.31
64.29
64.29
64.29
64.35
64.42
64.24
64.23
64.23
64.2
65.35
65.83
66.15
66.19
66.19
66.56
66.59
66.48
66.4
66.4
66.4
66.49
66.65
67.69
67.91
68.14
68.14
68.16
67.94
68
68.1
68.12
68.12
68.24
68.42
68.97
69.13
69.2




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'AstonUniversity' @ aston.wessa.net
R Engine error message
Error in if (del == 0 && to == 0) return(to) : 
  missing value where TRUE/FALSE needed
Calls: seq -> seq.default
Execution halted

\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 & 1 seconds \tabularnewline
R Server & 'AstonUniversity' @ aston.wessa.net \tabularnewline
R Engine error message & 
Error in if (del == 0 && to == 0) return(to) : 
  missing value where TRUE/FALSE needed
Calls: seq -> seq.default
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=163235&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'AstonUniversity' @ aston.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in if (del == 0 && to == 0) return(to) : 
  missing value where TRUE/FALSE needed
Calls: seq -> seq.default
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=163235&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=163235&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 time1 seconds
R Server'AstonUniversity' @ aston.wessa.net
R Engine error message
Error in if (del == 0 && to == 0) return(to) : 
  missing value where TRUE/FALSE needed
Calls: seq -> seq.default
Execution halted



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