knitr::opts_chunk$set( collapse = TRUE, warning = FALSE, message = FALSE, echo = FALSE, comment = "#>", fig.path = "../figures/", dpi = 300 )
Here, we reproduce figure 2b from [@nowicki2020] - "Gene expression correlates of social evolution in coral reef butterfly fishes". This study analyzed the neural correlates of pair-bonding vs. solitary mating behavior in several species of Chaetodontidae fishes. To accomplish this, they compared gene expression of socially-relevant neuropeptides in brain regions pertinent to social function. One of their main findings was that the males of pair-bonding species had elevated oxytocin receptor (OTR) expression in the supracommissural part of the ventral telencephalon (Vs) compared to the males of solitary species.
The authors of [@nowicki2020] used qRT-PCR to quantify relative mRNA expression of OTR in Vs of 5 Chaetodontidae species: C. vagabundus, C. lunatus, C. baronessa, C. trifascialis, and C. rainfordi. They employ a Markov Chain Monte Carlo approach to analyze differential gene expression. Below is a reproduction of figure 2b using our package buttRflyfish
.
library(buttRflyfish) data = load_data() # load data model = run_model(data) # run MCMC model p = otr_expression_plot(model,data) # plot differential expression results
::: {#refs} :::
\newpage
This report was generated on r Sys.time()
using the following computational environment and dependencies:
# which R packages and versions? if ("devtools" %in% installed.packages()) devtools::session_info()
The current Git commit details are:
# what commit is this file at? if ("git2r" %in% installed.packages() & git2r::in_repository(path = ".")) git2r::repository(here::here())
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.