Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_pairs.wasp
Title produced by softwareKendall tau Correlation Matrix
Date of computationFri, 24 Dec 2010 11:03:04 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/24/t1293196274pqv4trtndio3k67.htm/, Retrieved Tue, 30 Apr 2024 04:59:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=114902, Retrieved Tue, 30 Apr 2024 04:59:16 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact144
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Recursive Partitioning (Regression Trees)] [] [2010-12-05 18:59:57] [b98453cac15ba1066b407e146608df68]
- RMPD    [Kendall tau Correlation Matrix] [KCM] [2010-12-24 11:03:04] [fd751bc40fbbb4c72222c10190589d42] [Current]
-    D      [Kendall tau Correlation Matrix] [workshop 10 link 2] [2010-12-25 09:11:52] [cc4c09289ddf8962388fdbedfd8171c3]
-    D      [Kendall tau Correlation Matrix] [workshop 10 link 2] [2010-12-25 09:11:52] [cc4c09289ddf8962388fdbedfd8171c3]
Feedback Forum

Post a new message
Dataseries X:
1,8	0,8	2,9	1,8	2,3	0,8	2,6
1,7	-0,1	2,9	1,7	2,2	1	2,2
1,4	-1,5	2,9	1,6	2,1	0,6	2,3
1,2	-4,4	1,4	1,8	2,4	0,9	2,4
1	-4,2	1,1	1,6	2,5	0,6	2,1
1,7	3,5	1,9	1,5	2,4	0,6	1,9
2,4	10	2,8	1,5	2,3	0,4	2,2
2	8,6	1,4	1,3	2,1	0,3	1,9
2,1	9,5	0,7	1,4	2,3	0	2,3
2	9,9	-0,8	1,4	2,2	0,3	2,1
1,8	10,4	-3,1	1,3	2,1	0,1	2,2
2,7	16	0,1	1,3	2	0	2,3
2,3	12,7	1	1,2	2,1	0	1,9
1,9	10,2	1,9	1,1	2,1	0	1,7
2	8,9	-0,5	1,4	2,5	-0,2	2,5
2,3	12,6	1,5	1,2	2,2	-0,3	2,1
2,8	13,6	3,9	1,5	2,3	0,1	2,4
2,4	14,8	1,9	1,1	2,3	0,1	1,5
2,3	9,5	2,6	1,3	2,2	0,4	1,9
2,7	13,7	1,7	1,5	2,2	0,4	2,1
2,7	17	1,4	1,1	1,6	-0,5	2,2
2,9	14,7	2,8	1,4	1,8	0,5	2
3	17,4	0,5	1,3	1,7	0,4	2
2,2	9	1	1,5	1,9	0,7	2,2
2,3	9,1	1,5	1,6	1,8	0,8	2,3
2,8	12,2	1,8	1,7	1,9	0,8	2,3
2,8	15,9	2,7	1,1	1,5	0	2
2,8	12,9	3	1,6	1	1,1	2,2
2,2	10,9	-0,3	1,3	0,8	0,9	1,9
2,6	10,6	1,1	1,7	1,1	1,1	2,3
2,8	13,2	1,7	1,6	1,5	1	2,2
2,5	9,6	1,6	1,7	1,7	1,1	2,3
2,4	6,4	3	1,9	2,3	1,5	2,1
2,3	5,8	3,3	1,8	2,4	1	2,4
1,9	-1	6,7	1,9	3	1	2,3
1,7	-0,2	5,6	1,6	3	0,9	1,9
2	2,7	6	1,5	3,2	0,8	1,6
2,1	3,6	4,8	1,6	3,2	0,8	1,8
1,7	-0,9	5,9	1,6	3,2	0,8	1,8
1,8	0,3	4,3	1,7	3,5	0,8	2
1,8	-1,1	3,7	2	4	0,9	2,3
1,8	-2,5	5,6	2	4,3	0,8	2,2
1,3	-3,4	1,7	1,9	4,1	0,7	2,2
1,3	-3,5	3,2	1,7	4	0,6	2
1,3	-3,9	3,6	1,8	4,1	0,6	2
1,2	-4,6	1,7	1,9	4,2	1	1,9
1,4	-0,1	0,5	1,7	4,5	1	1,5
2,2	4,3	2,1	2	5,6	1	1,6
2,9	10,2	1,5	2,1	6,5	1,1	1,5
3,1	8,7	2,7	2,4	7,6	1,1	2
3,5	13,3	1,4	2,5	8,5	1,4	1,5
3,6	15	1,2	2,5	8,7	1,2	1,5
4,4	20,7	2,3	2,6	8,3	1,2	1,9
4,1	20,7	1,6	2,2	8,3	1,3	1,1
5,1	26,4	4,7	2,5	8,5	1,4	1,5
5,8	31,2	3,5	2,8	8,7	1,4	2,1
5,9	31,4	4,4	2,8	8,7	1,1	2,3
5,4	26,6	3,9	2,9	8,5	1,1	2,6
5,5	26,6	3,5	3	7,9	1,3	2,9
4,8	19,2	3	3,1	7	1,5	3,2
3,2	6,5	1,6	2,9	5,8	1,5	3,2




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 3 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=114902&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=114902&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=114902&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 time3 seconds
R Server'George Udny Yule' @ 72.249.76.132







