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 computationThu, 15 Dec 2011 10:42:48 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Dec/15/t1323963807jkfkh1ad6avcync.htm/, Retrieved Wed, 08 May 2024 07:25:08 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=155518, Retrieved Wed, 08 May 2024 07:25:08 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact87
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Kendall tau Correlation Matrix] [] [2011-12-15 15:42:48] [05d3841c0e91f0207133db830e88168b] [Current]
Feedback Forum

Post a new message
Dataseries X:
67	96	38	116	3	140824
63	67	34	127	4	110459
69	70	42	106	16	105079
103	134	38	133	2	112098
49	59	27	64	1	43929
28	8	35	89	3	76173
113	145	33	122	0	187326
19	1	18	22	0	22807
57	71	34	117	7	144408
43	82	33	82	0	66485
102	92	42	136	0	79089
110	106	55	184	7	81625
65	50	35	106	7	68788
74	113	51	159	4	103297
79	70	42	86	10	69446
174	168	59	199	0	114948
66	111	36	139	4	167949
154	96	39	92	4	125081
52	102	29	85	3	125818
82	135	46	174	8	136588
68	122	45	148	0	112431
102	86	39	144	1	103037
39	50	25	84	5	82317
54	97	52	208	9	118906
110	127	41	144	0	83515
112	86	38	139	0	104581
126	99	41	127	5	103129
84	117	39	136	0	83243
51	57	32	99	0	37110
63	125	41	135	0	113344
73	120	45	165	3	139165
72	44	46	135	5	86652
83	133	48	178	1	112302
35	43	37	137	4	69652
90	117	39	148	3	119442
56	83	42	127	0	69867
118	105	41	141	0	101629
79	79	36	89	2	70168
32	33	17	46	1	31081
180	116	39	143	2	103925
78	121	37	116	10	92622
62	67	38	103	8	79011
72	73	36	108	5	93487
56	68	42	126	6	64520
82	50	45	45	1	93473
146	101	38	122	2	114360
42	20	26	66	2	33032
75	101	52	180	0	96125
113	137	47	165	10	151911
54	99	45	146	3	89256
72	94	40	137	0	95671
24	8	4	7	0	5950
303	85	44	157	8	149695
17	21	18	61	5	32551
64	30	14	41	3	31701
56	96	37	120	1	100087
82	122	56	208	5	169707
171	115	36	127	5	150491
131	139	41	147	0	120192
82	89	36	127	12	95893
136	147	46	161	10	151715
113	135	28	73	12	176225
102	77	42	94	10	59900
86	72	38	129	8	104767
64	47	37	125	2	114799
65	96	30	87	0	72128
125	79	35	128	6	143592
139	85	44	148	9	89626
77	135	36	116	2	131072
66	143	28	89	5	126817
67	99	45	154	13	81351
32	22	23	67	6	22618
80	78	45	171	7	88977
52	77	38	90	2	92059
59	110	38	133	1	81897
76	132	42	137	4	108146
89	112	36	133	3	126372
106	78	41	125	6	249771
60	126	38	134	2	71154
60	73	37	110	0	71571
46	62	28	89	1	55918
111	143	45	138	0	160141
68	30	26	99	5	38692
103	117	44	92	2	102812
25	49	8	27	0	56622
53	26	27	77	0	15986
53	71	35	127	5	123534
175	59	37	137	1	108535
110	114	57	122	0	93879
102	161	41	143	1	144551
88	74	37	85	1	56750
73	151	38	131	3	127654
61	41	31	90	6	65594
72	121	36	135	1	59938
76	66	36	132	4	146975
36	83	36	139	3	143372
50	94	35	127	5	168553
74	154	39	104	0	183500
144	151	58	221	12	165986
105	164	30	106	13	184923
121	116	45	153	8	140358
62	140	41	130	0	149959
175	73	36	59	0	57224
14	13	19	64	4	43750
79	89	23	36	4	48029
130	90	40	88	0	104978
46	128	40	125	0	100046
87	169	40	124	0	101047
64	28	30	83	0	197426
86	116	41	127	0	160902
67	76	40	143	4	147172
85	145	45	115	0	109432
11	12	1	0	0	1168
70	120	36	94	0	83248
25	23	11	30	4	25162
48	83	45	119	0	45724
114	131	38	102	1	110529
16	4	0	0	0	855
52	81	30	77	5	101382
22	18	8	9	0	14116
110	103	39	137	3	89506
63	76	44	150	7	135356
83	55	44	137	13	116066
51	43	29	84	3	144244
34	16	8	21	0	8773
39	66	39	139	2	102153
80	137	47	168	0	117440
57	50	48	155	0	104128
77	134	46	161	4	134238
96	152	48	145	0	134047
121	137	50	175	3	279488
35	71	40	137	0	79756
42	42	36	100	0	66089
319	84	40	150	4	102070
164	103	46	163	4	146760
50	55	39	137	15	154771
127	127	42	149	0	165933
76	55	39	112	4	64593
46	104	41	135	1	92280
87	95	42	114	1	67150
111	35	32	45	0	128692
115	95	39	120	9	124089
83	121	35	111	1	125386
63	41	21	78	3	37238
98	143	45	136	11	140015
57	147	50	179	5	150047
81	97	36	118	2	154451
100	170	44	147	1	156349
0	0	0	0	9	0
10	4	0	0	0	6023
1	0	0	0	0	0
2	0	0	0	0	0
0	0	0	0	1	0
0	0	0	0	0	0
82	61	37	88	2	84601
139	130	47	115	3	68946
0	0	0	0	0	0
4	0	0	0	0	0
5	7	0	0	0	1644
20	12	5	13	0	6179
5	0	1	4	0	3926
42	37	43	76	0	52789
2	0	0	0	0	0
63	48	31	63	2	100350




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net

\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 & 'Gwilym Jenkins' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155518&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]'Gwilym Jenkins' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155518&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=155518&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'Gwilym Jenkins' @ jenkins.wessa.net







Correlations for all pairs of data series (method=pearson)
loginsblogged_computationsreviewed_compendiumslong_feedback_messagesshared_compendiumsnumber_characters
logins10.5520.5890.5560.1960.527
blogged_computations0.55210.7430.750.1230.709
reviewed_compendiums0.5890.74310.9130.2490.684
long_feedback_messages0.5560.750.91310.2880.709
shared_compendiums0.1960.1230.2490.28810.288
number_characters0.5270.7090.6840.7090.2881

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series (method=pearson) \tabularnewline
  & logins & blogged_computations & reviewed_compendiums & long_feedback_messages & shared_compendiums & number_characters \tabularnewline
logins & 1 & 0.552 & 0.589 & 0.556 & 0.196 & 0.527 \tabularnewline
blogged_computations & 0.552 & 1 & 0.743 & 0.75 & 0.123 & 0.709 \tabularnewline
reviewed_compendiums & 0.589 & 0.743 & 1 & 0.913 & 0.249 & 0.684 \tabularnewline
long_feedback_messages & 0.556 & 0.75 & 0.913 & 1 & 0.288 & 0.709 \tabularnewline
shared_compendiums & 0.196 & 0.123 & 0.249 & 0.288 & 1 & 0.288 \tabularnewline
number_characters & 0.527 & 0.709 & 0.684 & 0.709 & 0.288 & 1 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155518&T=1

[TABLE]
[ROW][C]Correlations for all pairs of data series (method=pearson)[/C][/ROW]
[ROW][C] [/C][C]logins[/C][C]blogged_computations[/C][C]reviewed_compendiums[/C][C]long_feedback_messages[/C][C]shared_compendiums[/C][C]number_characters[/C][/ROW]
[ROW][C]logins[/C][C]1[/C][C]0.552[/C][C]0.589[/C][C]0.556[/C][C]0.196[/C][C]0.527[/C][/ROW]
[ROW][C]blogged_computations[/C][C]0.552[/C][C]1[/C][C]0.743[/C][C]0.75[/C][C]0.123[/C][C]0.709[/C][/ROW]
[ROW][C]reviewed_compendiums[/C][C]0.589[/C][C]0.743[/C][C]1[/C][C]0.913[/C][C]0.249[/C][C]0.684[/C][/ROW]
[ROW][C]long_feedback_messages[/C][C]0.556[/C][C]0.75[/C][C]0.913[/C][C]1[/C][C]0.288[/C][C]0.709[/C][/ROW]
[ROW][C]shared_compendiums[/C][C]0.196[/C][C]0.123[/C][C]0.249[/C][C]0.288[/C][C]1[/C][C]0.288[/C][/ROW]
[ROW][C]number_characters[/C][C]0.527[/C][C]0.709[/C][C]0.684[/C][C]0.709[/C][C]0.288[/C][C]1[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155518&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=155518&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)
loginsblogged_computationsreviewed_compendiumslong_feedback_messagesshared_compendiumsnumber_characters
logins10.5520.5890.5560.1960.527
blogged_computations0.55210.7430.750.1230.709
reviewed_compendiums0.5890.74310.9130.2490.684
long_feedback_messages0.5560.750.91310.2880.709
shared_compendiums0.1960.1230.2490.28810.288
number_characters0.5270.7090.6840.7090.2881







