Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_boxcoxnorm.wasp
Title produced by softwareBox-Cox Normality Plot
Date of computationSun, 04 Nov 2007 04:48:04 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Nov/04/bc2i8mg0ej6voj11194176810.htm/, Retrieved Sun, 05 May 2024 10:20:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=293, Retrieved Sun, 05 May 2024 10:20:19 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordstaak 4, vraag 4, Box-Cox Normailty plot
Estimated Impact194
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Box-Cox Normality Plot] [Taak 4 Vraag 4] [2007-11-04 11:48:04] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message
Dataseries X:
63,1
63
74,9
64,7
68,5
73,9
57,3
62,5
75,7
73,5
71
76,1
65
66,3
74,4
71,1
71,9
75,7
64,3
62,6
75,3
80,3
72,2
80,3
71,3
70,7
77,5
79,2
75,4
83,4
72,9
68,4
85
86,2
76,7
91,2
76
77
87,3
83,1
81,1
88,4
76,7
70,9
85,9
84,8
79,6
89,9
74
74,8
91,4
84,1
80
92,2
77,3
78,1
93,6
92,9
89,8
101,8
84,8
90,8
106,8
92,2
101,5
101,6
86
89,7
106,3
107,8
106,7
118,1
100,3
97,5
112,9
101,3
107,8
110,3
94,6
93,9
104,8
109,4
107,6
112,3
99
99,9
111,1
111,2
109,5
111,6
102,2
94,5
114,4
116,8
110,5
121,1
106,3
104,6
117,3
112,4
107,5
116,9
104,9
95,9
121,5
124,1
110,7
128,1
108,9
111,2
131,9
120,5
115,4
130,8
112,6
111,1
129,5
131,8
127,7
142,6
115,8
119,5
139,8
130,6
126,2
141,7
118,9
120,3
141,2
138,9
141,3
153,7
131
136,5
160,1
139,8
143,8




Summary of compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\begin{tabular}{lllllllll}
\hline
Summary of compuational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 4 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=293&T=0

[TABLE]
[ROW][C]Summary of compuational 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]4 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=293&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=293&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 compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001







Box-Cox Normality Plot
# observations x137
maximum correlation0.928124787185849
optimal lambda0.29

\begin{tabular}{lllllllll}
\hline
Box-Cox Normality Plot \tabularnewline
# observations x & 137 \tabularnewline
maximum correlation & 0.928124787185849 \tabularnewline
optimal lambda & 0.29 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=293&T=1

[TABLE]
[ROW][C]Box-Cox Normality Plot[/C][/ROW]
[ROW][C]# observations x[/C][C]137[/C][/ROW]
[ROW][C]maximum correlation[/C][C]0.928124787185849[/C][/ROW]
[ROW][C]optimal lambda[/C][C]0.29[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=293&T=1

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

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Box-Cox Normality Plot
# observations x137
maximum correlation0.928124787185849
optimal lambda0.29



Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
n <- length(x)
c <- array(NA,dim=c(401))
l <- array(NA,dim=c(401))
mx <- 0
mxli <- -999
for (i in 1:401)
{
l[i] <- (i-201)/100
if (l[i] != 0)
{
x1 <- (x^l[i] - 1) / l[i]
} else {
x1 <- log(x)
}
c[i] <- cor(qnorm(ppoints(x), mean=0, sd=1),x1)
if (mx < c[i])
{
mx <- c[i]
mxli <- l[i]
}
}
c
mx
mxli
if (mxli != 0)
{
x1 <- (x^mxli - 1) / mxli
} else {
x1 <- log(x)
}
bitmap(file='test1.png')
plot(l,c,main='Box-Cox Normality Plot',xlab='Lambda',ylab='correlation')
mtext(paste('Optimal Lambda =',mxli))
grid()
dev.off()
bitmap(file='test2.png')
hist(x,main='Histogram of Original Data',xlab='X',ylab='frequency')
grid()
dev.off()
bitmap(file='test3.png')
hist(x1,main='Histogram of Transformed Data',xlab='X',ylab='frequency')
grid()
dev.off()
bitmap(file='test4.png')
qqnorm(x)
grid()
mtext('Original Data')
dev.off()
bitmap(file='test5.png')
qqnorm(x1)
grid()
mtext('Transformed Data')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Box-Cox Normality Plot',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'# observations x',header=TRUE)
a<-table.element(a,n)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'maximum correlation',header=TRUE)
a<-table.element(a,mx)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'optimal lambda',header=TRUE)
a<-table.element(a,mxli)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')