# Executive Summary

# Maternal Preconception Study

## Fecal Testosterone

# Introduction

The results of the analysis are presented using the following structure:

1. Searchable High‑Level Summary: Provides a brief overview of the key findings from our analysis with a table of effects
2. Visualisation, Descriptive Statistics and Diagnostics by Variable: Exploratory data analysis and model diagnostics for each behavioral measure
3. Searchable ROPE Analysis Table: Model diagnostics to check the fit of the models

# Effects of Interest

- All groups except for CRV have higher fecal testosterone following resident intruder behavioural testing.
- CRA has higher fecal testosterone, particularly prior to resident intruder behavioural testing.

# Searchable High‑Level Summary

### Enter search terms in the column headers to filter the table below

Note: TD in the Higher column indicates a Trend Difference for the continuous variable.

# Visualisation, Descriptive Statistics and Diagnostics by Variable

### Bar Plot

### Density Plot

Fecal Testosterone by Grouping Variables
Diet | Time | n | Mean | SD | Min | Max | % of Zeros | Zero Flag
--- | --- | --- | --- | --- | --- | --- | --- | ---
Control | 0 | 11.0 | 20.5 | 13.7 | 0.5 | 41.9 | 0.0 |
Control | 1 | 11.0 | 36.7 | 13.4 | 14.0 | 59.2 | 0.0 |
CR25 | 0 | 11.0 | 17.7 | 14.6 | 5.7 | 55.7 | 0.0 |
CR25 | 1 | 11.0 | 48.4 | 29.5 | 17.4 | 99.7 | 0.0 |
CRA | 0 | 10.0 | 31.9 | 16.1 | 9.6 | 65.9 | 0.0 |
CRA | 1 | 10.0 | 50.2 | 21.3 | 19.5 | 85.0 | 0.0 |
CRV | 0 | 9.0 | 29.9 | 14.9 | 15.8 | 58.0 | 0.0 |
CRV | 1 | 9.0 | 35.8 | 15.2 | 15.9 | 58.8 | 0.0 |
Warning symbols: * = >40% zeros in this group, ! = >20% zeros in entire sample Either flag triggers the use of a binary choice model.

### Summary Table

Model Summary
Family: gamma
Parameter | Estimate | Std. Error | 2.5% CI | 97.5% CI | Type
--- | --- | --- | --- | --- | ---
Fixed Effects
Intercept | 2.987 | 0.184 | 2.624 | 3.355 | Fixed Effects
DietCR25 | −0.156 | 0.251 | −0.649 | 0.340 | Fixed Effects
DietCRA | 0.446 | 0.256 | −0.050 | 0.954 | Fixed Effects
DietCRV | 0.368 | 0.261 | −0.138 | 0.884 | Fixed Effects
Time | 0.603 | 0.214 | 0.177 | 1.011 | Fixed Effects
DietCR25:Time | 0.412 | 0.294 | −0.160 | 0.994 | Fixed Effects
DietCRA:Time | −0.117 | 0.309 | −0.710 | 0.501 | Fixed Effects
DietCRV:Time | −0.378 | 0.311 | −0.975 | 0.249 | Fixed Effects

### Posterior Check

$data # A tibble: 902 × 7 y_id y_name rep_id rep_label is_y is_y_label value 1 1 1 1 italic(y)[rep] ( 1 ) FALSE italic(y)[rep] 21.2 2 1 1 2 italic(y)[rep] ( 2 ) FALSE italic(y)[rep] 21.5 3 1 1 3 italic(y)[rep] ( 3 ) FALSE italic(y)[rep] 20.9 4 1 1 4 italic(y)[rep] ( 4 ) FALSE italic(y)[rep] 11.8 5 1 1 5 italic(y)[rep] ( 5 ) FALSE italic(y)[rep] 9.10 6 1 1 6 italic(y)[rep] ( 6 ) FALSE italic(y)[rep] 12.2 7 1 1 7 italic(y)[rep] ( 7 ) FALSE italic(y)[rep] 13.9 8 1 1 8 italic(y)[rep] ( 8 ) FALSE italic(y)[rep] 16.1 9 1 1 9 italic(y)[rep] ( 9 ) FALSE italic(y)[rep] 11.2 10 1 1 10 italic(y)[rep] ( 10 ) FALSE italic(y)[rep] 24.0 # ℹ 892 more rows

