*C files: Change encoding to utf-8 so my name in comments is read correctly.
ltx2: Fix bugs in output.
GB_zeroin: C files removed.
Add logLik methods to coxreg and phreg.
coxreg: drop1 method added, handles missing values.
ltx2: New generic function. Similar to ltx, but prints confidence limits instead of standard errors. Prints LRT p-values, automatic use of drop1.
hazards.coxreg: Add start at zero to cumulative hazards plotting.
make.communal: Deal with zero length interval output and allow input for 'com.dat' with more than one variable.
Vignettes updated.
swepop, swedeaths: Data for the year 2020 included.
plot.tpchreg: Error with arguments 'xlab' and 'ylab' fixed.
coxreg: In case of request for bootstrap, 'boot.sd' is no longer calculated and returned: Caused too frequent errors due to ill balanced bootstrap samples, and was not really helpful anyway.
regtable: New function for creating simple regression output tables.
pchreg: New function, to replace phreg(..., dist = "pch").
phreg: 'dist = "pch"' deprecated.
age.window: No warning for an empty result, returns NULL in that case.
compHaz: New function, graphically comparing cumulative hazards.
plotHaz: New argument 'interval'.
hazards: Add method for aftreg objects.
hazards: Bug with 'pfixed = TRUE' fixed.
print.summary.aftreg: Error in heading fixed.
tpchreg: try for singular hessian, continue without error.
ltx: Move caption to top of table (not really a bug fix).
ltx.tpchreg: Bug fix (really!)
ltx.phreg: Remove double printing of last few rows.
init.c: Fix internal bug in 'breslow_rs2'.
ltx: Introduce \footnotesize for tables.
hazards: New generic function calculating baseline hazards from fits. Methods for regression outputs.
logrank: New function performing a log-rank test. A plot method comes with it.
coxfunk: New function for calculating partial likelihood with derivatives. Warning: Under development, do not trust results yet!
summary method to tpchreg with a print method.
Two new data sets, swepop and swedeaths: Tabular statistics, Sweden 1968-2019.
The data set swepop07 removed, replaced by swepop and swedeaths.
coxreg: Add log likelihood and df to null models.
The data set child is expanded with info about mother.
print.summary.coxreg: New argument 'short' added.
plot.phreg: Argument 'new.data' removed, 'score' added.
Gompertz: Argument 'rate' added to all functions.
risksets: New argument 'collate_sets', a logical. Defaults to FALSE, which conserves old behaviour. If TRUE, output is a list of lists, where each level 2 list describes one risk set.
coxreg: New argument 'coxph' (logical) is introduced. If TRUE, survival::agreg.fit (or survival::coxph.fit) is called for in the optimization. If FALSE, 'own' optimization is called for, can be slow for (even moderately) large data sets (many distinct event time points). Default value is TRUE.
coxreg: Arguments 'hazards' and 'center' deprecated.
phreg: Argument 'center' deprecated.
coxreg, phreg, aftreg: Fixed bug in calculating weighted means, which affects dates as covariates.
plot.hazdata: Fix y limits.
coxreg: Add stratum names when no covariates.
survival moved from Depends to Imports.
summary: Failed for null models with tpchreg, pchreg, now fixed.
summary.phreg: Bug fixed. It was introduced by differences between the distribution 'pch' and the other distributions. The piecewise constant hazards regression will eventually break out of phreg and form a separate function, named 'pchreg'. Have the same thoughts about the Gompertz distribution, to be put in 'gompreg'. We'll see... Main reason is that these guys need special treatment in the programming.
phreg: Fix issue with 'shape = ' something positive.
gompstartRate: Simpler starting values introduced (internal).
getsums.f: Replace 'dcopy' in three instances to please gfortran 10.
risksets: Bug fix: Now works when members = FALSE.
Rewriting tpchreg.
Add tpchreg methods: print, plot, extractAIC.
Add new function: oe (for occurrence/exposure data).
Submit to CRAN.
tpchreg: Preparations for pretty-printing.
print.tpchreg: Printing method for tpchreg objects introduced.
More fixes to work around the gfortran bug.
Add C wrapper, cdgemv, for Fortran call to BLAS subroutine dgemv.
risksets: Warning about sample_fraction fixed.
Errors in plot.aftreg fixed (thanks to Mark Clements).
In eha_fun.c, changed 'finite' to 'isfinite' (C99).
Fix tpchreg with extractAIC and nobs.
New functions tpchreg and toTpch.
Bug in phreg, Gompertz distribution, 'param ="rate"', fixed.
Bug in phreg, Gompertz distribution, 'param ="rate"', with strata, fixed.
Error in plot.phreg, Gompertz with 'param = "rate"' fixed.
New function: mpch gives the mean of a piecewise constant hazards distribution.
The function plot.Surv is removed, since it now is a plot method from the survival package, which eha depends on. Note though that this implies a new set of arguments. See the help page of survival:::plot.survfit.
New versioning system (Hadley Wickham's suggestion).
Documentation 'roxygenized'.
Development moved to gitHub.
The package glmmML was merged into eha in version 2.0. This has turned out to be a very bad idea ("double entry"), and all relevant functions are now Defunct and removed. Use glmmML instead.
glmmML(), glmmboot(), ghq(), glmmbootFit(), glmmML.fit() are all Defunct and removed. No reverse dependency used any of them. Those functions are all found in the package glmmML, on CRAN.
phreg(): Start values by 'coxreg' (bad idea) removed.
aftreg(): in 'addMeans' and 'aftreg1', an error
(spotted by Jingchunzi Shi) with fixed scale corrected.
init.c: Added after R-3.4.0.
coxreg(): "Geometric bootstrap" removed (it never worked and was
never used).
cal.window(), age.window(): Works now with 'tibbles'.
rpch: Bug fixed (spotted by Brigitte Mueller).
plot.phreg: parameters 'changeLog' and '...' enhanced. Documentation updated.
plot.coxreg: Bug concerning 'printLegend' fixed.
weibreg: Note added in documentation: "Use phreg instead". weibreg will soon be deprecated.
coxreg: Changed the default for 'hazards' to TRUE.
getHaz: Modified to handle missing strats and score.
print.risksets: Added '...' to argument list.
coxreg: Fitting a NULL model didn't produce enough information like log likelihood ('loglik'), for comparison of models. Fixed.
glmmboot: 'frail' set to Inf or -Inf for appropriate clusters.
plot.hazdata: Updated to acknowledge the argument 'hazards' that was included in 'coxreg' in version 2.4-2.
plot.coxreg: Updated to acknowledge the argument 'hazards' that was included in 'coxreg' in version 2.4-2. Errors in the manual page corrected.
risksets: Risk set sampling did not always work (memory leak in C code); fixed.
risksets: New value in the Value list: sample_fraction.
plot.aftreg: Now registered as an S3 method.
print.aftreg: Scale and shape parameters are only printed on the log scale.
src: 'abs' corrected to 'fabs' in two places in C code.
coxreg: Argument 'hazards' added, with default = FALSE.
plot.Surv: Argument 'printLegend' added.
coxreg: Changed labels in case of ordered factors.
Foreign calls: 'DUP = FALSE' removed everywhere.
print.pch: Bug preventing printing(!) fixed.
plot.pch: Bug preventing plotting(!) fixed.
coxreg: Call to survival::[agreg,coxph].fit in "standard" cases reintroduced; Thanks to Terry Therneau, who exports these functions from survival_2.37-6. So this version of 'eha' requires survival (>= 2.37-6). Argument 'center' reintroduced, with slightly different meaning.
plot.coxreg: argument 'newdata' not used.
phreg: Stratified analysis in the piecewise constant hazards model introduced. Argument 'center' reintroduced, with slightly different meaning.
plot.phreg: Argument 'newdata' not used.
plot.coxreg: added the argument col, for color line(s).
plot.hazdata: Ditto.
coxreg: Restoring equity.
glmmML: Removed spurious printing in 'glmmml.c'.
coxreg: 'Downgraded' due to limitations in using non-exported function (agreg) from 'survival'. This will be taken care of ASAP. Meanwhile, live with slightly slower 'coxreg' and choking on huge data sets.
Coxreg: Removed for the same reason as the downgrade of coxreg.
phreg: The lognormal and loglogistic distributions have an intercept added to the linear predictor; in effect extending these distributions to three-parameter families. A bug (found by Claire Williams) affecting the covariance matrix of the parameter estimates has been eliminated.
aftreg: The same bug as for phreg eliminated. The available parametrizations are changed to 'lifeAcc' and 'lifeExp'. For the Gompertz distribution, the 'canonical' representation is used regardless.
hlnorm, Hlnorm, hllogis, Hllogis: New parameter, 'prop', added. It simply multiplies the previous definitions by its value.
age.window, cal.window: speed up and handling of empty intervals
perstat: speed up.
plot.phreg, plot.coxreg: legend added for stratified fits.
Coxreg: New function designed to handle huge data sets. coxreg chokes on 'integer overflow', and memory overflow, when the collection of risk sets is created and saved. Coxreg is essentially a wrapper for 'coxph' in the 'survival' package, but with methods for nice(r) printing.
coxreg: Error in the calculation of 'total time at risk' fixed (gave negative values when no left truncation!).
phreg: The gompertz 'scale' parameter is now a real scale parameter through a reparametrization of the baseline gomperz distribution. It follows now the 'canonical' parametrization. There is also a new argument, 'param', which can take the values "canonical" (default) and "rate". Better start values for the Gompertz regressions.
check.dist : Has failed due to different scaling in coxreg and phreg since 2.2-2. Now fixed.
risksets: Added an argument, 'members', a logical with default value 'TRUE'. If TRUE, the members of all risk sets are returned. With huge data sets, this may be too much, and not wanted. In such cases, set 'members = FALSE'.
coxreg: Memory leaks fixed in getHaz.f, sizes.c and risk.c.
Gompertz distribution: new parameter 'param'; values 'default' (default, nothing changes), and 'canonical', which gives meaning 'PH' and 'AFT' to the shape and scale parameters, respectively.
aftreg: New input argument: 'param', with three values, 'default' (which is the default;), 'survreg', and 'canonical'. The value 'default' uses the standard parametrization, while 'survreg' uses the same parametrization as the survreg function in the survival package. 'canonical' is especially useful in the Gompertz distribution for separating PH and AFT parameters.
C code cleaned (avoiding some Warnings at CRAN).
check.dist: New argument, col, added. Default is now black-and-white plotting.
ltx: New function for printing coxreg fits in latex format.
plot.Surv: New arguments, col and lty, added. The default is now black-and-white plotting.
data: Two new data sets, "logrye" and "scania".
oldmort: Data set updated.
aftreg.fit: Errors in documentation corrected.
SurvSplit: Bug fixed. It works now.
glmmboot: Bug in calculation of 'frail' fixed.
plot.phreg: Bug in argument 'ylim' fixed. Bug with stratified fits fixed.
aftreg, phreg: summary method added (just printing).
glmmML: Merged into eha. This adds binary and count data regression for clustered data to eha. glmmML is still a standalone package on CRAN, for how long remains to be seen.
Pch: New class of distributions; Piecewise constant hazards (Pch).
phreg: piecewise constant baseline hazard function introduced.
phreg: 'center' reinstalled as a logical, default = TRUE, so nothing changes if 'center' is left untouched.
age.window: If the selected window is empty, an error was previously thrown. Now NULL is returned, with a warning.
cal.window: Now returns NULL, with a warning, if the selected window is empty.
aftreg: Gompertz distribution added (back again). Better start values.
phreg: Better start values for Gompertz distribution.
eha: Depends (R >= 2.13.0)
phreg, aftreg, weibreg: Errors in documentation corrected.
ChangeLog: Errors in 1.3-6 corrected.
coxph, phreg, aftreg: Bug fix: Logical covariates now OK. (converted to factor).
phreg: Removed "coxph" from "class".
aftreg: Ditto.
weibreg: Ditto.
oldmort: New data set.
coxreg: added return values 'df' and 'n'.
coxreg, phreg, aftreg: Added extractAIC and nobs methods.
toBinary: Now includes also risksets without survivors (if any).
fert: Row names changed to row numbers.
swe07: New data set.
Spelling errors (some!) corrected in documentation.
male.mortality: Old name (male.mortality) of data frame mort reinstalled. Now known under two names. Overlapping spells no longer exist.
fert: Documentation corrected.
toDate: Now works as intended.
dllogis: Bug fixed.
phreg: Fixed error and (some) numerical instability in the estimation of the Gompertz proportional hazards model (seems to be very dependent on good start values!).
plot.phreg: Fixed 'ylim' for cumulative hazards plot.
male.mortality: Data set renamed to
mort: The new name. (Just for convenience!)
coxreg: Major revision, calls 'coxph' in 'survival' (thanks to Terry Therneau!) for ordinary requests. Reason: speed. The argument 'center' is deprecated, centering of covariates is now routinely made. Affects the estimation of the baseline hazard, which is done at the means of the columns of the original design matrix.
phreg: Argument 'center' is deprecated. Centering is routinely performed, so the baseline hazard is estimated at the means of the columns of the original design matrix.
aftreg: Corrected error in documentation and made the argument 'center' deprecated.
aftreg: Estimation of intercept ("scale") is now comparable with corresponding estimate of Intercept in 'survreg'. 'shape' in aftreg is still '1 / scale' in survreg. Regression coefficients still have different signs but the same absolute values.
join.spells: a few bugs fixed, now allows for factor covariates.
male.mortality: data set "joined".
fert: is a new data set. Contains data on 19th century marital fertility from a parish in northern Sweden.
aftreg.fit: Check of data consistency in case of more than one record per individual added.
aftreg.fit: Fixed memory allocation bug in C code. This bug was not fatal (too much memory was allocated), but caused a severe performance problem, especially on Windows.
aftreg: 'dist = "gompertz"' should now work (realizing that it is a special case of 'dist = "EV").
phreg: Fixed bugs for the cases with fixed shape (includes 'dist = "gompertz"'.
aftreg: 'offset' works now.
phreg: Ditto.
check.dist: Bug fixed for 'dist = "gompertz"' (again!)
plot.phreg: Ditto.
#eha 1.2-15
aftreg: argument 'id' now works as intended, i.e., is looked for in 'data' (or 'parent.frame()')
print: functions return invisibly 'x'.
coxreg: argument 'weights' now works as intended, i.e., is looked for in 'data' (or 'parent.frame()')
check.dist: Bad 'ylim' fixed.
aftreg: Gompertz distribution had a bug; now fixed.
coxreg: Added (time-varying) weighted Cox regression. Returns error code 1 if failure, with a warning.
check.dist: Can now handle stratified models.
Bugs in qEV, rEV, all Gompertz functions, plus documentation of EV and Gompertz fixed. Note that these bugs did not affect aftreg or phreg! However,
check.dist had to be revised due to the above bugs.
Bug in coxreg, method = "ml" fixed.
Bug in coxreg, calculation of hazard, fixed.
Several errors in the documentation fixed.
Added full profiling to the "ml" and "mppl" methods in coxreg.
mlreg is now deprecated. Use coxreg with method = ml or = mppl
table.events: New parameter: 'strict = TRUE'.
make.communal: In case of 'communal=FALSE', spells are no longer truncated. Instead a value of NA is given if birthdate is out of range. Previous behavior was unpredictable when birthdate was out of range.
mlreg: Geometric distribution (i.e., constant baseline discrete hazard) added. Not for frailty models, yet (on the TODO list).
mlreg: New argument, 'n.points', added to 'control'. Controls the number of points used in the Gauss-Hermite quadrature.
mlreg: Stricter control of numerical problems, especially in the frailty fit.
clean: Replaced by check.spells and join.spells.
Return values changed to conform with R-1.8.0 (and later).