R/psi_1.R

In rodd: Optimal Discriminating Designs

Psi <- function(x, p, eta, theta.fix, theta.var)
{
    res <- 0
    for(i in 1:length(eta))
        for(j in 1:length(eta))
            if(p[i,j] != 0)
                res <- res + p[i,j] * psi(x, eta[[i]], eta[[j]], theta.fix[[i]], theta.var[[i,j]])
    res
}