aboutsummaryrefslogtreecommitdiff
path: root/R
AgeCommit message (Collapse)AuthorFilesLines
2019-10-26parms and confint methodsJohannes Ranke3-22/+190
The confint method can do profile likelihood based confidence intervals!
2019-10-25Use roxygen for functions and methodsJohannes Ranke33-2087/+3364
2019-10-22Go back to d_3 for errmod = "obs"Johannes Ranke1-1/+1
in the hope that this makes plotting cross-platform also for this error model
2019-10-21Improve some plotting routines, more testsJohannes Ranke2-26/+36
Static documentation rebuilt by pkgdown
2019-10-21Improve tests, remove geometric_meanJohannes Ranke1-4/+0
2019-10-21Refactor mkinfit, infrastructure workJohannes Ranke2-157/+178
mkinfit objects now include an ll() function to calculate the log-likelihood. Part of the code was refactored, hopefully making it easier to read and maintain. IRLS is currently the default algorithm for the error model "obs", for no particular reason. This may be subject to change when I get around to investigate. Slow tests are now in a separate subdirectory and will probably only be run by my own Makefile target. Formatting of test logs is improved. Roundtripping error model parameters works with a precision of 10% when we use lots of replicates in the synthetic data (see slow tests). This is not new in this commit, but as I think it is reasonable this closes #7.
2019-09-18Correct and rebuild docsJohannes Ranke1-9/+21
- Reconcile docs and code for max_twa_parent - Correct links to docs in twa vignette - Static documentation rebuilt by pkgdown
2019-09-02Some changes to improve plots on beamer slidesJohannes Ranke3-15/+22
2019-07-09Adaptations for gmkinJohannes Ranke1-1/+3
Address winbuilder check problems, update check log, update of static docs
2019-07-08Support summarizing older fitsJohannes Ranke1-8/+10
generated with mkin < 0.9.49.5
2019-07-05Static documentation rebuilt by pkgdownJohannes Ranke1-0/+1
2019-06-05Add error model algorithm to outputJohannes Ranke1-8/+19
2019-06-04Fix a bug introduced in the last commitJohannes Ranke1-2/+4
2019-06-04For the d_3 algorithm, report which was better, if anyJohannes Ranke1-2/+13
2019-06-04Additional algorithm "d_c", more tests, docsJohannes Ranke1-46/+71
The new algorithm tries direct optimization of the likelihood, as well as a three step procedure. In this way, we consistently get the model with the highest likelihood for SFO, DFOP and HS for all 12 new test datasets.
2019-06-04Algorithms direct, two-, three-, fourstep, IRLSJohannes Ranke1-9/+41
All of them are working now and allow for comparison Based on SFO, DFOP and HS fits to twelve test datasets, only the combination of direct and threestep is needed to find the lowest AIC
2019-06-03Status von Samstag morgen - untestedJohannes Ranke1-37/+63
2019-05-23Fix a bug that could lead to wrong chi2 error levelsJohannes Ranke1-0/+1
2019-05-20Only use colors for error model "obs"Johannes Ranke1-3/+3
2019-05-08Add functionality to plot the error modelJohannes Ranke4-15/+116
by plotting squared residuals against predicted values, and showing the variance function used in the fitted error model. Rebuild docs
2019-05-07Fix hessian calculationsJohannes Ranke1-7/+8
Static documentation rebuilt by pkgdown
2019-05-07After the OLS step, use OLS parameter estimatesJohannes Ranke1-24/+24
- Fix the respective error in the code - Static documentation rebuilt by pkgdown
2019-05-03Improve error model fittingJohannes Ranke2-15/+36
Now we have a three stage fitting process for nonconstant error models: - Unweighted least squares - Only optimize the error model - Optimize both Static documentation rebuilt by pkgdown
2019-05-02Better initials for error model parametersJohannes Ranke1-8/+21
- Also make it possible to specify initial values for error model parameters. - Run tests - Rebuild docs
2019-04-24Improve output and update testsJohannes Ranke1-5/+9
Remove skipped tests as I do not intend to reactivate them
2019-04-10Address failure on travisJohannes Ranke1-1/+1
Also reduce the digits in the representative half-live given by nafta()
2019-04-10Adapt tests, vignettes and examplesJohannes Ranke3-68/+90
- Write the NEWS - Static documentation rebuilt by pkgdown - Adapt mkinerrmin - Fix (hopefully all) remaining problems in mkinfit
2019-04-08Remove zero observations, adapt logLikJohannes Ranke2-24/+11
Zero observations at time zero made fitting the two-component error model fail. A concentration of exactly zero does not make sense anyways, as we generally have a limit of detection
2019-04-05Report negative log-likelihood during fitsJohannes Ranke1-4/+10
If not quiet = TRUE
2019-04-04Direct error model fitting worksJohannes Ranke1-366/+133
- No IRLS required - Removed optimization algorithms other than Port - Removed the dependency on FME - Fitting the error model 'obs' is much faster for the FOCUS_2006_D dataset and the FOMC_SFO model (1 second versus 3.4 seconds) - Vignettes build slower. Compiled models needs 3 minutes instead of 1.5 - For other vignettes, the trend is less clear. Some fits are faster, even for error_model = "const". FOCUS_Z is faster (34.9 s versus 44.1 s) - Standard errors and confidence intervals are slightly smaller - Removed code for plotting during the fit, as I hardly ever used it - Merged the two cost functions (using transformed and untransformed parameters) into one log-likelihood function
2019-03-04signif for DTx when printing nafta resultsJohannes Ranke1-1/+1
In order to be able to test cross-platform
2019-03-04Review tests and examples based on winbuilder resultsJohannes Ranke1-3/+3
Static documentation rebuilt by pkgdown
2019-02-27Comment out unused codeJohannes Ranke2-27/+27
2019-02-27Increase test coverageJohannes Ranke1-8/+17
2019-02-26NAFTA SOP Attachment vignetteJohannes Ranke2-3/+7
Also: - Change rounding in print.nafta - Add dots argument to nafta() - Use cores=1 in examples - Restrict N in IORE model to values > 0
2019-02-26Add another NAFTA dataset + testJohannes Ranke1-1/+2
2019-02-25Static documentation rebuilt by current pkgdownJohannes Ranke1-3/+7
2019-02-22Implement the NAFTA evaluation schemeJohannes Ranke1-0/+105
2019-02-22Further adaptations for the logistic modelJohannes Ranke3-4/+9
- Make mmkin compatible - Return DT50 values corresponding to k0 and kmax - Turn incompatible parameter names in parms.ini from an error to a warning, in order to make it possible to use this argument in calls to mmkin
2019-02-21Add the logistic modelJohannes Ranke6-17/+43
2019-01-31Restore NAMESPACEJohannes Ranke1-16/+1
which was accidentally overwritten by pkgdown -> roxygen
2019-01-31Export to CAKE study fileJohannes Ranke1-0/+80
2018-11-26Add AIC method for mmkin column objectsJohannes Ranke1-0/+30
Static documentation rebuilt by pkgdown
2018-11-26Support logLik and AIC for manual weightsJohannes Ranke1-2/+11
Static documentation rebuilt by pkgdown
2018-11-23Add logLik method to enable AIC() on mkinfit modelsJohannes Ranke2-1/+36
Further relax two tests to pass build on Travis
2018-11-23Improve initialisation of mkindsJohannes Ranke1-10/+13
Check in work from the beginning of the week
2018-09-21Improve fitting the two-component error modelJohannes Ranke1-47/+83
with respect to accuracy and robustness.
2018-09-20Use the median absolute deviance for fittingJohannes Ranke1-5/+16
as the absolute value is a biased estimator for the standard deviation
2018-09-18Make mkinpredict generic and write a method for fitsJohannes Ranke1-25/+52
2018-09-18Support maximum twa for HS, test, docsJohannes Ranke1-2/+16

Contact - Imprint