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, 08 Dec 2014 20:30:38 +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/08/t1418070699duomkbyu9zl6aen.htm/, Retrieved Sun, 19 May 2024 08:46:34 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=264220, Retrieved Sun, 19 May 2024 08:46:34 +0000
QR Codes:

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

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-08 20:30:38] [003c997d057e54927bd887526d955d96] [Current]
Feedback Forum

Post a new message

Dataset

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

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=264220&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 (paired)
Difference: Mean1 - Mean2-12.9383983572895
t-stat-27.172974708868
df486
p-value1.43687306611334e-99
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-13.873964120481,-12.0028325940981]
F-test to compare two variances
F-stat1.55434704450497
df486
p-value1.30065010983671e-06
H0 value1
Alternativetwo.sided
CI Level0.95
CI[1.30087257510972,1.85721090673125]

\begin{tabular}{lllllllll}
\hline
Two Sample t-test (paired) \tabularnewline
Difference: Mean1 - Mean2 & -12.9383983572895 \tabularnewline
t-stat & -27.172974708868 \tabularnewline
df & 486 \tabularnewline
p-value & 1.43687306611334e-99 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-13.873964120481,-12.0028325940981] \tabularnewline
F-test to compare two variances \tabularnewline
F-stat & 1.55434704450497 \tabularnewline
df & 486 \tabularnewline
p-value & 1.30065010983671e-06 \tabularnewline
H0 value & 1 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [1.30087257510972,1.85721090673125] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=264220&T=1

[TABLE]
[ROW][C]Two Sample t-test (paired)[/C][/ROW]
[ROW][C]Difference: Mean1 - Mean2[/C][C]-12.9383983572895[/C][/ROW]
[ROW][C]t-stat[/C][C]-27.172974708868[/C][/ROW]
[ROW][C]df[/C][C]486[/C][/ROW]
[ROW][C]p-value[/C][C]1.43687306611334e-99[/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][-13.873964120481,-12.0028325940981][/C][/ROW]
[ROW][C]F-test to compare two variances[/C][/ROW]
[ROW][C]F-stat[/C][C]1.55434704450497[/C][/ROW]
[ROW][C]df[/C][C]486[/C][/ROW]
[ROW][C]p-value[/C][C]1.30065010983671e-06[/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][1.30087257510972,1.85721090673125][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=264220&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=264220&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 (paired)
Difference: Mean1 - Mean2-12.9383983572895
t-stat-27.172974708868
df486
p-value1.43687306611334e-99
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-13.873964120481,-12.0028325940981]
F-test to compare two variances
F-stat1.55434704450497
df486
p-value1.30065010983671e-06
H0 value1
Alternativetwo.sided
CI Level0.95
CI[1.30087257510972,1.85721090673125]






Welch Two Sample t-test (paired)
Difference: Mean1 - Mean2-12.9383983572895
t-stat-27.172974708868
df486
p-value1.43687306611334e-99
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-13.873964120481,-12.0028325940981]

\begin{tabular}{lllllllll}
\hline
Welch Two Sample t-test (paired) \tabularnewline
Difference: Mean1 - Mean2 & -12.9383983572895 \tabularnewline
t-stat & -27.172974708868 \tabularnewline
df & 486 \tabularnewline
p-value & 1.43687306611334e-99 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
CI Level & 0.95 \tabularnewline
CI & [-13.873964120481,-12.0028325940981] \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=264220&T=2

[TABLE]
[ROW][C]Welch Two Sample t-test (paired)[/C][/ROW]
[ROW][C]Difference: Mean1 - Mean2[/C][C]-12.9383983572895[/C][/ROW]
[ROW][C]t-stat[/C][C]-27.172974708868[/C][/ROW]
[ROW][C]df[/C][C]486[/C][/ROW]
[ROW][C]p-value[/C][C]1.43687306611334e-99[/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][-13.873964120481,-12.0028325940981][/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=264220&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=264220&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 (paired)
Difference: Mean1 - Mean2-12.9383983572895
t-stat-27.172974708868
df486
p-value1.43687306611334e-99
H0 value0
Alternativetwo.sided
CI Level0.95
CI[-13.873964120481,-12.0028325940981]






Wicoxon rank sum test with continuity correction (paired)
W4693
p-value1.64779629675312e-68
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.5564681724846
p-value0
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0780287474332649
p-value0.103097280477303

\begin{tabular}{lllllllll}
\hline
Wicoxon rank sum test with continuity correction (paired) \tabularnewline
W & 4693 \tabularnewline
p-value & 1.64779629675312e-68 \tabularnewline
H0 value & 0 \tabularnewline
Alternative & two.sided \tabularnewline
Kolmogorov-Smirnov Test to compare Distributions of two Samples \tabularnewline
KS Statistic & 0.5564681724846 \tabularnewline
p-value & 0 \tabularnewline
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples \tabularnewline
KS Statistic & 0.0780287474332649 \tabularnewline
p-value & 0.103097280477303 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=264220&T=3

[TABLE]
[ROW][C]Wicoxon rank sum test with continuity correction (paired)[/C][/ROW]
[ROW][C]W[/C][C]4693[/C][/ROW]
[ROW][C]p-value[/C][C]1.64779629675312e-68[/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.5564681724846[/C][/ROW]
[ROW][C]p-value[/C][C]0[/C][/ROW]
[ROW][C]Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples[/C][/ROW]
[ROW][C]KS Statistic[/C][C]0.0780287474332649[/C][/ROW]
[ROW][C]p-value[/C][C]0.103097280477303[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=264220&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=264220&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 (paired)
W4693
p-value1.64779629675312e-68
H0 value0
Alternativetwo.sided
Kolmogorov-Smirnov Test to compare Distributions of two Samples
KS Statistic0.5564681724846
p-value0
Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples
KS Statistic0.0780287474332649
p-value0.103097280477303


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 1 ; par2 = 2 ; par3 = 0.95 ; par4 = two.sided ; par5 = paired ; par6 = 0.0 ;
Parameters (R input):
par1 = 1 ; par2 = 2 ; par3 = 0.95 ; par4 = two.sided ; par5 = paired ; 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')