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 computationThu, 18 Dec 2014 16:35:15 +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/2014/Dec/18/t14189205747hhsfo0m62p01ff.htm/, Retrieved Fri, 17 May 2024 15:28:24 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=271126, Retrieved Fri, 17 May 2024 15:28:24 +0000
QR Codes:

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

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] [] [2014-12-18 16:35:15] [18673d63f90870b9c004059cd6229007] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
149	NA
139	NA
148	NA
158	NA
128	NA
224	NA
159	NA
105	NA
159	NA
167	NA
165	NA
159	NA
119	NA
176	NA
54	NA
91	NA
163	NA
124	NA
137	NA
121	NA
153	NA
148	NA
221	NA
188	NA
NA	149
244	NA
NA	148
92	NA
150	NA
153	NA
94	NA
156	NA
NA	132
161	NA
105	NA
97	NA
151	NA
131	NA
166	NA
NA	157
111	NA
145	NA
162	NA
163	NA
59	NA
187	NA
109	NA
90	NA
105	NA
NA	83
116	NA
42	NA
148	NA
155	NA
125	NA
NA	116
128	NA
138	NA
49	NA
NA	96
NA	164
162	NA
99	NA
202	NA
186	NA
66	NA
183	NA
214	NA
188	NA
104	NA
NA	177
126	NA
76	NA
99	NA
139	NA
162	NA
108	NA
159	NA
74	NA
110	NA
96	NA
116	NA
87	NA
97	NA
127	NA
106	NA
80	NA
74	NA
91	NA
133	NA
74	NA
114	NA
140	NA
95	NA
98	NA
121	NA
NA	126
98	NA
95	NA
110	NA
70	NA
102	NA
86	NA
130	NA
96	NA
102	NA
100	NA
94	NA
52	NA
98	NA
118	NA
99	NA
48	NA
50	NA
150	NA
NA	154
NA	109
68	NA
NA	194
158	NA
159	NA
67	NA
147	NA
39	NA
100	NA
111	NA
138	NA
101	NA
131	NA
NA	101
NA	114
165	NA
114	NA
111	NA
75	NA
82	NA
121	NA
32	NA
150	NA
117	NA
71	NA
165	NA
154	NA
126	NA
149	NA
145	NA
120	NA
109	NA
132	NA
172	NA
169	NA
114	NA
156	NA
172	NA
NA	68
89	NA
NA	167
113	NA
NA	115
78	NA
NA	118
87	NA
173	NA
2	NA
NA	162
49	NA
NA	122
NA	96
100	NA
82	NA
NA	100
115	NA
NA	141
NA	165
NA	165
110	NA
118	NA
158	NA
146	NA
49	NA
90	NA
121	NA
NA	155
104	NA
147	NA
NA	110
NA	108
NA	113
115	NA
61	NA
NA	60
109	NA
68	NA
111	NA
77	NA
NA	73
NA	151
89	NA
78	NA
110	NA
220	NA
65	NA
141	NA
117	NA
122	NA
63	NA
44	NA
52	NA
131	NA
101	NA
NA	42
152	NA
107	NA
77	NA
154	NA
103	NA
96	NA
175	NA
57	NA
NA	112
143	NA
49	NA
NA	110
131	NA
167	NA
56	NA
137	NA
86	NA
121	NA
149	NA
168	NA
140	NA
88	NA
168	NA
94	NA
51	NA
48	NA
145	NA
66	NA
85	NA
109	NA
63	NA
102	NA
NA	162
NA	86
114	NA
164	NA
119	NA
126	NA
132	NA
142	NA
83	NA
94	NA
81	NA
NA	166
NA	110
64	NA
93	NA
104	NA
105	NA
49	NA
NA	88
95	NA
102	NA
99	NA
63	NA
76	NA
109	NA
NA	117
57	NA
120	NA
NA	73
91	NA
NA	108
105	NA
117	NA
119	NA
31	NA

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'Gertrude Mary Cox' @ cox.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 & 3 seconds \tabularnewline
R Server & 'Gertrude Mary Cox' @ cox.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271126&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]'Gertrude Mary Cox' @ cox.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271126&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271126&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 time3 seconds
R Server'Gertrude Mary Cox' @ cox.wessa.net






