Assignment 2 - 2013

Due date(s):	30 January 2013
(PDF)	Assignment questions
(PDF)	Assignment solutions

<rst> <rst-options: 'toc' = False/> <rst-options: 'reset-figures' = False/>

.. rubric:: Assignment objectives: Univariate data analysis

.. question:: :grading: 5

Compute the mean, median, standard deviation and MAD for salt content for the various soy sauces given `in this report <http://beta.images.theglobeandmail.com/archive/00245/Read_the_report_245543a.pdf>`_ (page 41) as described in `this article from the Globe and Mail <http://www.theglobeandmail.com/life/health/salt-variation-between-brands-raises-call-for-cuts/article1299117/>`_ on 24 September 2009. Plot a box plot of the data and report the interquartile range (IQR). Comment on the 3 measures of spread you have calculated: standard deviation, MAD, and interquartile range.

.. [460, 520, 580, 700, 760, 770, 890, 910, 920, 940, 960, 1060, 1100]

.. question:: :grading: 3

Give a reason why Statistics Canada reports the median income when reporting income by geographic area. Where would you expect the mean to lie, relative to the median? Use `this table <http://www.statcan.gc.ca/tables-tableaux/sum-som/l01/cst01/famil107a-eng.htm>`_ to look up the income for Hamilton. How does it compare to Toronto? And all of Canada?

Can you locate the data that shows the income for each riding in Hamilton? How are these data reported? [A "riding" is a political division of a city or region, used for election purposes].

.. question:: :grading: 6

Use the data set on "`Raw material properties <http://datasets.connectmv.com/info/raw-material-properties>`_".

- How many variables in the data set? - How many observations? - The data are properties of a plastic pellets. Plot each variable, one at a time, and locate any outliers. See how to use the ``identify`` function in `part 9 of the tutorial <http://learnche.mcmaster.ca/4C3/Software_tutorial>`_. - Compare values of the mean and median for a column containing strong outliers. - Also compare the standard deviation against the MAD for the same column.

.. question:: :grading: 8

A new wastewater treatment plant is being commissioned and part of the commissioning report requires a statement of the confidence interval of the `biochemical oxygen demand (BOD) <http://en.wikipedia.org/wiki/Biochemical_oxygen_demand>`__. How many samples must you send to the lab to be sure the true BOD is within a range of 2 mg/L, centered about the sample average? If there isn't enough information given here, specify your own numbers and assumptions and work with them to answer the question.

.. question:: :grading: 8

One of the questions we posed at the start of this chapter was: "`Here are the yields <http://datasets.connectmv.com/info/batch-yields>`_ from a batch bioreactor system for the last 3 years (300 data points; we run a new batch about every 3 to 4 days).

#. What sort of distribution do the yield data have? #. A recorded yield value today was less than 60%, what are the chances of that occurring? Express your answer as: *there's a 1 in x chance* of it occurring. #. Which assumptions do you have to make for the second part of this question?

.. From assignment 2, 2011

.. question:: :grading: 8 (required for 600-level students; extra credit for 400-level students)

#. At the 95% confidence level, for a sample size of 7, compare and comment on the upper and lower bounds of the confidence interval that you would calculate if:

a) you know the population standard deviation b) you have to estimate it for the sample.

Assume that the calculated standard deviation from the sample, :math:`s` matches the population :math:`\sigma = 4.19`.

#. As a follow up, overlay the probability distribution curves for the normal and :math:`t`-distribution that you would use for a sample of data of size :math:`n=7`.

#. Repeat part of this question, using larger sample sizes. At which point does the difference between the :math:`t`- and normal distributions become *practically* indistinguishable?

#. What is the implication of this?

.. question:: :grading: 0

For additional R practice, and data interpretation practice, review question 4 from `last year's Assignment 2 <http://stats4eng.connectmv.com/wiki/Assignment_2_-_2012>`_ (attempt it without looking at the solutions).

.. question:: :grading: 600-level: 0 points

*The solution appears as question 27 in* `PID <http://learnche.mcmaster.ca/pid/>`_

The paper by PJ Rousseeuw, "`Tutorial to Robust Statistics <http://dx.doi.org/10.1002/cem.1180050103>`_", *Journal of Chemometrics*, **5**, 1-20, 1991 discusses the breakdown point of a statistic.

#. Describe what the breakdown point is, and give two examples: one with a low breakdown point, and one with a high breakdown point. Use a vector of numbers to help illustrate your answer.

#. What is an advantage of using robust methods over their "classical" counterparts?

</rst>