02-data-design.qmd

# Study design {#sec-data-design}

```{r}
#| include: false
source("_common.R")
```

::: {.chapterintro data-latex=""}
Before digging into the details of working with data, we stop to think about how data come to be.
That is, if the data are to be used to make broad and complete conclusions, then it is important to understand who or what the data represent.
One important aspect of data provenance is sampling.
Knowing how the observational units were selected from a larger entity will allow for generalizations back to the population from which the data were randomly selected.
Additionally, by understanding the structure of the study, causal relationships can be separated from those relationships which are only associated.
A good question to ask oneself before working with the data at all is, "How were these observations collected?".
You will learn a lot about the data by understanding its source.
:::

## Sampling principles and strategies {#sec-sampling-principles-strategies}

The first step in conducting research is to identify topics or questions that are to be investigated.
A clearly laid out research question is helpful in identifying what subjects or cases should be studied and what variables are important.
It is also important to consider *how* data are collected so that the data are reliable and help achieve the research goals.

### Populations and samples

Consider the following three research questions:

1.  What is the average mercury content in swordfish in the Atlantic Ocean?
2.  Over the last five years, what is the average time to complete a degree for Duke undergrads?
3.  Does a new drug reduce the number of deaths in patients with severe heart disease?

Each research question refers to a target **population**.
In the first question, the target population is all swordfish in the Atlantic Ocean, and each fish represents a case.
Oftentimes, it is not feasible to collect data for every case in a population.
Collecting data for an entire population is called a **census**.
A census is difficult because it is too expensive to collect data for the entire population, but it might also be because it is difficult or impossible to identify the entire population of interest!
Instead, a sample is taken.
A **sample** is the data you have.
Ideally, a sample is a small fraction of the population.
For instance, 60 swordfish (or some other number) in the population might be selected, and this sample data may be used to provide an estimate of the population average and to answer the research question.

```{r}
#| include: false
terms_chp_2 <- c("population", "census", "sample")
```

::: {.guidedpractice data-latex=""}
For the second and third questions above, identify the target population and what represents an individual case.[^02-data-design-1]
:::

[^02-data-design-1]: The question *"Over the last five years, what is the average time to complete a degree for Duke undergrads?"* is only relevant to students who complete their degree; the average cannot be computed using a student who never finished their degree.
    Thus, only Duke undergrads who graduated in the last five years represent cases in the population under consideration.
    Each such student is an individual case.
    For the question *"Does a new drug reduce the number of deaths in patients with severe heart disease?"*, a person with severe heart disease represents a case.
    The population includes all people with severe heart disease.

### Parameters and statistics

In most statistical analysis procedures, the research question at hand boils down to understanding a numerical summary.
The number (or set of numbers) may be a quantity you are already familiar with (like the average) or it may be something you learn through this text (like the slope and intercept from a least squares model, provided in @sec-least-squares-regression).

A numerical summary can be calculated on either the sample of observations or the entire population.
However, measuring every unit in the population is usually prohibitive (so the parameter is very rarely calculated).
So, a "typical" numerical summary is calculated from a sample.
Yet, we can still conceptualize calculating the average income of all adults in Argentina.

We use specific terms in order to differentiate when a number is being calculated on a sample of data (**statistic**) and when it is being calculated or considered for calculation on the entire population (**parameter**).
The terms statistic and parameter are useful for communicating claims and models and will be used extensively in later chapters which delve into making inference on populations.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "parameter", "statistic")
```

### Anecdotal evidence

\index{bias}

Consider the following possible responses to the three research questions:

1.  A man on the news got mercury poisoning from eating swordfish, so the average mercury concentration in swordfish must be dangerously high.
2.  I met two students who took more than 7 years to graduate from Duke, so it must take longer to graduate at Duke than at many other colleges.
3.  My friend's dad had a heart attack and died after they gave him a new heart disease drug, so the drug must not work.

Each conclusion is based on data.
However, there are two problems.
First, the data only represent one or two cases.
Second, and more importantly, it is unclear whether these cases are actually representative of the population.
Data collected in this haphazard fashion are called **anecdotal evidence**.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "anecdotal evidence")
```

::: {.important data-latex=""}
**Anecdotal evidence.**

Be careful of data collected in a haphazard fashion.
Such evidence may be true and verifiable, but it may only represent extraordinary cases and therefore not be a good representation of the population.
:::

```{r}
#| label: mn-winter
#| out.width: 35%
#| fig-cap: In February 2010, some media pundits cited one large snowstorm as evidence
#|   against global warming. As comedian Jon Stewart pointed out, "It is one storm, in
#|   one region of one country."
#| fig-alt: Photograph shows a city street covered in snow and ice.
include_graphics("images/mn-winter/mn-winter.jpg")
```

