You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
When I try running the Lphy scripts seems always appear warnings like "WARNING: branch length used to get equilibrium distribution was not long enough!"
The script I use:
data {
L = 1000;
}
model {
N0 ~ LogNormal(meanlog=8,sdlog=0.5);
f0 ~ Beta(alpha=40, beta=2);
b ~ LogNormal(meanlog=-0.95,sdlog=0.2);
gompertzPopFunc = gompertz(N0=N0,f0=f0, b=b);
When I try running the Lphy scripts seems always appear warnings like "WARNING: branch length used to get equilibrium distribution was not long enough!"
The script I use:
data {
L = 1000;
}
model {
N0 ~ LogNormal(meanlog=8,sdlog=0.5);
f0 ~ Beta(alpha=40, beta=2);
b ~ LogNormal(meanlog=-0.95,sdlog=0.2);
gompertzPopFunc = gompertz(N0=N0,f0=f0, b=b);
gompertzPopFunc= gompertz(N0=N0,f0=f0, b=b);
ψ ~ CoalescentPopFunc(n=40, popFunc=gompertzPopFunc);
π ~ Dirichlet(conc=[30.0,3.0,3.0,3.0,3.0,30.0,3.0,3.0,3.0,3.0,30.0,3.0,3.0,3.0,3.0,30.0]);
rates ~ Dirichlet(conc=[1.0, 2.0, 1.0, 1.0, 2.0, 1.0]);
Q = gt16(freq=π, rates=rates); // construct the GT16 instantaneous rate matrix
A ~ PhyloCTMC(L=1000, Q=Q, tree=ψ, dataType=phasedGenotype(), mu=0.01);
// epsilon ~ Beta(alpha=1, beta=50);
// delta ~ Beta(alpha=1, beta=2);
// E ~ GT16ErrorModel(delta=delta, epsilon=epsilon, alignment=A);
D = unphase(A);
}
The text was updated successfully, but these errors were encountered: