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_twosampletests_mean.wasp
Title produced by softwarePaired and Unpaired Two Samples Tests about the Mean
Date of computationMon, 11 Dec 2017 14:51:26 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2017/Dec/11/t1513000652e9e9nsidodkrufc.htm/, Retrieved Wed, 15 May 2024 21:49:15 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=308985, Retrieved Wed, 15 May 2024 21:49:15 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Paired and Unpaired Two Samples Tests about the Mean] [] [2017-12-11 13:51:26] [f1ade19563a25eb31edff11eb9af1158] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
22	39
38	40
31	34
37	38
38	39
29	39
35	34
30	32
30	36
40	33
31	34
42	45
36	33
33	34
33	27
34	33
30	42
29	21
40	43
41	34
33	32
34	28
30	27
30	29
34	42
25	39
39	35
31	33
28	33
35	44
28	35
39	39
25	39
35	33
41	34
32	36
34	34
25	35
39	34
31	36
40	40
34	31
33	33
32	42
28	38
32	34
34	35
31	25
31	32
28	34
26	33
27	38
31	37
35	38
30	36
30	33
34	32
35	36
27	38
31	36
33	27
30	28
31	30
33	29
35	29
31	31
38	35
27	42
31	38
36	34
33	28
32	30
36	33
36	34
36	28
31	29
31	32
32	34
36	35
37	40
34	34
35	28
36	31
31	33
35	36
36	30
38	30
28	25
28	39
36	36
32	32
38	43
37	36
30	42
37	26
33	27
33	32
34	36
36	25
36	40
36	35
36	31
39	31
35	36
40	37
22	26
21	35
38	33
19	31
37	38
24	32
36	28
34	33
30	31
39	34
37	33
37	36
32	36
37	29
29	35
37	31
30	35
33	28
28	34
34	31
29	44
32	36
33	34
27	36
34	28
26	32
39	36
29	38
41	43
43	26
31	33
33	35
34	25
23	26
29	16
27	14
30	38
27	27
29	40
33	40
32	29
33	27
30	26
27	29
27	26
35	27
36	35
32	39
25	38
36	36
31	37
37	32
35	33
32	39
37	39
37	36
31	33
38	35
32	36
28	36
25	41
28	36
26	39
30	32
32	36
31	43
28	28
39	30
43	39
31	34
34	35
27	38
34	40
28	34
32	34
39	39
32	26
39	30
23	34
25	34
32	30
39	35
32	40
34	36
32	34
31	45
38	22
36	24
36	25
31	26
41	35
35	35
37	36
31	37
NA	33
NA	35
NA	37
NA	35
NA	29
NA	35
NA	30
NA	36
NA	34
NA	36
NA	39
NA	40
NA	35
NA	38
NA	41
NA	40
NA	37
NA	37
NA	40
NA	34
NA	33
NA	37
NA	31
NA	32
NA	39
NA	28
NA	36
NA	31
NA	32
NA	36
NA	31
NA	28
NA	28
NA	38
NA	35
NA	26
NA	32
NA	28
NA	33
NA	38
NA	31
NA	43
NA	37
NA	28
NA	35
NA	34
NA	40
NA	38

Tables (Output of Computation)





Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time2 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=308985&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]2 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=308985&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=308985&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center






Two Sample t-test (unpaired)
Mean of Sample 132.7386934673367
Mean of Sample 233.8866396761134
t-stat-2.53159889842498
df444
p-value0.0116983381367876
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.03911546904119,-0.256776948512174]
F-test to compare two variances
F-stat0.837211093799009
df198
p-value0.192053285006006
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.643337583855182,1.09371753174562]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 32.7386934673367 \tabularnewline
Mean of Sample 2 & 33.8866396761134 \tabularnewline
t-stat & -2.53159889842498 \tabularnewline
df & 444 \tabularnewline
p-value & 0.0116983381367876 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-2.03911546904119,-0.256776948512174] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 0.837211093799009 \tabularnewline
df & 198 \tabularnewline
p-value & 0.192053285006006 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [0.643337583855182,1.09371753174562] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=308985&T=1

