modello (version 0.1.1)

ldbeta.modello_number: Beta Distribution: log-density

Description

Considering a Beta Distribution calculates the log-density $$ld = (\alpha_1 - 1) * log(y) + (\alpha_2 - 1) * log(1 - y) - logBeta(\alpha_1, \alpha_2); y \in (0, 1)$$

Usage

# S3 method for modello_number
ldbeta(y, a1, a2)
# S3 method for default
ldbeta(y, a1, a2)
ldbeta(y, a1, a2)

Arguments

y

observations, numeric or reference object of class 'number'

a1

shape parameter, numeric of reference object of class 'number'

a2

shape parameter, numeric of reference object of class 'number'

Value

Returns a 'numeric' or a reference object of class 'number'

Examples

# NOT RUN {
modello.init(10, 10, 10, 10)
## For modello_numbers
y = number(rbeta(10, 1, 2))
ld = ldbeta(y, .k(1), .k(2))
print(ld)
print(ld$v)
modello.close()
## For numerics
y = rbeta(10, 1, 2)
ld = ldbeta(y, 1, 2)
print(ld) 
# }