Anecdotal evidence typically is composed of unusual cases that we recall based on their striking characteristics.
For instance, we are more likely to remember the two people we met who took 7 years to graduate than the six others who graduated in four years.
Instead, of looking at the most unusual cases, we should examine a sample of many cases that better represent the population.

### Sampling from a population

\index{random sample} \index{bias}

We might try to estimate the time to graduation for Duke undergraduates in the last five years by collecting a sample of graduates.
All graduates in the last five years represent the *population*\index{population}, and graduates who are selected for review are collectively called the *sample*\index{sample}.
In general, we always seek to *randomly* select a sample from a population.
The most basic type of random selection is equivalent to how raffles are conducted.
For example, in selecting graduates, we could write each graduate's name on a raffle ticket and draw 10 tickets.
The selected names would represent a random sample of 10 graduates.

```{r}
#| label: pop-to-sample
#| fig-cap: In this graphic, 10 graduates are randomly selected from the population to
#|   be included in the sample.
#| fig-asp: 0.5
#| fig-width: 8.0
#| fig-alt: A large circle contains many dots which indicate all the graduates. A smaller
#|   circle contains a few of the dots (i.e., graduates) which have been randomly selected
#|   from the larger circle.
#| out.width: 70%
set.seed(1234)
par_og <- par(no.readonly = TRUE) # save original par
par(mar = rep(0, 4))
plot(c(0, 2), c(0, 1.1), type = "n", axes = FALSE, xlab = "", ylab = "")
temp <- seq(0, 2 * pi, 2 * pi / 100)
x <- 0.5 + 0.5 * cos(temp)
y <- 0.5 + 0.5 * sin(temp)
lines(x, y)

s <- matrix(runif(1000), ncol = 2)
S <- matrix(NA, 350, 2)
j <- 0
for (i in 1:nrow(s)) {
  if (sum((s[i, ] - 0.5)^2) < 0.23) {
    j <- j + 1
    S[j, ] <- s[i, ]
  }
}
points(S, col = IMSCOL["blue", "f2"], pch = 20)
text(0.5, 1, "all graduates", pos = 3, cex = 1.3)

set.seed(50)
N <- sample(j, 25)
lines((x - 0.5) / 2 + 1.5, (y - 0.5) / 2 + 0.5, pch = 20)

SS <- (S[N, ] - 0.5) / 2 + 0.5
these <- c(2, 5, 10, 12, 20, 21, 22, 23, 1, 8)
points(SS[these, 1] + 1, SS[these, 2], col = IMSCOL["red", "f1"], pch = 20, cex = 1.5)
text(1.5, 0.75, "sample", pos = 3, cex = 1.3)

for (i in these) {
  arrows(S[N[i], 1], S[N[i], 2],
    SS[i, 1] + 1 - 0.03, SS[i, 2],
    length = 0.08, col = IMSCOL["black", "full"], lwd = 1.5
  )
}
par(par_og) # restore original par
```

::: {.workedexample data-latex=""}
Suppose we ask a student who happens to be majoring in nutrition to select several graduates for the study.
What kind of students do you think they might collect?
Do you think their sample would be representative of all graduates?

------------------------------------------------------------------------

Perhaps they would pick a disproportionate number of graduates from health-related fields.
Or perhaps their selection would be a good representation of the population.
When selecting samples by hand, we run the risk of picking a **biased** sample, even if our bias is unintended.
:::

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "sample bias")
```

```{r}
#| label: pop-to-sub-sample-graduates
#| fig-cap: Asked to pick a sample of graduates, a nutrition major might inadvertently
#|   pick a disproportionate number of graduates from health-related majors.
#| fig-asp: 0.5
#| fig-width: 8.0
#| fig-alt: A large circle contains many dots which indicate all the graduates, but some
#|   of the dots have been greyed out where others are dark dots from which the sample
#|   is taken.  A smaller circle contains a few of the dots (i.e., graduates) which have
#|   been selected from the biased group of dark dots in the large circle.
#| out.width: 70%
par_og <- par(no.readonly = TRUE) # save original par
par(mar = rep(0, 4))
plot(c(0, 2), c(0, 1.1), type = "n", axes = FALSE, xlab = "", ylab = "")
temp <- seq(0, 2 * pi, 2 * pi / 100)
x <- 0.5 + 0.5 * cos(temp)
y <- 0.5 + 0.5 * sin(temp)
lines(x, y)

