lrstat 0.3.2

• removed the pivoting step for pmvnormr and qmvnormr to take advantage of independent increment structure of score process
• removed the fast parameter for pmvnormr and qmvnormr so that the default is to use fast approximation
• updated the point estimate for getRCI
• renamed tsssd to seamless in the function names for two-stage seamless sequential design for treatment/dose selection
• updated the critical values when typeAlphaSpending = “none” in getDesign_seamless
• added integratedTrial to adaptDesign
• updated the backward image confidence interval method to account for the case when astar > conditional_power
• updated point estimate for repeated confidence interval
• added functions for multi-arm multi-stage (MAMS) design and anlysis, including exitprob_mams, getBound_mams, getDesign_mams, adaptDesign_mams, getCI_mams, getADCI_mams
• ensured that if typeAlphaSpending is “OF”, “P”, “WT”, or “none”, then informationRates, efficacyStopping, and spendingTime must be of full length kMax, and informationRates and spendingTime must end with 1 for getBound, getBound_seamless, and getBound_mams
• ensured that if typeAlphaSpending is “OF”, “P”, “WT”, or “none”, then informationRates, efficacyStopping, and spendingTime must be of full length kMax, and informationRates and spendingTime must end with 1 for getCI, getRCI
• ensured that if typeAlphaSpendingNew is “OF”, “P”, “WT”, or “none”, then informationRatesNew, efficacyStoppingNew, and spendingTimeNew must be of full length kNew, and informationRatesNew and spendingTimeNew must end with 1 for getADCI, getADRCI
• moved up the position of parameter MullerSchafer in getADCI, getADRCI
• added getCP_seamless for conditional power calculation for two-stage seamless sequential design for treatment/dose selection
• added getCP_mams for conditional power calculation for multi-arm multi-stage design
• added adaptDesign_seamless for power and sample size calculation of adaptive two-stage seamless sequential design for treatment/dose selection
• added getCI_seamless for confidence interval calculation for two-stage seamless sequential design for treatment/dose selection
• added getADCI_seamless for confidence interval calculation using the backward image method for adaptive two-stage seamless sequential design for treatment/dose selection
• moved the f_pvalue function from confidence_interval.cpp to generic_design.cpp
• updated infoRatesNew when informationRatesNew is not missing in adaptDesign_seamless_cpp and adaptDesign_mams_cpp
• added default value of 1 for parameter r of adaptDesign_seamless
• updated pmvnormmccpp to use s and T instead of Ivec
• added byLevelBounds in the output of getDesign_mams
• updated the condition for cpu0 calculation in adaptDesign_mams_cpp
• updated the description for adaptDesign_mams output
• updated the default value of r and rNew to 1 in adaptDesign_mams and getADCI_mams
• updated stage‑1 covariance used for bound calculation in getBound_seamless to account for corr_known when using alpha‑spending functions
• used memoization to avoid recalculating p0 and p1 in the lambda h for preject_by_arm (exitprob_seamless)
• added stopStage to output data sets and added bestArm and reject to sumdata2 output data set for lrsim_seamless
• used NULL as the default value for lambdas and gammas for lrsim_seamless
• added lrsim_mams for simulation of multi-arm multi-stage design
• added futilityBounds, futilityCP, futilityTheta, futilityBoundsInt, futilityCPInt, futilityThetaInt to adaptDesign, adaptDesign_seamless, and adaptDesign_mams
• removed typeBetaSpending and parameterBetaSpending from adaptDesign, adaptDesign_seamless, and adaptDesign_mams

lrstat 0.3.1

• calculated eventsPerStage, dropoutPerStage, subjectsPerStage, and analysisTimePerstage in the lrsim output based on trials having the stage (i.e., did not stop before the stage).
• added mvnormr for multivariate normal probability and quantile calculations
• added tsssd_design: analytical calculation for two-stage seamless sequential designs for treatment/dose selection.
• added lrstim_tsssd: simulation function for two-stage seamless sequential designs for treatment/dose selection.
• printed power with 4 instead of 3 decimal places in the output
• used std::size_t instead of size_t in all header files
• added pnorm_fast, qnorm_acklam, expand_stratified_to_slice, and invchol utility functions

lrstat 0.3.0

• replaced accrualDuration with n for lrsim, lrsim2e, lrsim3a, lrsim2e3a, lrsimsub, and binary_tte_sim to allow the simulation to run with a fixed number of subjects
• added time2 and weight parameters to kmest
• added time2, weight, and weight_readj parameters to lrtest
• added assess_phregr and zph_phregr functions
• updated findInterval3 to add rightmost_closed, all_inside, and left_open parameters
• added dtpwexp for density function of truncated piecewise exponential distribution
• updated hazard_pd to match probabilities at quantiles in addition to cut points
• added corr_pfs_os to obtain correlation between PFS and OS given correlation between PD and OS
• added hazard_pd2 to obtain the hazards for PD to match distribution for PFS and correlation between PFS and OS
• renamed rho to rho_pd_os for added clarity for lrsim2e and lrsim2e3a
• updated hazard_sub to match probabilities at quantiles in addition to cut points
• moved logparams to logistic_regression.cpp and removed logparams, f_der_0, f_ressco_0, logisregloop, and logisregplloop from logistic_regression.h
• updated lrsim2e and lrsim2e3a to account for additional quantile-based cut points for pd
• updated lrsimsub to account for additional quantile-based cut points for biomarker negative sub population
• moved aftparams and coxparams to survival_analysis.cpp and removed aftparams, f_der_1, f_ressco_1, f_jj_1, liferegloop, liferegplloop, f_ld_1, coxparams, f_der_2, phregloop, phregplloop, f_basehaz, f_ressco_2, f_resmart, f_ressch, and f_der_i_2 from survival_analysis.h
• added ptpwexpcpp1 and updated ptpwexpcpp accordingly
• added mtpwexp for the mean and variance of a truncated piecewise exponential distribution
• added quad2d for numerical integration in two dimensions
• removed hd, pd, ad, house, row_house, col_house, givens, row_rot, col_rot, house_bidiag, zero_diagonal, and svd_step from utilities.h
• renamed the parameters in lrsim2e and lrsim2e3a, replacing e1 by pfs and e2 by os
• run the simulations in serial mode to pass CRAN checks

lrstat 0.2.15

• added na.action = na.pass for model frame construction involving covariates for all methods and updated lrstat-package.R to add na.omit and na.pass to reflect this change
• added pbvnorm for the distribution function for standard bivariate normal
• added hazard_pd for the hazard function for progressive disease given the piecewise hazard functions for PFS and OS and the correlation between PD and OS
• added hazard_sub for the hazard function for the biomarker negative subpopulation given the hazard function of the ITT population, the hazard function and the prevalence of the biomarker positive subpopulation
• added the init parameter and the fail flag to the output to logisregr, liferegr, and phregr
• added the nullVariance parameter to getDesignOneProportion
• added details section to the documentation of getDesignPairedPropMcNemar, getDesignRiskDiff, getDesignRiskRatio, getDesignRiskRatioFM, getDesignOddsRatio, getDesignOneMultinom, getDesignTwoMultinom, getDesignTwoOrdinal, getDesignOrderedBinom, getDesignUnorderedBinom, getDesignUnorderedMultinom, pwexploglik, pwexpcuts, and lrschoenfeld
• added references to the documentation of ClopperPearsonCI, BOINTable, mTPI2Table
• added the tol parameter to pwexpcuts
• added the lrsimsub function for log-rank test simulation for enrichment design
• updated lrsim2e and lrsim2e3a to use the hazard_pd function when generating PFS time
• updated kmstat.cpp and rmstat.cpp to report an error message if the information at the first interim is less than the information at the milestone time
• updated lrstat.cpp to avoid the redundant calculation of the final analysis time in the lrpower, lrsamplesize, lrpowerequiv, and lrsamplesizeequiv functions
• updated getDesignMeanDiffCarryover and getDesignMeanDiffCarryoverEquiv to fix the calculation of p and v
• updated getDesignMeanDiff, getDesignMeanDiffXO, and getDesignSlopeDiff to add the parameter value at the null hypothesis to the boundaries for the parameter of interest for noninferiority trials
• updated getDesignPairedMeanDiffEquiv, getDesignMeanDiffEquiv, getDesignMeanDiffXOEquiv, and getDesignSlopeDiffMMRM to simplify the calculation of attainedAlpha for equivalence trials
• updated getDesignWilcoxon, getDesignMeanDiffMMRM, getDesignOneSlope, getDesignSlopeDiff to add details to the documentation
• updated getDesignMeanDiffMMRM to add numberOfCompleters to the output and a reference to the documentation
• updated getDesignMeanDiffCarryover and getDesignMeanDiffCarryoverEquiv to add half_width and nu to the output
• updated getDesignTwoWayANOVA to ensure the sample sizes are multiples of 4
• updated getDesignSlopeDiffMMRM to add more detailed description for the parameters w and N
• updated getDesignSlopeDiffMMRM to add the fixedFollowup parameter to allow fixed follow-up designs and modify the function for computing information and add number of completers accordingly
• updated getDesignSlopeDiffMMRM to add more content to the details section of the documentation
• updated the definition of the Z statistic in the details section of the documentation for getDesignEquiv
• updated print.designMeanDiffMMRM to add numberOfCompleters and analysisTime to the output
• updated print.designMeanDiffCarryover and print.designMeanDiffCarryoverEquiv to add half_width to the output
• updated print.designSlopeDiffMMRM to add numberOfCompleters to the output for fixed follow-up designs
• updated the documentation of the output of lrschoenfeld
• replaced the survreg initial value method with the OLS initial value method for liferegr

lrstat 0.2.14

• updated survival_analysis to ignore intervals not at risk within each stratum without creating non overlapping times across strata
• updated documentation for the survsplit utility function
• removed bc from logisregr
• added residuals_liferegr for residuals from parameteric regression models for failure time data
• added maxiter and eps to logisregr, liferegr, and phregr
• add initial values for liferegcpp
• moved survQuantile from getDesignSurvival to survival_analysis.cpp
• moved the detailed implementation from residuals_phregr to residuals_phregcpp
• add detail description of the spending functions for errorSpent

lrstat 0.2.13

• check whether the rounded up n is different from the original n0 before updating accrual intensity or accrual duration
• change accrualIntensity to accrualIntensity1 in lrsamplesizeequiv when obtaing the timeing of interim analysis
• update kmpower1s and rmpower1s to remove drift and inflation factor from the output while adding number of events, number of dropouts, and number of milestone subjects to the output
• simplify the calculation of underlying survival probability for each treatment group in kmstat1
• update kmpowersamplesize, kmsamplesize1s, kmsamplesizeequiv, rmpowersamplesize, rmsamplesize1s, rmsamplesizeequiv to remove the impossible case when followupTime is missing for fixed follow-up
• update kmpowersamplesize, kmsamplesize1s, kmsamplesizeequiv, rmpowersamplesize, rmsamplesize1s, rmsamplesizeequiv to add a check for solving accrualIntensity when study duration is shorter than or equal to milestone
• update kmpowersamplesize, kmsamplesize1s, kmsamplesizeequiv, rmpowersamplesize, rmsamplesize1s, rmsamplesizeequiv to add a case to decrease accrual duration for fixed follow-up under H0
• update getDesignMeanDiffCarryover to allow specification of the treatment comparison of interest and whether to account for carryover effects in power calculation
• add getDesignMeanDiffCarryoverEquiv to test for equivalence in mean difference for the direct treatment effect of interest in crossover designs
• add the keep_censor parameter to the kmest function to specify whether to retain censoring time points in the output data frame
• add interval to survsplit output
• use the same stopping criteria for SAS PROC LOGISTIC profile likelihood
• add the float_to_fraction utility function
• rename runShinyApp to runShinyApp_lrstat

lrstat 0.2.12

• allow time to take zero values in kmest, lrtest, rmest, and rmdiff
• add the bc parameter for bias correction for noncanonical parameterization and weighted logistic regression
• add loop to look for appropriate bracketing interval for the brent algorithm in lrstat, kmstat, rmstat, and nbstat
• update f_info nbstat to calculate the information directly
• update the mean exposure calculations in nbstat for starting values for reml
• add svdcpp for singular value decomposition of a rectangular matrix
• add number of events, number of dropouts and number of milestone subjects to kmstat and kmpower
• add number of events, number of dropouts and number of milestone subjects to rmstat and rmpower
• remove variance ratios for equivalence trials assuming Wald test statistics with variance of parameter estimator evaluated under the alternative hypothesis
• remove the interval argument from getDurationFromNevents
• add lrschoenfeld for sample size calculation using the Schoenfeld formula under proportional hazards

lrstat 0.2.11

• use numerical integration for lrstat for cases including unequal dropouts between treatment groups, hence removed the numSubintervals argument from lrstat1, lrstat, lrpower, lrsamplesize, lrpowerequiv, and lrsamplesizeequiv
• update shiny app accordingly
• add noncanonical parametrizations for logisregr

lrstat 0.2.10

• use newton-raphson instead of the vmmin algorithm for liferegr and phregr
• removed rpsft with all treatment switching functions moved to trtswitch package
• add the survfit_phregr function to obtain the survivor curve estimates
• add the residuals_phregr function to obtain the martingale, deviance, score, and schoenfeld residuals
• add the shilong dataset
• replace std::ceil(x) with std::ceil(x - 1.0e-12) to handle superficial precision in binary presentation of a double number
• allow rep and stratum parameters to represent more than one character variables in the source data
• add the logisregr function for logistic regression of a binary endpoint.

lrstat 0.2.9

• fix the runtime error: nan is outside the range of representable values of type ‘int’
• add the binary_tte_sim function to simulate two endpoints with one being a binary endpoint and the other being a time-to-event endpoint
• add the rpsft function to estimate hazard ratio using rank-preserving structured failure time model

lrstat 0.2.8

• update std::sort handling for Rcpp character vectors
• use static cast to explicitly cast “long”” of the member function size() to “int” and cast “double” of the floor, ceil, round functions to “int”

lrstat 0.2.7

• add the aml, heart, and tobin datasets
• in basket.cpp, replace y[k] = R::rbinom(1, p[j]) with u = R::runif(0,1); y[k] = u < p[j] ? 1 : 0;
• in misc.cpp, replace const in x = NA_REAL with const int x = NA_INTEGER for x = n, n1, y1, n2, y2
• add delete[] for iwork and work pointers in the quad function
• add simonBayesAnalysis and simonBayesSim for Simon’s Bayesian analysis of basket discovery trials
• update the algorithm for finding the backward image (J, zJ) in getADCI
• update the initial values for finding the confidence interval in getADCI
• add bmini for optimization and invsympd for inverting symmetric and positive defined matrices to utilities.cpp
• add phregr to estimate hazard ratios from the Cox model for right-censored or counting process data with robust variance
• add liferegr for parametric regression of failure time data with uncensored, right censored, left censored, and interval censored data

lrstat 0.2.6

• replace “R_isnancpp(x)” with “x == NA_INTEGER” for integer parameter x
• replace the check for equality of spendigTime and informationRates in settings in prints.R

lrstat 0.2.5

