## Free Statistics

of Irreproducible Research!

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, 11 Nov 2011 09:55:17 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2011/Nov/11/t13210233901b1q33pt5j4c7k2.htm/, Retrieved Sun, 16 Jun 2024 22:30:05 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=141356, Retrieved Sun, 16 Jun 2024 22:30:05 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact128
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Aston University Statistical Software] [Morning Sickness ...] [2009-11-16 16:26:06] [74be16979710d4c4e7c6647856088456]
- R     [Aston University Statistical Software] [Morning Sickness ...] [2009-11-16 17:22:16] [74be16979710d4c4e7c6647856088456]
-   P     [T-Tests] [Morning Sickness ...] [2010-11-09 11:12:43] [3fdd735c61ad38cbc9b3393dc997cdb7]
- RM        [T-Tests] [Morning Sickness ...] [2011-11-07 09:34:35] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R             [T-Tests] [week 6 scenario 2] [2011-11-11 14:55:17] [dc64b3d56ef70e149fd0278b1d06fde8] [Current]
Feedback Forum

Post a new message
Dataseries X:
68.95	70.18
70.84	69.94
70.26	69.61
70.93	71.39
70.71	72.54
70.73	70.48
70.20	70.48
69.76	71.68
71.46	70.59
70.31	71.61
69.24	70.56
70.94	72.74
68.98	72.63
71.16	71.58
69.96	70.14
69.66	71.22
70.20	70.72
69.99	71.37
71.97	69.73
68.29	71.29
67.67	73.19
70.07	70.64
70.89	69.30
69.66	71.76
71.12	70.28
70.60	69.77
70.09	70.37
70.55	71.70
71.63	70.85
69.30	71.07


 Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 1 seconds R Server 'George Udny Yule' @ yule.wessa.net R Engine error message Error in dimnames(data) <- dimnames : length of 'dimnames' [1] not equal to array extent Calls: array Execution halted 

\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 & 1 seconds \tabularnewline
R Server & 'George Udny Yule' @ yule.wessa.net \tabularnewline
R Engine error message & Error in dimnames(data) <- dimnames :
length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted
\tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=141356&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ yule.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]Error in dimnames(data) <- dimnames :
length of 'dimnames' [1] not equal to array extent
Calls: array
Execution halted
[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=141356&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=141356&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

 Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 1 seconds R Server 'George Udny Yule' @ yule.wessa.net R Engine error message Error in dimnames(data) <- dimnames : length of 'dimnames' [1] not equal to array extent Calls: array Execution halted 

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 <- TRUEx <- 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]), 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]), digits=5) )a <- table.row.end(a)a<-table.end(a)table.save(a,file='mytable1.tab')