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); }