“…# Reading data (Wang et al, 2015) t <-c (3,7,11,18,22,25,28,32,34,35,35,36,40,40,41,54,66,76,84,88,92) d <-c(1,1,0,1,0,1,1,0,0,1,0,0,0,0,1,0,0,0,0,0,0) n <-length(t) # the sample size K <-2 # number of parameters # Loading the maxLik package library(maxLik) # The likelihood function log.f <-function(parms…”