s <- matrix(runif(1000), ncol = 2)
S <- matrix(NA, 350, 2)
j <- 0
sub <- rep(FALSE, 1000)
for (i in 1:nrow(s)) {
  if (sum((s[i, ] - 0.5)^2) < 0.23) {
    j <- j + 1
    S[j, ] <- s[i, ]
  }
  if (sum((s[i, ] - c(0.05, 0.18) - 0.5)^2) < 0.07) {
    sub[j] <- TRUE
  }
}
points(S, col = IMSCOL["blue", 4 - 2 * sub], pch = 20)
text(0.5, 1, "all graduates", pos = 3, cex = 1.3)
lines(
  (x - 0.5) * 2 * sqrt(0.07) + 0.55,
  (y - 0.5) * 2 * sqrt(0.07) + 0.68
)

set.seed(7)
N <- sample((1:j)[sub], 25)
lines((x - 0.5) / 2 + 1.5,
  (y - 0.5) / 2 + 0.5,
  pch = 20
)

SS <- (S[N, ] - 0.5) / 2 + 0.5
these <- c(2, 5, 10, 12, 20, 21, 22, 23, 1, 8)
points(SS[these, 1] + 1, SS[these, 2], col = IMSCOL["red", "f1"], pch = 20, cex = 1.5)
text(1.5, 0.75, "sample", pos = 3, cex = 1.3)

for (i in these) {
  arrows(S[N[i], 1], S[N[i], 2],
    SS[i, 1] + 1 - 0.03, SS[i, 2],
    length = 0.08,
    col = IMSCOL["black", "full"],
    lwd = 1.5
  )
}
rect(0.143, 0.2, 0.952, 0.301,
  border = "#00000000",
  col = "#FFFFFF88"
)
rect(0.236, 0.301, 0.858, 0.403,
  border = "#00000000",
  col = "#FFFFFF88"
)
text(0.55, 0.5 + 0.18 - sqrt(0.07),
  "graduates from\nhealth-related fields",
  pos = 1, cex = 1.3
)
par(par_og) # restore original par
```

If someone was permitted to pick and choose exactly which graduates were included in the sample, it is entirely possible that the sample would overrepresent that person's interests, which may be entirely unintentional.
This introduces **bias** into a sample.
Sampling randomly helps address this problem.
The most basic random sample is called a **simple random sample** and is equivalent to drawing names out of a hat to select cases.
This means that each case in the population has an equal chance of being included and the cases in the sample are not related to each other.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "bias", "simple random sample")
```

The act of taking a simple random sample helps minimize bias.
However, bias can crop up in other ways.
Even when people are picked at random, e.g., for surveys, caution must be exercised if the **non-response rate**\index{non-response rate} is high.
For instance, if only 30% of the people randomly sampled for a survey actually respond, then it is unclear whether the results are **representative** of the entire population.
This **non-response bias**\index{non-response bias} can skew results.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "non-response rate", "representative", "non-response bias")
```

```{r}
#| label: survey-sample
#| fig-cap: Due to the possibility of non-response, survey studies may only reach a certain
#|   group within the population. It is difficult, and oftentimes impossible, to completely
#|   fix this problem.
#| fig-asp: 0.5
#| fig-width: 8.0
#| fig-alt: A large circle contains many dots which indicate the population of interest,
#|   but some of the dots have been greyed out where others are dark dots from which
#|   the sample is taken (where the grey dots are potentially due to non-response bias).  A
#|   smaller circle contains a few of the dots which have been selected from the group
#|   of dark dots in the large circle who were individuals willing to respond to the
#|   survey.
#| out.width: 70%
par_og <- par(no.readonly = TRUE) # save original par
par(mar = rep(0, 4))
plot(c(0, 2),
     c(0, 1.1),
     type='n',
     axes=FALSE,
     xlab="",
     ylab="")
temp <- seq(0, 2 * pi, 2 * pi / 100)
x <- 0.5 + 0.5 * cos(temp)
y <- 0.5 + 0.5 * sin(temp)
lines(x, y)

s <- matrix(runif(700), ncol = 2)
S <- matrix(NA, 350, 2)
j <- 0
sub <- rep(FALSE, 1000)
for (i in 1:nrow(s)) {
  if (sum((s[i,] - 0.5)^2) < 0.23) {
    j <- j + 1
    S[j, ] <- s[i, ]
  }
  if (sum((s[i, ] - c(-0.15, 0.05) - 0.5)^2) < 0.115) {
    sub[j] <- TRUE
  }
}
points(S, col = IMSCOL["blue", 4 - 2 * sub], pch = 20)
text(0.5, 1, 'population of interest', pos = 3, cex = 1.3)
lines((x - 0.5) * 2 * sqrt(0.115) + 0.35,
      (y - 0.5) * 2 * sqrt(0.115) + 0.55)

set.seed(7)
N <- sample((1:j)[sub], 25)
lines((x - 0.5) / 2 + 1.5,
      (y - 0.5) / 2 + 0.5,
      pch=20)

SS <- (S[N, ] - 0.5) / 2 + 0.5
these <- c(2, 5, 6, 7, 15)
points(SS[these, 1] + 1,
       SS[these, 2],
       col = IMSCOL["red", "f1"],
       pch = 20,
       cex = 1.5)
text(1.5, 0.75, 'sample', pos=3, cex = 1.3)

for(i in these){
	arrows(S[N[i], 1],
	       S[N[i], 2],
	       SS[i, 1] + 1 - 0.03,
	       SS[i, 2],
	       length=0.08,
	       col=IMSCOL["black", "full"],
	       lwd=1.5)
}
rect(0.145, 0.195, 0.775, 0.11,
     border="#00000000",
     col="#FFFFFF88")
rect(0.31, 0.018, 0.605, 0.11,
     border="#00000000",
     col="#FFFFFF88")
text(0.46, 0.5 + 0.06 - sqrt(0.115),
     'population actually\nsampled',
     pos=1,
     cex=1)
par(par_og) # restore original par
```

Another common downfall is a **convenience sample**\index{convenience sample}, where individuals who are easily accessible are more likely to be included in the sample.
For instance, if a political survey is done by stopping people walking in the Bronx, this will not represent all of New York City.
It is often difficult to discern what sub-population a convenience sample represents.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "convenience sample")
```

::: {.guidedpractice data-latex=""}
We can easily access ratings for products, sellers, and companies through websites.
These ratings are based only on those people who go out of their way to provide a rating.
If 50% of online reviews for a product are negative, do you think this means that 50% of buyers are dissatisfied with the product?
Why or why not?[^02-data-design-2]
:::

[^02-data-design-2]: Answers will vary.
    From our own anecdotal experiences, we believe people tend to rant more about products that fell below expectations than rave about those that perform as expected.
    For this reason, we suspect there is a negative bias in product ratings on sites like Amazon.
    However, since our experiences may not be representative, we also keep an open mind.

\index{random sample} \index{bias} \index{population} \index{sample}

### Four sampling methods {#sec-samp-methods}

Almost all statistical methods are based on the notion of implied randomness.
If data are not collected in a random framework from a population, these statistical methods -- the estimates and errors associated with the estimates -- are not reliable.
Here we consider four random sampling techniques: simple, stratified, cluster, and multistage sampling.
Figures @fig-simple-stratified and @fig-cluster-multistage provide graphical representations of these techniques.

\index{simple random sampling} \index{stratified sampling}

```{r}
#| label: fig-simple-stratified
#| fig-cap: 'Examples of simple random and stratified sampling. In the top panel, simple
#|   random sampling was used to randomly select the 18 cases (denoted in red). In the
#|   bottom panel, stratified sampling was used: cases were first grouped into strata,
#|   then simple random sampling was employed to randomly select 3 cases within each
#|   stratum.'
#| fig-width: 10.0
#| out.width: 100%
#| fig-alt: The top box shows a population of dots (i.e., individuals) where a handful
#|   of the dots have been sampled randomly.  The bottom box shows the same population
#|   of dots but grouped in such a way that there are six strata.  From each stratum
#|   three dots (i.e., individuals) are randomly selected.

source("helpers/helper-sampling.R")   # source helper
par_og <- par(no.readonly = TRUE)     # save original par
par(mar = rep(0.5,4), mfrow = c(2,1)) # no margin, 2 figures 
build_srs(n = 18, N = 108)            # build figure
build_stratified(N = 108)             # build figure
par(par_og)                           # restore original par
```

**Simple random sampling** is probably the most intuitive form of random sampling.
Consider the salaries of Major League Baseball (MLB) players, where each player is a member of one of the league's 30 teams.
To take a simple random sample of 120 baseball players and their salaries, we could write the names of that season's several hundreds of players onto slips of paper, drop the slips into a bucket, shake the bucket around until we are sure the names are all mixed up, then draw out slips until we have the sample of 120 players.
In general, a sample is referred to as "simple random" if each case in the population has an equal chance of being included in the final sample *and* knowing that a case is included in a sample does not provide useful information about which other cases are included.

\index{strata}

**Stratified sampling** is a divide-and-conquer sampling strategy.
The population is divided into groups called **strata**.
The strata are chosen so that similar cases are grouped together, then a second sampling method, usually simple random sampling, is employed within each stratum.
In the baseball salary example, each of the 30 teams could represent a stratum, since some teams have a lot more money (up to 4 times as much!).
Then we might randomly sample 4 players from each team for our sample of 120 players.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "simple random sampling", "stratified sampling", "strata")
```

**Stratified sampling** is especially useful when the cases in each stratum are very similar with respect to the outcome of interest.
The downside is that analyzing data from a stratified sample is a more complex task than analyzing data from a simple random sample.
The analysis methods introduced in this book would need to be extended to analyze data collected using stratified sampling.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "stratified sampling")
```

