Independent Study to prepare for workshop
Data Analysis 4: Immunobiology - Your data presentation
The independent study is to check your progress and understanding. We have covered the steps listed below. If you have got behind, start from wherever you got to.
📖 Read what the sample data are, get an overview of the analysis and tools
💻 Go through the sample data analysis and make sure you understand the steps and the code
You should have a single csv file with data from all 6 samples. The data has been cleaned with
flowAI
to remove cells ‘bad’ signals. The file has 19 columns with meaningful names. The TNFa_APC_Lin and E_coli_FITC_Lin columns have been logicle transformed. The file names were added to a columnsample
and columns indicating the treatment and antibody were created by extracting patterns from the file names. The code that does that relies on the names being in the format Treatment_Antibody. You can import this data as a dataframe and use familiar tidyverse tools.The treatments will be plotted in alphabetical order by default. We can change that order to “MEDIA”, “LPS”, “ECOLIGreen” using
fct_relevel()
. The code that does that relies on the treatment names being exactly “MEDIA”, “LPS”, “ECOLIGreen”. You can name the treaments differently but you will need to change the levels given infct_relevel()
accordingly.You should be able to plot side scatter against forward scatter for the six samples and add a rectangular gate for the live cells. You should be able to use this gate to filter the data to get a dataframe of live cells.
You should be able calculate the number and percentage of live cells in each sample and annotate the figure in iii. with those percentages. These were in a dataframe I named
clean_trans_nondebris_n
.You should be able to plot the logicle transformed TNFa_APC_Lin against the E_coli_FITC_Lin for the live cells and add a gate for the TNF-α positive cells and label the rows (cells) in the dataframe as TNF-α positive or TNF-α negative. You should be able to calculate the number and percentage of of TNF-α positive cells in each sample. These were in a dataframe I named
clean_trans_nondebris_tfna_pos
You should be able add a gate for the FITC positive cells and label the rows (cells) in the dataframe as FITC positive or FITC negative. You should be able to calculate the number and percentage of FITC positive cells in each sample using the same logic as in v.
💻 Prepare to analyse your own data / the model data
Make a new RStudio Project, copy in the script and tidy it up
Save your data. Your data are on googledrive. If you do not have any data you can use the Flow cytometry MODEL DATA
💻 Analyse your own data / the model data
Make your life easier renaming your data files to match the sample data.
Note that your data has 16 columns. You need a different metadata file than used for the sample data. You will not need to Drop the unused channels
💻 Enter data from your analysis in to BIO00066I Biomedical Sciences class data
The columns you must add are:
apc_mfi
: Mean fluorescence intensity of the logicle transformed TNFa_APC_Lin in the TNF-α positive cellsperc_tfna_pos
: % non debris cells that are TNF-α positive cells
The other columns are calculations you make along the way and may help you get to the
apc_mfi
andperc_tfna_pos
values. The column names are the same as those used in the Data Analysis 2: Immunobiology - Sample data analysis workshop