emvolz-phylodynamics / mlesky

MIT License
10 stars 2 forks source link

what is tau meaning here, thanks. #8

Closed liamxg closed 2 years ago

liamxg commented 2 years ago

tau = 18.4332061932401

liamxg commented 2 years ago
image

is this normal or not?

xavierdidelot commented 2 years ago

Maybe, it depends what your data looks like. It would be better to plot the effective population size on a log scale, and also to look at the growth rate.

liamxg commented 2 years ago
image
liamxg commented 2 years ago

plot(fit,logy=TRUE) like this, right?

liamxg commented 2 years ago
image
liamxg commented 2 years ago

could you please tell me how to look at the growth rate, thanks.

xavierdidelot commented 2 years ago

Yes logy=TRUE is indeed what you need to use to have the y-axis on a log-scale. To plot the growth rate you need to use growth=TRUE It looks like your analysis is completely uninformative.

liamxg commented 2 years ago

thanks for your help, like the following?

plot(fit,logy=TRUE,growth=TRUE) Warning message: In xy.coords(x, y, xlabel, ylabel, log) : 17 y values <= 0 omitted from logarithmic plot

liamxg commented 2 years ago
image
xavierdidelot commented 2 years ago

plot(fit,logy=FALSE,growth=TRUE) would be better, for the growth rate you usually don't want the logged y axis.

liamxg commented 2 years ago

thanks. like this?

image
liamxg commented 2 years ago

no 95% interval

xavierdidelot commented 2 years ago

Yes but this looks hard to believe. You should check your data is used correctly. Email me if you need help.

liamxg commented 2 years ago

thanks, is this your email (xavier.didelot@gmail.com), thanks.

xavierdidelot commented 2 years ago

Yes it is.

liamxg commented 2 years ago

thanks.