User-provided, replication-wise confidence intervals for coverage

[ellessenne]

Idea from the simulation course, t-test example: passing the t-based confidence intervals as columns in data:

simsum(data, ..., ci.limits = c("column.lower", "column.upper"), ...)

This is compatible with the current functionality, in which case the column should have the same values for each replication.