function(level = 0, slope = 1) { #brass.general.mltf # logit.lx <- data.matrix(read.table("d:\\aidsmodel\\data\\brass.general.logits.txt"))[, 2] rows <- 1:100 cols <- c("Age", "logit.lx", "lx", "Lx", "Tx", "ex") logit.table <- matrix(data = 0, nrow = 100, ncol = 6, dimnames = list(rows, cols)) logit.table[, 1] <- 0:99 logit.table[, 2] <- logit.lx lx <- 1/(1 + exp(2 * (level + slope * logit.lx))) lx[1] <- 1 logit.table[, 3] <- lx Lx <- 0.5 * (lx + c(lx[2:length(lx)], 0)) Lx[1] <- 0.3 + 0.7 * lx[2] logit.table[, 4] <- Lx Tx <- rev(cumsum(rev(Lx))) logit.table[, 5] <- Tx ex <- Tx/lx logit.table[, 6] <- ex write.ascii(brass.general.mltf) logit.table }