Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_boxcoxnorm.wasp
Title produced by softwareBox-Cox Normality Plot
Date of computationMon, 09 Nov 2009 06:02:38 -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/2009/Nov/09/t12577718690vraiweldkl5e2a.htm/, Retrieved Wed, 24 Apr 2024 00:46:42 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=54797, Retrieved Wed, 24 Apr 2024 00:46:42 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact182
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Box-Cox Normality Plot] [3/11/2009] [2009-11-02 22:22:24] [b98453cac15ba1066b407e146608df68]
-    D    [Box-Cox Normality Plot] [Box Cox Norm Plot...] [2009-11-09 13:02:38] [4f297b039e1043ebee7ff7a83b1eaaaa] [Current]
Feedback Forum

Post a new message
Dataseries X:
124.06
124.58
122.00
124.02
124.16
124.29
123.93
124.62
121.81
124.14
124.31
125.15
125.35
125.48
124.17
125.33
124.46
123.39
123.14
122.24
119.31
120.87
120.43
119.41
118.85
119.08
117.25
118.51
118.42
118.56
117.97
117.98
115.25
117.23
117.08
116.83
117.17
117.73
115.74
116.99
116.90
116.49
115.84
115.92
113.32
114.84
114.75
114.84
115.03
115.03
112.99
114.15
113.77
113.57
113.38
112.71
110.27
111.73
112.12
112.31
111.73
111.83
109.99
111.15
111.25
110.87
110.27
110.18
108.15
109.60
109.60
109.41
109.80
109.60
107.76
109.02
108.62
109.02
109.22
108.92
106.69
107.76
107.66
107.85
107.95
107.85
106.30
107.37
107.66
107.46
107.37
107.18
105.43
106.39
106.50
106.50
106.69
106.50
105.14
106.50
106.20
105.72
104.76
104.55
102.71
104.36
104.65
104.46
104.65
103.88
102.32
103.39
103.00
102.71
102.51
102.04
100.00




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 2 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=54797&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=54797&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=54797&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 time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135







Box-Cox Normality Plot
# observations x117
maximum correlation0
optimal lambda-999

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=54797&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 x117
maximum correlation0
optimal lambda-999



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)
qqline(x)
grid()
mtext('Original Data')
dev.off()
bitmap(file='test5.png')
qqnorm(x1)
qqline(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')