Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_Tests to Compare Two Means.wasp
Title produced by softwareT-Tests
Date of computationWed, 02 Jun 2010 09:51:16 +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/2010/Jun/02/t1275472314xgxh339ina8dsih.htm/, Retrieved Thu, 28 Mar 2024 19:34:03 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=77083, Retrieved Thu, 28 Mar 2024 19:34:03 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact134
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Two-Way ANOVA] [two-way anova wit...] [2010-05-26 17:02:24] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R PD  [Two-Way ANOVA] [ANOVA with good l...] [2010-05-28 23:09:47] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R       [Variability] [ANOVA with better...] [2010-05-29 09:47:12] [98fd0e87c3eb04e0cc2efde01dbafab6]
- R         [Variability] [ANOVA with better...] [2010-05-29 09:54:40] [98fd0e87c3eb04e0cc2efde01dbafab6]
- RMPD          [T-Tests] [] [2010-06-02 09:51:16] [e819fbd06db9b5cbe8452a0c7030b310] [Current]
Feedback Forum

Post a new message
Dataseries X:
Weight at week 0	Weight at week 8	Treatment group	Weight loss
91.3	90.2	1	1.1
88.8	87.0	2	1.8
87.0	81.8	2	5.2
114.0	110.5	1	3.5
107.4	103.5	2	3.9
89.2	90.2	1	-1.0
84.1	84.0	1	0.1
90.0	87.6	2	2.4
95.4	95.0	1	0.4
97.3	93.6	1	3.7
108.1	100.9	2	7.2
101.5	101.9	1	-0.4
83.4	75.0	2	8.4
102.2	96.3	1	5.9
84.0	82.6	1	1.4
103.7	95.7	2	8.0
99.2	99.2	2	0.0
126.0	123.2	2	2.8
103.7	95.5	2	8.2
117.9	117.0	1	0.9
112.4	111.8	1	0.6
85.0	80.0	2	5.0
83.8	77.9	2	5.9
73.8	74.8	1	-1.0
67.7	66.1	1	1.6
106.9	103.7	2	3.2
102.5	102.0	1	0.5
81.5	78.9	2	2.6
79.5	82.2	1	-2.7
103.0	97.2	2	5.8
127.5	124.7	2	2.8
101.6	100.3	1	1.3
106.0	101.0	1	5.0
112.8	111.2	1	1.6
98.4	95.0	1	3.4
114.9	105.3	2	9.6
103.4	96.0	2	7.4




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server184.73.214.54 @ 184.73.214.54
R Engine error message
Error in as.vector(data) : object 'Weight' not found
Calls: array -> as.vector
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 & 0 seconds \tabularnewline
R Server & 184.73.214.54 @ 184.73.214.54 \tabularnewline
R Engine error message & 
Error in as.vector(data) : object 'Weight' not found
Calls: array -> as.vector
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=77083&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]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]184.73.214.54 @ 184.73.214.54[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in as.vector(data) : object 'Weight' not found
Calls: array -> as.vector
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=77083&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=77083&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server184.73.214.54 @ 184.73.214.54
R Engine error message
Error in as.vector(data) : object 'Weight' not found
Calls: array -> as.vector
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 <- 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]), 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')