\(layers\)layers[[1]] mapping: group = ~.data$rep_id, colour = yrep geom_line: na.rm = FALSE, orientation = NA stat_density: bw = nrd0, adjust = 1, kernel = gaussian, n = 1024, trim = FALSE, na.rm = FALSE, bounds = c(-Inf, Inf), orientation = NA position_identity

$layers[[2]] mapping: colour = y geom_line: na.rm = FALSE, orientation = NA, lineend = round stat_density: bw = nrd0, adjust = 1, kernel = gaussian, n = 1024, trim = FALSE, na.rm = FALSE, bounds = c(-Inf, Inf), orientation = NA position_identity

$scales <ggproto object: Class ScalesList, gg> add: function add_defaults: function add_missing: function backtransform_df: function clone: function find: function get_scales: function has_scale: function input: function map_df: function n: function non_position_scales: function scales: list set_palettes: function train_df: function transform_df: function super: <ggproto object: Class ScalesList, gg>

$guides <Guides[0] ggproto object>

\(mapping\)x expr: ^.data$value env: 0x000001e4e4e1e200

attr(,“class”) [1] “uneval”

\(theme\)line $colour [1] “black”

$linewidth [1] 0.5

$linetype [1] 1

$lineend [1] “butt”

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

\(rect\)fill [1] “white”

$colour [1] “black”

$linewidth [1] 0.5

$linetype [1] 1

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_rect” “element”

\(text\)family [1] “”

$face [1] “plain”

$colour [1] “black”

$size [1] 11

$hjust [1] 0.5

$vjust [1] 0.5

$angle [1] 0

$lineheight [1] 0.9

$margin [1] 0points 0points 0points 0points

