Get a formatted data.frame of inputs from which metabolism can be modeled. These test data were provided by Bob Hall.
Arguments
- num_days
the number of days to include in the data. character format because only certain numbers of days are permitted (see defaults in Usage for the accepted options)
- res
character specifying the desired resolution of the data in minutes (character; see defaults in Usage for the accepted options)
- flaws
character specifying one or more flaws to include in the data, or empty (
c()) for no flaws. default is no flaws.- day_start
start time (inclusive) of a day's data in number of hours from the midnight that begins the date. For example, day_start=-1.5 indicates that data describing 2006-06-26 begin at 2006-06-25 22:30, or at the first observation time that occurs after that time if day_start doesn't fall exactly on an observation time. For metabolism models working with single days of input data, it is conventional/useful to begin the day the evening before, e.g., -1.5, and to end just before the next sunrise, e.g., 30. For multiple consecutive days, it may make the most sense to start just before sunrise (e.g., 4) and to end 24 hours later. For nighttime regression, the date assigned to a chunk of data should be the date whose evening contains the data. The default is therefore 12 to 36 for metab_night, of which the times of darkness will be used.
- day_end
end time (exclusive) of a day's data in number of hours from the midnight that begins the date. For example, day_end=30 indicates that data describing 2006-06-26 end at the last observation time that occurs before 2006-06-27 06:00. See day_start for recommended start and end times.
- attach.units
(deprecated, effectively FALSE in future) logical, default TRUE for backward compatibility. Should units be attached to the data.frame?
Examples
head(data_metab())
#> solar.time DO.obs DO.sat depth temp.water light
#> 5688 2012-09-18 04:00:58 8.39 9.076040 0.16 3.63 0
#> 5689 2012-09-18 04:05:58 8.41 9.083329 0.16 3.60 0
#> 5690 2012-09-18 04:10:58 8.40 9.088194 0.16 3.58 0
#> 5691 2012-09-18 04:15:58 8.41 9.090628 0.16 3.57 0
#> 5692 2012-09-18 04:20:58 8.42 9.093063 0.16 3.56 0
#> 5693 2012-09-18 04:25:58 8.42 9.097936 0.16 3.54 0
head(data_metab(res='30'))
#> solar.time DO.obs DO.sat depth temp.water light
#> 5689 2012-09-18 04:05:58 8.41 9.083329 0.16 3.60 0.0000
#> 5695 2012-09-18 04:35:58 8.42 9.105254 0.16 3.51 0.0000
#> 5701 2012-09-18 05:05:58 8.45 9.127267 0.16 3.42 0.0000
#> 5707 2012-09-18 05:35:58 8.47 9.149369 0.16 3.33 0.0000
#> 5713 2012-09-18 06:05:58 8.55 9.166621 0.16 3.26 134.3296
#> 5719 2012-09-18 06:35:58 8.71 9.174032 0.16 3.23 361.1082