Simulates Supcluster Function

Share:

Description

Produces summary statistics from a simulation of supcluster

Usage

1
2
3
4
tab1(ratio=4,reps=100,n=1000,start=500,fbeta=FALSE,
               maxclusters=5,chains=1,clusts=c(15,15,20),
               sig=1,gamma=1,npats=80,beta=seq(-5,5,5),
               plot=FALSE)

Arguments

ratio

The ratio of tau to sigma

reps

The number of runs

n

The number of MCMC iterations

start

The first MCMC iteration used

fbeta

If TRUE the outcome is not used

maxclusters

The maximum number of clusters for the estimation step

chains

The number of chains to run

clusts

A list of the number of genes in each cluster

sig,gamma,beta

The parameters sigma,gamma,beta

npats

The number of experimental units(patients)

plot

Plots the first run

Value

A data frame is returned with the mean parameter value, it's standard error and the mean of it's standard error calculated from the MCMC

Author(s)

David A. Schoenfeld, Jessie Hsu

See Also

supcluster,,compare.chains,concordmap

Examples

1
2
3
#very few iterations done so that this runs in less than 5 seconds.
#You need to change reps=100,start=2000,n=3000 to get enough iterations
tab1(ratio=2,reps=5,n=10,start=1,maxclusters=5)