# compareICL: Compare corrected ICL values after data transformation In coseq: Co-Expression Analysis of Sequencing Data

## Description

Compare the corrected ICL values after applying the arcsin, logit, and logMedianRef transformations in a coseq analysis

## Usage

 `1` ```compareICL(x) ```

## Arguments

 `x` A list made up of `coseq` objects. At the current time, this function only supports the comparison of `coseq` objects using `model="Normal"` and `transformation = c("arcsin", "logit", "logMedianRef")`

## Value

A plot of corrected ICL values for the models included in `x` (the list of `coseq` objects)

## Author(s)

Andrea Rau, Cathy Maugis-Rabusseau

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19``` ```## Simulate toy data, n = 300 observations set.seed(12345) countmat <- matrix(runif(300*4, min=0, max=500), nrow=300, ncol=4) countmat <- countmat[which(rowSums(countmat) > 0),] conds <- rep(c("A","B","C","D"), each=2) ## Run the Normal mixture model for K = 2,3,4 run_arcsin <- coseq(y=countmat, K=2:4, iter=5, transformation="arcsin") ## Plot and summarize results plot(run_arcsin) summary(run_arcsin) ## Compare ARI values for all models (no plot generated here) ARI <- compareARI(run_arcsin, plot=FALSE) ## Compare ICL values for models with arcsin and logit transformations run_logit <- coseq(y=countmat, K=2:4, iter=5, transformation="logit") compareICL(list(run_arcsin, run_logit)) ```

coseq documentation built on May 31, 2017, 2:38 a.m.