Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_tukeylambda.wasp
Title produced by softwareTukey lambda PPCC Plot
Date of computationSun, 26 Oct 2008 02:33:26 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Oct/26/t1225010039gunb2r82hj4wcoo.htm/, Retrieved Sat, 18 May 2024 13:59:00 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=18814, Retrieved Sat, 18 May 2024 13:59:00 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsInvestigating Distributions
Estimated Impact168
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Univariate Data Series] [Werkloosheid, bru...] [2008-10-13 18:17:44] [6743688719638b0cb1c0a6e0bf433315]
-   PD  [Univariate Data Series] [Reeks 1 Werkloosheid] [2008-10-20 05:40:56] [6743688719638b0cb1c0a6e0bf433315]
F RMP       [Tukey lambda PPCC Plot] [Investigating Dis...] [2008-10-26 08:33:26] [9b05d7ef5dbcfba4217d280d9092f628] [Current]
Feedback Forum
2008-12-01 18:57:47 [8e2cc0b2ef568da46d009b2f601285b2] [reply
De student heeft de juiste techniek toegepast maar komt tot de conclusie dat de hoogste correlatie zich niet voordoet bij lambda = 0,14. De student bespreekt dit spijtig genoeg niet verder wat de mogelijke gevolgen hiervan zijn. Een normaal verdeling doet zich voor bij lambda = 0,14. Vermits de gevonde correlatie zich voordoet bij 0.5 kunnen we al stellen dat deze reeks niet normaal verdeeld is.

Post a new message
Dataseries X:
569
580
578
565
547
555
562
561
555
544
537
543
594
611
613
611
594
595
591
589
584
573
567
569
621
629
628
612
595
597
593
590
580
574
573
573
620
626
620
588
566
557
561
549
532
526
511
499
555
565
542
527
510
514
517
508
493
490
469
478
528




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=18814&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







Tukey Lambda - Key Values
Distribution (lambda)Correlation
Approx. Cauchy (lambda=-1)0.658237974131817
Exact Logistic (lambda=0)0.980929859116424
Approx. Normal (lambda=0.14)0.988126664959712
U-shaped (lambda=0.5)0.988363027443174
Exactly Uniform (lambda=1)0.980562496534658

\begin{tabular}{lllllllll}
\hline
Tukey Lambda - Key Values \tabularnewline
Distribution (lambda) & Correlation \tabularnewline
Approx. Cauchy (lambda=-1) & 0.658237974131817 \tabularnewline
Exact Logistic (lambda=0) & 0.980929859116424 \tabularnewline
Approx. Normal (lambda=0.14) & 0.988126664959712 \tabularnewline
U-shaped (lambda=0.5) & 0.988363027443174 \tabularnewline
Exactly Uniform (lambda=1) & 0.980562496534658 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=18814&T=1

[TABLE]
[ROW][C]Tukey Lambda - Key Values[/C][/ROW]
[ROW][C]Distribution (lambda)[/C][C]Correlation[/C][/ROW]
[ROW][C]Approx. Cauchy (lambda=-1)[/C][C]0.658237974131817[/C][/ROW]
[ROW][C]Exact Logistic (lambda=0)[/C][C]0.980929859116424[/C][/ROW]
[ROW][C]Approx. Normal (lambda=0.14)[/C][C]0.988126664959712[/C][/ROW]
[ROW][C]U-shaped (lambda=0.5)[/C][C]0.988363027443174[/C][/ROW]
[ROW][C]Exactly Uniform (lambda=1)[/C][C]0.980562496534658[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=18814&T=1

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

As an alternative you can also use a QR Code:  

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

Tukey Lambda - Key Values
Distribution (lambda)Correlation
Approx. Cauchy (lambda=-1)0.658237974131817
Exact Logistic (lambda=0)0.980929859116424
Approx. Normal (lambda=0.14)0.988126664959712
U-shaped (lambda=0.5)0.988363027443174
Exactly Uniform (lambda=1)0.980562496534658



Parameters (Session):
Parameters (R input):
R code (references can be found in the software module):
gp <- function(lambda, p)
{
(p^lambda-(1-p)^lambda)/lambda
}
sortx <- sort(x)
c <- array(NA,dim=c(201))
for (i in 1:201)
{
if (i != 101) c[i] <- cor(gp(ppoints(x), lambda=(i-101)/100),sortx)
}
bitmap(file='test1.png')
plot((-100:100)/100,c[1:201],xlab='lambda',ylab='correlation',main='PPCC Plot - Tukey lambda')
grid()
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Tukey Lambda - Key Values',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Distribution (lambda)',1,TRUE)
a<-table.element(a,'Correlation',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Approx. Cauchy (lambda=-1)',header=TRUE)
a<-table.element(a,c[1])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Exact Logistic (lambda=0)',header=TRUE)
a<-table.element(a,(c[100]+c[102])/2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Approx. Normal (lambda=0.14)',header=TRUE)
a<-table.element(a,c[115])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'U-shaped (lambda=0.5)',header=TRUE)
a<-table.element(a,c[151])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Exactly Uniform (lambda=1)',header=TRUE)
a<-table.element(a,c[201])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')