llyuleall {degreenet}

Calculate the log-likelihood for Count Distributions

Description

Functions to Estimate the Log-likelihood for Discrete Probability Distributions.

Usage

llyuleall(v, x, cutoff = 2, cutabove = 1000,  np=1)

Arguments

v

A vector of parameters for the Yule (a 1-vector - the scaling exponent).

x

A vector of counts (one per observation).

cutoff

Calculate estimates conditional on exceeding this value.

cutabove

Calculate estimates conditional on not exceeding this value.

np

wnumber of parameters in the model. For the Yule this is 1.

Value

the log-likelihood for the data x at parameter value v.

Note

See the papers on https://handcock.github.io/?q=Holland for details

References

Jones, J. H. and Handcock, M. S. "An assessment of preferential attachment as a mechanism for human sexual network formation," Proceedings of the Royal Society, B, 2003, 270, 1123-1128.

See Also

ayulemle, llyule, dyule, llwarall

Examples

# Simulate a Yule distribution over 100
# observations with rho=4.0

set.seed(1)
s4 <- simyule(n=100, rho=4)
table(s4)

#
# Calculate the MLE and an asymptotic confidence
# interval for rho
#

s4est <- ayulemle(s4)
s4est

# Calculate the MLE and an asymptotic confidence
# interval for rho under the Waring model (i.e., rho=4, p=2/3)
#

s4warest <- awarmle(s4)
s4warest

#
# Compare the AICC and BIC for the two models
#

llyuleall(v=s4est$theta,x=s4)
llwarall(v=s4warest$theta,x=s4)