BPS651 Exercise V

BPS651 Research Methodology
Laboratory Exercise V
Course BPS651 Research Methodology
R.S. Rajput
Assistant Professor Computer Science
T Test
A t-test is a statistical hypothesis test, in which the test statistic follows a Student's t-distribution
if the null hypothesis is supported. It can be used to determine if two sets of data are
significantly different from each other, and is most commonly applied when the test statistic
would follow a normal distribution if the value of a scaling term in the test statistic were known.
t-test Function in R
The R function t.test() can be used to perform both one and two sample t-tests on vectors of
data. The function contains a variety of options and can be called as follows:
t.test(x, y = NULL, alternative = c("two.sided", "less", "greater"), mu = 0, paired = FALSE,

var.equal = FALSE, conf.level = 0.95)
Here
x is a numeric vector of data values
y is an optional numeric vector of data values. If y is excluded, the function performs a onesample t-test on the data contained in x, if it is included it performs a two-sample t-tests
using both x and y.
The option mu provides a number indicating the true value of the mean (or difference in
means if you are performing a two sample test) under the null hypothesis.
The option alternative is a character string specifying the alternative hypothesis, and must
be one of the following: "two. sided" (which is the default), "greater" or "less" depending
on whether the alternative hypothesis is that the mean is different than, greater than or less
than mu, respectively. For example the following call:
> t.test(x, alternative = "less", mu = 10)
Department of Mathematics, Statistics & Computer Science, GBPUAT
Page 1
performs a one sample t-test on the data contained in x where the null hypothesis is that =10
and the alternative is that <10.
The option paired indicates whether or not you want a paired t-test (TRUE = yes and FALSE
= no). If you leave this option out it defaults to FALSE.
The option var.equal is a logical variable indicating whether or not to assume the two
variances as being equal when performing a two-sample t-test. If TRUE then the pooled
variance is used to estimate the variance otherwise the Welch (or Satterthwaite)
approximation to the degrees of freedom is used. If you leave this option out it defaults to
FALSE.
The option conf.level determines the confidence level of the reported confidence interval
One-sample t-test
Example: An outbreak of Salmonella-related illness was attributed to ice cream produced at a
certain factory. Scientists measured the level of Salmonella in 9 randomly sampled batches of
ice cream. The levels (in MPN/g) were: 0.593, 0.142, 0.329, 0.691, 0.231, 0.793, 0.519, 0.392,
0.418 is there evidence that the mean level of Salmonella in the ice cream is greater than 0.3
MPN/g?
Solution: Let be the mean level of Salmonella in all batches of ice cream. Here the hypothesis of
interest can be expressed as:
H0: = 0.3
H1: > 0.3
Hence, we will need to include the options alternative="greater", mu=0.3. Below is the relevant
R-code:
>x=c(0.593,0.142,0.329,0.691,0.231,0.793,0.519,0.392,0.418)
>t.test(x, alternative="greater", mu=0.3)
One Sample t-test

data: x
t = 2.2051, df = 8, p-value = 0.02927
alternative hypothesis: true mean is greater than 0.3
Page 2

95 percent confidence interval:
0.3245133
Inf
sample estimates:
mean of x
0.4564444
From the output we see that the p-value = 0.029. Hence, there is moderately strong evidence
that the mean Salmonella level in the ice cream is above 0.3 MPN/g.
Exercise-17
Ten pieces of cloth each of 100 square meters were selected and number of weaving defects
were counted as given below. Test whether average number of weaving defects on such a cloth
is less than 5. (Given that table value= t(9,0.05)=1.833).
4, 0, 3, 3, 2, 3, 5, 8, 6, 6.
Exercise-18
Daily protein intake ( in gm) by an adult during a fortnight was reported as 48.1, 52.1, 52.0, 48.4,
49.4, 52.0, 46.6, 49.8, 46.4, 50.2, 48.8, 46.5, 50.0, 51.4 and 48.0 .Can we say that on an average
daily protein intake by that adult is 50 gm ?(Given that table value= t(14, 0.05/2)=2.145).
Exercise-19
Over the years an instructor has observed that on an average students get 11.5 marks in first
pre-final examination. Marks secured by 16 students of a particular batch in the first pre-final
examination were 9.5, 10.0, 14.5, 13.5, 14.0, 15.0, 10.5, 12.2, 11.0, 10.5, 14.0, 12.0, 13.5, 12.0,
14.0 and 15.0 . Can we conclude that this is a superior batch ? (Given that table value= t(15,
0.05) = 1.753).
Exercise-20
The weight of a canned food product is specified as 500 gm. For a random sample of 20 cans the
weights were observed as 480, 475, 510, 500, 505, 495, 498, 504, 490, 485, 505, 492, 490, 495,
515, 504, 486, 494, 496 and 491. Test whether on an average the weight is as per specification.
(Given that table value= t(19,0.05/2)=2.093) .
Two-sample t-tests
Example: 6 subjects were given a drug (treatment group) and an additional 6 subjects a placebo
(control group). Their reaction time to a stimulus was measured (in ms). We want to perform a
two-sample t-test for comparing the means of the treatment and control groups.
Let 1 be the mean of the population taking medicine and 2 the mean of the untreated
population. Here the hypothesis of interest can be expressed as:
H0: 1-2=0
Ha: 1-2<0
Page 3

