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_boxcoxlin.wasp
Title produced by softwareBox-Cox Linearity Plot
Date of computationThu, 12 Nov 2009 10:19:04 -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/12/t125804638146t9c9l7yx3jdsx.htm/, Retrieved Fri, 03 May 2024 19:08:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56266, Retrieved Fri, 03 May 2024 19:08:40 +0000
QR Codes:

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

Tree of Dependent Computations

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]
-    D    [Box-Cox Linearity Plot] [WS6.4] [2009-11-12 17:19:04] [dd4f17965cad1d38de7a1c062d32d75d] [Current]
- RMP       [Bagplot] [WS6.5] [2009-11-12 17:25:25] [626f1d98f4a7f05bcb9f17666b672c60]
- RMP       [Bivariate Kernel Density Estimation] [Ws6.6] [2009-11-12 17:28:51] [626f1d98f4a7f05bcb9f17666b672c60]
- RM        [Kendall tau Rank Correlation] [WS6.7] [2009-11-12 17:33:05] [626f1d98f4a7f05bcb9f17666b672c60]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
269285
269829
270911
266844
271244
269907
271296
270157
271322
267179
264101
265518
269419
268714
272482
268351
268175
270674
272764
272599
270333
270846
270491
269160
274027
273784
276663
274525
271344
271115
270798
273911
273985
271917
273338
270601
273547
275363
281229
277793
279913
282500
280041
282166
290304
283519
287816
285226
287595
289741
289148
288301
290155
289648
288225
289351
294735
305333
309030
310215
321935
325734
320846
323023
319753
321753
320757
324479
Dataseries Y:
Download CSV
Histogram 
258596
259056
264193
260325
261890
260683
257941
258151
262434
261577
262188
261092
263571
265031
270388
265458
266218
266386
263486
263620
267755
266554
266981
264133
265980
267183
272113
267261
269117
269034
266609
267261
271406
269529
270282
268663
269847
270998
277068
273529
275307
276488
274455
274507
279528
277673
278102
275131
277162
278799
285502
280672
281342
281132
278286
279120
289131
294453
295733
302233
308859
311054
318130
315823
316517
316907
314969
316107

Tables (Output of Computation)





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=56266&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=56266&T=0

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






Box-Cox Linearity Plot
# observations x68
maximum correlation0.982856319795715
optimal lambda(x)2
Residual SD (orginial)3158.18323896516
Residual SD (transformed)3089.82044069362

\begin{tabular}{lllllllll}
\hline
Box-Cox Linearity Plot \tabularnewline
# observations x & 68 \tabularnewline
maximum correlation & 0.982856319795715 \tabularnewline
optimal lambda(x) & 2 \tabularnewline
Residual SD (orginial) & 3158.18323896516 \tabularnewline
Residual SD (transformed) & 3089.82044069362 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=56266&T=1

[TABLE]
[ROW][C]Box-Cox Linearity Plot[/C][/ROW]
[ROW][C]# observations x[/C][C]68[/C][/ROW]
[ROW][C]maximum correlation[/C][C]0.982856319795715[/C][/ROW]
[ROW][C]optimal lambda(x)[/C][C]2[/C][/ROW]
[ROW][C]Residual SD (orginial)[/C][C]3158.18323896516[/C][/ROW]
[ROW][C]Residual SD (transformed)[/C][C]3089.82044069362[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=56266&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=56266&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 Linearity Plot
# observations x68
maximum correlation0.982856319795715
optimal lambda(x)2
Residual SD (orginial)3158.18323896516
Residual SD (transformed)3089.82044069362


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