Correlations for all pairs of data series with p-values
pairPearson rSpearman rhoKendall tau
logins;blogged_computations0.55160.65930.4842
p-value(0)(0)(0)
logins;reviewed_compendiums0.58910.58690.4355
p-value(0)(0)(0)
logins;long_feedback_messages0.55590.55070.407
p-value(0)(0)(0)
logins;shared_compendiums0.19650.21960.1579
p-value(0.0117)(0.0047)(0.005)
logins;number_characters0.5270.58980.4296
p-value(0)(0)(0)
blogged_computations;reviewed_compendiums0.74270.64690.4932
p-value(0)(0)(0)
blogged_computations;long_feedback_messages0.75020.66750.4976
p-value(0)(0)(0)
blogged_computations;shared_compendiums0.12310.0940.0687
p-value(0.1162)(0.2311)(0.2221)
blogged_computations;number_characters0.70890.68380.5228
p-value(0)(0)(0)
reviewed_compendiums;long_feedback_messages0.91340.82150.6828
p-value(0)(0)(0)
reviewed_compendiums;shared_compendiums0.24880.18530.1326
p-value(0.0013)(0.0175)(0.0202)
reviewed_compendiums;number_characters0.68410.50990.375
p-value(0)(0)(0)
long_feedback_messages;shared_compendiums0.28750.2650.1931
p-value(2e-04)(6e-04)(6e-04)
long_feedback_messages;number_characters0.70920.63290.4736
p-value(0)(0)(0)
shared_compendiums;number_characters0.28810.30370.2273
p-value(2e-04)(1e-04)(1e-04)