Here we will need to include the data for the treatment group in x and the data for the control
group in y. We will also need to include the options alternative="less", mu=0.

Finally, we need to decide whether or not the standard deviations are the same in both groups.
Below is the relevant R-code when assuming equal standard deviation:
> Control = c(91, 87, 99, 77, 88, 91)

> Treat = c(101, 110, 103, 93, 99, 104)
> t.test(Control,Treat,alternative="less", var.equal=TRUE)
Two Sample t-test
data: Control and Treat
t = -3.4456, df = 10, p-value = 0.003136
alternative hypothesis: true difference in means is less than 0
-Inf -6.082744
sample estimates:
mean of x mean of y
88.83333 101.66667
Below is the relevant R-code when not assuming equal standard deviation:
> t.test(Control,Treat,alternative="less")
Welch Two Sample t-test
data: Control and Treat
t = -3.4456, df = 9.48, p-value = 0.003391
alternative hypothesis: true difference in means is less than 0
-Inf -6.044949
sample estimates:
mean of x mean of y
88.83333 101.66667
Here the pooled t-test and the Welsh t-test give roughly the same results (p-value = 0.00313 and
0.00339, respectively).
Exercise-21
Tensile strength of threads manufactured by firm A for a random sample of 12 pieces of threads
were measured as 2.21, 2.16, 2.17, 2.34, 2.22, 2.20, 2.18, 2.06, 2.26,2.14, 2.28, 2.24 while that
for threads of firm B for a random sample of 10 pieces were measured as 2.46, 2.42, 2.38, 2.26,
2.46, 2.38, 2.37, 2.35, 2.43, 2.41. Test whether average tensile strength of threads of these firms
is same. (Given that t(20,0.05/2)=2.086).
Exercise-22
Page 4

In a survey expenditure( in Rs.) on bakery products by a number of families from rural and urban
areas were recorded as given below. Can we conclude from the data that on an average urban
families tend to spend more than rural families? (Given that table value = t(28, 0.05) = 1.701).
Rural families 152, 144, 162, 182, 176, 175, 120, 128, 140, 210, 160, 180, 190, 195,182, 195,
175, 180, 162, 152.
Urban families 210, 220, 190, 310, 342, 360, 296, 314, 308, 300.
Exercise-23
It is felt that students from Public schools fair better than those from Government schools. To
verify this students were selected from both the schools and their marks were obtained as given
below. Draw the conclusion.(Given that table value=t(24,0.05)=1.711).
Public school 88 75 64 58 85 70 86 58 66 74 70 58
Government school- 66 76 58 54 80 68 58 84 70 56 52 60 58 48
Page 5

BPS651 Exercise V

Uploaded by

Copyright:

Available Formats

BPS651 Exercise V

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

BPS651 Exercise V

Uploaded by

Copyright:

Available Formats

BPS651 Research Methodology

t.test(x, y = NULL, alternative = c("two.sided", "less", "greater"), mu = 0, paired = FALSE,

x is a numeric vector of data values

> t.test(x, alternative = "less", mu = 10)

Department of Mathematics, Statistics & Computer Science, GBPUAT

BPS651 Research Methodology

One Sample t-test

Department of Mathematics, Statistics & Computer Science, GBPUAT

BPS651 Research Methodology

Department of Mathematics, Statistics & Computer Science, GBPUAT

BPS651 Research Methodology

group in y. We will also need to include the options alternative="less", mu=0.

> Control = c(91, 87, 99, 77, 88, 91)

Department of Mathematics, Statistics & Computer Science, GBPUAT

BPS651 Research Methodology

Department of Mathematics, Statistics & Computer Science, GBPUAT

You might also like