Draw Disproportionate Sample from Data Frame in R (Example)

 

This article demonstrates how to draw a sample with different probabilities by group in R programming.

Table of contents:

Let’s dig in…

 

Creation of Example Data

First of all, we’ll need to define some data that we can use in the following examples:

data <- data.frame(value = 1:50,                             # Create example data frame
                   group = rep(letters[1:5], each = 10))
head(data)                                                   # Head of example data frame

 

table 1 data frame draw disproportionate sample from data frame r

 

Table 1 shows that our example data is composed of two columns called “value” and “group”. The variable value is an integer and the column group is a character.

 

Example: Create Random Sample of Data Frame with Multiple Probabilities

In this example, I’ll explain how to create a random subsample of a data frame with different sampling probabilities by group.

To achieve this, we first have to specify a vector of probabilities that has the same length as the number of rows in our data frame:

my_prob <- rep(NA, 50)                                       # Create vector of probabilities
my_prob[data$group == "a" | data$group == "b"] <- 0.05
my_prob[data$group == "c" | data$group == "d"] <- 0.1
my_prob[data$group == "e"] <- 0.7
my_prob                                                      # Print vector of probabilities
#  [1] 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05
# [16] 0.05 0.05 0.05 0.05 0.05 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10
# [31] 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.70 0.70 0.70 0.70 0.70
# [46] 0.70 0.70 0.70 0.70 0.70

As you can see based on the previous output of the RStudio console, our vector of probabilities contains 50 values.

We have specified that we want to draw rows of the groups a and b with a probability of 5%, rows of the groups c and d with a likelihood of 10%, and rows of the group e with a probability of 70%.

In other words: the probabilities to draw a group are disproportionate.

Next, we should specify a random seed to make the following data sampling process reproducible:

set.seed(239678564)                                          # Set random seed

Now, we can apply the sample function to generate a random subsample of our data frame. Note that we are specifying the prob argument to be equal to the vector of probabilities that we have created before:

data_samp <- data[sample(nrow(data), 10, prob = my_prob), ]  # Draw sample of data frame
data_samp                                                    # Print sample of data frame

 

table 2 data frame draw disproportionate sample from data frame r

 

After executing the previous code the data frame subsample shown in Table 2 has been created.

As you can see, the group e was drawn the most often, since it had the highest chance of being drawn (i.e. 70%). In contrast, the groups a and c have not been selected at all.

Note that we have generated our sample without replacement. In case you want to draw a sample with replacement, you may specify the replace argument within the sample function to be equal to TRUE.

 

Video, Further Resources & Summary

I have recently released a video on my YouTube channel, which illustrates the contents of this tutorial. You can find the video instruction below:

 

The YouTube video will be added soon.

 

Additionally, you might want to have a look at some of the other tutorials on my website. A selection of articles that are related to the creation of a sample with multiple probabilities by group can be found below:

 

In this R tutorial you have learned how to take a random sample with multiple disproportionate probabilities. Let me know in the comments section, if you have any additional questions or comments.

 

Subscribe to the Statistics Globe Newsletter

Get regular updates on the latest tutorials, offers & news at Statistics Globe.
I hate spam & you may opt out anytime: Privacy Policy.


2 Comments. Leave new

Leave a Reply

Your email address will not be published.

Fill out this field
Fill out this field
Please enter a valid email address.

Menu
Top