• update survQuantile in getDesignMeans.R to handle new censoring and inestimable quantiles
• update simon2stage to increase the search space of sample size
• add zstatRiskDiff, zstatRiskRatio, zstatOddsRatio, zstatRateDiff, zstatRateRatio for the score test of two-sample crude rates and exposure-adjusted incidence rates
• add mini for minimization of a univariate function over a finite interval
• update powerRiskDiffExact, samplesizeRiskDiffExact, powerRiskRatioExact, samplesizeRiskRatioExact, powerRiskDiffExactEquiv, samplesizeRiskDiffExactEquiv, powerRiskRatioExactEquiv, and samplesizeRiskRatioExactEquiv to divide the interval of the response rate in the control group into 100 subintervals and obtain the minimum value within each subinterval in order to obtain the global minimum. This replaces the random draw of 500 values of the control group response rate for global minimum. In addition, the attained alpha calculation is updated to use the global maximum
• add riskDiffExactPValue, riskDiffExactCI, riskRatioExactPValue, and riskRatioExactCI for exact unconditional test and confidence limits based on the Miettinen & Nurminen score statistics
• add rawdata to be used as an input data frame to test survival analysis functions
• add kmest and kmdiff for Kaplan-Meier estimates of survival curves and milestone survival difference
• add lrtest for log-rank test using the Fleming-Harrington family of weights
• add rmest and rmdiff for estimates of restricted mean survival times and restricted mean survival time difference
• add lower.tail and log.p arguments to ptpwexp and qtpwexp functions
• add the fquantile function to obtain the quantiles of a survival function
• add the pwexploglik function to obtain the profile log-likelihood function for the change points in the piecewise exponential approximation to a survival function
• add the pwexpcuts function to obtain a piecewise exponential distribution that approximates a survival distribution
• update getDesignMeanDiffMMRM and getDesignSlopeDiffMMRM to add accrual and dropout information

lrstat 0.2.4

• update the check for pconfigs of getDesignLogistic to avoid an error on macOS
• replace which_min with which_max in lrsim2e and lrsim2e3a
• replace qromb with quad for numerical integration
• use numerical integration for lrstat
• rename kmest to kmstat and use numerical integration
• add kmpower and kmsamplesize for power and sample size calculation for difference in milestone survival probability
• add kmpower1s and kmsamplesize1s for power and sample size calculation for one-sample milestone survival probability
• add rmst, covrmst, and rmstat for restricted mean survival time analysis
• add rmpower and rmsamplesize for power and sample size calculation for difference in restricted mean survival time
• add rmpower1s and rmsamplesize1s for power and sample size calculation for one-sample restricted mean survival time
• update rmstat, rmpower, rmsamplesize, rmpower1s, and rmsamplesize1s to account for stratification
• add lrpowerequiv and lrsamplesizeequiv for power and sample size calculation for equivalence in hazard ratio
• add kmpowerequiv and kmsamplesizeequiv for power and sample size calculation for equivalence in milestone survival probability difference
• add rmpowerequiv and rmsamplesizeequiv for power and sample size calculation for equivalence in restricted mean survival time difference

lrstat 0.2.3

• issue a warning message when unequal spacing is used with O’Brien-Fleming, Pocock, or Wang-Tsiatis boundaries in getBound
• rename maxInformation to information in the overallResults data frame of the getDesign function output
• add simon2stage for Simon’s two-stage design
• add nbstat to calculate the number of events and information for negative binomial rate ratio
• add nbpower to calculate the power for negative binomial rate ratio test
• add nbpower1s to calculate the power for one-sample negative binomial rate
• add nbpowerequiv to calculate the power for equivalence in negative binomial rate ratio
• add nbsamplesize to calculate the sample size for negative binomial rate ratio
• add nbsamplesize1s to calculate the sample size for one-sample negative binomial rate
• add nbsamplesizeequiv to calculate the sample size for equivalence in negative binomial rate ratio
• add runShinyApp to run a Shiny app for power and sample size calculation for log-rank tests
• add getDesignOneMean for group sequential design for one-sample mean
• add getDesignPairedMeanDiff for group sequential design for paired mean difference
• add getDesignPairedMeanRatio for group sequential design for paired mean ratio
• add getDesignMeanDiff for group sequential design for two-sample mean difference
• add getDesignMeanRatio for group sequential design for two-sample mean ratio
• add getDesignMeanDiffXO for group sequential design for mean difference in 2x2 crossover
• add getDesignMeanRatioXO for group sequential design for mean ratio in 2x2 crossover
• add getDesignPairedMeanDiffEquiv for group sequential design for equivalence in paired mean difference
• add getDesignPairedMeanRatioEquiv for group sequential design for equivalence in paired mean ratio
• add getDesignMeanDiffEquiv for group sequential design for equivalence in two-sample mean difference
• add getDesignMeanRatioEquiv for group sequential design for equivalence in two-sample mean ratio
• add getDesignMeanDiffXOEquiv for group sequential design for equivalence in mean difference in 2x2 crossover
• add getDesignMeanRatioXOEquiv for group sequential design for equivalence in mean ratio in 2x2 crossover
• add getDesignWilcoxon for group sequential design for two-sample Wilcoxon test
• add getDesignMeanDiffMMRM for two-sample mean difference at the last time point from the MMRM model
• add getDesignMeanDiffCarryover for direct treatment effects in crossover trials accounting for the carryover effects
• add getDesignANOVA for one-way analysis of variance
• add getDesignANOVAContrast for one-way analysis of variance contrast
• add getDesignRepeatedANOVA for one-way repeated analysis of variance
• add getDesignRepeatedANOVAContrast for one-way repeated analysis of variance contrast
• add getDesignTwoWayANOVA for two-way analysis of variance
• add getDesignOneSlope for group sequential design for one-sample slope
• add getDesignSlopeDiff for group sequential design for two-sample slope difference
• add getDesignSlopeDiffMMRM for two-sample slope difference from the MMRM model
• add getDesignOneProportion for group sequential design for one-sample proportion
• add getDesignPairedPropMcNemar for group sequential design for McNemar’s test for paired proportions
• add getDesignRiskDiff for group sequential design for two-sample risk difference
• add getDesignRiskDiffExact for exact unconditional test for risk difference
• add getDesignRiskRatio for group sequential design for two-sample risk ratio
• add getDesignRiskRatioFM for the Farrington-Manning score test for risk ratio
• add getDesignRiskRatioExact for exact unconditional test for risk ratio
• add getDesignOddsRatio for group sequential design for two-sample odds ratio
• add getDesignRiskDiffEquiv for group sequential design for equivalence in two-sample risk difference
• add getDesignRiskDiffExactEquiv for exact unconditional test for equivalence in risk difference
• add getDesignRiskRatioEquiv for group sequential design for equivalence in two-sample risk ratio
• add getDesignRiskRatioExactEquiv for exact unconditional test for equivalence in risk ratio
• add getDesignOddsRatioEquiv for group sequential design for equivalence in two-sample odds ratio
• add getDesignFisherExact for Fisher’s exact conditional test for two proportions
• add ClopperPearsonCI for Clopper-Pearson confidence interval for one-sample proportion
• add survQuantile for Brookmeyer-Crowley confidence interval of quantiles of right-censored time-to-event data
• add mTPI2Table for mTPI-2 decision table
• add BOINTable for BOIN decision table
• add mnRiskDiffCI for the Miettinen-Nurminen score confidence interval for two-sample risk difference
• add mnRiskRatioCI for the Miettinen-Nurminen score confidence interval for two-sample risk ratio
• add mnOddsRatioCI for the Miettinen-Nurminen score confidence interval for two-sample odds ratio
• add mnRateDiffCI for the Miettinen-Nurminen score confidence interval for two-sample rate difference
• add mnRateRatioCI for the Miettinen-Nurminen score confidence interval for two-sample rate ratio
• add getDesignOneMultinom for one-sample multinomial response
• add getDesignTwoMultinom for difference in two-sample multinomial response
• add getDesignTwoOrdinal for Wilcoxon test for two-sample ordinal response
• add getDesignOrderedBinom for Cochran-Armitage trend test for ordered multi-sample binomial response
• add getDesignUnorderedBinom for unordered multi-sample binomial response
• add getDesignUnorderedMultinom for unordered multi-sample multinomial response
• add getDesignLogistic for logistic regression
• add getDesignAgreement for Cohen’s kappa agreement coefficient
• add getDesignOneRateExact for exact test of one-sample Poisson rate
• add ptpwexp for distribution function of truncated piecewise exponential distribution
• add rtpwexp for random number generation of truncated piecewise exponential distribution
• add hedgesg for Hedges’ g effect size estimate and confidence interval
• add getDesignEquiv for a generic group sequential equivalence design
• add remlRiskDiff for REML estimates of individual proportions with specified risk difference
• add remlRiskRatio for REML estimates of individual proportions with specified risk ratio
• add remlOddsRatio for REML estimates of individual proportions with specified odds ratio
• add remlRateDiff for REML estimates of individual rates with specified rate difference
• add remlRateRatio for REML estimates of individual rates with specified rate ratio

