R: Interval Estimation of the Population Mean

Interval estimation of the population mean can be computed from functions of the following R packages:

stats – contains the t.test;

TeachingDemos – contains the z.test; and,

BSDA – contains the zsum.test and tsum.test.

The t.test of the stats package is a student’s t test, and is use when raw dataset is given. The same case for z.test, but this function is specifically for z-test of known population standard deviation. When dataset is not given and only the summary statistics (mean, and standard deviation) are presented, then the appropriate functions are zsum.test or tsum.test. Note that, t.test and tsum.test are functions of the same statistical test, and that of z.test and zsum.test. Consider the example below,

From the data, we obtain the following information: (i) the sample size is more than 30, and (ii) the population standard deviation is known. Therefore, the appropriate test is z-test. And the function to use is z.test, that is

Interpretation: The true mean of all SASE scores in the school year 2013-2014 from CSM is likely between 88.01327 and 92.65340 (95% CI). And the true mean of all SASE scores for the said college and school year is likely between 87.28425 and 93.38241 (99% CI).Aside from the confidence interval, the function returns also the computed z-statistics with p-value, and as well as the point estimate of the mean. To get rid of this, one can add a suffix $conf.int to the function to extract the confidence interval only.

The appropriate test for this is t-test since the sample size is small, n < 30, and the population variance is unknown. And thus,

Interpretation: The true mean of the population of the given data above is likely between 285.5911 and 356.1423 (90% CI).

Often in textbooks, however, we are presented with summary statistics of the data like the next example below from Simplified Biostatistics by Abubakar S. Asaad.

Example 3. The biostatistician took a random sample of 49 patients from a list of all patients ever admitted to the hospital within a three-month period and the number of drugs prescribed per admission was determined for each. The average drug per case was found to be 7.5 with standard deviation of 2.5. Calculate and interpret the 95% confidence interval for true mean of all the patients ever admitted to the hospital.

In this example, no dataset is given, but we have the computed mean = 7.5 of this dataset, standard deviation = 2.5, and sample size = 49. Thus, to compute for the interval estimate of the population mean in R, we use the zsum.test

Interpretation: The true mean of all the patients ever admitted to the hospital is likely between 6.800013 and 8.199987 (95% CI).

The tsum.test function is used in situation like in Example 3, but this time the population variance should be unknown and the sample size should be less than 30.