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 07:10:42 -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/t1257948823d9ot8eowuixejo8.htm/, Retrieved Fri, 29 Mar 2024 12:08:08 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=55608, Retrieved Fri, 29 Mar 2024 12:08:08 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsshwws6vr8
Estimated Impact159
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Box-Cox Linearity Plot] [3/11/2009] [2009-11-02 21:47:57] [b98453cac15ba1066b407e146608df68]
- R  D    [Box-Cox Linearity Plot] [] [2009-11-11 14:10:42] [4407d6264e55b051ec65750e6dca2820] [Current]
Feedback Forum

Post a new message
Dataseries X:
100
99,31916948
114,7106298
114,9941367
114,6085397
116,3633855
108,0492516
112,8502449
127,2518452
116,4999655
113,8511416
126,3771815
103,8290681
108,296889
125,0969166
115,6701386
130,2545354
131,2064565
124,8665241
122,6191626
148,2368766
117,7160792
128,9232255
130,7380837
113,1192661
116,5454922
127,9782024
132,9240533
133,9732358
121,2623301
128,559702
124,692695
140,9850314
124,6823481
134,9451611
137,2642616
132,6308891
122,7122853
132,5846727
151,5410085
145,764641
134,4043595
152,9192247
143,9387459
152,9709595
162,3356557
148,1658274
167,9499207
156,9400566
140,2241843
163,13996
160,969166
135,3031662
124,4554046
119,2143202
116,21163
125,2638477
114,0987791
110,0413879
126,6317169
113,0875354
Dataseries Y:
100
87,14054095
112,0054296
112,312101
109,474134
104,9746116
100,4926851
104,2154743
120,1768388
112,1028355
108,1481575
116,802197
102,1699512
95,15358705
120,6707808
111,5234277
119,9669448
113,3697401
110,0717661
111,5567342
132,424212
107,900558
122,1626615
124,3992258
110,4450505
101,5874013
122,3203962
125,2582826
125,4411543
108,9902468
118,9243879
116,7242723
134,1724901
116,8530994
124,5732995
130,9914031
123,4239103
111,4536725
124,5135991
139,2589613
129,8596099
112,3460359
131,381655
133,0004776
134,3220552
144,2379719
134,1278719
150,1891559
140,722563
114,8389975
143,1973003
140,2738676
112,1248303
102,8951536
100,5090242
103,3513901
111,4134533
104,5887587
101,7840983
114,7007441
108,7426474




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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=55608&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=55608&T=0

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







Box-Cox Linearity Plot
# observations x61
maximum correlation0.935100267940797
optimal lambda(x)0.88
Residual SD (orginial)4.78879405578631
Residual SD (transformed)4.78749157900284

\begin{tabular}{lllllllll}
\hline
Box-Cox Linearity Plot \tabularnewline
# observations x & 61 \tabularnewline
maximum correlation & 0.935100267940797 \tabularnewline
optimal lambda(x) & 0.88 \tabularnewline
Residual SD (orginial) & 4.78879405578631 \tabularnewline
Residual SD (transformed) & 4.78749157900284 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=55608&T=1

[TABLE]
[ROW][C]Box-Cox Linearity Plot[/C][/ROW]
[ROW][C]# observations x[/C][C]61[/C][/ROW]
[ROW][C]maximum correlation[/C][C]0.935100267940797[/C][/ROW]
[ROW][C]optimal lambda(x)[/C][C]0.88[/C][/ROW]
[ROW][C]Residual SD (orginial)[/C][C]4.78879405578631[/C][/ROW]
[ROW][C]Residual SD (transformed)[/C][C]4.78749157900284[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=55608&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=55608&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 x61
maximum correlation0.935100267940797
optimal lambda(x)0.88
Residual SD (orginial)4.78879405578631
Residual SD (transformed)4.78749157900284



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