tor_fit: Mixture model aimed at assigning metabolic rate measurements...

In vullioud/toRpoR: Assigning metabolic rate measurements to torpor and euthermia in heterothermic endotherms

Mixture model aimed at assigning metabolic rate measurements (M) to torpor and euthermia

Description

The function tor_fit() considers the relation between metabolic rate (M) and ambient temperature (Ta) assumed by the Scholander-Irving model and its later extensions.

Usage

tor_fit(
  Ta,
  M,
  Mtnz = NULL,
  Tlc = NULL,
  fitting_options = list(ni = 50000, nt = 10, nb = 30000, nc = 3, parallel = TRUE),
  confidence = 0.5
)

Arguments

Ta	A vector of ambient temperature.
M	A vector of methabolic measure.
Mtnz	Mtnz value if not estimated
Tlc	Tlc value if not estimated
fitting_options	a list of fitting option to pass to jags.
confidence	Assignment confidence threshold

Details

Resting M measured within the thermoneutral zone (TNZ) is independent of Ta. This rate is hereafter referred to as Mtnz, although it would correspond to the basal metabolic rate (BMR) provided that the specific criteria for the BMR are met (see Fasel et al. in prep.). Below the lower critical temperature of TNZ (Tlc), M of euthermic animals increases linearly with decreasing Ta. M of torpid animals increases linearly with decreasing Ta to maintain a minimal body temperature below some threshold ambient temperature (Tt). This state is usually referred to as "regulated torpor". Between Tt and Tlc, M of torpid animals follows an exponential curve. In this Ta range, torpor is referred to as "conforming torpor".

Value

A list of class tor_obj

Examples

## Not run: 
test_mod <- tor_fit(Ta = test_data2$Ta, M = test_data2$VO2ms,
                   fitting_options = list(parallel = TRUE, ni = 500, nt = 2,
                   nb = 200), confidence = 0.9)
test_mod2 <- tor_fit(Ta = test_data2$Ta, M = test_data2$VO2ms, Mtnz = 1.8, Tlc = 29.2,
                   fitting_options = list(parallel = TRUE))
test_mod3 <- tor_fit(Ta = test_data2$Ta, M = test_data2$VO2ms, Mtnz = 1.8,
                   fitting_options = list(parallel = TRUE))

test_mod4 <- tor_fit(Ta = test_data2$Ta, M = test_data2$VO2ms, Tlc = 29.2,
                   fitting_options = list(parallel = TRUE))

## End(Not run)

vullioud/toRpoR documentation built on Nov. 30, 2024, 11:13 a.m.