::: {.workedexample data-latex=""}
Why would it be good for cases within each stratum to be very similar?

------------------------------------------------------------------------

We might get a more stable estimate for the subpopulation in a stratum if the cases are very similar, leading to more precise estimates within each group.
When we combine these estimates into a single estimate for the full population, that population estimate will tend to be more precise since each individual group estimate is itself more precise.
:::

In a **cluster sample**, we break up the population into many groups, called **clusters**.
Then we sample a fixed number of clusters and include all observations from each of those clusters in the sample.
A **multistage sample** is like a cluster sample, but rather than keeping all observations in each cluster, we would collect a random sample within each selected cluster.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "cluster sampling", "cluster", "multistage sample")
```

\index{cluster sampling}

```{r}
#| label: fig-cluster-multistage
#| fig-cap: 'Examples of cluster and multistage sampling. In the top panel, cluster sampling
#|   was used: data were binned into nine clusters, three of these clusters were sampled,
#|   and all observations within these three clusters were included in the sample. In
#|   the bottom panel, multistage sampling was used, which differs from cluster sampling
#|   only in that we randomly select a subset of each cluster to be included in the sample
#|   rather than measuring every case in each sampled cluster.'
#| fig-width: 10.0
#| out.width: 100%
#| fig-alt: In the top figure, dots are grouped into clusters, three clusters are selected,
#|   and every dot (i.e., all individuals) from each of the three clusters are sampled.  In
#|   the bottom figure, dots are again grouped into clusters and three clusters are selected.  However,
#|   random sampling is applied so that a random sample from each of the three selected
#|   clusters is taken.

source("helpers/helper-sampling.R")   # source helper
par_og <- par(no.readonly = TRUE)     # save original par
par(mar = rep(0.5,4), mfrow = c(2,1)) # no margin, 2 figures 
build_cluster()                       # build figure
build_multistage()                    # build figure
par(par_og)                           # restore original par
```

Sometimes cluster or multistage sampling can be more economical than the alternative sampling techniques.
Also, unlike stratified sampling, these approaches are most helpful when there is a lot of case-to-case variability within a cluster but the clusters themselves do not look very different from one another.
For example, if neighborhoods represented clusters, then cluster or multistage sampling work best when the populations inside each neighborhood are very diverse.
A downside of these methods is that more advanced techniques are typically required to analyze the data, though the methods in this book can be extended to handle such data.

::: {.workedexample data-latex=""}
Suppose we are interested in estimating the malaria rate in a densely tropical portion of rural Indonesia.
We learn that there are 30 villages in that part of the Indonesian jungle, each more or less like the next, but the distances between the villages are substantial.
Our goal is to test 150 individuals for malaria.
What sampling method should be employed?

------------------------------------------------------------------------

A simple random sample would likely draw individuals from all 30 villages, which could make data collection extremely expensive.
Stratified sampling would be a challenge since it is unclear how we would build strata of similar individuals.
However, cluster sampling or multistage sampling seem like very good ideas.
If we decided to use multistage sampling, we might randomly select half of the villages, then randomly select 10 people from each.
This would probably reduce our data collection costs substantially in comparison to a simple random sample, and the cluster sample would still give us reliable information, even if we would need to analyze the data with slightly more advanced methods than we discuss in this book.
:::

\clearpage

## Experiments {#sec-experiments}

Studies where the researchers assign treatments to cases are called **experiments**.
When this assignment includes randomization, e.g., using a coin flip to decide which treatment a patient receives, it is called a **randomized experiment**.
Randomized experiments are fundamentally important when trying to show a causal connection between two variables.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "experiment", "randomized experiment")
```

### Principles of experimental design {#sec-principles-experimental-design}

1.  **Controlling.** Researchers assign treatments to cases, and they do their best to **control** any other differences in the groups[^02-data-design-3]. For example, when patients take a drug in pill form, some patients take the pill with only a sip of water while others may have it with an entire glass of water. To control for the effect of water consumption, a doctor may instruct every patient to drink a 12-ounce glass of water with the pill.

