Interpreting Gas Production Models

Introduction

Fitting a model is only the first step in the analysis of rumen gas production data.

Researchers must also interpret:

This vignette summarizes the most common interpretations used in rumen gas production studies.

library(rumenGP)

Understanding Common Parameters

Although different models use different equations, many share similar biological concepts.

Asymptotic Gas Production

Common parameter names:

A
VF
Vf
V1F
V2F

These parameters represent the maximum gas production that the model predicts after long incubation times.

Example:

A = 120 mL

Interpretation:

The model predicts approximately
120 mL of gas at fermentation completion.

Higher values generally indicate:

However, interpretation should always be made within the context of the substrate being studied.

Fermentation Rate

Common parameter names:

k
k1
k2
mu

These parameters describe how rapidly gas production approaches the asymptote.

Example:

Treatment A
k = 0.08

Treatment B
k = 0.04

Interpretation:

Treatment A ferments more rapidly
than Treatment B.

Higher rates generally suggest:

Lag Time

Common parameter name:

lambda

or:

\[ \lambda \]

Lag time represents the delay before substantial fermentation begins.

Example:

lambda = 2 h

Interpretation:

Approximately two hours are required
before active fermentation starts.

Large lag values often occur with:

Half-Time Parameters

Common parameter names:

b
K

Used in:

These parameters determine the time required to achieve approximately half of the asymptotic gas production.

Example:

K = 12 h

Interpretation:

Approximately 50% of total gas production
is achieved after 12 hours.

Smaller values indicate faster fermentation.

Shape Parameters

Common parameter names:

c
d
m

Shape parameters modify the curvature of the fermentation profile.

Interpretation:

Shape parameters control how fermentation
accelerates and decelerates through time.

Unlike asymptotes or rates, shape parameters often have no simple biological interpretation.

They are usually considered:

Empirical flexibility parameters.

Interpreting Dual-Pool Models

Dual-pool models separate fermentation into:

Rapid fraction
Slow fraction

Parameters:

V1F
V2F
k1
k2

Rapid Fraction

V1F
k1

Typically associated with:

Slow Fraction

V2F
k2

Typically associated with:

Example:

V1F = 30 mL

V2F = 90 mL

Interpretation:

Most fermentation derives from
the slowly degradable fraction.

Understanding Goodness-of-Fit Metrics

Model fit should never be evaluated using a single statistic.

R-Squared

\[ R^2 \]

Measures the proportion of observed variation explained by the model.

Example:

R² = 0.99

Interpretation:

99% of variation is explained by
the fitted model.

RMSE

Root Mean Squared Error:

\[ RMSE \]

Measures average prediction error.

Example:

RMSE = 1.5 mL

Interpretation:

Predictions differ from observations
by approximately 1.5 mL on average.

Smaller values are preferred.

RSS

Residual Sum of Squares:

\[ RSS \]

Represents total unexplained variation.

Smaller values indicate better fit.

AIC

Akaike Information Criterion:

\[ AIC \]

Balances:

Fit quality
+
Model complexity

Smaller values are preferred.

BIC

Bayesian Information Criterion:

\[ BIC \]

Similar to AIC but applies a stronger penalty for additional parameters.

Smaller values are preferred.

Why Higher R² Does Not Always Mean a Better Model

Consider:

Model Parameters R² AIC
Groot 3 0.9992 33
Richards 4 0.9994 35

The Richards model explains slightly more variation.

However:

Additional complexity

may not justify:

Minimal improvement

AIC correctly penalizes the extra parameter.

Therefore:

Higher R² alone should not determine
model selection.

Model Selection Strategy

Recommended workflow:

1. Fit multiple models

2. Evaluate convergence

3. Compare RMSE

4. Compare AIC and BIC

5. Examine residual plots

6. Consider biological interpretation

7. Select the most appropriate model

Interpreting Failed Fits

Common reasons include:

Poor starting values

Too many parameters

Insufficient observations

Parameter redundancy

Inappropriate model structure

When convergence problems occur:

Biological Reality Matters

The statistically best model is not always the biologically most meaningful model.

Researchers should consider:

alongside fit statistics.

Practical Recommendations

Use Simple Models When

Examples:

Use Lag Models When

Examples:

Use Flexible Sigmoidal Models When

Examples:

Use Dual-Pool Models When

Example:

Summary

A successful analysis combines:

Researchers are encouraged to fit multiple models and evaluate both statistical and biological performance before selecting a final model.