aboutsummaryrefslogtreecommitdiff
path: root/R
AgeCommit message (Collapse)AuthorFilesLines
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
2018-09-15Remove correlation test as gmkin on windows still crashesJohannes Ranke1-18/+10
Clean up the code a bit
2018-09-15Respect reps argument in add_errJohannes Ranke1-2/+2
2018-09-14Update NEWSJohannes Ranke1-2/+4
Static documentation rebuilt by pkgdown
2018-09-14Avoid nls crashes by checking for a correlation firstJohannes Ranke1-0/+7
2018-09-13Enable new error model in gmkinJohannes Ranke1-7/+20
Also improve the irls fitting of the error model and add a test for FOCUS_2006_C where the second component of the error model is zero
2018-07-18Correct references to the Rocke and Lorenzato modelJohannes Ranke3-6/+6
Rename 'sigma_rl' to 'sigma_twocomp' as the Rocke and Lorenzato model assumes lognormal distribution for large y. Rebuild static documentation.
2018-07-04Improve IORE fitting, use 1.1 as N starting valueJohannes Ranke1-1/+4
2018-02-02mkinfit: Work around a bug in the current FME versionJohannes Ranke1-6/+15
that prevented the convergence message to be returned in the case of non-convergence.
2018-02-02Show R and mkin version used for fitting in the summaryJohannes Ranke1-2/+15
For fits generated with previous version, the version numbers used for fitting were not stored in the fit object. Therefore, the versions used for fitting can only be shown for fits generated with mkin containing this commit.
2018-01-30Reweighting with two-component error modelJohannes Ranke2-20/+92
Static documentation except articles rebuilt by pkgdown
2017-11-16Add the code used for generating synthetic_data_for_UBAJohannes Ranke1-3/+3
Static documentation except articles rebuilt by pkgdown
2017-09-06Use full fitted curve(s) for y axis scaling ...Johannes Ranke1-2/+1
... when using plot_sep() or plot.mkinfit(..., sep_obs = TRUE)
2017-09-06Fix comment that was copied without checkingJohannes Ranke1-1/+1
2017-09-06Fix scaling of residual plots ...Johannes Ranke1-2/+6
... in the case of separate plots for each observed variable as obtained with plot_sep() or plot.mkinfit(..., sep_obs = TRUE)
2017-07-21Rename twa to max_twa_parent, update docsJohannes Ranke1-2/+2

Contact - Imprint