Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Q: please explain the layout of the model shown in .model #81

Open
IanCodes opened this issue Dec 15, 2020 · 1 comment
Open

Q: please explain the layout of the model shown in .model #81

IanCodes opened this issue Dec 15, 2020 · 1 comment

Comments

@IanCodes
Copy link

IanCodes commented Dec 15, 2020

How does the explanation of how decide whether the model (in .model) is correct relates to the statement in the readme.

"The open state should have the highest emission parameters for all 4 signal components, the nucleosome state should have the second highest emission parameters for all 4 signal components and the background state should have the lowest emission parameters for all 4 signal components."

State 0
  Pi: 0.3333333333333333
  Aij: 0.962 0.031 0.007
  Opdf: Multi-variate Gaussian distribution --- Mean: [ 0 0 0 0 ]

State 1
  Pi: 0.3333333333333333
  Aij: 0.045 0.939 0.015
  Opdf: Multi-variate Gaussian distribution --- Mean: [ 0.075 0.125 0.009 0 ]

State 2
  Pi: 0.3333333333333333
  Aij: 0.005 0.007 0.988
  Opdf: Multi-variate Gaussian distribution --- Mean: [ 0.134 0.217 0.188 0.097 ]

How does 'open', 'nuclear' and 'background' states relate to the above model? Should the 'aij' values or 'mean' values be used?
Thank you.

@EvanTarbell
Copy link
Collaborator

State 0 is the background State, State 1 is the nucleosome state and state 2 is the open state (assuming default -k value)
Pi is the starting probability, ie the chance that the sequence starts with State i
Aij is the transition probability of moving from state i to state j
And the Opdf means are the mean values for the separated signal components. So for this model, the signal values for the open state are:
Mean: [ 0.134 0.217 0.188 0.097 ]
which satisfies the recommendation

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants