LightFitR (version 1.0.0)

makeRegime: Create a regime (matrix) to program lights to achieve intended irradiances

Description

This is a wrapper function that carries out multiple steps:

1. Calculate closest intensities

2. Predict the intensities to use to achieve the target irradiance (via a system of linear equations or non-negative least squares)

3. Tidy the intensities (rounding to integer, keep within the range of intensities that the lights can be set to)

4. Format the intensities and timestamps into a human-readable regime matrix

Usage

makeRegime(
  timeVector_POSIXct,
  irradiance_matrix,
  calibration_leds,
  calibration_wavelengths,
  calibration_intensities,
  calibration_irradiances,
  peaks = LightFitR::helio.dyna.leds$wavelength,
  method = "nnls"
)

Value

Matrix with light regime needed to program the lights

Arguments

timeVector_POSIXct

Vector of schedule timepoints in POSICxt format

irradiance_matrix

Matrix of intended irradiances. rows = leds and columns = events

calibration_leds

A numeric vector of LED values from calibration, mapping to intensities and irradiances (i.e. the same length)

calibration_wavelengths

A numeric vector of wavelengths from calibration, corresponding to intensities and irradiances

calibration_intensities

A numeric vector of intensities (heliospectra units) from calibration

calibration_irradiances

A numeric vector of measured irradiances (any units, as long as it is consistently used) from calibration

peaks

Vector of length 8 or 9. Containing wavelengths at which each LED peaks.

method

Use 'nnls' (non-negative least squares) or 'sle' (system of linear equations)

NNLS vs SLE

NNLS and SLE largely predict the same intensities, except in outlier cases. The default is NNLS, but if your predicted intensities end up being very far off, try SLE.

Examples

# Prep variables
calib <- LightFitR::calibration
times <- LightFitR::time_vector
target_irradiance <- LightFitR::target_irradiance

# Run function
makeRegime(times, target_irradiance, calib$led, calib$wavelength, calib$intensity, calib$irradiance)