Two Sample t-test (unpaired)
Mean of Sample 1115.286324786325
Mean of Sample 2122.340909090909
t-stat-1.07808421862461
df276
p-value0.281937496310305
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-19.9363505273513,5.82718191818268]
F-test to compare two variances
F-stat1.35510194750297
df233
p-value0.232951264289959
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.819976065117057,2.07798658248498]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 115.286324786325 \tabularnewline
Mean of Sample 2 & 122.340909090909 \tabularnewline
t-stat & -1.07808421862461 \tabularnewline
df & 276 \tabularnewline
p-value & 0.281937496310305 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-19.9363505273513,5.82718191818268] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 1.35510194750297 \tabularnewline
df & 233 \tabularnewline
p-value & 0.232951264289959 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [0.819976065117057,2.07798658248498] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271126&T=1

[TABLE]
[ROW][C]Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]115.286324786325[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]122.340909090909[/C][/ROW]
[ROW][C]t-stat[/C][C]-1.07808421862461[/C][/ROW]
[ROW][C]df[/C][C]276[/C][/ROW]
[ROW][C]p-value[/C][C]0.281937496310305[/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][-19.9363505273513,5.82718191818268][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]1.35510194750297[/C][/ROW]
[ROW][C]df[/C][C]233[/C][/ROW]
[ROW][C]p-value[/C][C]0.232951264289959[/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.819976065117057,2.07798658248498][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271126&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271126&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 1115.286324786325
Mean of Sample 2122.340909090909
t-stat-1.07808421862461
df276
p-value0.281937496310305
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-19.9363505273513,5.82718191818268]
F-test to compare two variances
F-stat1.35510194750297
df233
p-value0.232951264289959
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.819976065117057,2.07798658248498]






Welch Two Sample t-test (unpaired)
Mean of Sample 1115.286324786325
Mean of Sample 2122.340909090909
t-stat-1.19595147469897
df66.9034450887021
p-value0.235936995107758
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-18.8287934294087,4.71962482024012]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 115.286324786325 \tabularnewline
Mean of Sample 2 & 122.340909090909 \tabularnewline
t-stat & -1.19595147469897 \tabularnewline
df & 66.9034450887021 \tabularnewline
p-value & 0.235936995107758 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-18.8287934294087,4.71962482024012] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271126&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]115.286324786325[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]122.340909090909[/C][/ROW]
[ROW][C]t-stat[/C][C]-1.19595147469897[/C][/ROW]
[ROW][C]df[/C][C]66.9034450887021[/C][/ROW]
[ROW][C]p-value[/C][C]0.235936995107758[/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][-18.8287934294087,4.71962482024012][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271126&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=271126&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 1115.286324786325
Mean of Sample 2122.340909090909
t-stat-1.19595147469897
df66.9034450887021
p-value0.235936995107758
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-18.8287934294087,4.71962482024012]






Wicoxon rank sum test with continuity correction (unpaired)
W4522
p-value0.201063685006759
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.171717171717172
p-value0.224827002572769
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0990675990675991
p-value0.860422974820543

\begin{tabular}{lllllllll}
\hline
Wicoxon rank sum test with continuity correction (unpaired) \tabularnewline
W & 4522 \tabularnewline
p-value & 0.201063685006759 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
Kolmogorov-Smirnov Test to compare Distributions of two Samples \tabularnewline
KS Statistic & 0.171717171717172 \tabularnewline
p-value & 0.224827002572769 \tabularnewline
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples \tabularnewline
KS Statistic & 0.0990675990675991 \tabularnewline
p-value & 0.860422974820543 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=271126&T=3

[TABLE]
[ROW][C]Wicoxon rank sum test with continuity correction (unpaired)[/C][/ROW]
[ROW][C]W[/C][C]4522[/C][/ROW]
[ROW][C]p-value[/C][C]0.201063685006759[/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.171717171717172[/C][/ROW]
[ROW][C]p-value[/C][C]0.224827002572769[/C][/ROW]
[ROW][C]Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples[/C][/ROW]
[ROW][C]KS Statistic[/C][C]0.0990675990675991[/C][/ROW]
[ROW][C]p-value[/C][C]0.860422974820543[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=271126&T=3

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

Wicoxon rank sum test with continuity correction (unpaired)
W4522
p-value0.201063685006759
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.171717171717172
p-value0.224827002572769
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0990675990675991
p-value0.860422974820543


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):
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)
a<-table.element(a,paste('Wicoxon rank sum test 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')