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_Tests to Compare Two Means.wasp
Title produced by softwareT-Tests
Date of computationFri, 01 Jun 2012 06:11:55 -0400
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Jun/01/t1338545539vybwqaacgg01r8y.htm/, Retrieved Wed, 01 May 2024 23:15:45 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=168594, Retrieved Wed, 01 May 2024 23:15:45 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Histogram, QQplot and Density] [Women measured we...] [2012-06-01 09:52:57] [b687f6e38eec57c340962174507f7a8f]
- R  D  [Histogram, QQplot and Density] [Women reported w] [2012-06-01 09:54:56] [b687f6e38eec57c340962174507f7a8f]
-    D    [Histogram, QQplot and Density] [women - measured ...] [2012-06-01 09:58:23] [b687f6e38eec57c340962174507f7a8f]
-    D      [Histogram, QQplot and Density] [women reported he...] [2012-06-01 09:59:53] [b687f6e38eec57c340962174507f7a8f]
-    D        [Histogram, QQplot and Density] [male measured ] [2012-06-01 10:01:59] [b687f6e38eec57c340962174507f7a8f]
-    D          [Histogram, QQplot and Density] [male height measure] [2012-06-01 10:04:16] [b687f6e38eec57c340962174507f7a8f]
-    D            [Histogram, QQplot and Density] [men reported weight] [2012-06-01 10:06:02] [b687f6e38eec57c340962174507f7a8f]
- RMPD                [T-Tests] [Women - measured ...] [2012-06-01 10:11:55] [088dec10f3cee2e79fb056e2eecc70f7] [Current]
- R  D                  [T-Tests] [men measured vs r...] [2012-06-01 10:15:43] [b687f6e38eec57c340962174507f7a8f]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
58	51
53	54
59	59
57	56
51	52
64	64
52	57
65	66
62	62
61	61
61	61
54	59
50	50
63	61
58	60
39	41
71	71
52	52
68	63
56	54
54	53
63	59
54	55
49	NA
54	56
75	75
56	57
66	65
78	75
60	59
64	63
64	62
52	51
62	61
55	54
56	57
50	50
50	NA
50	55
63	64
61	60
53	52
60	55
56	56
53	53
57	59
57	56
56	56
56	57
50	50
52	52
55	NA
55	55
47	47
45	45
62	63
53	51
52	51
57	55
64	64
59	55
55	57
76	77
62	62
68	68
55	55
52	56
47	NA
45	45
68	68
44	44
62	61
56	53
50	47
53	53
64	62
62	NA
52	53
53	55
54	55
64	66
55	55
55	55
59	55
70	67
57	58
47	47
47	45
55	NA
48	44
59	NA
58	NA
57	56
51	50
54	54
53	52
59	58
59	59
63	62
66	66
53	50
54	NA
60	NA
43	NA
63	59
56	56
60	55
58	54
50	49
59	59
51	51
62	61

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=168594&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'Gwilym Jenkins' @ jenkins.wessa.net






T-Test
Mean156.89286
Mean256.74257
T Statistic0.16116
P-value0.87212
Lower Confidence Limit-1.68799
Upper Confidence Limit1.98855

\begin{tabular}{lllllllll}
\hline
T-Test \tabularnewline
Mean1 & 56.89286 \tabularnewline
Mean2 & 56.74257 \tabularnewline
T Statistic & 0.16116 \tabularnewline
P-value & 0.87212 \tabularnewline
Lower Confidence Limit & -1.68799 \tabularnewline
Upper Confidence Limit & 1.98855 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168594&T=1