[^02-data-design-3]: This is a different concept than a *control group*, which we discuss in the second principle and in @sec-reducing-bias-human-experiments.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "control")
```

2.  **Randomization.** Researchers randomize patients into treatment groups to account for variables that cannot be controlled.
    For example, some patients may be more susceptible to a disease than others due to their dietary habits.
    In this example dietary habit is a **confounding variable**[^02-data-design-4], which is defined as a variable that is associated with both the explanatory and response variables.
    Randomizing patients into the treatment or control group helps even out such differences.

    ```{r}
    #| include: false
    terms_chp_2 <- c(terms_chp_2, "confounding variable")
    ```

3.  **Replication.** The more cases researchers observe, the more accurately they can estimate the effect of the explanatory variable on the response.
    In a single study, we **replicate** by collecting a sufficiently large sample.
    What is considered sufficiently large varies from experiment to experiment, but at a minimum we want to have multiple subjects (experimental units) per treatment group.
    Another way of achieving replication is replicating an entire study to verify an earlier finding.
    The term **replication crisis** refers to the ongoing methodological crisis in which past findings from scientific studies in several disciplines have failed to be replicated.
    **Pseudoreplication** occurs when individual observations under different treatments are heavily dependent on each other.
    For example, suppose you have 50 subjects in an experiment where you're taking blood pressure measurements at 10 time points throughout the course of the study.
    By the end, you will have 50 $\times$ 10 = 500 measurements.
    Reporting that you have 500 observations would be considered pseudoreplication, as the blood pressure measurements of a given individual are not independent of each other.
    Pseudoreplication often happens when the wrong entity is replicated, and the reported sample sizes are exaggerated.

[^02-data-design-4]: Also called a **lurking variable**, **confounding factor**, or a **confounder**.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "replication", "pseudoreplication", "replication crisis")
```

4.  **Blocking.** Researchers sometimes know or suspect that variables, other than the treatment, influence the response. Under these circumstances, they may first group individuals based on this variable into **blocks** and then randomize cases within each block to the treatment groups. This strategy is often referred to as **blocking**. For instance, if we are looking at the effect of a drug on heart attacks, we might first split patients in the study into low-risk and high-risk blocks, then randomly assign half the patients from each block to the control group and the other half to the treatment group, as shown in @fig-blocking. This strategy ensures that each treatment group has the same number of low-risk patients and the same number of high-risk patients.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "blocking")
```

```{r}
#| label: fig-blocking
#| fig-cap: Blocking using a variable depicting patient risk. Patients are first divided
#|   into low-risk and high-risk blocks, then each block is evenly separated into the
#|   treatment groups using randomization. This strategy ensures an equal representation
#|   of patients in each treatment group from both the low-risk and high-risk categories.
#| warning: false
#| fig-width: 6.0
#| fig-asp: 1.0
#| fig-alt: Before randomly allocating, the red low risk patients and blue high risk
#|   patients are split into two separate groups. Subsequently, half of the red low risk
#|   patients are randomly chosen to receive the treatment, and half of the blue high
#|   risk patients are randomly chosen to receive the treatment.
#| out.width: 100%

set.seed(12345)
source("helpers/helper-blocking.R") # source helper
par_og <- par(no.readonly = TRUE)   # save original par
par(mar = rep(0,4))                 # no margins
build_blocking()                    # build figure
par(par_og)                         # restore original par
```

It is important to incorporate the first three experimental design principles into any study, and this book describes applicable methods for analyzing data from such experiments.
Blocking is a slightly more advanced technique, and statistical methods in this book may be extended to analyze data collected using blocking.

\clearpage

### Reducing bias in human experiments {#sec-reducing-bias-human-experiments}

Randomized experiments have long been considered to be the gold standard for data collection, but they do not ensure an unbiased perspective into the cause-and-effect relationship in all cases.
Human studies are perfect examples where bias can unintentionally arise.
Here we reconsider a study where a new drug was used to treat heart attack patients.
In particular, researchers wanted to know if the drug reduced deaths in patients.

These researchers designed a randomized experiment because they wanted to draw causal conclusions about the drug's effect.
Study volunteers[^02-data-design-5] were randomly placed into two study groups.
One group, the **treatment group**, received the drug.
The other group, called the **control group**, did not receive any drug treatment.

[^02-data-design-5]: Human subjects are often called **patients**, **volunteers**, or **study participants**.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "treatment group", "control group")
```

Put yourself in the place of a person in the study.
If you are in the treatment group, you are given a fancy new drug that you anticipate will help you.
On the other hand, a person in the other group does not receive the drug and sits idly, hoping her participation does not increase her risk of death.
These perspectives suggest there are actually two effects in this study: the one of interest is the effectiveness of the drug, and the second is an emotional effect of (not) taking the drug, which is difficult to quantify.

