MatrixArrange: Arrange Vectorised Matrix Samples using Rcpp

In seabbs/SpeedyMarkov: Speed up Discrete Markov Model Cost Effectiveness Simulations

Description

A convenience function used to arrange vectorised matrix samples into the correct matrix format for several functions used to specify Markov model. See example_two_state_markov for an example use case. Implemented using Rcpp.

Usage

MatrixArrange(samples)

Arguments

samples

A list of vectorised matrix samples

Value

A list of matrices with each matrix representing a single sample.

Examples

matrix_samples <- list(VGAM::rdiric(1:5, c(88, 12)),
                      VGAM::rdiric(1:5, c(8, 92)))

# R implementation
samples_r <- matrix_arrange_inner(matrix_samples)


# Rcpp implementation
samples_rcpp <- MatrixArrange(matrix_samples)

all.equal(samples_r, samples_rcpp)

# Benchmark
library(microbenchmark)
microbenchmark(matrix_arrange_inner(matrix_samples), 
               MatrixArrange(matrix_samples), 
               times = 1000)

seabbs/SpeedyMarkov documentation built on Dec. 26, 2019, 4:38 a.m.