[TABLE]
[ROW][C]Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]32.7386934673367[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]33.8866396761134[/C][/ROW]
[ROW][C]t-stat[/C][C]-2.53159889842498[/C][/ROW]
[ROW][C]df[/C][C]444[/C][/ROW]
[ROW][C]p-value[/C][C]0.0116983381367876[/C][/ROW]
[ROW][C]H0 value[/C][C]0[/C][/ROW]
[ROW][C]Alternative[/C][C]two.sided[/C][/ROW]
[ROW][C]CI Level[/C][C]0.95[/C][/ROW]
[ROW][C]CI[/C][C][-2.03911546904119,-0.256776948512174][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]0.837211093799009[/C][/ROW]
[ROW][C]df[/C][C]198[/C][/ROW]
[ROW][C]p-value[/C][C]0.192053285006006[/C][/ROW]
[ROW][C]H0 value[/C][C]1[/C][/ROW]
[ROW][C]Alternative[/C][C]two.sided[/C][/ROW]
[ROW][C]CI Level[/C][C]0.95[/C][/ROW]
[ROW][C]CI[/C][C][0.643337583855182,1.09371753174562][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=308985&T=1

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

Two Sample t-test (unpaired)
Mean of Sample 132.7386934673367
Mean of Sample 233.8866396761134
t-stat-2.53159889842498
df444
p-value0.0116983381367876
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.03911546904119,-0.256776948512174]
F-test to compare two variances
F-stat0.837211093799009
df198
p-value0.192053285006006
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.643337583855182,1.09371753174562]






Welch Two Sample t-test (unpaired)
Mean of Sample 132.7386934673367
Mean of Sample 233.8866396761134
t-stat-2.55591114207235
df436.836328015048
p-value0.0109284914193781
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.03067808764367,-0.265214329909691]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 32.7386934673367 \tabularnewline
Mean of Sample 2 & 33.8866396761134 \tabularnewline
t-stat & -2.55591114207235 \tabularnewline
df & 436.836328015048 \tabularnewline
p-value & 0.0109284914193781 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-2.03067808764367,-0.265214329909691] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=308985&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]32.7386934673367[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]33.8866396761134[/C][/ROW]
[ROW][C]t-stat[/C][C]-2.55591114207235[/C][/ROW]
[ROW][C]df[/C][C]436.836328015048[/C][/ROW]
[ROW][C]p-value[/C][C]0.0109284914193781[/C][/ROW]
[ROW][C]H0 value[/C][C]0[/C][/ROW]
[ROW][C]Alternative[/C][C]two.sided[/C][/ROW]
[ROW][C]CI Level[/C][C]0.95[/C][/ROW]
[ROW][C]CI[/C][C][-2.03067808764367,-0.265214329909691][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=308985&T=2

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

As an alternative you can also use a QR Code:  
QR Code


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

Welch Two Sample t-test (unpaired)
Mean of Sample 132.7386934673367
Mean of Sample 233.8866396761134
t-stat-2.55591114207235
df436.836328015048
p-value0.0109284914193781
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.03067808764367,-0.265214329909691]






Wilcoxon Rank-Sum Test (Mann–Whitney U test) with continuity correction (unpaired)
W20919
p-value0.00675102316141609
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.151404797265681
p-value0.0127839621775114
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0905539844973857
p-value0.326715027113208

\begin{tabular}{lllllllll}
\hline
Wilcoxon Rank-Sum Test (Mann–Whitney U test) with continuity correction (unpaired) \tabularnewline
W & 20919 \tabularnewline
p-value & 0.00675102316141609 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
Kolmogorov-Smirnov Test to compare Distributions of two Samples \tabularnewline
KS Statistic & 0.151404797265681 \tabularnewline
p-value & 0.0127839621775114 \tabularnewline
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples \tabularnewline
KS Statistic & 0.0905539844973857 \tabularnewline
p-value & 0.326715027113208 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=308985&T=3

