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 12:40:11 +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/t1418910633r1m2zry3cbfkkqr.htm/, Retrieved Fri, 17 May 2024 01:10:59 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=270931, Retrieved Fri, 17 May 2024 01:10:59 +0000
QR Codes:

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

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 12:40:11] [6ac057e9f6255a74ae39891d7e02481c] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
96	NA
75	NA
NA	70
88	NA
NA	114
NA	69
NA	176
114	NA
NA	121
NA	110
NA	158
NA	116
NA	181
NA	77
141	NA
35	NA
80	NA
NA	152
97	NA
NA	99
84	NA
NA	68
NA	101
107	NA
NA	88
NA	112
NA	171
NA	137
77	NA
NA	66
93	NA
105	NA
131	NA
NA	89
NA	102
NA	161
NA	120
NA	127
77	NA
NA	108
NA	85
168	NA
NA	48
NA	152
NA	75
NA	107
NA	62
121	NA
NA	124
NA	72
40	NA
NA	58
NA	97
NA	88
NA	126
NA	104
NA	148
NA	146
80	NA
NA	97
25	NA
NA	99
NA	118
58	NA
63	NA
NA	139
50	NA
NA	60
152	NA
NA	142
NA	94
66	NA
127	NA
67	NA
90	NA
NA	75
NA	96
128	NA
NA	41
146	NA
NA	69
186	NA
81	NA
NA	85
54	NA
46	NA
106	NA
NA	34
60	NA
NA	95
NA	57
62	NA
36	NA
56	NA
NA	54
NA	64
NA	76
98	NA
NA	88
35	NA
NA	102
NA	61
NA	80
NA	49
NA	78
90	NA
NA	45
NA	55
NA	96
43	NA
52	NA
60	NA
54	NA
51	NA
51	NA
NA	38
NA	41
NA	146
NA	182
NA	192
263	NA
NA	35
NA	439
214	NA
NA	341
58	NA
292	NA
NA	85
NA	200
NA	158
NA	199
NA	297
NA	227
NA	108
NA	86
302	NA
NA	148
NA	178
NA	120
NA	207
NA	157
NA	128
296	NA
NA	323
NA	79
NA	70
NA	146
NA	246
145	NA
196	NA
199	NA
NA	127
91	NA
153	NA
299	NA
NA	228
190	NA
NA	180
NA	212
269	NA
NA	130
NA	179
NA	243
190	NA
299	NA
121	NA
137	NA
NA	305
157	NA
NA	96
183	NA
NA	52
238	NA
NA	40
226	NA
190	NA
NA	214
145	NA
NA	119
NA	222
NA	222
NA	159
NA	165
249	NA
NA	125
122	NA
186	NA
148	NA
NA	274
172	NA
NA	84
168	NA
102	NA
106	NA
2	NA
NA	139
NA	95
NA	130
NA	72
141	NA
113	NA
NA	206
268	NA
175	NA
77	NA
125	NA
NA	255
NA	111
132	NA
211	NA
NA	92
76	NA
NA	171
NA	83
NA	119

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270931&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'George Udny Yule' @ yule.wessa.net






Two Sample t-test (unpaired)
Mean of Sample 1126.977777777778
Mean of Sample 2127.672
t-stat-0.0704354631710778
df213
p-value0.943913151353689
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-20.1223016487487,18.7338572043042]
F-test to compare two variances
F-stat1.10301214845436
df89
p-value0.610148049466359
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.753545114386745,1.63480486238868]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 126.977777777778 \tabularnewline
Mean of Sample 2 & 127.672 \tabularnewline
t-stat & -0.0704354631710778 \tabularnewline
df & 213 \tabularnewline
p-value & 0.943913151353689 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-20.1223016487487,18.7338572043042] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 1.10301214845436 \tabularnewline
df & 89 \tabularnewline
p-value & 0.610148049466359 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [0.753545114386745,1.63480486238868] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270931&T=1

[TABLE]
[ROW][C]Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]126.977777777778[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]127.672[/C][/ROW]
[ROW][C]t-stat[/C][C]-0.0704354631710778[/C][/ROW]
[ROW][C]df[/C][C]213[/C][/ROW]
[ROW][C]p-value[/C][C]0.943913151353689[/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][-20.1223016487487,18.7338572043042][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]1.10301214845436[/C][/ROW]
[ROW][C]df[/C][C]89[/C][/ROW]
[ROW][C]p-value[/C][C]0.610148049466359[/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.753545114386745,1.63480486238868][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270931&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270931&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 1126.977777777778
Mean of Sample 2127.672
t-stat-0.0704354631710778
df213
p-value0.943913151353689
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-20.1223016487487,18.7338572043042]
F-test to compare two variances
F-stat1.10301214845436
df89
p-value0.610148049466359
H0 value1
Alternativetwo.sided
CI Level0.95
CI[0.753545114386745,1.63480486238868]






Welch Two Sample t-test (unpaired)
Mean of Sample 1126.977777777778
Mean of Sample 2127.672
t-stat-0.0698733953247592
df186.17596906584
p-value0.944369423194088
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-20.2947181790081,18.9062737345637]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (unpaired) \tabularnewline
Mean of Sample 1 & 126.977777777778 \tabularnewline
Mean of Sample 2 & 127.672 \tabularnewline
t-stat & -0.0698733953247592 \tabularnewline
df & 186.17596906584 \tabularnewline
p-value & 0.944369423194088 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-20.2947181790081,18.9062737345637] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=270931&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (unpaired)[/C][/ROW]
[ROW][C]Mean of Sample 1[/C][C]126.977777777778[/C][/ROW]
[ROW][C]Mean of Sample 2[/C][C]127.672[/C][/ROW]
[ROW][C]t-stat[/C][C]-0.0698733953247592[/C][/ROW]
[ROW][C]df[/C][C]186.17596906584[/C][/ROW]
[ROW][C]p-value[/C][C]0.944369423194088[/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][-20.2947181790081,18.9062737345637][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=270931&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=270931&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 1126.977777777778
Mean of Sample 2127.672
t-stat-0.0698733953247592
df186.17596906584
p-value0.944369423194088
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-20.2947181790081,18.9062737345637]






Wicoxon rank sum test with continuity correction (unpaired)
W5503
p-value0.787151976311066
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.103555555555556
p-value0.628709540697964
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.103555555555556
p-value0.628709540697964

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

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

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


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
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')