Log-likelihood values for the NB estimator
This will return return the log-likelihood value given a value of N.
NB.likelihood(N, infile, alleles, sample.interval)
The effective population size dor diploid individuals.
Your input data file in a plain text format. This data contains the allele counts at each locus and from each sampling time point.
a vector containing the number of alleles at each locus. For example, c(4, 4, 4) would mean that 3 loci are sampled, with 4 alleles each.
a vector stating at which generations the samples were taken. For example, c(0, 8) would indicate that two samples were collected from the 0th and 8th generation.
More details please see
This function returns one single element, the log-likelihood of the model given the effective population size N and your inputs.
This function allows you to use your own optimisation algorithms or customise the parameters using
nlm. Otherwise please use
NB.estimator which has the internal
optim wrapped inside.
1 2 3 4 5 6 7 8 9 10 11
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.