[TABLE]
[ROW][C]Wilcoxon Rank-Sum Test (Mann–Whitney U test) with continuity correction (unpaired)[/C][/ROW]
[ROW][C]W[/C][C]20919[/C][/ROW]
[ROW][C]p-value[/C][C]0.00675102316141609[/C][/ROW]
[ROW][C]H0 value[/C][C]0[/C][/ROW]
[ROW][C]Alternative[/C][C]two.sided[/C][/ROW]
[ROW][C]Kolmogorov-Smirnov Test to compare Distributions of two Samples[/C][/ROW]
[ROW][C]KS Statistic[/C][C]0.151404797265681[/C][/ROW]
[ROW][C]p-value[/C][C]0.0127839621775114[/C][/ROW]
[ROW][C]Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples[/C][/ROW]
[ROW][C]KS Statistic[/C][C]0.0905539844973857[/C][/ROW]
[ROW][C]p-value[/C][C]0.326715027113208[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=308985&T=3

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

As an alternative you can also use a QR Code:  
QR Code


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

Wilcoxon Rank-Sum Test (Mann–Whitney U test) with continuity correction (unpaired)
W20919
p-value0.00675102316141609
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.151404797265681
p-value0.0127839621775114
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0905539844973857
p-value0.326715027113208


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1 ; par2 = 2 ; par3 = 0.95 ; par4 = two.sided ; par5 = unpaired ; par6 = 0.0 ;
Parameters (R input):
par1 = 1 ; par2 = 2 ; par3 = 0.95 ; par4 = two.sided ; par5 = unpaired ; par6 = 0.0 ;
R code (references can be found in the software module):
par6 <- '0.0'
par5 <- 'unpaired'
par4 <- 'greater'
par3 <- '0.95'
par2 <- '2'
par1 <- '1'
par1 <- as.numeric(par1) #column number of first sample
par2 <- as.numeric(par2) #column number of second sample
par3 <- as.numeric(par3) #confidence (= 1 - alpha)
if (par5 == 'unpaired') paired <- FALSE else paired <- TRUE
par6 <- as.numeric(par6) #H0
z <- t(y)
if (par1 == par2) stop('Please, select two different column numbers')
if (par1 < 1) stop('Please, select a column number greater than zero for the first sample')
if (par2 < 1) stop('Please, select a column number greater than zero for the second sample')
if (par1 > length(z[1,])) stop('The column number for the first sample should be smaller')
if (par2 > length(z[1,])) stop('The column number for the second sample should be smaller')
if (par3 <= 0) stop('The confidence level should be larger than zero')
if (par3 >= 1) stop('The confidence level should be smaller than zero')
(r.t <- t.test(z[,par1],z[,par2],var.equal=TRUE,alternative=par4,paired=paired,mu=par6,conf.level=par3))
(v.t <- var.test(z[,par1],z[,par2],conf.level=par3))
(r.w <- t.test(z[,par1],z[,par2],var.equal=FALSE,alternative=par4,paired=paired,mu=par6,conf.level=par3))
(w.t <- wilcox.test(z[,par1],z[,par2],alternative=par4,paired=paired,mu=par6,conf.level=par3))
(ks.t <- ks.test(z[,par1],z[,par2],alternative=par4))
m1 <- mean(z[,par1],na.rm=T)
m2 <- mean(z[,par2],na.rm=T)
mdiff <- m1 - m2
newsam1 <- z[!is.na(z[,par1]),par1]
newsam2 <- z[,par2]+mdiff
newsam2 <- newsam2[!is.na(newsam2)]
(ks1.t <- ks.test(newsam1,newsam2,alternative=par4))
mydf <- data.frame(cbind(z[,par1],z[,par2]))
colnames(mydf) <- c('Variable 1','Variable 2')
bitmap(file='test1.png')
boxplot(mydf, notch=TRUE, ylab='value',main=main)
dev.off()
bitmap(file='test2.png')
qqnorm(z[,par1],main='Normal QQplot - Variable 1')
qqline(z[,par1])
dev.off()
bitmap(file='test3.png')
qqnorm(z[,par2],main='Normal QQplot - Variable 2')
qqline(z[,par2])
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,paste('Two Sample t-test (',par5,')',sep=''),2,TRUE)
a<-table.row.end(a)
if(!paired){
a<-table.row.start(a)
a<-table.element(a,'Mean of Sample 1',header=TRUE)
a<-table.element(a,r.t$estimate[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Mean of Sample 2',header=TRUE)
a<-table.element(a,r.t$estimate[[2]])
a<-table.row.end(a)
} else {
a<-table.row.start(a)
a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE)
a<-table.element(a,r.t$estimate)
a<-table.row.end(a)
}
a<-table.row.start(a)
a<-table.element(a,'t-stat',header=TRUE)
a<-table.element(a,r.t$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'df',header=TRUE)
a<-table.element(a,r.t$parameter[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,r.t$p.value)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'H0 value',header=TRUE)
a<-table.element(a,r.t$null.value[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Alternative',header=TRUE)
a<-table.element(a,r.t$alternative)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI Level',header=TRUE)
a<-table.element(a,attr(r.t$conf.int,'conf.level'))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI',header=TRUE)
a<-table.element(a,paste('[',r.t$conf.int[1],',',r.t$conf.int[2],']',sep=''))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'F-test to compare two variances',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'F-stat',header=TRUE)
a<-table.element(a,v.t$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'df',header=TRUE)
a<-table.element(a,v.t$parameter[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,v.t$p.value)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'H0 value',header=TRUE)
a<-table.element(a,v.t$null.value[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Alternative',header=TRUE)
a<-table.element(a,v.t$alternative)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI Level',header=TRUE)
a<-table.element(a,attr(v.t$conf.int,'conf.level'))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI',header=TRUE)
a<-table.element(a,paste('[',v.t$conf.int[1],',',v.t$conf.int[2],']',sep=''))
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,paste('Welch Two Sample t-test (',par5,')',sep=''),2,TRUE)
a<-table.row.end(a)
if(!paired){
a<-table.row.start(a)
a<-table.element(a,'Mean of Sample 1',header=TRUE)
a<-table.element(a,r.w$estimate[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Mean of Sample 2',header=TRUE)
a<-table.element(a,r.w$estimate[[2]])
a<-table.row.end(a)
} else {
a<-table.row.start(a)
a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE)
a<-table.element(a,r.w$estimate)
a<-table.row.end(a)
}
a<-table.row.start(a)
a<-table.element(a,'t-stat',header=TRUE)
a<-table.element(a,r.w$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'df',header=TRUE)
a<-table.element(a,r.w$parameter[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,r.w$p.value)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'H0 value',header=TRUE)
a<-table.element(a,r.w$null.value[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Alternative',header=TRUE)
a<-table.element(a,r.w$alternative)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI Level',header=TRUE)
a<-table.element(a,attr(r.w$conf.int,'conf.level'))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'CI',header=TRUE)
a<-table.element(a,paste('[',r.w$conf.int[1],',',r.w$conf.int[2],']',sep=''))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
myWlabel <- 'Wilcoxon Signed-Rank Test'
if (par5=='unpaired') myWlabel = 'Wilcoxon Rank-Sum Test (Mann–Whitney U test)'
a<-table.element(a,paste(myWlabel,' with continuity correction (',par5,')',sep=''),2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'W',header=TRUE)
a<-table.element(a,w.t$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,w.t$p.value)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'H0 value',header=TRUE)
a<-table.element(a,w.t$null.value[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Alternative',header=TRUE)
a<-table.element(a,w.t$alternative)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Kolmogorov-Smirnov Test to compare Distributions of two Samples',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'KS Statistic',header=TRUE)
a<-table.element(a,ks.t$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,ks.t$p.value)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'KS Statistic',header=TRUE)
a<-table.element(a,ks1.t$statistic[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,ks1.t$p.value)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable2.tab')