negllc_expo: Function to calculate the loglikelihood of a exponential...

In xiangfstats/GenieR: Inference of demographic history from phylogeny

Description

Function to calculate the loglikelihood of a exponential coalescent model

Usage

negllc_expo(parr, t, A)

Arguments

parr	A parameter vector for population size and growth rate.
t	A sampling and coalescent time vector.
A	A lineages vector

Value

loglikelihood of a exponential coalescent model.

Author(s)

Fei Xiang ([e-mail])

Examples

library(ape)
t1=rcoal(20)
x=att(t1)
negllc_expo(c(1,1), x$t, x$A)
library(minqa)
bobyqa(c(1,2),negllc_expo,lower=0,upper=Inf,t=x$t,A=x$A)
Geniefit(t1,Model="expo",start=c(10,2),upper=Inf,lower=0)

xiangfstats/GenieR documentation built on May 4, 2019, 1:06 p.m.