\begin{tabular}{lllllllll}
\hline
Correlations for all pairs of data series with p-values \tabularnewline
pair & Pearson r & Spearman rho & Kendall tau \tabularnewline
logins;blogged_computations & 0.5516 & 0.6593 & 0.4842 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
logins;reviewed_compendiums & 0.5891 & 0.5869 & 0.4355 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
logins;long_feedback_messages & 0.5559 & 0.5507 & 0.407 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
logins;shared_compendiums & 0.1965 & 0.2196 & 0.1579 \tabularnewline
p-value & (0.0117) & (0.0047) & (0.005) \tabularnewline
logins;number_characters & 0.527 & 0.5898 & 0.4296 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
blogged_computations;reviewed_compendiums & 0.7427 & 0.6469 & 0.4932 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
blogged_computations;long_feedback_messages & 0.7502 & 0.6675 & 0.4976 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
blogged_computations;shared_compendiums & 0.1231 & 0.094 & 0.0687 \tabularnewline
p-value & (0.1162) & (0.2311) & (0.2221) \tabularnewline
blogged_computations;number_characters & 0.7089 & 0.6838 & 0.5228 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
reviewed_compendiums;long_feedback_messages & 0.9134 & 0.8215 & 0.6828 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
reviewed_compendiums;shared_compendiums & 0.2488 & 0.1853 & 0.1326 \tabularnewline
p-value & (0.0013) & (0.0175) & (0.0202) \tabularnewline
reviewed_compendiums;number_characters & 0.6841 & 0.5099 & 0.375 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
long_feedback_messages;shared_compendiums & 0.2875 & 0.265 & 0.1931 \tabularnewline
p-value & (2e-04) & (6e-04) & (6e-04) \tabularnewline
long_feedback_messages;number_characters & 0.7092 & 0.6329 & 0.4736 \tabularnewline
p-value & (0) & (0) & (0) \tabularnewline
shared_compendiums;number_characters & 0.2881 & 0.3037 & 0.2273 \tabularnewline
p-value & (2e-04) & (1e-04) & (1e-04) \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=155518&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]logins;blogged_computations[/C][C]0.5516[/C][C]0.6593[/C][C]0.4842[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]logins;reviewed_compendiums[/C][C]0.5891[/C][C]0.5869[/C][C]0.4355[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]logins;long_feedback_messages[/C][C]0.5559[/C][C]0.5507[/C][C]0.407[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]logins;shared_compendiums[/C][C]0.1965[/C][C]0.2196[/C][C]0.1579[/C][/ROW]
[ROW][C]p-value[/C][C](0.0117)[/C][C](0.0047)[/C][C](0.005)[/C][/ROW]
[ROW][C]logins;number_characters[/C][C]0.527[/C][C]0.5898[/C][C]0.4296[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]blogged_computations;reviewed_compendiums[/C][C]0.7427[/C][C]0.6469[/C][C]0.4932[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]blogged_computations;long_feedback_messages[/C][C]0.7502[/C][C]0.6675[/C][C]0.4976[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]blogged_computations;shared_compendiums[/C][C]0.1231[/C][C]0.094[/C][C]0.0687[/C][/ROW]
[ROW][C]p-value[/C][C](0.1162)[/C][C](0.2311)[/C][C](0.2221)[/C][/ROW]
[ROW][C]blogged_computations;number_characters[/C][C]0.7089[/C][C]0.6838[/C][C]0.5228[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]reviewed_compendiums;long_feedback_messages[/C][C]0.9134[/C][C]0.8215[/C][C]0.6828[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]reviewed_compendiums;shared_compendiums[/C][C]0.2488[/C][C]0.1853[/C][C]0.1326[/C][/ROW]
[ROW][C]p-value[/C][C](0.0013)[/C][C](0.0175)[/C][C](0.0202)[/C][/ROW]
[ROW][C]reviewed_compendiums;number_characters[/C][C]0.6841[/C][C]0.5099[/C][C]0.375[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]long_feedback_messages;shared_compendiums[/C][C]0.2875[/C][C]0.265[/C][C]0.1931[/C][/ROW]
[ROW][C]p-value[/C][C](2e-04)[/C][C](6e-04)[/C][C](6e-04)[/C][/ROW]
[ROW][C]long_feedback_messages;number_characters[/C][C]0.7092[/C][C]0.6329[/C][C]0.4736[/C][/ROW]
[ROW][C]p-value[/C][C](0)[/C][C](0)[/C][C](0)[/C][/ROW]
[ROW][C]shared_compendiums;number_characters[/C][C]0.2881[/C][C]0.3037[/C][C]0.2273[/C][/ROW]
[ROW][C]p-value[/C][C](2e-04)[/C][C](1e-04)[/C][C](1e-04)[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=155518&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=155518&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
logins;blogged_computations0.55160.65930.4842
p-value(0)(0)(0)
logins;reviewed_compendiums0.58910.58690.4355
p-value(0)(0)(0)
logins;long_feedback_messages0.55590.55070.407
p-value(0)(0)(0)
logins;shared_compendiums0.19650.21960.1579
p-value(0.0117)(0.0047)(0.005)
logins;number_characters0.5270.58980.4296
p-value(0)(0)(0)
blogged_computations;reviewed_compendiums0.74270.64690.4932
p-value(0)(0)(0)
blogged_computations;long_feedback_messages0.75020.66750.4976
p-value(0)(0)(0)
blogged_computations;shared_compendiums0.12310.0940.0687
p-value(0.1162)(0.2311)(0.2221)
blogged_computations;number_characters0.70890.68380.5228
p-value(0)(0)(0)
reviewed_compendiums;long_feedback_messages0.91340.82150.6828
p-value(0)(0)(0)
reviewed_compendiums;shared_compendiums0.24880.18530.1326
p-value(0.0013)(0.0175)(0.0202)
reviewed_compendiums;number_characters0.68410.50990.375
p-value(0)(0)(0)
long_feedback_messages;shared_compendiums0.28750.2650.1931
p-value(2e-04)(6e-04)(6e-04)
long_feedback_messages;number_characters0.70920.63290.4736
p-value(0)(0)(0)
shared_compendiums;number_characters0.28810.30370.2273
p-value(2e-04)(1e-04)(1e-04)



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