Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_boxcoxlin.wasp
Title produced by softwareBox-Cox Linearity Plot
Date of computationWed, 11 Nov 2009 08:06:58 -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/11/t1257952073dxnrbcvxe6h7cke.htm/, Retrieved Fri, 29 Mar 2024 09:19:07 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=55662, Retrieved Fri, 29 Mar 2024 09:19:07 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsshwws6v1
Estimated Impact133
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Box-Cox Linearity Plot] [] [2009-11-11 15:06:58] [efdfe680cd785c4af09f858b30f777ec] [Current]
Feedback Forum

Post a new message
Dataseries X:
6539
6699
6962
6981
7024
6940
6774
6671
6965
6969
6822
6878
6691
6837
7018
7167
7076
7171
7093
6971
7142
7047
6999
6650
6475
6437
6639
6422
6272
6232
6003
5673
6050
5977
5796
5752
5609
5839
6069
6006
5809
5797
5502
5568
5864
5764
5615
5615
5681
5915
6334
6494
6620
6578
6495
6538
6737
6651
6530
6563
Dataseries Y:
2605
2682
2755
2760
2735
2659
2654
2670
2785
2845
2723
2746
2767
2940
2977
2993
2892
2824
2771
2686
2738
2723
2731
2632
2606
2605
2646
2627
2535
2456
2404
2319
2519
2504
2382
2394
2381
2501
2532
2515
2429
2389
2261
2272
2439
2373
2327
2364
2388
2553
2663
2694
2679
2611
2580
2627
2732
2707
2633
2683




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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 & 2 seconds \tabularnewline
R Server & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=55662&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=55662&T=0

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







Box-Cox Linearity Plot
# observations x60
maximum correlation0.931541927756756
optimal lambda(x)0.03
Residual SD (orginial)63.7604043862645
Residual SD (transformed)63.57574164573

\begin{tabular}{lllllllll}
\hline
Box-Cox Linearity Plot \tabularnewline
# observations x & 60 \tabularnewline
maximum correlation & 0.931541927756756 \tabularnewline
optimal lambda(x) & 0.03 \tabularnewline
Residual SD (orginial) & 63.7604043862645 \tabularnewline
Residual SD (transformed) & 63.57574164573 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=55662&T=1

[TABLE]
[ROW][C]Box-Cox Linearity Plot[/C][/ROW]
[ROW][C]# observations x[/C][C]60[/C][/ROW]
[ROW][C]maximum correlation[/C][C]0.931541927756756[/C][/ROW]
[ROW][C]optimal lambda(x)[/C][C]0.03[/C][/ROW]
[ROW][C]Residual SD (orginial)[/C][C]63.7604043862645[/C][/ROW]
[ROW][C]Residual SD (transformed)[/C][C]63.57574164573[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=55662&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=55662&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 Linearity Plot
# observations x60
maximum correlation0.931541927756756
optimal lambda(x)0.03
Residual SD (orginial)63.7604043862645
Residual SD (transformed)63.57574164573



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(x1,y)
if (mx < abs(c[i]))
{
mx <- abs(c[i])
mxli <- l[i]
}
}
c
mx
mxli
if (mxli != 0)
{
x1 <- (x^mxli - 1) / mxli
} else {
x1 <- log(x)
}
r<-lm(y~x)
se <- sqrt(var(r$residuals))
r1 <- lm(y~x1)
se1 <- sqrt(var(r1$residuals))
bitmap(file='test1.png')
plot(l,c,main='Box-Cox Linearity Plot',xlab='Lambda',ylab='correlation')
grid()
dev.off()
bitmap(file='test2.png')
plot(x,y,main='Linear Fit of Original Data',xlab='x',ylab='y')
abline(r)
grid()
mtext(paste('Residual Standard Deviation = ',se))
dev.off()
bitmap(file='test3.png')
plot(x1,y,main='Linear Fit of Transformed Data',xlab='x',ylab='y')
abline(r1)
grid()
mtext(paste('Residual Standard Deviation = ',se1))
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Box-Cox Linearity 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(x)',header=TRUE)
a<-table.element(a,mxli)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Residual SD (orginial)',header=TRUE)
a<-table.element(a,se)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Residual SD (transformed)',header=TRUE)
a<-table.element(a,se1)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')