lrstat 0.2.2

• add the intnorm utility function to integrate a function with respect to a normal density
• add predictive power calculation to adaptDesign
• add the ftrunc function to calculate the adjusted p-values for truncated Holm, Hochberg, or Hommel procedures
• reuse the efficacy and futility stopping boundaries calculated under H1 for H0 in lrsamplesize
• add capabilities to calculate Haybittle & Peto boundaries in getDesign, lrpower, and lrsamplesize
• use informationRates as event fractions for conventional log-rank test and information fractions for weighted log-rank tests in lrpower and lrsamplesize
• match the number of events under H0 with the number of events under H1 for conventional log-rank test and match the information under H0 with the information under H1 for weighted log-rank tests
• remove getCriticalValues and getCumAlphaSpent function in lrstat.cpp
• adjust test-f_lrpower and test-f_lrsamplesize to reflect changes to the definition of informationRates
• rename informationTime to informationRates in lrsim for consistency

lrstat 0.2.1

• use markdown for Roxygen documentations
• rename getAccrualDuration to getAccrualDurationFromN
• replace predictEventOnly with predictTarget for the lrstat function
• add number of subjects reaching the maximum follow-up for fixed follow-up design for the lrstat function
• add efficacyStopping to the getBound function to improve coding efficiency
• add the getPower utility function to improve coding efficiency
• apply only equal spacing of looks for typeAlphaSpending of “OF”, “P”, or “WT” in the getBound function
• replace the drift parameter with Imax and theta parameters in the getDesign function
• calculate alpha when critical values are not missing for the getDesign and lrpower functions
• add expected information under H0 to the getDesign function output
• add rejectPerStageH0, futilityPerStageH0, cumulativeRejectionH0, cumulativeFutilityH0, and attainedAlpha to the output of the getDesign function
• add the getCI function for parameter estimation after termination of a group sequential trial
• add the getRCI function to calculate repeated confidence intervals of a group sequential trial
• add the adaptDesign function for sample size re-estimation and conditional power calculation
• add the getADCI function for parameter estimation using the backward image method after termination of an adaptive group sequential trial
• add the getADRCI function to calculate repeated confidence intervals for an adpaptive group sequential trial
• add the getCP function to calculate the conditional power when the parameter value may vary over time

lrstat 0.2.0

• add fadjpdun to calculate the adjusted p-values for Dunnett-based graphical approaches.

lrstat 0.1.15

• add fstp2seq for stepwise gatekeeping procedures with or without retesting for multiplicity problems involving two sequences of hypotheses.
• add fstdmix to obtain adjusted p-values for standard mixture gatekeeping procedures
• add fmodmix to obtain adjusted p-values for modified mixture gatekeeping procedures

lrstat 0.1.14

• add the getAccrualDuration function to obtain the accrual duration to enroll the target number of subjects.
• add the getDurationFromNevents function to obtain a range of accrual duration to reach the target number of events.
• add the getNeventsFromHazardRatio function to obtain the required number of events given the hazard ratios under the null and alternative hypotheses for a group sequential design.
• allow studyDuration < accrualDuration + followupTime for fixed follow-up in lrpower
• update the handling of rounding for fixed follow-up design in lrsamplesize
• update the handling of null hypothesis for fixed follow-up design

lrstat 0.1.13

• add a rounding argument to lrsamplesize to round up the total sample size and events at each stage.
• add by treatment counts of events, counts, and subjects to lrpower output.
• add results under H0 to lrsamplesize output.

lrstat 0.1.12

• use tolower to make typeAlphaSpending and typeBetaSpending into case insensitive inputs.

lrstat 0.1.11

• Add Kaplan-Meier estimate of milestone survival, Greendwood variance estimate, difference in milestone survival, and Z test statistic for survival difference.

lrstat 0.1.10