$debug [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$title NULL

$aspect.ratio NULL

$axis.title NULL

\(axis.title.x\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 2.75points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.title.x.top\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 0

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 2.75points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.title.x.bottom NULL

\(axis.title.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle [1] 90

$lineheight NULL

$margin [1] 0points 2.75points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.title.y.left NULL

\(axis.title.y.right\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle [1] -90

$lineheight NULL

$margin [1] 0points 0points 0points 2.75points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text\)family NULL

$face NULL

$colour [1] “grey30”

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text.x\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 2.2points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(axis.text.x.top\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust [1] 0

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 2.2points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.x.bottom NULL

\(axis.text.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 1

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 2.2points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.y.left NULL

\(axis.text.y.right\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 0points 2.2points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.text.theta NULL

\(axis.text.r\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0.5

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 0points 2.2points 0points 2.2points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$axis.ticks list() attr(,“class”) [1] “element_blank” “element”

$axis.ticks.x NULL

$axis.ticks.x.top NULL

$axis.ticks.x.bottom NULL

$axis.ticks.y NULL

$axis.ticks.y.left NULL

$axis.ticks.y.right NULL

$axis.ticks.theta NULL

$axis.ticks.r NULL

$axis.minor.ticks.x.top NULL

$axis.minor.ticks.x.bottom NULL

$axis.minor.ticks.y.left NULL

$axis.minor.ticks.y.right NULL

$axis.minor.ticks.theta NULL

$axis.minor.ticks.r NULL

$axis.ticks.length [1] 2.75points

$axis.ticks.length.x NULL

$axis.ticks.length.x.top NULL

$axis.ticks.length.x.bottom NULL

$axis.ticks.length.y NULL

$axis.ticks.length.y.left NULL

$axis.ticks.length.y.right NULL

$axis.ticks.length.theta NULL

$axis.ticks.length.r NULL

$axis.minor.ticks.length [1] 0.75 *

$axis.minor.ticks.length.x NULL

$axis.minor.ticks.length.x.top NULL

$axis.minor.ticks.length.x.bottom NULL

$axis.minor.ticks.length.y NULL

$axis.minor.ticks.length.y.left NULL

$axis.minor.ticks.length.y.right NULL

$axis.minor.ticks.length.theta NULL

$axis.minor.ticks.length.r NULL

$axis.line list() attr(,“class”) [1] “element_blank” “element”

$axis.line.x NULL

$axis.line.x.top NULL

$axis.line.x.bottom NULL

$axis.line.y NULL

$axis.line.y.left NULL

$axis.line.y.right NULL

$axis.line.theta NULL

$axis.line.r NULL

$legend.background list() attr(,“class”) [1] “element_blank” “element”

$legend.margin [1] 5.5points 5.5points 5.5points 5.5points

$legend.spacing [1] 11points

$legend.spacing.x NULL

$legend.spacing.y NULL

$legend.key list() attr(,“class”) [1] “element_blank” “element”

$legend.key.size [1] 1.2lines

$legend.key.height NULL

$legend.key.width NULL

$legend.key.spacing [1] 5.5points

$legend.key.spacing.x NULL

$legend.key.spacing.y NULL

$legend.frame NULL

$legend.ticks NULL

$legend.ticks.length [1] 0.2 *

$legend.axis.line NULL

\(legend.text\)family NULL

$face NULL

$colour NULL

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$legend.text.position NULL

\(legend.title\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust NULL

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$legend.title.position NULL

$legend.position [1] “right”

$legend.position.inside NULL

$legend.direction NULL

$legend.byrow NULL

$legend.justification [1] “center”

$legend.justification.top NULL

$legend.justification.bottom NULL

$legend.justification.left NULL

$legend.justification.right NULL

$legend.justification.inside NULL

$legend.location NULL

$legend.box NULL

$legend.box.just NULL

$legend.box.margin [1] 0cm 0cm 0cm 0cm

$legend.box.background list() attr(,“class”) [1] “element_blank” “element”

$legend.box.spacing [1] 11points

$panel.background list() attr(,“class”) [1] “element_blank” “element”

$panel.border list() attr(,“class”) [1] “element_blank” “element”

$panel.spacing [1] 5.5points

$panel.spacing.x NULL

$panel.spacing.y NULL

\(panel.grid\)colour [1] “grey92”

$linewidth NULL

$linetype NULL

$lineend NULL

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

$panel.grid.major NULL

\(panel.grid.minor\)colour NULL

$linewidth [1] 0.5 *

$linetype NULL

$lineend NULL

$arrow [1] FALSE

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_line” “element”

$panel.grid.major.x NULL

$panel.grid.major.y NULL

$panel.grid.minor.x NULL

$panel.grid.minor.y NULL

$panel.ontop [1] FALSE

$plot.background list() attr(,“class”) [1] “element_blank” “element”

\(plot.title\)family NULL

$face NULL

$colour NULL

$size [1] 1.2 *

$hjust [1] 0.5

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 5.5points 0points

$debug NULL

$inherit.blank [1] FALSE

attr(,“class”) [1] “element_text” “element”

$plot.title.position [1] “panel”

\(plot.subtitle\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust [1] 0

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 0points 0points 5.5points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(plot.caption\)family NULL

$face NULL

$colour NULL

$size [1] 0.8 *

$hjust [1] 1

$vjust [1] 1

$angle NULL

$lineheight NULL

$margin [1] 5.5points 0points 0points 0points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$plot.caption.position [1] “panel”

\(plot.tag\)family NULL

$face NULL

$colour NULL

$size [1] 1.2 *

$hjust [1] 0.5

$vjust [1] 0.5

$angle NULL

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$plot.tag.position [1] “topleft”

$plot.tag.location NULL

$plot.margin [1] 5.5points 5.5points 5.5points 5.5points

$strip.background list() attr(,“class”) [1] “element_blank” “element”

$strip.background.x NULL

$strip.background.y NULL

$strip.clip [1] “inherit”

$strip.placement [1] “inside”

\(strip.text\)family NULL

$face NULL

$colour [1] “grey10”

$size [1] 0.8 *

$hjust NULL

$vjust NULL

$angle NULL

$lineheight NULL

$margin [1] 4.4points 4.4points 4.4points 4.4points

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$strip.text.x NULL

$strip.text.x.bottom NULL

$strip.text.x.top NULL

\(strip.text.y\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust NULL

$angle [1] -90

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

\(strip.text.y.left\)family NULL

$face NULL

$colour NULL

$size NULL

$hjust NULL

$vjust NULL

$angle [1] 90

$lineheight NULL

$margin NULL

$debug NULL

$inherit.blank [1] TRUE

attr(,“class”) [1] “element_text” “element”

$strip.text.y.right NULL

$strip.switch.pad.grid [1] 2.75points

$strip.switch.pad.wrap [1] 2.75points

attr(,“class”) [1] “theme” “gg”
attr(,“complete”) [1] TRUE attr(,“validate”) [1] TRUE

$coordinates <ggproto object: Class CoordCartesian, Coord, gg> aspect: function backtransform_range: function clip: on default: FALSE distance: function draw_panel: function expand: FALSE is_free: function is_linear: function labels: function limits: list modify_scales: function range: function render_axis_h: function render_axis_v: function render_bg: function render_fg: function reverse: none setup_data: function setup_layout: function setup_panel_guides: function setup_panel_params: function setup_params: function train_panel_guides: function transform: function super: <ggproto object: Class CoordCartesian, Coord, gg>

$facet <ggproto object: Class FacetNull, Facet, gg> attach_axes: function attach_strips: function compute_layout: function draw_back: function draw_front: function draw_labels: function draw_panel_content: function draw_panels: function finish_data: function format_strip_labels: function init_gtable: function init_scales: function map_data: function params: list set_panel_size: function setup_data: function setup_panel_params: function setup_params: function shrink: TRUE train_scales: function vars: function super: <ggproto object: Class FacetNull, Facet, gg>

$plot_env <environment: 0x000001e4e4e1e200>

$layout <ggproto object: Class Layout, gg> coord: NULL coord_params: list facet: NULL facet_params: list finish_data: function get_scales: function layout: NULL map_position: function panel_params: NULL panel_scales_x: NULL panel_scales_y: NULL render: function render_labels: function reset_scales: function resolve_label: function setup: function setup_panel_guides: function setup_panel_params: function train_position: function super: <ggproto object: Class Layout, gg>

\(labels\)labels$title [1] “Posterior Predictive Check”

\(labels\)x NULL

\(labels\)y NULL

\(labels\)group [1] “rep_id”

\(labels\)colour [1] “colour”

\(labels\)fill [1] “fill” attr(,“fallback”) [1] TRUE

\(labels\)weight [1] “weight” attr(,“fallback”) [1] TRUE

attr(,“class”) [1] “gg” “ggplot”

### Convergence

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Diet: Control vs CR25
Diet: Control vs CR25 | -0.05 | [-0.45, 0.38] | 61.00 | 28.25
Diet: Control vs CRA
Diet: Control vs CRA | -0.38 | [-0.80, 0.03] | 96.90 | 4.85
Diet: Control vs CRV
Diet: Control vs CRV | -0.18 | [-0.64, 0.23] | 79.75 | 20.85
Diet: CR25 vs CRA
Diet: CR25 vs CRA | -0.34 | [-0.74, 0.10] | 94.50 | 8.85
Diet: CR25 vs CRV
Diet: CR25 vs CRV | -0.13 | [-0.57, 0.30] | 72.20 | 23.55
Diet: CRA vs CRV
Diet: CRA vs CRV | 0.21 | [-0.21, 0.67] | 82.55 | 19.65

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time: Overall Change
Time: Overall Change 0 - 1 | 0.58 | [0.35, 0.82] | 100 | 0

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time × Diet: Control - CR25
Time × Diet: Control - CR25 0 | 0.15 | [-0.37, 0.60] | 72.60 | 22.55
Time × Diet: Control - CR25 1 | -0.26 | [-0.81, 0.24] | 82.15 | 9.75
Time × Diet: Control - CRA
Time × Diet: Control - CRA 0 | -0.44 | [-0.91, 0.08] | 96.00 | 5.15
Time × Diet: Control - CRA 1 | -0.33 | [-0.87, 0.16] | 90.70 | 5.30
Time × Diet: Control - CRV
Time × Diet: Control - CRV 0 | -0.38 | [-0.88, 0.13] | 93.25 | 8.65
Time × Diet: Control - CRV 1 | 0.01 | [-0.54, 0.52] | 50.80 | 15.50
Time × Diet: CR25 - CRA
Time × Diet: CR25 - CRA 0 | -0.60 | [-1.10, -0.06] | 98.55 | 2.20
Time × Diet: CR25 - CRA 1 | -0.08 | [-0.61, 0.44] | 61.05 | 14.00
Time × Diet: CR25 - CRV
Time × Diet: CR25 - CRV 0 | -0.53 | [-1.11, -0.06] | 97.60 | 3.80
Time × Diet: CR25 - CRV 1 | 0.25 | [-0.28, 0.76] | 82.50 | 9.50
Time × Diet: CRA - CRV
Time × Diet: CRA - CRV 0 | 0.07 | [-0.49, 0.58] | 59.85 | 25.30
Time × Diet: CRA - CRV 1 | 0.34 | [-0.21, 0.90] | 88.40 | 6.70

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Time: Control
Time: Control 0 - 1 | 0.60 | [0.18, 1.06] | 99.70 | 0.60
Time: CR25
Time: CR25 0 - 1 | 1.01 | [0.55, 1.44] | 100.00 | 0.00
Time: CRA
Time: CRA 0 - 1 | 0.49 | [0.00, 0.90] | 98.35 | 2.30
Time: CRV
Time: CRV 0 - 1 | 0.22 | [-0.24, 0.70] | 82.00 | 10.45

Effect | Median | [95% HDI] | Direction Prob. | ROPE %
--- | --- | --- | --- | ---
Diet: Control vs CR25 | Time = 0
Diet: Control vs CR25 | Time = 0 0 | 0.15 | [-0.37, 0.60] | 72.60 | 22.55
Diet: Control vs CRA | Time = 0
Diet: Control vs CRA | Time = 0 0 | -0.44 | [-0.91, 0.08] | 96.00 | 5.15
Diet: Control vs CRV | Time = 0
Diet: Control vs CRV | Time = 0 0 | -0.38 | [-0.88, 0.13] | 93.25 | 8.65
Diet: CR25 vs CRA | Time = 0
Diet: CR25 vs CRA | Time = 0 0 | -0.60 | [-1.10, -0.06] | 98.55 | 2.20
Diet: CR25 vs CRV | Time = 0
Diet: CR25 vs CRV | Time = 0 0 | -0.53 | [-1.11, -0.06] | 97.60 | 3.80
Diet: CRA vs CRV | Time = 0
Diet: CRA vs CRV | Time = 0 0 | 0.07 | [-0.49, 0.58] | 59.85 | 25.30
Diet: Control vs CR25 | Time = 1
Diet: Control vs CR25 | Time = 1 1 | -0.26 | [-0.81, 0.24] | 82.15 | 9.75
Diet: Control vs CRA | Time = 1
Diet: Control vs CRA | Time = 1 1 | -0.33 | [-0.87, 0.16] | 90.70 | 5.30
Diet: Control vs CRV | Time = 1
Diet: Control vs CRV | Time = 1 1 | 0.01 | [-0.54, 0.52] | 50.80 | 15.50
Diet: CR25 vs CRA | Time = 1
Diet: CR25 vs CRA | Time = 1 1 | -0.08 | [-0.61, 0.44] | 61.05 | 14.00
Diet: CR25 vs CRV | Time = 1
Diet: CR25 vs CRV | Time = 1 1 | 0.25 | [-0.28, 0.76] | 82.50 | 9.50
Diet: CRA vs CRV | Time = 1
Diet: CRA vs CRV | Time = 1 1 | 0.34 | [-0.21, 0.90] | 88.40 | 6.70

# End of Document