many thanks for help. Your code geos excellent! ]]>

Dear Ula, beta0 is the exponential parameter estimated by the model. If you look in older comments by Massub (and my answers) you will see that you are not the first experiencing this error. It may be due to the presence of NAs in your data or to some other problems. Try to follow those suggestions and see if they work. If they don’t we will have to investigate more in detail.

]]>I’m new in that issue.

What’s mean ‘beta0’ parameter. What can I get it?

I have an error :

Error in numericDeriv(form[[3L]], names(ind), env) :

Missing value or an infinity produced when evaluating the mode

thanx in advance

]]>Yes!

]]>Number of sample used of genotype : 189

Total number of chromosomes in my organism : 18

Total number of SNPs used for LD calculation : 16986

In this case ‘n’ for half decay calculation is 189*2=378.

Is all right Fabio ?

I hope it will help others in future.

Grazie Fabio.

]]>Sorry, I never manage to explain this correctly! If it is diploid, then it is 2*number_of_genotyped_individuals. So, for example, if you genotyped 50 diploid individuals, your n will be 2*50=100. I hope this clarifies the issue.

]]>Thanks for your quick reply. This organism is having 18 chromosomes and diploid.Then ‘n’ is 36. Is it ?

]]>None of them. By “number of chromosomes” I meant the number of each homologue chromosome. I change that to “number of alleles” which I hope is more immediate. So, if your organism is diploid and the number of individuals is “i” it will b n=2*i.

]]>This is again regarding ‘n’ (half-decay distance).

I have used ‘16986’ SNPs from ‘199’ scaffolds for LD calculation (it is a non model organism).In this case what will be the value of ‘n’ for half-decay distance calculation ?

n=16986*2 or n=199*2.

please could you clarify my request.

]]>