search

fgvieira / ngsF

Estimation of per-individual inbreeding coefficients under a probabilistic framework
Other
20 stars 8 forks source link

When should I set --Freq_fixed=TRUE in calculating inbreeding Coefficient? #21

Open wul88 opened 1 year ago

wul88 commented 1 year ago

Hi, Filipe,

I am using ngsF to estimate the inbreeding coefficient for a population of giraffes. I am not sure about the parameter "--freq_fixed". The readMe stated "assume initial MAF as fixed parameters (only estimates F)". However, the examples shown in the ngsTools tutorials didn't mention this parameter in either initial values estimate (with --approx_EM) or the final run.

I tried running ngsF with and without setting "--freq_fixed=TRUE". The results are quite different. I would like to learn more about this parameter.

fgvieira commented 1 year ago

This option is only useful in few cases, when you don't want to estimate allele frequencies from your data, usually when you already have the freqs from another set of individuals of the same population.

wul88 commented 1 year ago

Thank you very much for responding quickly!!

We try to assess the inbreeding situation of this population of giraffes. So it is important we understand the method very well. I went through the original (2013) paper. The allele frequency and inbreeding coefficient per individual per site are both estimated in each iteration. If using "--freq_fixed=TRUE", does that mean the allele frequency per site is fixed and taken out of the iteration updating?

Thank you again!

fgvieira commented 1 year ago

Yep! With --freq_fixed=TRUE per site freqs are fixed, and only per individual inbreeding coef are estimated.