estimate_state_probs: Estimate conditional/marginal state probabilities

Description Usage Arguments Value Examples

View source: R/estimate_state_probs.R

Description

Estimates P(S = s_k; \mathbf{W}), k = 1, …, K, the probability of being in state s_k using the weight matrix \mathbf{W}.

These probabilites can be marginal (P(S = s_k; \mathbf{W})) or conditional (P(S = s_k \mid \ell^{-}, \ell^{+}; \mathbf{W})), depending on the provided information (pdfs$PLC and pdfs$FLC).

Usage

1
2
estimate_state_probs(weight.matrix = NULL, states = NULL, pdfs = list(FLC = NULL, 
    PLC = NULL), num.states = NULL)

Arguments

weight.matrix

N \times K weight matrix

states

vector of length N with entry i being the label k = 1, …, K of PLC i

pdfs

a list with estimated pdfs for PLC and/or FLC evaluated at each PLC, i=1, …, N and/or FLC, i=1, …, N

num.states

number of states in total. If NULL (default) then it sets it to max(states) or ncol(weight.matrix) - depending on which one is provided.

Value

A vector of length K or a N \times K matrix.

Examples

1
2
3
WW <- matrix(runif(10000), ncol = 10)
WW <- normalize(WW)
estimate_state_probs(WW)

LICORS documentation built on May 29, 2017, 1:02 p.m.