Comments (4)
You are correct, this was supposed to include T2 rather than T1. This mistake has been fixed in the updated code. Thanks for taking the time to read through the models carefully enough to catch this issue!
from dadi_pipeline.
@hemstrow I just wanted to follow up on this and report that there were other issues I found with the original island model set, including mis-assigned and nonsense parameters. These models were incorrectly specified in addition to the issue you raised.
I discussed some better approaches to these models with Ryan (threads here and here). Afterwards, I completely re-wrote the original island model set. I released the revised island models as part of dadi_pipeline v3.1.1 in December 2019.
I'm not sure if you have seen these updates, but wanted to make you aware of them. If you plan to use the revised island models, please let me know if you find any other issues or discrepancies. These resources are a work in progress, and I always appreciate the feedback.
Thanks!
from dadi_pipeline.
Dear dportik
To continue this question, I want to build 1D model that my population experienced first bottleneck and then a founder effect, I draft the code, while it shows error:
code:
def three_epoch_founder(params, ns, pts):
"""
params = (nuB,nuF,nuFd,TB,TF,TFe)
ns = (n1,)
nuB: Ratio of bottleneck population size to ancient pop size
nuF: Ratio of contemporary to ancient pop size ##Ratio of founder to ancient
nuFd: Ratio of contemporary to founder
TB: Length of bottleneck (in units of 2*Na generations)
TF: Time since bottleneck recovery (in units of 2*Na generations)
TFe: Time of founder effect
n1: Number of samples in resulting Spectrum
pts: Number of grid points to use in integration.
"""
nuB,nuF,nuFd,TB,TF,TFe = params
xx = Numerics.default_grid(pts)
phi = PhiManip.phi_1D(xx)
phi = Integration.one_pop(phi, xx, TB, nuB)
phi = Integration.one_pop(phi, xx, TF, nuF)
nuFd_func = lambda t: nuF * (nuFd/nuF)**(t/TFe)
phi = Integration.one_pop(phi, xx, TFe, nuFd=nuFd_func)
fs = Spectrum.from_phi(phi, ns, (xx,))
return fs
three_epoch_founder.param_names = ['nuB', 'nuF','nuFd', 'TB', 'TF', 'TFe']
Looking forward for your reply, thankyou
Best
Xue
from dadi_pipeline.
Hi @xuefenfei712 ,
For questions related to model specification (and not the pipeline workflow), please post on the dadi user group:
https://groups.google.com/g/dadi-user
from dadi_pipeline.
Related Issues (20)
- Question about the simulations HOT 1
- Fixing theta HOT 1
- AIC HOT 4
- Optimization is too slow HOT 1
- Question about input file HOT 2
- does the pipeline incorporates correction for ancestral state misidentification? HOT 4
- Goodness of Fit test suggests model overfit HOT 2
- nan AIC on edited models HOT 3
- VCF2snp HOT 2
- How to transform the parameter to true parameter HOT 1
- about asymmetric migration and adjiacent asymmetric migration models HOT 2
- nu1,nu2,nu3 HOT 1
- pop1 or pop2 HOT 2
- dadi plot error HOT 1
- Summarize_Outputs.py HOT 6
- Folded spectrum vs. unfolded spectrum HOT 2
- Question about theta, Nref, and nuA HOT 4
- Default parameter boundaries HOT 2
- running time HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from dadi_pipeline.