FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_CARE Data Boxplot.wasp
Title produced by softwareCARE Data - Boxplots and Scatterplot Matrix
Date of computationTue, 27 Nov 2012 10:45:30 -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/Nov/27/t1354031155qyvehm6yf1n2e3c.htm/, Retrieved Sat, 04 May 2024 07:24:59 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=193938, Retrieved Sat, 04 May 2024 07:24:59 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact109

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [CARE Data - Boxplots and Scatterplot Matrix] [IQ on Average Age] [2012-11-27 15:45:30] [01d189c67e410f55b92337d59f7551a8] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
111	52
102	55
108	80
109	45
118	60
79	1
88	45
102	68
105	26
92	70
11	85
104	54
83	55
84	40
85	55
110	50
121	71
120	55
100	70
94	55
89	60
93	65
128	66
84	55
127	90
106	55
129	60
82	1
106	55
109	26
91	14
111	45
105	1
118	65
103	1
101	60
101	60
95	60
108	65
95	45
98	20
82	50
100	60
100	48
107	40
95	55
97	54
93	40
81	40
89	1
111	60
95	1
106	75
83	24
81	1
115	80
112	60
92	46
85	1
95	60
115	75
91	54
107	78
102	20
86	45
96	60
114	70
105	1
82	20
120	60
88	20
90	50
85	50
106	75
109	70
75	20
91	45
96	20
108	50
86	55
98	15
99	26
95	25
88	1
111	60
103	40
107	40
118	50

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Engine error message
Error in if (par3 > 45) { : missing value where TRUE/FALSE needed
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 & 0 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
R Engine error message & 
Error in if (par3 > 45) { : missing value where TRUE/FALSE needed
Execution halted

 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=193938&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in if (par3 > 45) { : missing value where TRUE/FALSE needed
Execution halted

[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=193938&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=193938&T=0

As an alternative you can also use a QR Code:  
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 time0 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Engine error message
Error in if (par3 > 45) { : missing value where TRUE/FALSE needed
Execution halted


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 3 ; par2 = TRUE ; par3 = 10% ;
Parameters (R input):
par1 = 3 ; par2 = TRUE ; par3 = 10% ;
R code (references can be found in the software module):
par1 <- as.numeric(par1) #colour
par2<- as.logical(par2) # Notches ?
par3<-as.numeric(par3) # % trim
if(par3>45){par3<-45;warning('trim limited to 45%')}
if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')}
lotrm<-as.integer(length(y[1,])*par3/100)+1
hitrm<-as.integer(length(y[1,])*(100-par3)/100)
y1<-array(dim=c(dim(y)[1], hitrm-lotrm+1), dimnames=list(dimnames(y)[[1]], 1:(hitrm-lotrm+1) ))
for(i in 1:dim(y)[1]){
tmp<-order(y[i,])
y1[i,]<- y[i, tmp[lotrm:hitrm] ]
}
bitmap(file='test2.png')
pairs(t(y))
dev.off()
y<-y1
z <- as.data.frame(t(y))
bitmap(file='test1.png')
(r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=par2,col=par1))
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('overview.htm','Boxplot statistics','Boxplot overview'),6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,hyperlink('lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE)
a<-table.element(a,hyperlink('lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE)
a<-table.element(a,hyperlink('central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE)
a<-table.element(a,hyperlink('upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE)
a<-table.element(a,hyperlink('upper_whisker.htm','upper whisker','definition of upper whisker'),1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
for (j in 1:5)
{
a<-table.element(a,signif(r$stats[j,i], digits=4))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
if (par2){
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Notches',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,'lower bound',1,TRUE)
a<-table.element(a,'median',1,TRUE)
a<-table.element(a,'upper bound',1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
a<-table.element(a,signif(r$conf[1,i], digits=4))
a<-table.element(a, signif(r$stats[3,i], digits=4))
a<-table.element(a,signif(r$conf[2,i], digits=4))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
}
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Means',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,hyperlink('trimmed_mean.htm','trimmed mean','definition of trimmed mean'),1,TRUE)
a<-table.element(a,hyperlink('unbiased1.htm','unbiased SD','definition of unbiased SD'),1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
a<-table.element(a,signif(mean(z[i], trim=par3/100, na.rm=TRUE), digits=4))
a<-table.element(a,signif(sd(z[i], na.rm=TRUE), digits=4))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')