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_boxcoxnorm.wasp
Title produced by softwareBox-Cox Normality Plot
Date of computationWed, 29 Dec 2010 12:42:30 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/29/t12936264128und47rj729zhqk.htm/, Retrieved Fri, 03 May 2024 10:40:37 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=116763, Retrieved Fri, 03 May 2024 10:40:37 +0000
QR Codes:

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

Tree of Dependent Computations

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]
- R  D  [Box-Cox Normality Plot] [] [2009-12-19 15:37:16] [69c775ce4d55db2aa75a88e773e8d700]
-   P     [Box-Cox Normality Plot] [] [2010-12-26 18:31:19] [69c775ce4d55db2aa75a88e773e8d700]
-   PD        [Box-Cox Normality Plot] [] [2010-12-29 12:42:30] [e180d4cd19004beeddc12e67012247dc] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
09,166456
07,970589
07,104091
06,621064
07,529215
08,170938
08,157450
07,378962
07,921496
08,156740
08,856365
08,817177
08,734347
09,345927
08,992970
10,785120
08,886867
08,818847
08,823744
09,165298
08,652657
08,173054
07,563416
07,595809
08,381467
07,216432
06,540178
06,238914
05,487288
05,759462
05,993215
07,474726
07,348907
07,303379
07,119314
06,993780
06,958153
07,595706
08,088153
07,555753
07,315433
07,893427
08,858794
08,839367
08,014733
07,873465
08,930377
10,500550
12,611440
11,417870
11,872490
11,060820
12,043310
09,776299
09,557194
09,202590
10,224020
09,350807
08,300913
08,365779
08,133595
07,660470
08,074839
07,848597
07,998220
07,396895
07,900419
08,100500
07,899453
07,599783
08,100929
09,002175
10,298900
10,101520
10,699150
09,698140
09,800951
10,900470
10,697850
09,297252
10,397440
10,900720
12,901270
13,099060
11,698280
11,099870
11,301570
10,702110
10,099310
09,591119

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk

\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 & 3 seconds \tabularnewline
R Server & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=116763&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=116763&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116763&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 time3 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk






Box-Cox Normality Plot
# observations x90
maximum correlation0.503391276354205
optimal lambda1.65

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

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

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

As an alternative you can also use a QR Code:  
QR Code


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

Box-Cox Normality Plot
# observations x90
maximum correlation0.503391276354205
optimal lambda1.65


Figures (Output of Computation)

Input Parameters & R Code

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')