For most organisms, metabolic rate does not scale linearly, but rather according to a power function: *MO2 = b0 * M^b*. This function estimates the scaling coefficient, `b`

, given MO2s from different sized individuals.

1 |

`mass` |
a vector of animal masses. |

`MO2` |
a vector of metabolic rates. |

`plot` |
a string defining what kind of plot to display. "linear" for linear axes, "log" for log10-scale axes, and "none" for no plot. Default is "linear". |

*MO2 = b0 * M^b*

where `b0`

is species-specific normalization constant, `M`

is mass and `b`

is the scaling coefficient.

Matthew A. Birk, matthewabirk@gmail.com

1 2 3 4 5 6 7 8 9 | ```
# Simple example
mass <- c(1, 10, 100, 1000, 40, 4, 400, 60, 2, 742, 266, 983) # made up values
MO2 <- mass ^ 0.65 + rnorm(n = length(mass)) # make up some data
calc_b(mass = mass, MO2 = MO2)
# How about some mass-specific MO2s?
msMO2 <- mass ^ -0.25 + rnorm(n = length(mass), sd = 0.05)
calc_b(mass = mass, MO2 = msMO2)
calc_b(mass = mass, MO2 = msMO2, plot = "log")
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.