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, 01 Dec 2014 18:25:59 +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/01/t1417458394hnsj5q3s8e5lomt.htm/, Retrieved Thu, 16 May 2024 13:20:33 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=262167, Retrieved Thu, 16 May 2024 13:20:33 +0000
QR Codes:

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

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 intr motiv M-V] [2014-12-01 18:25:59] [7919944b2c0818d4401807e8f8057775] [Current]
Feedback Forum

Post a new message

Dataset

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

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=262167&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=262167&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=262167&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.075
Mean of Sample 253.4285714285714
t-stat-0.40633876993642
df558
p-value0.684649491313667
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.06272077614445,1.35557791900159]
F-test to compare two variances
F-stat0.935561567949137
df279
p-value0.578388873641405
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.739509244652338,1.18358959506306]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 53.075 \tabularnewline
Mean of Sample 2 & 53.4285714285714 \tabularnewline
t-stat & -0.40633876993642 \tabularnewline
df & 558 \tabularnewline
p-value & 0.684649491313667 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-2.06272077614445,1.35557791900159] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 0.935561567949137 \tabularnewline
df & 279 \tabularnewline
p-value & 0.578388873641405 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [0.739509244652338,1.18358959506306] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=262167&T=1

[TABLE]
[ROW][C]Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]53.075[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]53.4285714285714[/C][/ROW]
[ROW][C]t-stat[/C][C]-0.40633876993642[/C][/ROW]
[ROW][C]df[/C][C]558[/C][/ROW]
[ROW][C]p-value[/C][C]0.684649491313667[/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.06272077614445,1.35557791900159][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]0.935561567949137[/C][/ROW]
[ROW][C]df[/C][C]279[/C][/ROW]
[ROW][C]p-value[/C][C]0.578388873641405[/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.739509244652338,1.18358959506306][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=262167&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=262167&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.075
Mean of Sample 253.4285714285714
t-stat-0.40633876993642
df558
p-value0.684649491313667
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.06272077614445,1.35557791900159]
F-test to compare two variances
F-stat0.935561567949137
df279
p-value0.578388873641405
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.739509244652338,1.18358959506306]






Welch Two Sample t-test (unpaired)
Mean of Sample 153.075
Mean of Sample 253.4285714285714
t-stat-0.40633876993642
df557.382226799207
p-value0.684649663937371
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.06272489378672,1.35558203664387]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 53.075 \tabularnewline
Mean of Sample 2 & 53.4285714285714 \tabularnewline
t-stat & -0.40633876993642 \tabularnewline
df & 557.382226799207 \tabularnewline
p-value & 0.684649663937371 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-2.06272489378672,1.35558203664387] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=262167&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]53.075[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]53.4285714285714[/C][/ROW]
[ROW][C]t-stat[/C][C]-0.40633876993642[/C][/ROW]
[ROW][C]df[/C][C]557.382226799207[/C][/ROW]
[ROW][C]p-value[/C][C]0.684649663937371[/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.06272489378672,1.35558203664387][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=262167&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=262167&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.075
Mean of Sample 253.4285714285714
t-stat-0.40633876993642
df557.382226799207
p-value0.684649663937371
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-2.06272489378672,1.35558203664387]






Wicoxon rank sum test with continuity correction (unpaired)
W39269
p-value0.971437479364924
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.0571428571428571
p-value0.750526292621831
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.1
p-value0.121592776881046

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=262167&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)
W39269
p-value0.971437479364924
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.0571428571428571
p-value0.750526292621831
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.1
p-value0.121592776881046


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