Getting Started with Experimental Design in R

This quick blog is designed to help you get off to the races quickly in world of data science; and here specifically, Experimental design. Enjoy!

When it comes to experiemental design there are three main streps it can be broken down to:

  • Planning
  • Design
  • Analysis

Planning & Design

Planning should always begin with a well formed hypothesis.

Some of the major considerations you want to make in this process are as follows:

  • What is the question you want answered?
  • What is the population in question?
  • What are your dependent and independent variables?

Analysis

When conducting an experiment, there are 3 key components to consider. These three aspects of an experiment allows us to assess our population’s variability.

  • Randomization
  • Replication
  • Blocking

Randomization

The purpose of randomization is to make sure that if there is variation in outcomes that is related to outside factors, then it is distributed across treatment groups.

Replication

When conducting an experiment, we seek to understand the variability of outcomes. For instance; if I were to run a given experiment but once, and I was depending on an outcome that may have occurred due to random chance. The point here is to understand the broad spectrum of possibilities or outcomes, it’s important that we replicate the experiment accordingly.

Statistical Power

The concept of statistical power means that if your experiment concludes such that you reject the NULL hypothesis and accept the alternative hypothesis, it is the likelihood that it would not be due to random chance. Best practice is 80% statistical power.

So to simplify this even further; If your hypothesis appears to be correct, what’s the likelihood that you didn’t get that outcome just due to random chance.

Blocking

Blocking is used to help control variability by making treatment groups more alike. Inside of a given group, you might see that differences are minimal, however across other groups that could be much larger. One example of this might be blocking an experiment by gender.

  • Blocking by variables – use aov for the sake of blocking
  • Randomized Complete Block Design (RCBD) experiment

t-test

After accumulating data from your experiment; one quick and east test of statistical significance you might run is called a t-test.

  1. Consider your hypothesis or central research question:
    1. NULL hypothesis – lets keep this simple, the null hypothesis is pretty much when you’re wrong. For the mtcars dataset, the null hypothesis might be something like a vehicles horse power has no effect on miles per gallon.
    2. Alternative hypothesis – conversely the alternative hypothesis means that there was a difference. If you are able to determine with statistical significance the effect of the independent variable on the dependent variable, you would say that you reject the null hypothesis and accept the alternative hypothesis.
  2. Is this a one or two sided test?
    1. One sided test – when you are testing whether a given variable is greater than another then it’s a one-sided test; if you’re testing whether it’s less than another… still one-sided.
    2. Two sided test – when you are testing that a given variable is not equal to another, then that is two sided. Greater or less than in a single test.
  3. Were your results statistically significant?
    1. People use the term statistically significant left and right with little consideration of what it actually means. What that is saying is that if you run your test and your data is suggesting a that your hypothesis is correct, statistical significance is effectively knowing that it’s not likely due to random chance.
    2. The standard here is 95% confidence, or a less than or equal to likelihood of 5%.
  4. What is statistical power?
    1. Similar to statistical significance; given that the alternative hypothesis is true, power represents the likelihood that the null hypothesis will be rejected.
    2. The standard for power is 80%.

Sample Size

For a given experiment one thing to consider is the sample size. In order to arrive at a required number for this is requires a handful of other variables including targeted statistical power & significance.

Another measure is that of effect size. Effect size represents the difference between the average of two groups divided by the standard deviation of both groups combined.

The greater the distance between groups the less of a sample to validate it. The smaller the difference the greater the likelihood that the observed distance is only due to chance.

In order to calculate any of these values including effect size, statistical power, p-value, etc. you need all buy one of them.

Load up the package pwr and use the pwr.anova.test to identify the odd variable out here.

k – number of groups
n – sample size per group
f = effect size
sig.level = significance level
power = statistical power

library(pwr)
 pwr.anova.test(k = 2,
                n = NULL,
                f = .1,
                sig.level = 0.05,
                power = .8)

As you can see above, given the demands that we passed to the pwr.anova.test function, we can see that for those things to hold true, we at least 394 samples per group.

Conclusion

I hope this helps to get you started in experimental design!

Join the Conversation

5 Comments

  1. I bought online courses that have your explanation on it! They simple copied it and are now selling it by 250UDS !!! How could I send you an email with the plagiarism content ?

    Like

Leave a comment

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: