lognorm_rp {marp} | R Documentation |
A function to fit Log-Normal renewal model
lognorm_rp(data, t, y)
data |
as input inter-event times |
t |
as user-specified time intervals (used to compute hazard rate) |
y |
as user-specified time point (used to compute time-to-event probability) |
returns list of estimates after fitting Log-Normal renewal model
Estimated mean (on the log scale) of the Log-Normal model
Estimated standard deviation (on the log scale)of the Log-Normal model
Negative log-likelihood
Akaike information criterion (AIC)
Bayesian information criterion (BIC)
Estimated mean
Estimated (logit) probabilities
Estimated (log) hazard rates
set.seed(42)
data <- rgamma(100,3,0.01)
# set some parameters
t = seq(100, 200, by=10) # time intervals
y = 304 # cut-off year for estimating probablity
# fit Log-Normal renewal model
result <- marp::lognorm_rp(data, t, y)
# print result
cat("par1 = ", result$par1, "\n")
cat("par2 = ", result$par2, "\n")
cat("logL = ", result$logL, "\n")
cat("AIC = ", result$AIC, "\n")
cat("BIC = ", result$BIC, "\n")
cat("mu_hat = ", result$mu_hat, "\n")
cat("pr_hat = ", result$pr_hat, "\n")