Researchers aren't usually interested in the emotional effect, which might bias the study.
To circumvent this problem, researchers do not want patients to know which group they are in.
When researchers keep the patients uninformed about their treatment, the study is said to be **blind**.
But there is one problem: if a patient does not receive a treatment, they will know they're in the control group.
A solution to this problem is to give a fake treatment to patients in the control group.
This is called a **placebo**, and an effective placebo is the key to making a study truly blind.
A classic example of a placebo is a sugar pill that is made to look like the actual treatment pill.
However, offering such a fake treatment may not be ethical in certain experiments.
For example, in medical experiments, typically the control group must get the current standard of care.
Oftentimes, a placebo results in a slight but real improvement in patients.
This effect has been dubbed the **placebo effect**.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "blind", "placebo", "placebo effect")
```

The patients are not the only ones who should be blinded: doctors and researchers can unintentionally bias a study.
When a doctor knows a patient has been given the real treatment, they might inadvertently give that patient more attention or care than a patient that they know is on the placebo.
To guard against this bias, which again has been found to have a measurable effect in some instances, most modern studies employ a **double-blind** setup where doctors or researchers who interact with patients are, just like the patients, unaware of who is or is not receiving the treatment.[^02-data-design-6]

[^02-data-design-6]: There are always some researchers involved in the study who do know which patients are receiving which treatment.
    However, they do not interact with the study's patients and do not tell the blinded health care professionals who is receiving which treatment.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "double-blind")
```

::: {.guidedpractice data-latex=""}
Look back to the study in @sec-case-study-stents-strokes where researchers were testing whether stents were effective at reducing strokes in at-risk patients.
Is this an experiment?
Was the study blinded?
Was it double-blinded?[^02-data-design-7]
:::

[^02-data-design-7]: The researchers assigned the patients into their treatment groups, so this study was an experiment.
    However, the patients could distinguish what treatment they received because a stent is a surgical procedure.
    There is no equivalent surgical placebo, so this study was not blind.
    The study could not be double-blind since it was not blind.

::: {.guidedpractice data-latex=""}
For the study in @sec-case-study-stents-strokes, could the researchers have employed a placebo?
If so, what would that placebo have looked like?[^02-data-design-8]
:::

[^02-data-design-8]: Ultimately, can we make patients think they got treated from a surgery?
    In fact, we can, and some experiments use a **sham surgery**.
    In a sham surgery, the patient does undergo surgery, but the patient does not receive the full treatment, though they will still get a placebo effect.

You may have many questions about the ethics of sham surgeries to create a placebo.
These questions may have even arisen in your mind when in the general experiment context, where a possibly helpful treatment was withheld from individuals in the control group; the main difference is that a sham surgery tends to create additional risk, while withholding a treatment only maintains a person's risk.

There are always multiple viewpoints of experiments and placebos, and rarely is it obvious which is ethically "correct".
For instance, is it ethical to use a sham surgery when it creates a risk to the patient?
However, if we do not use sham surgeries, we may promote the use of a costly treatment that has no real effect; if this happens, money and other resources will be diverted away from other treatments that are known to be helpful.
Ultimately, this is a difficult situation where we cannot perfectly protect both the patients who have volunteered for the study and the patients who may benefit (or not) from the treatment in the future.

## Observational studies {#sec-observational-studies}

Data where no treatment has been explicitly applied (or explicitly withheld) is called **observational data**.
For instance, the loan data and county data described in @sec-data-basics are both examples of observational data.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "observational data")
```

Making causal conclusions based on experiments is often reasonable, since we can randomly assign the explanatory variable(s), i.e., the treatments.
However, making the same causal conclusions based on observational data can be treacherous and is not recommended.
Thus, observational studies are generally only sufficient to show associations or form hypotheses that can be later checked with experiments.

::: {.guidedpractice data-latex=""}
Suppose an observational study tracked sunscreen use and skin cancer, and it was found that the more sunscreen someone used, the more likely the person was to have skin cancer.
Does this mean sunscreen *causes* skin cancer?[^02-data-design-9]
:::

[^02-data-design-9]: No.
    See the paragraph following the question!

Some previous research tells us that using sunscreen actually reduces skin cancer risk, so maybe there is another variable that can explain this hypothetical association between sunscreen usage and skin cancer.
One important piece of information that is absent is sun exposure.
If someone is out in the sun all day, they are more likely to use sunscreen *and* more likely to get skin cancer.
Exposure to the sun is unaccounted for in the simple observational investigation.

```{r}
#| label: sun-causes-cancer
#| fig-asp: 0.4
#| fig-alt: 'Three boxes are shown in a triangle arrangement representing: sun exposure,
#|   using sunscreen, and skin cancer.  A solid arrow connects sun exposure as a causal
#|   mechanism to using sunscreen; a solid arrow also connects sun exposure as a causal
#|   mechanism to skin cancer.  A questioning arrow indicates that the causal effect
#|   of using sunscreen on skin cancer is unknown.'
#| out.width: 60%
par_og <- par(no.readonly = TRUE) # save original par
par(mar = rep(0, 4))
plot(c(-0.05, 1.2),
     c(0.39, 1),
     type = 'n',
     axes = FALSE)
