sim_1tm: Simulate from the one thermal mass model (1TM)

In paulnorthrop/walls: Statistical Inference for the Thermophysical Properties of Walls from In-Situ Surveys

Description

Simulates from the one thermal mass model (1TM)

Usage

sim_1tm(data, start = NULL, end = NULL, params = NULL, tau = 1/12,
  sigma_in = 1, sigma_out = 1, n = NULL, tc0 = NULL)

Arguments

data	A dataset in the format of either cwall_east or cwall_north.
tau	A numeric scalar. Time interval between successive measurements, in hours.
sigma_in	A numeric scalar. Error standard deviation inside.
n	An integer scalar.
sigma_in	A numeric scalar. Error standard deviation outside.
t0	A numeric scalar. Initial temperature of the thermal mass C0. If t0 is not supplied then the mean of the earliest

See Also

cwall_east, cwall_north.

Examples

set.seed(23072018)
fit_n <- one_tm(data = cwall_north, flux = "both", hetero = FALSE)
params_n <- coef(fit_n)
params_n2 <- coef(fit_n2)
sigma(fit_n)

# Estimates from Table 5.13 on page 173 of Gori (2017)
# Based on 3 days of data
pars <- c(R1 = 1.297, R2 = 0.069, C1 = 1.02e5 / 3600)
sim_n <- sim_1tm(data = cwall_north, params = pars, tc0 = 14.65,
                 start = "07/05/2015 14:00", end = "10/05/2015 13:55")
sim_n <- sim_1tm(data = cwall_north, params = pars, tc0 = 14.65,
                 start = "07/05/2015 14:00", end = "07/06/2015 13:55")
fit2 <- one_tm(data = cwall_north, flux = "both", sim_data = sim_n, hetero = FALSE)

pjn <- sim_1tm(data = cwall_north, params = params_n, n = 1400)
pjn2 <- sim_1tm(data = cwall_north, params = params_n2, n = 1400)

fit_n <- one_tm(data = cwall_north, flux = "both", hetero = FALSE)
params_n <- coef(fit_n)
sim_north <- sim_1tm(data = cwall_north, params = params_n, flux = "both")
sim_north <- sim_1tm(data = cwall_north, flux = "both")
fit_ns <- one_tm(data = cwall_north, flux = "both", sim_data = sim_north, hetero = FALSE)

fit_e <- one_tm(data = cwall_east, flux = "both")
params_e <- coef(fit_e)
sim_east <- sim_1tm(data = cwall_east, params = params_e, flux = "both")
fit_es <- one_tm(data = cwall_east, flux = "out", sim_data = sim_east)

paulnorthrop/walls documentation built on May 15, 2019, 10:02 p.m.