aboutsummaryrefslogtreecommitdiff
path: root/R
AgeCommit message (Collapse)AuthorFilesLines
2022-12-19Rename template folder, create formatJohannes Ranke2-1/+40
Instead of rmarkdown::pdf_document, mkin::hierarchical_kinetics is used as a document format in the template. In this way, the template file can be freed from some R code and yaml options that the average user does not have to be aware of.
2022-12-15List random effects correlations in output if anyJohannes Ranke2-4/+20
Update docs
2022-12-07Fix parplot for the case of failed multistart runsJohannes Ranke2-2/+2
2022-12-07Check slopes in saemix covariate modelsJohannes Ranke1-1/+13
2022-12-06Some parplot improvementsJohannes Ranke1-4/+12
llquant argument, improved legend text, tests
2022-12-02Avoid redundant warnings in summariesJohannes Ranke1-2/+3
2022-12-01Possibility to manually specify no_random_effects in mhmkinJohannes Ranke1-20/+71
2022-11-29Complete starting values in summary for saem.mmkin fitsJohannes Ranke1-1/+10
Also update tests to the changes in mhmkin (see NEWS)
2022-11-25mhmkin: Easy specification of ill-defined parmsJohannes Ranke1-18/+20
The argument 'no_random_effect' now accepts an illparms.mhmkin object
2022-11-18Make mixed model test data permanent to ensure reproducibilityJohannes Ranke1-0/+17
To ensure that tests on different platforms work on the same data, the mixed modelling test data previosly generated in tests/testthat/setup_script.R were generated once using the script in inst/dataset/generation/ds_mixed.R, and are now distributed with the package.
2022-11-18Pass error.init to saemix_model, show in parplotJohannes Ranke2-5/+17
Due to an oversight, error.init was not really passed to saemix_model in saem.mmkin. The new initial values were reverted to c(1, 1), in order to avoid changing the test results. Initial values for error model parameters are now shown in parplot.multistart.
2022-11-18Improve logLik.mkinfit to attach nobs attributeJohannes Ranke1-0/+1
The lack of that attribute lead to a failure to calculate the BIC in test_AIC.R on R-devel from yesterday.
2022-11-16Read in all data per defaultJohannes Ranke1-1/+1
2022-11-15Import saemix for easier profiling, update online docsJohannes Ranke1-1/+2
2022-11-14Automatic starting parameters for saem.mmkinJohannes Ranke1-1/+28
For the case of mkin transformations. This gives faster convergence, and appears to avoid problems with numeric ODE solutions
2022-11-08Same SFORB simplification for saem.mmkin fitsJohannes Ranke1-2/+2
2022-11-08Show DFOP g equivalent of SFORB in summary, testJohannes Ranke2-8/+11
2022-11-08Simplify SFORB also in the endpoint functionJohannes Ranke1-11/+10
2022-11-08Simplify SFORB analytical solution, whitespaceJohannes Ranke1-64/+64
I do not know why the formulae for b1 and b2 on page 64 of FOCUS kinetics (2014) were not simplified. Clearly, the term k12 * k21 - (k12 + k1output) * k21) can be simplified to - k1output * k21 The test for equivalence of DFOP and SFORB fits verifies that the change is OK. I also removed trailing whitespaces, substituted tab characters by two whitespaces and removed indenting of text in paragraphs describing parameters in roxygen comments to unify formatting.
2022-11-07Make print and summary for saem.mmkin more robustJohannes Ranke2-8/+18
If the likelihood computed by importance sampling, these methods failed. Now they report "Not available" or NA, respectively.
2022-11-04Further improvement of robustness of anova methodJohannes Ranke1-9/+9
2022-11-04Fix anova.saem.mmkin for unavailable likelihoodsJohannes Ranke1-4/+10
2022-11-04Attempt at automatic setting of random effectsJohannes Ranke1-4/+24
Based on parameters in the separate fits that fail the t-test.
2022-11-04Make it easy to override maxsteps for lsodaJohannes Ranke1-23/+26
The hope was to be more successful in saemix fits, but I did not notice an improvement (e.g. for a failing likelihood calculation with importance sampling, where the error message suggested to increase maxsteps).
2022-11-04Make illparms.mhmkin and anova.saem.mmkin more robustJohannes Ranke3-12/+14
2022-11-04Never fail with errors in saem.mmkinJohannes Ranke1-6/+0
I think it is always preferable to be able to investigate the saemix object returned.
2022-11-01Add tex_listing()Johannes Ranke1-0/+32
2022-11-01Add read_spreadsheetJohannes Ranke1-0/+113
2022-11-01Fix plotting saem fits with covariatesJohannes Ranke2-47/+72
2022-11-01Fix a bug, couldn't set class attribute on NULLJohannes Ranke1-1/+1
2022-11-01Fix illparms testsJohannes Ranke1-1/+1
2022-11-01Improved printing for illparms()Johannes Ranke1-4/+20
For the case of single fits and no ill-defined parameters found
2022-10-31Fix R CMD checkJohannes Ranke1-3/+5
2022-10-31Possibility to override the error model in update.saemJohannes Ranke2-11/+14
2022-10-31Add anova.mhmkin methodJohannes Ranke1-0/+16
2022-10-30Correct documentation of sorting in anova.saem.mmkin()Johannes Ranke1-3/+4
2022-10-30saem summary in case of non-invertible FIM partsJohannes Ranke1-7/+13
2022-10-29Update multistart docsJohannes Ranke1-2/+2
2022-10-28Rename parhist to parplot and make it genericJohannes Ranke4-8/+15
That parhist name was not the brightest idea, as it does not show histograms.
2022-10-28Export status method for mhmkin, move to status.RJohannes Ranke2-39/+43
2022-10-28Rename 'convergence' method to 'status'Johannes Ranke5-27/+25
The reason is that it is misleading in the case of saem.mmkin objects, because convergence is not really checked there.
2022-10-26Really test anova output, test covariatesJohannes Ranke1-0/+5
2022-10-26Update multistart example code and vignetteJohannes Ranke1-9/+3
2022-10-26Remove kernel density estimate from llhistJohannes Ranke1-11/+9
2022-10-26Make multistart work with testthat scopesJohannes Ranke1-3/+20
Also skip long-running test with saem using the analytical from mkin
2022-10-25Complete documentation and fix a bugJohannes Ranke3-2/+6
The bug was introduced by the changes in summary.saem.mmkin.R and surfaced in the tests when using saemix transformations.
2022-10-25Revert accidental changeJohannes Ranke1-21/+27
2022-10-21Explain chisq tests in anovaJohannes Ranke1-1/+2
Now that I understand how anova.lmer from which I ported the code produces mostly meaningful likelihood ratio tests in the table.
2022-10-21Adapt summary.saem.mmkin to covariate modelsJohannes Ranke1-7/+6
Dissipation times are calculated for covariate values of zero currently, which will usually not be very useful.
2022-10-21Add a simple anova method for model comparisonJohannes Ranke2-4/+116
Update docs

Contact - Imprint