{"action":"closed","issue":{"url":"https://api.github.com/repos/D3DEnergetic/FIDASIM/issues/28","repository_url":"https://api.github.com/repos/D3DEnergetic/FIDASIM","labels_url":"https://api.github.com/repos/D3DEnergetic/FIDASIM/issues/28/labels{/name}","comments_url":"https://api.github.com/repos/D3DEnergetic/FIDASIM/issues/28/comments","events_url":"https://api.github.com/repos/D3DEnergetic/FIDASIM/issues/28/events","html_url":"https://github.com/D3DEnergetic/FIDASIM/issues/28","id":48119938,"number":28,"title":"Error Estimation","user":{"login":"lstagner","id":4305891,"avatar_url":"https://avatars.githubusercontent.com/u/4305891?v=3","gravatar_id":"","url":"https://api.github.com/users/lstagner","html_url":"https://github.com/lstagner","followers_url":"https://api.github.com/users/lstagner/followers","following_url":"https://api.github.com/users/lstagner/following{/other_user}","gists_url":"https://api.github.com/users/lstagner/gists{/gist_id}","starred_url":"https://api.github.com/users/lstagner/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/lstagner/subscriptions","organizations_url":"https://api.github.com/users/lstagner/orgs","repos_url":"https://api.github.com/users/lstagner/repos","events_url":"https://api.github.com/users/lstagner/events{/privacy}","received_events_url":"https://api.github.com/users/lstagner/received_events","type":"User","site_admin":false},"labels":[{"url":"https://api.github.com/repos/D3DEnergetic/FIDASIM/labels/enhancement","name":"enhancement","color":"84b6eb"}],"state":"closed","locked":false,"assignee":null,"milestone":null,"comments":1,"created_at":"2014-11-07T19:19:47Z","updated_at":"2016-04-04T08:25:41Z","closed_at":"2016-04-04T08:25:41Z","body":"We need to be able to propagate errors in the plasma profiles to the calculated spectra. The way to do this is to do the simulation with plasma profiles that are sampled according to their errors. You do this many times and this would build up a posterior distribution. It would take to too long to do this with the Monte Carlo routines but it can be done with the weight functions.\r\n\r\nThis project should be rather straight forward. I believe this would be a good project for an undergraduate to undertake. The following changes will have to be made\r\n\r\n1. Input error estimates for plasma profiles into FIDASIM\r\n2. Change colrad routine to take dne,te,ti ... arguments instead of grid cell index\r\n3. User input switches e.g. calc_error, num_samples, error_type...\r\n4. Routines to handle different error types e.g. Gaussian, Poisson, systematic....(Gaussian probably only one needed)\r\n5. Required changes to the FIDA and NPA weight function routines\r\n6. Output. Should the W.F. be outputted for every sample or just the mean and standard deviation? Perhaps have the output style be chosen by the user. Should we output the sampled parameters as well? \r\n7. How will we propagate the errors in the weight functions to the spectra? If the posterior for each E-p grid point is Gaussian then we can just use standard error propagation. If not then we will have to repeat the above procedure for the spectra calculation.\r\n\r\nOnce this is in place we can then do a detailed sensitivity study\r\n"}}