R/make_constants.R

BEERS_INT_MOD_FIRST <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  1.0000,  0.0000,  0.0000,  0.0000,  0.0000,  0.0000,
  0.6668,  0.5270, -0.2640,  0.0820, -0.0140,  0.0022,
  0.4099,  0.8592, -0.3598,  0.1052, -0.0173,  0.0028,
  0.2196,  1.0279, -0.3236,  0.0874, -0.0136,  0.0023,
  0.0862,  1.0644, -0.1916,  0.0464, -0.0066,  0.0012))

BEERS_INT_MOD_SECOND <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  1.0000,  0.0000,  0.0000,  0.0000,  0.0000,
 -0.0486,  0.8655,  0.2160, -0.0350,  0.0030, -0.0009,
 -0.0689,  0.6903,  0.4238, -0.0442,  0.0003, -0.0013,
 -0.0697,  0.5018,  0.5938, -0.0152, -0.0097, -0.0010,
 -0.0589,  0.3233,  0.7038,  0.0578, -0.0257, -0.0003))

BEERS_INT_MOD_MID <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
 -0.0430,  0.1720,  0.7420,  0.1720, -0.0430,  0.0000,
 -0.0270,  0.0587,  0.7072,  0.3162, -0.0538, -0.0013,
 -0.0141, -0.0132,  0.6098,  0.4708, -0.0477, -0.0056,
 -0.0056, -0.0477,  0.4708,  0.6098, -0.0132, -0.0141,
 -0.0013, -0.0538,  0.3162,  0.7072,  0.0587, -0.0270))

BEERS_INT_MOD_PENULT <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000, -0.0430,  0.1720,  0.7420,  0.1720, -0.0430,
 -0.0003, -0.0257,  0.0578,  0.7038,  0.3233, -0.0589,
 -0.0010, -0.0097, -0.0152,  0.5938,  0.5018, -0.0697,
 -0.0013,  0.0003, -0.0442,  0.4238,  0.6903, -0.0689,
 -0.0009,  0.0030, -0.0350,  0.2160,  0.8655, -0.0486))

BEERS_INT_MOD_FINAL <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  0.0000,  0.0000,  0.0000,  1.0000,  0.0000,
  0.0012, -0.0066,  0.0464, -0.1916,  1.0644,  0.0862,
  0.0023, -0.0136,  0.0874, -0.3236,  1.0279,  0.2196,
  0.0028, -0.0173,  0.1052, -0.3598,  0.8592,  0.4099,
  0.0022, -0.0140,  0.0820, -0.2640,  0.5270,  0.6668))

BEERS_INT_ORD_FIRST <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  1.0000,  0.0000,  0.0000,  0.0000,  0.0000,  0.0000,
  0.6667,  0.4969, -0.1426, -0.1006,  0.1079, -0.0283,
  0.4072,  0.8344, -0.2336, -0.0976,  0.1224, -0.0328,
  0.2148,  1.0204, -0.2456, -0.0536,  0.0884, -0.0244,
  0.0819,  1.0689, -0.1666, -0.0126,  0.0399, -0.0115))

BEERS_INT_ORD_SECOND <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  1.0000,  0.0000,  0.0000,  0.0000,  0.0000,
 -0.0404,  0.8404,  0.2344, -0.0216, -0.0196,  0.0068,
 -0.0497,  0.6229,  0.5014, -0.0646, -0.0181,  0.0081,
 -0.0389,  0.3849,  0.7534, -0.1006, -0.0041,  0.0053,
 -0.0191,  0.1659,  0.9354, -0.0906,  0.0069,  0.0015))

BEERS_INT_ORD_MID <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  0.0000,  1.0000,  0.0000,  0.0000,  0.0000,
  0.0117, -0.0921,  0.9234,  0.1854, -0.0311,  0.0027,
  0.0137, -0.1101,  0.7194,  0.4454, -0.0771,  0.0087,
  0.0087, -0.0771,  0.4454,  0.7194, -0.1101,  0.0137,
  0.0027, -0.0311,  0.1854,  0.9234, -0.0921,  0.0117))

BEERS_INT_ORD_PENULT <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  0.0000,  0.0000,  1.0000,  0.0000,  0.0000,
  0.0015,  0.0069, -0.0906,  0.9354,  0.1659, -0.0191,
  0.0053, -0.0041, -0.1006,  0.7534,  0.3849, -0.0389,
  0.0081, -0.0181, -0.0646,  0.5014,  0.6229, -0.0497,
  0.0068, -0.0196, -0.0216,  0.2344,  0.8404, -0.0404))

BEERS_INT_ORD_FINAL <- matrix(nrow = 5, ncol = 6, byrow = TRUE, data = c(
  0.0000,  0.0000,  0.0000,  0.0000,  1.0000,  0.0000,
 -0.0115,  0.0399, -0.0126, -0.1666,  1.0689,  0.0819,
 -0.0244,  0.0884, -0.0536, -0.2456,  1.0204,  0.2148,
 -0.0328,  0.1224, -0.0976, -0.2336,  0.8344,  0.4072,
 -0.0283,  0.1079, -0.1006, -0.1426,  0.4969,  0.6667))

