Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_BoxplotAndScatterplotMatrix.wasp
Title produced by softwareBoxplot and Trimmed Means
Date of computationThu, 23 Oct 2014 16:42:38 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Oct/23/t1414078979gymp99owog4ogl9.htm/, Retrieved Fri, 10 May 2024 14:14:22 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=245907, Retrieved Fri, 10 May 2024 14:14:22 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact65
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Boxplot and Trimmed Means] [Care Age 10 Data] [2009-10-26 09:01:50] [98fd0e87c3eb04e0cc2efde01dbafab6]
-   PD  [Boxplot and Trimmed Means] [Care Age 7 Data] [2009-10-26 18:36:29] [98fd0e87c3eb04e0cc2efde01dbafab6]
-   P     [CARE Data - Boxplots and Scatterplot Matrix] [CARE Data] [2010-10-19 14:16:27] [3fdd735c61ad38cbc9b3393dc997cdb7]
- RM        [CARE Data - Boxplots and Scatterplot Matrix] [CARE data - works...] [2011-10-17 10:23:12] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMP         [Boxplot and Trimmed Means] [CARE Study Age 7 ] [2013-10-17 12:59:45] [34296d8f7657c52ed60d5bff9133afec]
- RMPD            [Boxplot and Trimmed Means] [Year 10 data boxplot] [2014-10-23 15:42:38] [0798657923d04bfb89a94c1a1eece86b] [Current]
Feedback Forum

