Learn R Programming

BANFF (version 1.1)

HyperPara.Select: Selecting Hyper Parameters by Bayesian Model Averaging

Description

This function is aiming at selecting appropriate hyper parameters by Bayesian Model Averaging. This process of this function has already embedded in function Networks.STD() and Networks.Fast(). The purpose of this function is for having the selected hyper-parameters independently.

Usage

HyperPara.Select(net,pvalue,piall,rhoall,n=30)

Arguments

net
an "n" by "n" binary adjacent matrix (0/1) for the network configuration with n=length(pvalue)
pvalue
a vector of p-values obtained from large scale statistical hypothesis testing
piall
a vector of possible choices for "pi0" in an increasing order
rhoall
a vector of possible choices of "rho0" and "rho1" in an increasing order
n
a number of iterations you set for Bayesian Model Averaging. The default setting is 30, which is accord with the embedded inner process of in function Networks.STD() and Networks.Fast(). More iterations could be set if you have a desire to have more accurate parameter

Value

a list of selected parameter
pi0
the selected pi0
rho0
the selected rho0
rho1
the selected rho1

References

Yize Zhao, Jian Kang, Tianwei Yu (2014) A Bayesian nonparametric model for selecting gene and gene sub network, Annals of Applied Statistics, in press.

Zhou Lan, Jian Kang, Tianwei Yu, Yize Zhao, BANFF: an R package for network identifications via Bayesian nonparametric mixture models, working paper.