Add evaluateKinetics.py script to testing folder #117
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Raise errors if depository isn't found. Remove averaging and adding of training reactions in the family, assumes you have done it outside the function. Compares reverse kinetics if the NIST reaction is in the reverse direction from the family template
retrieves data involving the family Now shows the top node names, and also the number of rules csv file has headers now. Does not average the families or add training data. Assumes any type of manipulation has been done outside the folder. Also, don't load any kineticsLibraries or thermoLibraries by default except the primaryThermoLibrary. But do load all the depositories
Add NIST folder to ketoenol. Fix atom labeling for several reactions
Is not quite the original leaveOneOut test as it was done before. The original test takes extremely long to run because it must average up the family after each exact rate rule. In this method, we remove the original rate rule, and see if using an average that (should incorporate the original) does a good job of estimating that rate rule. The option to use the original version of the algorithm can still be used by setting averaging=True Write data to csv file after each family rather than after entire set of families so that they are in alphabetical order
|
|
||
| The original family should not contained averaged nodes when starting out. The | ||
| leave one out test should be performed only for original exact matches. | ||
| """ |
There was a problem hiding this comment.
Need to update comment here explaining about the averaging boolean
|
Maybe |
|
I think |
This used to excise parity data, it's better to keep them for complete statistical analysis
…MG kinetics in parityData User gets to choose whether it wants all kinetics in the parity Data (both exact and non-exact), or only the exact ones
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This is a revised and updated script of Nathan's old kinetics evaluation script containing 3 main functionalities:
obtainKineticsFamilyStatistics: gathers and saves data about families into a .csv files containing info on the number of exact rate rules, number of groups per tree in the groupsleaveOneOut: performs a leave one out test by deleting an exact node (excludes averaged ones), averaging the tree and then trying to estimate that tree. Creates parity plots for each family and creates a .csv of the root mean square error.compareNIST: compares the family estimated kinetics against the exact kinetics from the NIST depository within each family. Also creates parity plots and .csv of error data.