Mendel's Genetics#16
Conversation
|
Once the build has completed, you can preview your PR at this URL: https://biojulia.dev/BiojuliaDocs/previews/PR16/ |
Just noting that the comment is being made, but the link doesn't actually work. Probably unrelated to the above, your pull request is for some reason requesting to merge into another branch, rather than into |
kescobo
left a comment
There was a problem hiding this comment.
Another solution would be to use StatsBase.jl and do a weighted probability.
One other thing that would be nice to include here is a bit more didactic discussion about how often times we make algorithms that are narrowly tailored, but then we either repeat ourselves or get more complicated as additional requirements get tacked on. Eg, for this problem, your solution works for the specific problem, but we'd have to derive a new equation if the question is something like "What's the probability of a heterozygous offspring?" It also doesn't scale up if we add another trait etc.
Nice thing about the StatsBase.jl solution and even a simulation is that they can be made generic and then can be used to ask more types of questions. I'm not necessarily demanding we add this to a first draft, but maybe open an issue as a potential enhancement.
|
I like the idea of a simulation, though it will generally not give a precisely correct answer for rosalind. I think that's fine if that's explained. |
| !!! warning "The Problem" | ||
|
|
||
| Probability is the mathematical study of randomly occurring phenomena. | ||
|
|
There was a problem hiding this comment.
These newlines in between lines need to be removed, else we'll end up with paragraph breaks in the rendered version. But these splits look good
|
|
||
| ```julia | ||
| function mendel(k,m,n) | ||
|
|
|
|
||
| This solution is generic and can be used to ask more types of questions. | ||
|
|
||
| <br> |
There was a problem hiding this comment.
| <br> |
Not needed


Making a draft PR here. There's multiple ways to solve the problem, and I added a first approach. I'm thinking that the second would be a more statistical/simulation approach. Basically, based on the values of k, m, n, we can make a vector containing all of the possible organisms (eg. [HH, Hh, hh, HH, etc.]). Then, we can calculate the percentage of dominant individuals/total individuals.
Wanted to run this by you first and see if you had any suggestions on packages to use.