Correlations for all pairs of data series (method=pearson)
HICPEDNBLITBLNEID
HICP10.9080.1480.6990.7040.4160.142
ED0.9081-0.1120.3570.4380.080.045
NBL0.148-0.11210.3020.2410.3450.008
IT0.6990.3570.30210.880.7950.266
BL0.7040.4380.2410.8810.602-0.124
NEI0.4160.080.3450.7950.60210.042
D0.1420.0450.0080.266-0.1240.0421

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & HICP & ED & NBL & IT & BL & NEI & D \tabularnewline
HICP & 1 & 0.908 & 0.148 & 0.699 & 0.704 & 0.416 & 0.142 \tabularnewline
ED & 0.908 & 1 & -0.112 & 0.357 & 0.438 & 0.08 & 0.045 \tabularnewline
NBL & 0.148 & -0.112 & 1 & 0.302 & 0.241 & 0.345 & 0.008 \tabularnewline
IT & 0.699 & 0.357 & 0.302 & 1 & 0.88 & 0.795 & 0.266 \tabularnewline
BL & 0.704 & 0.438 & 0.241 & 0.88 & 1 & 0.602 & -0.124 \tabularnewline
NEI & 0.416 & 0.08 & 0.345 & 0.795 & 0.602 & 1 & 0.042 \tabularnewline
D & 0.142 & 0.045 & 0.008 & 0.266 & -0.124 & 0.042 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=114902&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]HICP[/C][C]ED[/C][C]NBL[/C][C]IT[/C][C]BL[/C][C]NEI[/C][C]D[/C][/ROW]
[ROW][C]HICP[/C][C]1[/C][C]0.908[/C][C]0.148[/C][C]0.699[/C][C]0.704[/C][C]0.416[/C][C]0.142[/C][/ROW]
[ROW][C]ED[/C][C]0.908[/C][C]1[/C][C]-0.112[/C][C]0.357[/C][C]0.438[/C][C]0.08[/C][C]0.045[/C][/ROW]
[ROW][C]NBL[/C][C]0.148[/C][C]-0.112[/C][C]1[/C][C]0.302[/C][C]0.241[/C][C]0.345[/C][C]0.008[/C][/ROW]
[ROW][C]IT[/C][C]0.699[/C][C]0.357[/C][C]0.302[/C][C]1[/C][C]0.88[/C][C]0.795[/C][C]0.266[/C][/ROW]
[ROW][C]BL[/C][C]0.704[/C][C]0.438[/C][C]0.241[/C][C]0.88[/C][C]1[/C][C]0.602[/C][C]-0.124[/C][/ROW]
[ROW][C]NEI[/C][C]0.416[/C][C]0.08[/C][C]0.345[/C][C]0.795[/C][C]0.602[/C][C]1[/C][C]0.042[/C][/ROW]
[ROW][C]D[/C][C]0.142[/C][C]0.045[/C][C]0.008[/C][C]0.266[/C][C]-0.124[/C][C]0.042[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=114902&T=1

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

As an alternative you can also use a QR Code:  

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

Correlations for all pairs of data series (method=pearson)
HICPEDNBLITBLNEID
HICP10.9080.1480.6990.7040.4160.142
ED0.9081-0.1120.3570.4380.080.045
NBL0.148-0.11210.3020.2410.3450.008
IT0.6990.3570.30210.880.7950.266
BL0.7040.4380.2410.8810.602-0.124
NEI0.4160.080.3450.7950.60210.042
D0.1420.0450.0080.266-0.1240.0421







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
HICP;ED0.90820.89180.7531
p-value(0)(0)(0)
HICP;NBL0.14820.02590.0454
p-value(0.2543)(0.8427)(0.6133)
HICP;IT0.69870.30140.1892
p-value(0)(0.0183)(0.0389)
HICP;BL0.70420.1720.0592
p-value(0)(0.1851)(0.5121)
HICP;NEI0.41570.40640.2726
p-value(9e-04)(0.0012)(0.0029)
HICP;D0.14220.06030.0355
p-value(0.2743)(0.6446)(0.7016)
ED;NBL-0.1115-0.1405-0.0939
p-value(0.3921)(0.28)(0.2896)
ED;IT0.3574-0.0034-0.0497
p-value(0.0047)(0.9791)(0.5824)
ED;BL0.4377-0.0187-0.0634
p-value(4e-04)(0.8861)(0.4773)
ED;NEI0.08040.1120.0684
p-value(0.538)(0.3901)(0.4495)
ED;D0.0452-0.0039-0.0103
p-value(0.7296)(0.976)(0.9101)
NBL;IT0.30160.39490.2659
p-value(0.0182)(0.0016)(0.0035)
NBL;BL0.24060.38530.255
p-value(0.0618)(0.0022)(0.0045)
NBL;NEI0.34520.30750.2143
p-value(0.0064)(0.0159)(0.0187)
NBL;D0.00810.04430.0359
p-value(0.9504)(0.7346)(0.6972)
IT;BL0.880.7820.6209
p-value(0)(0)(0)
IT;NEI0.79520.87370.7323
p-value(0)(0)(0)
IT;D0.26570.15390.1242
p-value(0.0385)(0.2363)(0.1862)
BL;NEI0.60160.5570.422
p-value(0)(0)(0)
BL;D-0.1244-0.1859-0.1374
p-value(0.3395)(0.1514)(0.1383)
NEI;D0.04240.00390.0024
p-value(0.7453)(0.9764)(0.9799)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
HICP;ED & 0.9082 & 0.8918 & 0.7531 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
HICP;NBL & 0.1482 & 0.0259 & 0.0454 \tabularnewline
p-value & (0.2543) & (0.8427) & (0.6133) \tabularnewline
HICP;IT & 0.6987 & 0.3014 & 0.1892 \tabularnewline
p-value & (0) & (0.0183) & (0.0389) \tabularnewline
HICP;BL & 0.7042 & 0.172 & 0.0592 \tabularnewline
p-value & (0) & (0.1851) & (0.5121) \tabularnewline
HICP;NEI & 0.4157 & 0.4064 & 0.2726 \tabularnewline
p-value & (9e-04) & (0.0012) & (0.0029) \tabularnewline
HICP;D & 0.1422 & 0.0603 & 0.0355 \tabularnewline
p-value & (0.2743) & (0.6446) & (0.7016) \tabularnewline
ED;NBL & -0.1115 & -0.1405 & -0.0939 \tabularnewline
p-value & (0.3921) & (0.28) & (0.2896) \tabularnewline
ED;IT & 0.3574 & -0.0034 & -0.0497 \tabularnewline
p-value & (0.0047) & (0.9791) & (0.5824) \tabularnewline
ED;BL & 0.4377 & -0.0187 & -0.0634 \tabularnewline
p-value & (4e-04) & (0.8861) & (0.4773) \tabularnewline
ED;NEI & 0.0804 & 0.112 & 0.0684 \tabularnewline
p-value & (0.538) & (0.3901) & (0.4495) \tabularnewline
ED;D & 0.0452 & -0.0039 & -0.0103 \tabularnewline
p-value & (0.7296) & (0.976) & (0.9101) \tabularnewline
NBL;IT & 0.3016 & 0.3949 & 0.2659 \tabularnewline
p-value & (0.0182) & (0.0016) & (0.0035) \tabularnewline
NBL;BL & 0.2406 & 0.3853 & 0.255 \tabularnewline
p-value & (0.0618) & (0.0022) & (0.0045) \tabularnewline
NBL;NEI & 0.3452 & 0.3075 & 0.2143 \tabularnewline
p-value & (0.0064) & (0.0159) & (0.0187) \tabularnewline
NBL;D & 0.0081 & 0.0443 & 0.0359 \tabularnewline
p-value & (0.9504) & (0.7346) & (0.6972) \tabularnewline
IT;BL & 0.88 & 0.782 & 0.6209 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
IT;NEI & 0.7952 & 0.8737 & 0.7323 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
IT;D & 0.2657 & 0.1539 & 0.1242 \tabularnewline
p-value & (0.0385) & (0.2363) & (0.1862) \tabularnewline
BL;NEI & 0.6016 & 0.557 & 0.422 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
BL;D & -0.1244 & -0.1859 & -0.1374 \tabularnewline
p-value & (0.3395) & (0.1514) & (0.1383) \tabularnewline
NEI;D & 0.0424 & 0.0039 & 0.0024 \tabularnewline
p-value & (0.7453) & (0.9764) & (0.9799) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=114902&T=2

[TABLE]
[ROW][C]Correlations for all pairs of data series with p-values[/C][/ROW]
[ROW][C]pair[/C][C]Pearson r[/C][C]Spearman rho[/C][C]Kendall tau[/C][/ROW]
[ROW][C]HICP;ED[/C][C]0.9082[/C][C]0.8918[/C][C]0.7531[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]HICP;NBL[/C][C]0.1482[/C][C]0.0259[/C][C]0.0454[/C][/ROW]
[ROW][C]p-value[/C][C](0.2543)[/C][C](0.8427)[/C][C](0.6133)[/C][/ROW]
[ROW][C]HICP;IT[/C][C]0.6987[/C][C]0.3014[/C][C]0.1892[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0.0183)[/C][C](0.0389)[/C][/ROW]
[ROW][C]HICP;BL[/C][C]0.7042[/C][C]0.172[/C][C]0.0592[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0.1851)[/C][C](0.5121)[/C][/ROW]
[ROW][C]HICP;NEI[/C][C]0.4157[/C][C]0.4064[/C][C]0.2726[/C][/ROW]
[ROW][C]p-value[/C][C](9e-04)[/C][C](0.0012)[/C][C](0.0029)[/C][/ROW]
[ROW][C]HICP;D[/C][C]0.1422[/C][C]0.0603[/C][C]0.0355[/C][/ROW]
[ROW][C]p-value[/C][C](0.2743)[/C][C](0.6446)[/C][C](0.7016)[/C][/ROW]
[ROW][C]ED;NBL[/C][C]-0.1115[/C][C]-0.1405[/C][C]-0.0939[/C][/ROW]
[ROW][C]p-value[/C][C](0.3921)[/C][C](0.28)[/C][C](0.2896)[/C][/ROW]
[ROW][C]ED;IT[/C][C]0.3574[/C][C]-0.0034[/C][C]-0.0497[/C][/ROW]
[ROW][C]p-value[/C][C](0.0047)[/C][C](0.9791)[/C][C](0.5824)[/C][/ROW]
[ROW][C]ED;BL[/C][C]0.4377[/C][C]-0.0187[/C][C]-0.0634[/C][/ROW]
[ROW][C]p-value[/C][C](4e-04)[/C][C](0.8861)[/C][C](0.4773)[/C][/ROW]
[ROW][C]ED;NEI[/C][C]0.0804[/C][C]0.112[/C][C]0.0684[/C][/ROW]
[ROW][C]p-value[/C][C](0.538)[/C][C](0.3901)[/C][C](0.4495)[/C][/ROW]
[ROW][C]ED;D[/C][C]0.0452[/C][C]-0.0039[/C][C]-0.0103[/C][/ROW]
[ROW][C]p-value[/C][C](0.7296)[/C][C](0.976)[/C][C](0.9101)[/C][/ROW]
[ROW][C]NBL;IT[/C][C]0.3016[/C][C]0.3949[/C][C]0.2659[/C][/ROW]
[ROW][C]p-value[/C][C](0.0182)[/C][C](0.0016)[/C][C](0.0035)[/C][/ROW]
[ROW][C]NBL;BL[/C][C]0.2406[/C][C]0.3853[/C][C]0.255[/C][/ROW]
[ROW][C]p-value[/C][C](0.0618)[/C][C](0.0022)[/C][C](0.0045)[/C][/ROW]
[ROW][C]NBL;NEI[/C][C]0.3452[/C][C]0.3075[/C][C]0.2143[/C][/ROW]
[ROW][C]p-value[/C][C](0.0064)[/C][C](0.0159)[/C][C](0.0187)[/C][/ROW]
[ROW][C]NBL;D[/C][C]0.0081[/C][C]0.0443[/C][C]0.0359[/C][/ROW]
[ROW][C]p-value[/C][C](0.9504)[/C][C](0.7346)[/C][C](0.6972)[/C][/ROW]
[ROW][C]IT;BL[/C][C]0.88[/C][C]0.782[/C][C]0.6209[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]IT;NEI[/C][C]0.7952[/C][C]0.8737[/C][C]0.7323[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]IT;D[/C][C]0.2657[/C][C]0.1539[/C][C]0.1242[/C][/ROW]
[ROW][C]p-value[/C][C](0.0385)[/C][C](0.2363)[/C][C](0.1862)[/C][/ROW]
[ROW][C]BL;NEI[/C][C]0.6016[/C][C]0.557[/C][C]0.422[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]BL;D[/C][C]-0.1244[/C][C]-0.1859[/C][C]-0.1374[/C][/ROW]
[ROW][C]p-value[/C][C](0.3395)[/C][C](0.1514)[/C][C](0.1383)[/C][/ROW]
[ROW][C]NEI;D[/C][C]0.0424[/C][C]0.0039[/C][C]0.0024[/C][/ROW]
[ROW][C]p-value[/C][C](0.7453)[/C][C](0.9764)[/C][C](0.9799)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=114902&T=2

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

As an alternative you can also use a QR Code:  

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

Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
HICP;ED0.90820.89180.7531
p-value(0)(0)(0)
HICP;NBL0.14820.02590.0454
p-value(0.2543)(0.8427)(0.6133)
HICP;IT0.69870.30140.1892
p-value(0)(0.0183)(0.0389)
HICP;BL0.70420.1720.0592
p-value(0)(0.1851)(0.5121)
HICP;NEI0.41570.40640.2726
p-value(9e-04)(0.0012)(0.0029)
HICP;D0.14220.06030.0355
p-value(0.2743)(0.6446)(0.7016)
ED;NBL-0.1115-0.1405-0.0939
p-value(0.3921)(0.28)(0.2896)
ED;IT0.3574-0.0034-0.0497
p-value(0.0047)(0.9791)(0.5824)
ED;BL0.4377-0.0187-0.0634
p-value(4e-04)(0.8861)(0.4773)
ED;NEI0.08040.1120.0684
p-value(0.538)(0.3901)(0.4495)
ED;D0.0452-0.0039-0.0103
p-value(0.7296)(0.976)(0.9101)
NBL;IT0.30160.39490.2659
p-value(0.0182)(0.0016)(0.0035)
NBL;BL0.24060.38530.255
p-value(0.0618)(0.0022)(0.0045)
NBL;NEI0.34520.30750.2143
p-value(0.0064)(0.0159)(0.0187)
NBL;D0.00810.04430.0359
p-value(0.9504)(0.7346)(0.6972)
IT;BL0.880.7820.6209
p-value(0)(0)(0)
IT;NEI0.79520.87370.7323
p-value(0)(0)(0)
IT;D0.26570.15390.1242
p-value(0.0385)(0.2363)(0.1862)
BL;NEI0.60160.5570.422
p-value(0)(0)(0)
BL;D-0.1244-0.1859-0.1374
p-value(0.3395)(0.1514)(0.1383)
NEI;D0.04240.00390.0024
p-value(0.7453)(0.9764)(0.9799)



Parameters (Session):
par1 = pearson ;
Parameters (R input):
par1 = pearson ;
R code (references can be found in the software module):
panel.tau <- function(x, y, digits=2, prefix='', cex.cor)
{
usr <- par('usr'); on.exit(par(usr))
par(usr = c(0, 1, 0, 1))
rr <- cor.test(x, y, method=par1)
r <- round(rr$p.value,2)
txt <- format(c(r, 0.123456789), digits=digits)[1]
txt <- paste(prefix, txt, sep='')
if(missing(cex.cor)) cex <- 0.5/strwidth(txt)
text(0.5, 0.5, txt, cex = cex)
}
panel.hist <- function(x, ...)
{
usr <- par('usr'); on.exit(par(usr))
par(usr = c(usr[1:2], 0, 1.5) )
h <- hist(x, plot = FALSE)
breaks <- h$breaks; nB <- length(breaks)
y <- h$counts; y <- y/max(y)
rect(breaks[-nB], 0, breaks[-1], y, col='grey', ...)
}
bitmap(file='test1.png')
pairs(t(y),diag.panel=panel.hist, upper.panel=panel.smooth, lower.panel=panel.tau, main=main)
dev.off()
load(file='createtable')
n <- length(y[,1])
n
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,paste('Correlations for all pairs of data series (method=',par1,')',sep=''),n+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,' ',header=TRUE)
for (i in 1:n) {
a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
}
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],header=TRUE)
for (j in 1:n) {
r <- cor.test(y[i,],y[j,],method=par1)
a<-table.element(a,round(r$estimate,3))
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Correlations for all pairs of data series with p-values',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'pair',1,TRUE)
a<-table.element(a,'Pearson r',1,TRUE)
a<-table.element(a,'Spearman rho',1,TRUE)
a<-table.element(a,'Kendall tau',1,TRUE)
a<-table.row.end(a)
cor.test(y[1,],y[2,],method=par1)
for (i in 1:(n-1))
{
for (j in (i+1):n)
{
a<-table.row.start(a)
dum <- paste(dimnames(t(x))[[2]][i],';',dimnames(t(x))[[2]][j],sep='')
a<-table.element(a,dum,header=TRUE)
rp <- cor.test(y[i,],y[j,],method='pearson')
a<-table.element(a,round(rp$estimate,4))
rs <- cor.test(y[i,],y[j,],method='spearman')
a<-table.element(a,round(rs$estimate,4))
rk <- cor.test(y[i,],y[j,],method='kendall')
a<-table.element(a,round(rk$estimate,4))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=T)
a<-table.element(a,paste('(',round(rp$p.value,4),')',sep=''))
a<-table.element(a,paste('(',round(rs$p.value,4),')',sep=''))
a<-table.element(a,paste('(',round(rk$p.value,4),')',sep=''))
a<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable1.tab')