Post a new message
Dataseries X:
WJ10AFS	WJ10AVA	WJ10ARD	WJ10AMA	WJ10AKN			
88	95	95	81	84			
98	92	87	90	97			
89	95	91	91	88			
81	97	84	87	87			
88	92	92	89	92			
95	104	84	85	91			These 5 variables have been abstracted from the complete dataset for 
79	79	84	81	80			The Carolina Abecedarian Project and
91	78	69	80	73			the Carolina Approach to
80	79	84	75	73			Responsive Education (CARE),
91	86	84	84	77			1972-1992
85	92	85	91	86			
94	78	89	87	82			Bibliographic Citation: Ramey, Craig T., James J. Gallagher, Frances A. Campbell, Barbara
86	86	83	90	82			H. Wasik, and Joseph J. Sparling. CAROLINA ABECEDARIAN
82	71	86	72	73			PROJECT AND THE CAROLINA APPROACH TO RESPONSIVE
107	90	91	96	99			EDUCATION (CARE), 1972-1992 [Computer file]. ICPSR version.
97	93	90	88	90			Chapel Hill, NC: University of North Carolina [producer], 2002. Ann
78	80	77	73	78			Arbor, MI: Inter-university Consortium for Political and Social Research
91	90	88	86	100			[distributor], 2004.
74	77	88	67	73			
95	86	88	89	95			CAUTION:
97	89	102	94	93			For study use missing values have been removed. This excel file 
112	88	103	88	91			is not for use other than for educational study purposes.
80	86	84	77	81			
98	97	87	89	98			The full datset is available to interested professionals from ICPSR
102	81	88	104	88			Please register and retrieve the origianl dataset from
100	88	93	89	90			the Inter-UnivesrityConsortium on Political and Social research
97	82	101	107	85			http://www.icpsr.umich.edu/icpsrweb/ICPSR/
92	82	88	96	79			
114	100	113	111	107			
88	91	95	84	82			
96	111	99	83	98			
91	102	85	82	84			
114	105	105	105	99			
91	92	105	100	89			
114	100	96	111	105			
89	82	88	90	94			
103	91	93	87	92			
83	71	93	90	81			
74	71	69	67	NA			
88	89	83	84	86			
88	82	83	84	90			
98	105	93	96	101			
74	71	69	67	73			
86	91	86	90	88			
93	105	90	83	90			
85	82	87	78	84			
91	86	92	100	89			
75	72	75	67	73			
97	98	89	94	98			
88	80	88	89	82			
79	84	87	89	77			
106	109	113	111	107			
102	80	102	88	84			
98	112	101	94	93			
98	86	87	111	96			
114	112	105	111	107			
95	90	101	105	95			
98	96	79	97	89			
82	84	86	83	82			
97	94	95	105	92			
94	112	104	69	100			
84	81	79	80	75			
88	105	86	87	94			
83	77	81	67	83			
98	105	92	92	92			
106	112	95	102	102			
94	91	97	90	98			
91	92	87	98	94			
89	72	94	94	83			
102	100	105	101	107			
103	103	94	96	107			
96	89	88	94	91			
89	86	105	105	83			
112	112	113	105	100			
74	71	69	84	77			
74	73	69	72	74			
102	99	100	96	92			
106	85	113	94	100			
91	80	88	107	92			
94	88	104	103	93			
98	95	94	95	89			
85	86	91	73	86			
108	92	103	110	97			
90	96	95	84	82			
74	71	69	67	73			
81	88	86	73	93			
104	98	108	84	97			
90	108	83	79	93			
95	103	102	91	89			
85	81	94	92	86			
114	96	113	107	89			
114	112	113	111	107			
93	95	93	89	90			
104	92	105	110	104			
92	103	88	98	96			
100	95	90	95	93			
92	94	92	84	104			
77	79	81	73	83			
86	91	86	86	92			
87	87	85	98	83			
91	82	86	90	83			
102	97	97	97	96			
100	92	103	106	84			
113	112	102	106	107			
87	72	78	85	82
96	90	96	80	95
80	87	84	77	85
85	89	87	67	88
89	110	105	98	98
89	95	89	89	87
73	72	78	63	71
88	70	96	97	89
86	79	79	77	86
83	92	90	96	88
91	93	96	84	93
98	97	95	94	95
80	75	79	70	78
84	88	86	100	84
90	105	99	88	97
92	82	111	97	95
86	89	105	78	88
117	114	123	108	121
88	81	102	81	97
86	73	82	96	82
78	82	86	84	82
104	110	117	104	108
89	93	94	87	90
99	114	110	98	113
84	91	86	89	86
89	97	109	93	105
95	97	118	89	105
58	68	64	63	60
77	72	78	84	75
77	81	86	82	79
98	102	98	105	104
58	62	64	70	60
81	102	84	94	87
92	84	87	101	93
94	75	88	83	85
87	102	106	88	97
75	62	78	69	67
80	100	79	94	85
84	82	95	82	86
87	78	98	82	88
114	106	101	109	114
97	88	104	86	97
102	85	91	99	93
80	97	78	82	80
112	114	123	114	118
104	105	109	99	109
75	100	77	91	81
82	72	77	81	79
92	106	102	102	104
102	78	88	101	97
73	80	81	72	75
78	88	83	91	82
74	62	72	83	71
92	93	95	102	92
97	87	109	103	104
92	109	91	93	99
91	91	87	88	89
102	110	112	110	114
89	98	95	100	95
98	99	103	88	106
117	95	108	119	117
92	97	116	98	109
70	88	71	90	74
63	72	70	69	63
91	102	105	98	103
106	95	112	98	107
81	93	86	81	86
98	90	109	93	99
88	80	112	91	87
95	75	97	86	93
105	109	122	95	114
90	75	97	87	82
58	72	64	63	60
86	77	84	89	83
112	93	110	108	111
71	72	70	99	75
98	91	100	90	94
87	95	102	91	91
98	111	108	82	100
117	105	123	119	121
98	90	78	97	90
98	105	109	108	105
90	90	106	101	97
84	92	95	96	92
86	95	108	97	92
68	75	74	69	67
86	92	86	89	90
77	100	82	84	85
83	82	88	79	85
90	102	96	103	97
107	104	114	92	118
102	98	111	104	109
90	95	88	102	94
113	114	123	106	121
86	93	89	93	85
105	110	114	108	112
110	104	121	119	118
82	91	100	102	93
89	79	83	83	85
99	97	92	92	97
113	102	108	119	121
91	77	87	100	83
92	105	84	97	89
86	95	95	103	92
91	104	90	95	89
82	75	82	73	82
86	82	84	86	85
86	91	96	88	92
102	105	109	102	112
84	71	84	87	79
80	89	83	100	84
96	79	93	88	92
117	113	120	100	118
88	82	83	79	82
99	102	103	90	100
89	82	103	83	87
86	92	102	95	97
98	82	90	89	87
58	86	69	64	70
98	108	114	108	110
58	74	64	80	60
81	87	87	83	81
90	97	116	96	94
112	102	121	108	114
84	75	90	86	85
84	82	92	88	85
75	72	81	97	79
78	88	84	110	82
94	97	106	102	102
75	85	75	92	73
82	62	79	91	78
87	95	80	92	84
59	95	74	72	70
68	72	66	88	66
80	92	81	99	79
94	91	98	95	95
87	97	99	83	95
70	78	68	68	69
113	95	100	101	109
99	91	97	91	100
90	95	99	89	91
102	88	104	98	99
80	92	82	89	84
98	97	92	98	101
98	104	99	97	99
95	113	109	98	104
104	113	106	113	109
86	82	89	83	87
84	90	87	92	88
85	91	98	113	87
90	95	95	103	95
86	102	92	96	92
97	92	90	94	100




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net
R Framework error message
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.

\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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
R Framework error message & 
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=245907&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[ROW][C]R Framework error message[/C][C]
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=245907&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=245907&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 time0 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net
R Framework error message
Warning: there are blank lines in the 'Data X' field.
Please, use NA for missing data - blank lines are simply
 deleted and are NOT treated as missing values.



Parameters (Session):
par1 = 2 ; par2 = FALSE ; par3 = 10 ;
Parameters (R input):
par1 = 2 ; par2 = FALSE ; 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,round(r$stats[j,i], digits=2))
}
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,round(r$conf[1,i], digits=2))
a<-table.element(a, round(r$stats[3,i], digits=2))
a<-table.element(a,round(r$conf[2,i], digits=2))
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,round(mean(unlist(z[i]), trim=par3/100, na.rm=TRUE), digits=2))
a<-table.element(a,round(sd(unlist(z[i]), na.rm=TRUE), digits=2))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')