Add contractGrid and contractGridPremium functions to compare tariffs

DESCRIPTION

@@ -23,6 +23,7 @@ Collate:
     'ProfitParticipation.R'
     'InsuranceTarif.R'
     'InsuranceContract.R'
+    'contractGrid.R'
     'exportInsuranceContract_xlsx.R'
     'showVmGlgExamples.R'
     'exportInsuranceContractExample.R'

NAMESPACE

@@ -8,6 +8,8 @@ export(InsuranceContract.ParametersFill)
 export(InsuranceContract.Values)
 export(InsuranceTarif)
 export(ProfitParticipation)
+export(contractGrid)
+export(contractGridPremium)
 export(deathBenefit.annuityDecreasing)
 export(deathBenefit.linearDecreasing)
 export(exportInsuranceContract.xlsx)

R/contractGrid.R

+#' @export
+contractGrid = function(axes = list(ages = seq(20, 60, 10), policyPeriod = seq(5, 35, 5)), ...) {
+    # Create all combinations of the variables given for the axes:
+    gridByRow = expand.grid(axes);
+    # Apply InsuranceContract$new to each combination (and add the additional arguments)
+    vals = apply(gridByRow, 1, function(axisVals) {
+        do.call(
+            InsuranceContract$new,
+            c(as.list(axisVals), ...))
+    })
+    dimnames = makeContractGridDimnames(axes)
+    matrix(vals, nrow = length(axes[[1]]), ncol = length(axes[[2]]), byrow = FALSE, dimnames = dimnames)
+}
+
+makeContractGridDimname = function(value) { UseMethod("makeContractGridDimname", value) }
+makeContractGridDimname.InsuranceTarif = function(tarif) { tarif$name }
+makeContractGridDimname.mortalityTable = function(table) { table@name }
+makeContractGridDimname.default = function(value) { value }
+makeContractGridDimnames = function(axes) {
+    lapply(axes, function(axis) { lapply(axis, makeContractGridDimname); } )
+}
+
+
+#' @export
+contractGridPremium = function(contractGrid = NULL, premium="written", ...) {
+    if (missing(contractGrid) || is.null(contractGrid)) {
+        contractGrid = contractGrid(...)
+    }
+    apply(contractGrid, 1:2, function(c) { c[[1]]$Values$premiums[[premium]] })
+}
+
+if (FALSE) {
+    gg = contractGrid(tarif = Generali.U188, axes = list(age = seq(20, 60, 10), policyPeriod = seq(5, 35, 5)), sumInsured = 100000)
+
+    contractGridPremium(gg)
+    apply(gg, 1:2, function(c) { c[[1]]$Parameters$ContractData$sumInsured })
+
+    contractGridPremium(tarif = Generali.U188, axes = list(age = seq(20, 60, 10), policyPeriod = seq(5, 35, 5)), sumInsured = 100000)
+
+    contractGridPremium(age = 30, axes = list(tarif = c(Generali.U188, Generali.U180, Generali.U3.1_2015), policyPeriod = seq(5, 35, 5)), sumInsured = 100000)
+}
+
+