{"action":"created","issue":{"url":"https://api.github.com/repos/tylerhether/hmmspore/issues/16","labels_url":"https://api.github.com/repos/tylerhether/hmmspore/issues/16/labels{/name}","comments_url":"https://api.github.com/repos/tylerhether/hmmspore/issues/16/comments","events_url":"https://api.github.com/repos/tylerhether/hmmspore/issues/16/events","html_url":"https://github.com/tylerhether/hmmspore/issues/16","id":55424438,"number":16,"title":"Maximum likelihood estimation for p.assign and p.trans parameters needed","user":{"login":"tylerhether","id":5805450,"avatar_url":"https://avatars.githubusercontent.com/u/5805450?v=3","gravatar_id":"","url":"https://api.github.com/users/tylerhether","html_url":"https://github.com/tylerhether","followers_url":"https://api.github.com/users/tylerhether/followers","following_url":"https://api.github.com/users/tylerhether/following{/other_user}","gists_url":"https://api.github.com/users/tylerhether/gists{/gist_id}","starred_url":"https://api.github.com/users/tylerhether/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/tylerhether/subscriptions","organizations_url":"https://api.github.com/users/tylerhether/orgs","repos_url":"https://api.github.com/users/tylerhether/repos","events_url":"https://api.github.com/users/tylerhether/events{/privacy}","received_events_url":"https://api.github.com/users/tylerhether/received_events","type":"User","site_admin":false},"labels":[{"url":"https://api.github.com/repos/tylerhether/hmmspore/labels/question","name":"question","color":"cc317c"}],"state":"closed","locked":false,"assignee":{"login":"mwpennell","id":1786485,"avatar_url":"https://avatars.githubusercontent.com/u/1786485?v=3","gravatar_id":"","url":"https://api.github.com/users/mwpennell","html_url":"https://github.com/mwpennell","followers_url":"https://api.github.com/users/mwpennell/followers","following_url":"https://api.github.com/users/mwpennell/following{/other_user}","gists_url":"https://api.github.com/users/mwpennell/gists{/gist_id}","starred_url":"https://api.github.com/users/mwpennell/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/mwpennell/subscriptions","organizations_url":"https://api.github.com/users/mwpennell/orgs","repos_url":"https://api.github.com/users/mwpennell/repos","events_url":"https://api.github.com/users/mwpennell/events{/privacy}","received_events_url":"https://api.github.com/users/mwpennell/received_events","type":"User","site_admin":false},"milestone":null,"comments":4,"created_at":"2015-01-25T18:33:43Z","updated_at":"2015-05-22T05:19:49Z","closed_at":"2015-05-22T05:19:49Z","body":"We talked earlier about creating a (better) method for estimating jointly the two parameters of the forward-backward function.  Systematically going through parameter spaces can be really inefficient.  What methods do you have in mind? "},"comment":{"url":"https://api.github.com/repos/tylerhether/hmmspore/issues/comments/104513134","html_url":"https://github.com/tylerhether/hmmspore/issues/16#issuecomment-104513134","issue_url":"https://api.github.com/repos/tylerhether/hmmspore/issues/16","id":104513134,"user":{"login":"tylerhether","id":5805450,"avatar_url":"https://avatars.githubusercontent.com/u/5805450?v=3","gravatar_id":"","url":"https://api.github.com/users/tylerhether","html_url":"https://github.com/tylerhether","followers_url":"https://api.github.com/users/tylerhether/followers","following_url":"https://api.github.com/users/tylerhether/following{/other_user}","gists_url":"https://api.github.com/users/tylerhether/gists{/gist_id}","starred_url":"https://api.github.com/users/tylerhether/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/tylerhether/subscriptions","organizations_url":"https://api.github.com/users/tylerhether/orgs","repos_url":"https://api.github.com/users/tylerhether/repos","events_url":"https://api.github.com/users/tylerhether/events{/privacy}","received_events_url":"https://api.github.com/users/tylerhether/received_events","type":"User","site_admin":false},"created_at":"2015-05-22T05:19:49Z","updated_at":"2015-05-22T05:19:49Z","body":"Code has been implemented. The C++ variant is between 50 and 100 times faster than the R variant when comparing 100 runs using microbenchmark. The results are the same save the R variant provides 1 more significant digit than the C++ variant. \r\n\r\nDuring the updating we switched the algorithms (and simulations) so that p.trans is no longer used. Instead of p.trans (i.e., the per base pair transition probability) we have the \"scale\" parameter. scale is used to rescale the physical distance between snps to mapping distance. The mapping distance is used during the generation of the hmm transition matrix during both the forward and the backward algorithms; it is used in the haldane function. \r\n\r\nFor both the R variant and the C++ variant of the FB algorithm, the maximum lnL can be estimated for each parameter combination by summing up the lnL for each tetrad, spores, and chromosomes and picking the parameter combination with the maximum lnL. "}}