Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_Reddy-Moores Data Boxplot V2.0.wasp
Title produced by softwareBoxplot and Trimmed Means
Date of computationThu, 09 May 2013 05:34:11 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/May/09/t13680920658z29a2of3vb1hem.htm/, Retrieved Sat, 27 Apr 2024 16:43:39 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=208792, Retrieved Sat, 27 Apr 2024 16:43:39 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact114
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [T-Tests] [Salk Heatbeat Study] [2012-02-14 16:56:38] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMPD  [Histogram, QQplot and Density] [Salk Heatbeat Stu...] [2012-02-14 17:58:15] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMPD      [Boxplot and Trimmed Means] [Heart Beat Sample] [2013-05-09 09:34:11] [0d99488d7fbe7b1f0f6ad61929450c7c] [Current]
Feedback Forum

Post a new message
Dataseries X:
190	140
80	100
80	100
75	70
50	25
40	20
30	10
20	0
20	-10
10	-10
10	-25
10	-25
0	-25
0	-30
-10	-30
-25	-30
-30	-45
-45	-45
-60	-45
-85	-50
NA	-50
NA	-50
NA	-60
NA	-75
NA	-75
NA	-85
NA	-85
NA	-100
NA	-110
NA	-130
NA	-130
NA	-155
NA	-155
NA	-180
NA	-240
NA	-290




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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 & 1 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=208792&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]'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=208792&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=208792&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'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Engine error message
Error in if (par3 > 45) { : missing value where TRUE/FALSE needed
Execution halted



Parameters (Session):
par1 = grey ; par2 = TRUE ;
Parameters (R input):
par1 = grey ; par2 = TRUE ; par3 = ; par4 = ; par5 = ;
R code (references can be found in the software module):
bitmap(file='test1.png')
par2<- as.logical(par2) # notches
par3<-as.numeric(par3) # percentage 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.')}
par4 <- as.numeric(par4) #factor column
par5 <- as.numeric(par5) # response column
x <- t(x)
x1<-as.numeric(x[,par5]) # response
f1<-as.character(x[,par4]) # factor
x2<-x1[f1=='no']
f2 <- f1[f1=='no']
lotrm<-as.integer(length(x2)*par3/100)
hitrm<-as.integer(length(x2)*(100-par3)/100)
srt<-order(x2,f2)
trmx1<-x2[srt[lotrm:hitrm]]
trmf1<-f2[srt[lotrm:hitrm]]
x3<-x1[f1=='yes']
f3 <- f1[f1=='yes']
lotrm<-as.integer(length(x3)*par3/100)
hitrm<-as.integer(length(x3)*(100-par3)/100)
srt<-order(x3,f3)
trmx2<-x3[srt[lotrm:hitrm]]
trmf2<-f3[srt[lotrm:hitrm]]
xtrm<-c(trmx1,trmx2)
ftrm<-c(trmf1,trmf2)
xtrm[1:6]
ftrm[1:6]
r<-boxplot(xtrm~as.factor(ftrm), col=par1, notch=par2, main=main, xlab=xlab, ylab=ylab)
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,'Placement',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)
a<-table.row.start(a)
a<-table.element(a,'no',1,TRUE)
for (j in 1:5)
{
a<-table.element(a,r$stats[j,1])
}
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'yes',1,TRUE)
for (j in 1:5)
{
a<-table.element(a,r$stats[j,2])
}
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
tr.mns<-tapply(x1,f1,mean, trim=par3/100)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('trimmed_mean.htm','Trimmed Mean Equation','Trimmed Mean'),2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'No Placement')
a<-table.element(a,'Placement')
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,tr.mns[1])
a<-table.element(a,tr.mns[2])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')