BEERS_SUB_MOD_FIRST <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.3332, -0.1938, +0.0702, -0.0118, +0.0022,
 +0.2569, -0.0753, +0.0205, -0.0027, +0.0006,
 +0.1903, +0.0216, -0.0146, +0.0032, -0.0005,
 +0.1334, +0.0969, -0.0351, +0.0059, -0.0011,
 +0.0862, +0.1506, -0.0410, +0.0054, -0.0012))

BEERS_SUB_MOD_SECOND <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.0486, +0.1831, -0.0329, +0.0021, -0.0009,
 +0.0203, +0.1955, -0.0123, -0.0031, -0.0004,
 +0.0008, +0.1893, +0.0193, -0.0097, +0.0003,
 -0.0108, +0.1677, +0.0577, -0.0153, +0.0007,
 -0.0159, +0.1354, +0.0972, -0.0170, +0.0003))

BEERS_SUB_MOD_MID <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 -0.0160, +0.0973, +0.1321, -0.0121, -0.0013,
 -0.0129, +0.0590, +0.1564, +0.0018, -0.0043,
 -0.0085, +0.0260, +0.1650, +0.0260, -0.0085,
 -0.0043, +0.0018, +0.1564, +0.0590, -0.0129,
 -0.0013, -0.0121, +0.1321, +0.0973, -0.0160))

BEERS_SUB_MOD_PENULT <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.0003, -0.0170, +0.0972, +0.1354, -0.0159,
 +0.0007, -0.0153, +0.0577, +0.1677, -0.0108,
 +0.0003, -0.0097, +0.0193, +0.1893, +0.0008,
 -0.0004, -0.0031, -0.0123, +0.1955, +0.0203,
 -0.0009, +0.0021, -0.0329, +0.1831, +0.0486))

BEERS_SUB_MOD_FINAL <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 -0.0012, +0.0054, -0.0410, +0.1506, +0.0862,
 -0.0011, +0.0059, -0.0351, +0.0969, +0.1334,
 -0.0005, +0.0032, -0.0146, +0.0216, +0.1903,
 +0.0006, -0.0027, +0.0205, -0.0753, +0.2569,
 +0.0022, -0.0118, +0.0702, -0.1938, +0.3332))

BEERS_SUB_ORD_FIRST <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.3333, -0.1636, -0.0210, +0.0796, -0.0283,
 +0.2595, -0.0780, +0.0130, +0.0100, -0.0045,
 +0.1924, +0.0064, +0.0184, -0.0256, +0.0084,
 +0.1329, +0.0844, +0.0054, -0.0356, +0.0129,
 +0.0819, +0.1508, -0.0158, -0.0284, +0.0115))

BEERS_SUB_ORD_SECOND <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.0404, +0.2000, -0.0344, -0.0128, +0.0068,
 +0.0093, +0.2268, -0.0402, +0.0028, +0.0013,
 -0.0108, +0.2272, -0.0248, +0.0112, -0.0028,
 -0.0198, +0.1992, +0.0172, +0.0072, -0.0038,
 -0.0191, +0.1468, +0.0822, -0.0084, -0.0015))

BEERS_SUB_ORD_MID <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 -0.0117, +0.0804, +0.1570, -0.0284, +0.0027,
 -0.0020, +0.0160, +0.2200, -0.0400, +0.0060,
 +0.0050, -0.0280, +0.2460, -0.0280, +0.0050,
 +0.0060, -0.0400, +0.2200, +0.0160, -0.0020,
 +0.0027, -0.0284, +0.1570, +0.0804, -0.0117))

BEERS_SUB_ORD_PENULT <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 -0.0015, -0.0084, +0.0822, +0.1468, -0.0191,
 -0.0038, +0.0072, +0.0172, +0.1992, -0.0198,
 -0.0028, +0.0112, -0.0248, +0.2272, -0.0108,
 +0.0013, +0.0028, -0.0402, +0.2268, +0.0093,
 +0.0068, -0.0128, -0.0344, +0.2000, +0.0404))

BEERS_SUB_ORD_FINAL <- matrix(nrow = 5, ncol = 5, byrow = TRUE, data = c(
 +0.0115, -0.0284, -0.0158, +0.1508, +0.0819,
 +0.0129, -0.0356, +0.0054, +0.0844, +0.1329,
 +0.0084, -0.0256, +0.0184, +0.0064, +0.1924,
 -0.0045, +0.0100, +0.0130, -0.0780, +0.2595,
 -0.0283, +0.0796, -0.0210, -0.1636, +0.3333))

save(BEERS_INT_MOD_FIRST,  BEERS_INT_MOD_SECOND, BEERS_INT_MOD_MID, 
     BEERS_INT_MOD_PENULT, BEERS_INT_MOD_FINAL,
     BEERS_INT_ORD_FIRST,  BEERS_INT_ORD_SECOND, BEERS_INT_ORD_MID,
     BEERS_INT_ORD_PENULT, BEERS_INT_ORD_FINAL, 
     BEERS_SUB_MOD_FIRST,  BEERS_SUB_MOD_SECOND, BEERS_SUB_MOD_MID, 
     BEERS_SUB_MOD_PENULT, BEERS_SUB_MOD_FINAL,
     BEERS_SUB_ORD_FIRST,  BEERS_SUB_ORD_SECOND, BEERS_SUB_ORD_MID,
     BEERS_SUB_ORD_PENULT, BEERS_SUB_ORD_FINAL, 
  file = "R/sysdata.rda")
mrc-ide/beers documentation built on Dec. 10, 2020, 12:10 a.m.