Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_boxcoxnorm.wasp

Title produced by software

Box-Cox Normality Plot

Date of computation

Thu, 13 Nov 2008 01:38:29 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Nov/13/t1226565568n4oiupkoc86o5af.htm/, Retrieved Sun, 19 May 2024 11:40:09 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=24483, Retrieved Sun, 19 May 2024 11:40:09 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

natalie en evelyn

Estimated Impact

166

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Bivariate Kernel Density Estimation] [Bivariate Kernel ...] [2008-11-11 14:08:00] [adb6b6905cde49db36d59ca44433140d]
F RMPD  [Hierarchical Clustering] [Hierarchical Clus...] [2008-11-11 23:34:50] [b591abfa820a394aeb0c5ebd9cfa1091]
- RMPD      [Box-Cox Normality Plot] [box cox normal plot] [2008-11-13 08:38:29] [32a7b12f2bdf14b45f7a9a96ba1ab98d] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

3.253	11.836	345	6.06	519164
3.233	11.85	334	5.983	517009
3.196	11.897	345	6.11	509933
3.138	12.082	333	6.143	509127
3.091	11.936	336	6.093	500857
3.17	11.928	324	6.148	506971
3.378	12.646	320	6.464	569323
3.468	12.747	330	6.532	579714
3.33	12.447	313	6.321	577992
3.413	12.445	301	6.23	565464
3.356	12.257	288	6.176	547344
3.525	12.878	294	6.338	554788
3.633	13.69	302	6.462	562325
3.597	13.665	294	6.401	560854
3.6	13.78	293	6.46	555332
3.522	13.608	290	6.519	543599
3.503	13.375	283	6.542	536662
3.532	13.376	286	6.637	542722
3.686	13.918	293	7.114	593530
3.748	14.304	334	7.579	610763
3.672	13.877	329	7.408	612613
3.843	14.543	411	8.243	611324
3.905	14.291	416	8.243	594167
3.999	14.788	418	8.434	595454
4.07	15.241	408	8.576	590865
4.084	15.265	402	8.58	589379
4.042	15.322	401	8.645	584428
3.951	15.175	400	8.66	573100
3.933	14.817	389	8.72	567456
3.958	14.579	371	8.787	569028
4.147	15.247	364	9.162	620735
4.221	15.385	350	9.144	628884
4.058	14.891	332	8.806	628232
4.057	14.766	323	8.778	612117
4.089	14.42	316	8.66	595404
4.268	14.85	312	8.826	597141
4.309	15.117	315	8.609	593408
4.303	15.352	314	8.628	590072
4.177	15.099	313	8.619	579799
4.117	15.291	314	8.775	574205
4.065	15.208	317	8.84	572775
3.983	14.995	308	8.745	572942
4.091	15.454	312	9.092	619567
4.067	15.251	306	8.934	625809
4.024	14.975	304	8.749	619916
3.868	14.005	297	8.298	587625
3.8	13.55	284	8.067	565742
3.804	13.422	278	7.969	557274
3.862	13.848	273	7.999	560576
3.792	13.376	265	7.865	548854
3.674	13.038	259	7.746	531673
3.56	12.974	252	7.633	525919
3.489	12.554	245	7.458	511038
3.412	11.971	235	7.391	498662
3.674	12.916	232	7.856	555362
3.672	12.757	229	7.72	564591
3.463	11.924	219	7.297	541657
3.429	11.693	218	7.123	527070
3.4	11.382	215	7.004	509846
3.533	11.821	211	7.151	514258

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=24483&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=24483&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=24483&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 Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code