negllc_plog: Function to calculate the loglikelihood of a piecewise...

Description Usage Arguments Value Author(s) Examples

View source: R/RcppExports.R

Description

Function to calculate the loglikelihood of a piecewise logistic coalescent model

Usage

1
negllc_plog(parr, t, A)

Arguments

parr

A parameter vector for population size, exponential growth rate and change time.

t

A sampling and coalescent time vector.

A

A lineages vector

Value

loglikelihood of a piecewise logistic coalescent model.

Author(s)

Fei Xiang ([email protected])

Examples

1
2
3
4
5
6
7
8
9
library(ape)
t1=rcoal(20)
x=att(t1)
negllc_plog(c(1,1,.1), x$t, x$A)
library(minqa)
bobyqa(c(1,2,1),negllc_plog,lower=0,upper=Inf,t=x$t,A=x$A)
Geniefit(t1,Model="plog",start=c(10,2,1),upper=Inf,lower=0)
library(dfoptim)
nmkb(c(1,2,10),negllc_plog,lower=0,upper=Inf,t=x$t,A=x$A)

xiangfstats/GenieR documentation built on June 14, 2018, 2:44 a.m.