• Add drift parameter to the getDesign function to compute power given the drift parameter.
• Update the repeatedPValue function to respect the range of repeated p-values and to allow matrix input of raw p-values.
• Remove repeatedPValueFlag from the fseqbon function.
• Remove numSubintervals from the caltime function.
• Update the description of selected functions, parameters, and output.

lrstat 0.1.9

• Add fwgtmat and fadjpsim to calculate the adjusted p-values for Simes-based graphical approaches.
• update the print method for design, lrpower, and lrsim.

lrstat 0.1.8

• Add spendingTime to getDesign, lrpower, and lrsamplesize to allow the error spending time to be different from the information time.
• Rewrite lrsamplesize to simplify and accelerate the computation for typeOfComputation == “Schoenfeld”.
• Add getBound to obtain the efficacy stopping boundaries for a group sequential design allowing the error spending time to be different from the information time.
• Add fadjpbon to obtain the adjusted p-values for graphical approaches using weighted Bonferroni tests for fixed design.
• Add updateGraph to update the weights and transition matrix after removing a hypothesis from the set of indices of yet to be rejected null hypotheses.
• Add repeatedPValue to Obtain the repeated p-values for a group sequential design based on a given alpha spending function.
• Add fseqbon to obtain the test results for group sequential trials using graphical approaches based on weighted Bonferroni tests with the option to provide repeated p-values for each hypothesis over time.
• Add lrsim3a to perform simulation for three-arm group sequential trials based on weighted log-rank test. The looks are driven by the total number of events in Arm A and Arm C combined.
• Add lrsim2e to perform simulation for two-endpoint two-arm group sequential trials based on weighted log-rank test. The first few looks are driven by the total number of PFS events in two arms combined, and the subsequent looks are driven by the total number of OS events in two arms combined.
• Add lrsim2e3a to perform simulation for two-endpoint three-arm group sequential trials based on weighted log-rank test. The first few looks are driven by the total number of PFS events in Arm A and Arm C combined, and the subsequent looks are driven by the total number of OS events in Arm A and Arm C combined.

lrstat 0.1.7

• Add getDesign for creating a generic group sequential design with constant treatment effect over time.

lrstat 0.1.6

• Add capability for performing noninferiority tests in lrpower, lrsamplesize, and lrsim.
• Add capability for simulating analyses based on calendar times in lrsim.
• Adjust the critical value at the final look if the observed total number of events is less than the planned total number of events in lrsim.
• Retain summary statistics for all stages even after crossing the efficacy and futility boundaries in lrsim.
• Add number of dropouts to lrpower/lrsamplesize and lrsim output.
• Add Schoenfeld method for proportional hazards and conventional log-rank test in lrpower and lrsamplesize.

lrstat 0.1.5

• Replace Inf with 6 and -Inf with -6 for test statistic stopping boundaries to avoid potential memory issue.

lrstat 0.1.4

New features

• Add capability for lrstat to calculate hazard ratios from weighted Cox regression model.
• Add capability for lrsamplesize to calculate absolute accrual rate from relative accrual rate given power, accrual duration, and follow-up duration.

Bug fixes

• Use specified informationRates to calculate Wang-Tsiatis boundaries.
• Use hazard ratios from weighted Cox regression model to determine crossing boundaries on the hazard ratio scale for lrpower.
• Replace stratum-specific output with overall results for lrstat.
• Remove hazard ratio estimate from weighted log-rank test from lrsim output.

lrstat 0.1.3

• Add more statistics to lrpower output.

lrstat 0.1.2

New features

• Add capability for lrpower and lrsamplesize to use error spending functions.
• Add more statistics to lrstat, lrpower and lrsim output.
• Allow user to specify numSubintervals to control approximation.

Bug fixes

• Add parameter checking for lrpower, lrsamplesize, and lrsim.
• Add test files.
• Add print_lrpower.R to print lrpower objects.
• Use informationTime instead of informationRates in lrsim to differentiate information based on weighted log-rank tests score statistic variance from information based on number of events.
• Rename sumstat to overview in lrsim output.

lrstat 0.1.1

• Fix hyperlinks.

lrstat 0.1.0

• Initial release.