qch.fit: Infer Hconfig posteriors

In qch: Query Composed Hypotheses

Description

Infer Hconfig posteriors

Usage

qch.fit(pValMat, Hconfig, plotting = FALSE)

Arguments

pValMat	a matrix of p-values, each column corresponding to a p-value serie.
Hconfig	an Hconfig list as generated by the GetHinfo() function.
plotting	a boolean. Should some diagnostic graphs be plotted ? Default is FALSE.

Value

A list of 2 objects 'prior' and 'posterior'. Object 'prior' is a vector of estimated prior probabilities for each of the H-configurations. Object 'posterior' is a matrix providing for each item (in row) its posterior probability to belong to each of the H-configurations (in columns).

Examples

data(PvalSets)
PvalMat <- as.matrix(PvalSets[,-3])
## Build the Hconfig objects
Q <- 2
AtLeast <- 2
Hconfig <- GetHinfo(Q,AtLeast)$Hconfig

## Run the function
res.fit <- qch.fit(PvalMat,Hconfig)

## Display the prior of each class of items
res.fit$prior

## Display the first posteriors
head(res.fit$posterior)

qch documentation built on May 7, 2021, 5:06 p.m.