sum_A_mat: Calculate the a parameter in the conditional Hurdle density...

In sqyu/ZiDAG: Directed Graphical Models and Causal Discovery for Zero-Inflated Data

Description

Calculate the a parameter in the conditional Hurdle density of the child node given the parents for multiple samples, assuming a is linear in the parents and their indicators.

Usage

sum_A_mat(aa, Vo, Yo)

Arguments

aa	A numerical vector of length 2*ncol(Yo)+1. See details.
Vo	A numerical vector of the same dimension as Yo indicating if each entry in Yo is non-zero, i.e. Vo = (Yo != 0).
Yo	A numerical vector, a sample for the parent nodes (regressors).

Details

Matrix version of sum_a(). See examples.

Value

A number, the mean log probability.

Examples

set.seed(1)
n <- 100; p_others <- 10
Vo <- matrix(stats::rbinom(n*p_others, 1, 0.8), nrow=n, ncol=p_others)
Yo <- matrix(stats::rnorm(n*p_others) * Vo, nrow=n, ncol=p_others)
aa <- rnorm(2*p_others+1)
sum_A_mat(aa, Vo, Yo) - sapply(1:n, function(i){sum_a(aa, Vo[i,], Yo[i,])})

sqyu/ZiDAG documentation built on Jan. 19, 2021, 4:11 p.m.