fitACi: fitACi

View source: R/fitAci.R

fitACiR Documentation

fitACi

Description

Fit A/Ci curves using forced values, and bounds

Usage

fitACi(
  data,
  name_assimilation = "A",
  name_ci = c("pCi", "Ci"),
  gammastar = 3.52,
  O2 = 21,
  pressure = 101,
  tleaf = 25,
  initialGuess = NA,
  forceValues = c(NA, NA, NA, NA, NA, NA, NA),
  bound_l = c(1, 1, 1, 0.001, 0.001, 0, 0),
  bound_h = c(1000, 1000, 1000, 30, 30, 1, 0.75),
  ignoreTPU = F,
  maxiter = 250,
  Kc = exp(35.9774 - (80.99/(0.008314 * (273.15 + tleaf)))),
  Ko = exp(12.3772 - (23.72/(0.008314 * (273.15 + tleaf))))
)

Arguments

data

The A/Ci data to be used

name_assimilation

The name given to the assimilation column in "data"

name_ci

Name given to the Ci column in "data"

gammastar

The Cc compensation point. Default: Tobacco.

O2

The oxygen concentration in parts per hundred.

pressure

Atmospheric pressure in kPa

tleaf

The leaf temperature, in celsius

initialGuess

Initial guesses for fitting. Will generate ok guesses automatically if not provided

forceValues

A vector of forced values in order: VcMax, J, TPU, gm, Rd, aG, aS. Include NAs where not forced, if you modify this!

bound_l

The minimum values to be considered for each parameter

bound_h

the maximum value to be considered for each parameter

ignoreTPU

Whether to fit TPU or not. Leave false if you don't know what you're doing!

maxiter

How many iterations to allow for the fitting algorithm. If you have a lot of points, large values for maxiter will take a long time!

Kc

Michaelis-menten kinetic parameter for carboxylation

Ko

Michaelis-menten kinetic parameter for oxygenation


poales/msuRACiFit documentation built on Jan. 4, 2024, 3:57 p.m.