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 computationSun, 07 Dec 2014 19:44:35 +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/07/t14179816059t9gpaxf536oxan.htm/, Retrieved Thu, 16 May 2024 09:49:10 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=263890, Retrieved Thu, 16 May 2024 09:49:10 +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] [Paper BS intr mot...] [2014-12-07 19:44:35] [7919944b2c0818d4401807e8f8057775] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
26	62
51	56
57	57
37	51
67	56
43	30
52	61
52	47
43	56
84	50
67	67
49	41
70	45
52	48
58	44
68	37
43	56
56	66
74	38
65	34
63	49
58	55
57	49
63	59
53	40
64	58
53	60
29	63
54	56
51	54
58	52
43	34
51	69
53	32
54	48
61	67
47	58
39	57
48	42
50	64
35	58
68	66
49	26
67	61
43	52
62	51
57	55
54	50
61	60
56	56
41	63
43	61
53	52
66	55
58	72
46	33
51	66
51	66
45	64
37	40
59	46
42	58
66	51
53	50
52	52
16	54
46	66
56	61
50	80
59	51
60	56
52	53
44	47
67	50
52	39
55	58
37	35
54	58
51	60
48	62
60	63
50	53
63	46
33	67
67	59
46	64
54	38
59	50
61	48
47	47
69	66
52	63
55	44
55	43
41	38
73	56
51	45
52	50
50	54
51	55
60	37
56	46
56	51
29	64
73	47
55	62
43	67
61	56
56	65
56	50
47	57
25	47
46	47
51	57
48	50
47	22
58	59
51	56
55	53
57	42
60	52
56	54
49	44
43	62
59	53
58	50
53	36
48	76
51	66
59	62
62	59
51	47
64	55
52	58
50	60
54	57
58	45
63	49
31	62
71	56
54	60
43	67
41	52
63	52
63	53
56	45
51	47
41	41
66	53
44	34
58	45
51	44
57	60
30	53
46	53
51	51
56	65
58	51
44	49
14	58
53	62
42	52
44	50
30	53
46	62
50	66
54	50
48	58
55	53
35	59
55	58
41	52
59	58
54	71
66	58
55	46
45	64
51	67
47	44
42	69
53	64
53	38
41	59
55	47
55	57
46	51
63	67
43	43
65	41
59	58
39	64
44	50
57	59
69	55
46	59
46	58
40	41
70	77
54	58
77	62
60	60
50	56
66	43
60	54
51	54
69	56
60	65
58	66
39	62
51	67
52	53
49	49
63	56
51	76
52	33
52	72
31	51
61	42
54	69
72	51
65	51
56	67
63	64
45	58
52	
68	
45	
70	
69	
46	
39	
54	
41	
68	
63	
57	
61	
39	
59	
51	
51	
65	
50	
21	
47	
37	
58	
51	
40	
64	
58	
56	
63	
60	
64	
47	
46	
50	
46	
44	
58	
58	
25	
56	
56	
59	
46	
49	
53	
58	
54	
52	
59	
53

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263890&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 time2 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net






Two Sample t-test (unpaired)
Mean of Sample 153.6642335766423
Mean of Sample 252.9598540145985
t-stat0.804433827383204
df546
p-value0.421496589257732
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-1.0156197812433,2.4243789053309]
F-test to compare two variances
F-stat0.875919808339166
df273
p-value0.274363230436035
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.69058119378051,1.11099971668903]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 53.6642335766423 \tabularnewline
Mean of Sample 2 & 52.9598540145985 \tabularnewline
t-stat & 0.804433827383204 \tabularnewline
df & 546 \tabularnewline
p-value & 0.421496589257732 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-1.0156197812433,2.4243789053309] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 0.875919808339166 \tabularnewline
df & 273 \tabularnewline
p-value & 0.274363230436035 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [0.69058119378051,1.11099971668903] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=263890&T=1

[TABLE]
[ROW][C]Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]53.6642335766423[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]52.9598540145985[/C][/ROW]
[ROW][C]t-stat[/C][C]0.804433827383204[/C][/ROW]
[ROW][C]df[/C][C]546[/C][/ROW]
[ROW][C]p-value[/C][C]0.421496589257732[/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][-1.0156197812433,2.4243789053309][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]0.875919808339166[/C][/ROW]
[ROW][C]df[/C][C]273[/C][/ROW]
[ROW][C]p-value[/C][C]0.274363230436035[/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.69058119378051,1.11099971668903][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=263890&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263890&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 153.6642335766423
Mean of Sample 252.9598540145985
t-stat0.804433827383204
df546
p-value0.421496589257732
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-1.0156197812433,2.4243789053309]
F-test to compare two variances
F-stat0.875919808339166
df273
p-value0.274363230436035
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.69058119378051,1.11099971668903]






Welch Two Sample t-test (unpaired)
Mean of Sample 153.6642335766423
Mean of Sample 252.9598540145985
t-stat0.804433827383204
df543.621664443432
p-value0.421498120420493
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-1.01563649838253,2.42439562247012]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 53.6642335766423 \tabularnewline
Mean of Sample 2 & 52.9598540145985 \tabularnewline
t-stat & 0.804433827383204 \tabularnewline
df & 543.621664443432 \tabularnewline
p-value & 0.421498120420493 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-1.01563649838253,2.42439562247012] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=263890&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]53.6642335766423[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]52.9598540145985[/C][/ROW]
[ROW][C]t-stat[/C][C]0.804433827383204[/C][/ROW]
[ROW][C]df[/C][C]543.621664443432[/C][/ROW]
[ROW][C]p-value[/C][C]0.421498120420493[/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][-1.01563649838253,2.42439562247012][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=263890&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263890&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 153.6642335766423
Mean of Sample 252.9598540145985
t-stat0.804433827383204
df543.621664443432
p-value0.421498120420493
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-1.01563649838253,2.42439562247012]






Wicoxon rank sum test with continuity correction (unpaired)
W39192
p-value0.371966083965661
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.0693430656934307
p-value0.525326021901978
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0693430656934307
p-value0.525326021901978

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263890&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)
W39192
p-value0.371966083965661
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.0693430656934307
p-value0.525326021901978
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0693430656934307
p-value0.525326021901978


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
Parameters (R input):
par1 = 2 ; par2 = 1 ; 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')