Smoking habits of UK residents. A survey was conducted to study the smoking habits of UK residents. Below is a data matrix displaying a portion of the data collected in this survey. Note that “£” stands for British Pounds Sterling, “cig” stands for cigarettes, and “N/A” refers to a missing component of the data.
| ID | sex | age | marital | grossIncome | smoke | amtWeekends | amtWeekdays |
|---|---|---|---|---|---|---|---|
| 1 | Female | 42 | Single | Under£2,600 | Yes | 12 cig/day | 12 cig/day |
| 2 | Male | 44 | Single | £10,400 to £15,600 | No | N/A | N/A |
| 3 | Male | 53 | Married | Above£36,400 | Yes | 6 cig/day | 6 cig/day |
| … | … | … | … | … | … | … | … |
| 1691 | Male | 40 | Single | £2,600 to £5,200 | Yes | 8 cig/day | 8 cig/day |
| sex | age | marital | grossIncome | smoke | amtWeekends | amtWeekdays |
|---|---|---|---|---|---|---|
| categorical | continuous | categorical | ordinal | categorical | discrete | discrete |
Cheaters, scope of inference. Exercise 1.5 introduces a study where researchers studying the relationship between honesty, age, and self-control conducted an experiment on 160 children between the ages of 5 and 15. The researchers asked each child to toss a fair coin in private and to record the outcome (white or black) on a paper sheet, and said they would only reward children who report white. Half the students were explicitly told not to cheat and the others were not given any explicit instructions. Differences were observed in the cheating rates in the instruction and no instruction groups, as well as some differences across children’s characteristics within each group.
Reading the paper. Below are excerpts from two articles published in the NY Times:
“Researchers analyzed data from 23,123 health plan members who participated in a voluntary exam and health behavior survey from 1978 to 1985, when they were 50-60 years old. 23 years later, about 25% of the group had dementia, including 1,136 with Alzheimer’s disease and 416 with vascular dementia. After adjusting for other factors, the researchers concluded that pack-aday smokers were 37% more likely than nonsmokers to develop dementia, and the risks went up with increased smoking; 44% for one to two packs a day; and twice the risk for more than two packs.”
Based on this study, can we conclude that smoking causes dementia later in life? Explain your reasoning.
“The University of Michigan study, collected survey data from parents on each child’s sleep habits and asked both parents and teachers to assess behavioral concerns. About a third of the students studied were identified by parents or teachers as having problems with disruptive behavior or bullying. The researchers found that children who had behavioral issues and those who were identified as bullies were twice as likely to have shown symptoms of sleep disorders.”
A friend of yours who read the article says, “The study shows that sleep disorders lead to bullying in school children.” Is this statement justified? If not, how best can you describe the conclusion that can be drawn from this study?
Exercise and mental health. A researcher is interested in the effects of exercise on mental health and he proposes the following study: Use stratified random sampling to ensure representative proportions of 18-30, 31-40 and 41- 55 year old from the population. Next, randomly assign half the subjects from each age group to exercise twice a week, and instruct the rest not to exercise. Conduct a mental health exam at the beginning and at the end of the study, and compare the results.
Stats scores. Below are the final exam scores of twenty introductory statistics students.
57, 66, 69, 71, 72, 73, 74, 77, 78, 78, 79, 79, 81, 81, 82, 83, 83, 88, 89, 94
Create a box plot of the distribution of these scores. The five number summary provided below may be useful.
| Min | Q1 | Q2 | (Median) | Q3 | Max |
|---|---|---|---|---|---|
| 57 | 72.5 | 78.5 | 82.5 | 94 |
scores <- c(57, 66, 69, 71, 72, 73, 74, 77, 78, 78, 79, 79, 81, 81, 82, 83, 83, 88, 89, 94)
boxplot(scores, ylab="Final Scores", xlab="Scores")
title("Boxplot of Final Exam Scores")
Mix-and-match. Describe the distribution in the histograms below and match them to the box plots.
Histogram A has a distribution that appears to be normal, it is unimodal and symmetric, therefore it matches with boxplot 2 because the mean, median, and mode all see to occur at the same point.
Histogram B looks more like multimodel than uniform since it starts high, dips, high again before it dips once more and forms another smaller mode and it is symmetric, therefore it matches with boxplot 3 because the interquartile range is very wide.
Histogram C is unimodal and skewed to the right, therefore it matches with boxplot 1 because upper whisker is further than where most of the data is, with many outliers.
Distributions and appropriate statistics, Part II . For each of the following, state whether you expect the distribution to be symmetric, right skewed, or left skewed. Also specify whether the mean or median would best represent a typical observation in the data, and whether the variability of observations would be best represented using the standard deviation or IQR. Explain your reasoning.
This distribution is expected to be extremely left skewed since many houses are priced over $6,000,000. The median would be a better represention of a typical obsevation in this data. The IQR would depict the variability of observations because it is not affected by the outliers in the 6 million price range.
This distribution is expected to be symmetric since there is not many houses that are overpriced, and the quartiles fall evenly. The mean would be a better represention of a typical obsevation in this data. The standard deviation would depict the variability of observations because there is not many extreme prices, either too high or too low.
This distribution is expected to be extremely right skewed since many students don’t drink. The median would be a better represention of a typical obsevation in this data. The IQR would depict the variability of observations because it is not affected by those who do drink.
This distribution is expected to be right skewed since there is a few executives with a high salary. The median would be a better represention of a typical obsevation in this data. The IQR would depict the variability of observations because it is not affected by the high salary values.
Heart transplants. The Stanford University Heart Transplant Study was conducted to determine whether an experimental heart transplant program increased lifespan. Each patient entering the program was designated an official heart transplant candidate, meaning that he was gravely ill and would most likely benefit from a new heart. Some patients got a transplant and some did not. The variable transplant indicates which group the patients were in; patients in the treatment group got a transplant and those in the control group did not. Another variable called survived was used to indicate whether or not the patient was alive at the end of the study.
We write alive on [xxx] cards representing patients who were alive at the end of the study, and dead on [xxx] cards representing patients who were not. Then, we shuffe these cards and split them into two groups: one group of size [xxx] representing treatment, and another group of size [xxx] representing control. We calculate the difference between the proportion of dead cards in the treatment and control groups (treatment - control) and record this value. We repeat this 100 times to build a distribution centered at [xxx]. Lastly, we calculate the fraction of simulations where the simulated differences in proportions are [xxx]. If this fraction is low, we conclude that it is unlikely to have observed such an outcome by chance and that the null hypothesis should be rejected in favor of the alternative.
iii. What do the simulation results shown below suggest about the effectiveness of the transplant program?
library(openintro)
## Please visit openintro.org for free statistics materials
##
## Attaching package: 'openintro'
## The following objects are masked from 'package:datasets':
##
## cars, trees
data(heartTr)
table1<-table(heartTr$survived,heartTr$transplant)
table1[2,]/colSums(table1)
## control treatment
## 0.8823529 0.6521739
The proportion of patients in the control group that dies is 88.2%. The proportion of patients in the treatment group that dies is 65.2%