This is probably also the reason why I get an error with dimension mismatch using the package, as in the example max_t is set to 22 but china_cases.R and world_cases.R still contain 24 days – thus assertthat::assert_that(all(dim(y_prime) == dim(exp_incr_prime))) fails. While if I put 24 or reduce the dataset in the package it works.
Also, just to flag, in the example on github you have
However p_hat is a list containings also $value and $counts beside $par, so I believe this should be p_hat$par – as likely also in the bootstrap function, but I didn’t check that
From: Luca Brugnolini Luca.Brugnolini@rokoscapital.com Sent: 28 January 2020 13:27 To: Read, Jonathan jonathan.read@lancaster.ac.uk Subject: RE: [External] Paper data Novel coronavirus 2019-nCov
Thanks a lot for the answer,
This is probably also the reason why I get an error with dimension mismatch using the package, as in the example max_t is set to 22 but china_cases.R and world_cases.R still contain 24 days – thus assertthat::assert_that(all(dim(y_prime) == dim(exp_incr_prime))) fails. While if I put 24 or reduce the dataset in the package it works.
Also, just to flag, in the example on github you have
However p_hat is a list containings also
$value
and$counts
beside$par
, so I believe this should bep_hat$par
– as likely also in the bootstrap function, but I didn’t check thatHope this may help,
Best
Luca