[TABLE]
[ROW][C]T-Test[/C][/ROW]
[ROW][C]Mean1[/C][C]56.89286[/C][/ROW]
[ROW][C]Mean2[/C][C]56.74257[/C][/ROW]
[ROW][C]T Statistic[/C][C]0.16116[/C][/ROW]
[ROW][C]P-value[/C][C]0.87212[/C][/ROW]
[ROW][C]Lower Confidence Limit[/C][C]-1.68799[/C][/ROW]
[ROW][C]Upper Confidence Limit[/C][C]1.98855[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168594&T=1

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

T-Test
Mean156.89286
Mean256.74257
T Statistic0.16116
P-value0.87212
Lower Confidence Limit-1.68799
Upper Confidence Limit1.98855






Standard Deviations
Variable 16.85941
Variable 26.73744

\begin{tabular}{lllllllll}
\hline
Standard Deviations \tabularnewline
Variable 1 & 6.85941 \tabularnewline
Variable 2 & 6.73744 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=168594&T=2

[TABLE]
[ROW][C]Standard Deviations[/C][/ROW]
[ROW][C]Variable 1[/C][C]6.85941[/C][/ROW]
[ROW][C]Variable 2[/C][C]6.73744[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=168594&T=2

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

Standard Deviations
Variable 16.85941
Variable 26.73744


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = two.sided ; par2 = 1 ; par3 = 2 ; par4 = T-Test ; par5 = unpaired ; par6 = 0.0 ; par7 = 0.95 ; par8 = TRUE ;
Parameters (R input):
par1 = two.sided ; par2 = 1 ; par3 = 2 ; par4 = T-Test ; par5 = unpaired ; par6 = 0.0 ; par7 = 0.95 ; par8 = TRUE ;
R code (references can be found in the software module):
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.character(par4)
par5 <- as.character(par5)
par6 <- as.numeric(par6)
par7 <- as.numeric(par7)
par8 <- as.logical(par8)
if ( par5 == 'unpaired') paired <- FALSE else paired <- TRUE
x <- t(y)
if(par8){
bitmap(file='test1.png')
(r<-boxplot(x ,xlab=xlab,ylab=ylab,main=main,notch=FALSE,col=2))
dev.off()
}
load(file='createtable')
if( par4 == 'Wilcoxon-Mann_Whitney'){
a<-table.start()
a <- table.row.start(a)
a <- table.element(a,'Wilcoxon Test',3,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'',1,TRUE)
a <- table.element(a,'Statistic',1,TRUE)
a <- table.element(a,'P-value',1,TRUE)
a <- table.row.end(a)
W <- wilcox.test(x[,par2],x[,par3],alternative=par1, paired = paired)
a<-table.row.start(a)
a<-table.element(a,'Wilcoxon Test',1,TRUE)
a<-table.element(a,W$statistic[[1]])
a<-table.element(a,round(W$p.value, digits=5) )
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
}
if( par4 == 'T-Test')
{
T <- t.test(x[,par2],x[,par3],alternative=par1, paired=paired, mu=par6, conf.level=par7)
a<-table.start()
a <- table.row.start(a)
a <- table.element(a,'T-Test',3,TRUE)
a <- table.row.end(a)
if(paired){
a <- table.row.start(a)
a <- table.element(a,'Difference: Mean1 - Mean2',1,TRUE)
a<-table.element(a,round(T$estimate, digits=5) )
a <- table.row.end(a)
}
if(!paired){
a <- table.row.start(a)
a <- table.element(a,'Mean1',1,TRUE)
a<-table.element(a,round(T$estimate[1], digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Mean2',1,TRUE)
a<-table.element(a,round(T$estimate[2], digits=5) )
a <- table.row.end(a)
}
a <- table.row.start(a)
a <- table.element(a,'T Statistic',1,TRUE)
a<-table.element(a,round(T$statistic, digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'P-value',1,TRUE)
a<-table.element(a,round(T$p.value, digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Lower Confidence Limit',1,TRUE)
a<-table.element(a,round(T$conf.int[1], digits=5) )
a <- table.row.end(a)
a<-table.row.start(a)
a <- table.element(a,'Upper Confidence Limit',1,TRUE)
a<-table.element(a,round(T$conf.int[2], digits=5) )
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,'Standard Deviations',3,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Variable 1',1,TRUE)
a<-table.element(a,round(sd(x[,par2], na.rm=TRUE), digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Variable 2',1,TRUE)
a<-table.element(a,round(sd(x[,par3], na.rm=TRUE), digits=5) )
a <- table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')