text(0.59, 0.89, 'sun exposure', cex = 1)
rect(0.4, 0.8, 0.78, 1)
text(0.3, 0.49, 'use sunscreen', cex = 1)
rect(0.1, 0.4, 0.48, 0.6)
arrows(0.49, 0.78, 0.38, 0.62,
       length = 0.08, lwd = 1.5)
text(0.87, 0.5, 'skin cancer', cex = 1)
rect(0.71,0.4, 1.01, 0.6)
arrows(0.67, 0.78, 0.8, 0.62,
       length = 0.08, lwd = 1.5)

arrows(0.5, 0.5, 0.69, 0.5,
       length = 0.08, col = IMSCOL["gray", "f1"])
text(0.595, 0.565, "?",
     cex = 1.5, col = IMSCOL["gray", "full"])
par(par_og) # restore original par
```

In this example, sun exposure is a confounding variable.
The presence of confounding variables is what inhibits the ability for observational studies to make causal claims.
While one method to justify making causal conclusions from observational studies is to exhaust the search for confounding variables, there is no guarantee that all confounding variables can be examined or measured.

::: {.guidedpractice data-latex=""}
@fig-county-multi-unit-homeownership shows a negative association between the homeownership rate and the percentage of housing units that are in multi-unit structures in a county.
However, it is unreasonable to conclude that there is a causal relationship between the two variables.
Suggest a variable that might explain the negative relationship.[^02-data-design-10]
:::

[^02-data-design-10]: Answers will vary.
    Population density may be important.
    If a county is very dense, then this may require a larger percentage of residents to live in housing units that are in multi-unit structures.
    Additionally, the high density may contribute to increases in property value, making homeownership unfeasible for many residents.

Observational studies come in two forms: prospective and retrospective studies.
A **prospective study** identifies individuals and collects information as events unfold.
For instance, medical researchers may identify and follow a group of patients over many years to assess the possible influences of behavior on cancer risk.
One example of such a study is The Nurses' Health Study.
Started in 1976 and expanded in 1989, the Nurses' Health Study has collected data on over 275,000 nurses and is still enrolling participants.
This prospective study recruits registered nurses and then collects data from them using questionnaires.
**Retrospective studies** collect data after events have taken place, e.g., researchers may review past events in medical records.
Some datasets may contain both prospectively- and retrospectively collected variables, such as medical studies which gather information on participants' lives before they enter the study and subsequently collect data on participants throughout the study.

```{r}
#| include: false
terms_chp_2 <- c(terms_chp_2, "prospective study", "retrospective study")
```

\clearpage

## Chapter review {#sec-chp2-review}

### Summary

A strong analyst will have a good sense of the types of data they are working with and how to visualize the data in order to gain a complete understanding of the variables.
Equally important however, is an understanding of the data source.
In this chapter, we have discussed randomized experiments and taking good, random, representative samples from a population.
When we discuss inferential methods (starting in @sec-foundations-randomization), the conclusions that can be drawn will be dependent on how the data were collected.
@fig-randsampValloc summarizes the differences between random assignment of treatments and random samples.[^02-data-design-11]
Regularly revisiting @fig-randsampValloc will be important when making conclusions from a given data analysis.

[^02-data-design-11]: Derived from similar figures in @ISCAM and @sleuth.

```{r}
#| label: fig-randsampValloc
#| out.width: 100%
#| fig-cap: As we will see, analysis conclusions should be made carefully according to
#|   how the data were collected.  Note that very few datasets come from the top left
#|   box because usually ethics require that random assignment of treatments can only
#|   be given to volunteers. Both representative (ideally random) sampling and experiments
#|   (random assignment of treatments) are important for how statistical conclusions
#|   can be made on populations.
knitr::include_graphics("images/randsampValloc.png")
```

### Terms

We introduced the following terms in the chapter.
If you're not sure what some of these terms mean, we recommend you go back in the text and review their definitions.
We are purposefully presenting them in alphabetical order, instead of in order of appearance, so they will be a little more challenging to locate.
However, you should be able to easily spot them as **bolded text**.

```{r}
make_terms_table(terms_chp_2)
```

\clearpage

## Exercises {#sec-chp2-exercises}

Answers to odd-numbered exercises can be found in [Appendix -@sec-exercise-solutions-02].

::: {.exercises data-latex=""}
{{< include exercises/_02-ex-data-design.qmd >}}
:::