gmhta {puma} | R Documentation

## Compute gene and transcript expression values and standard deviatons from hta2.0 CEL Files

### Description

This function converts an object of `FeatureSet` into an object of class `exprReslt` using the gamma model for hta2.0 chips. This function obtains confidence of measures, standard deviation and 5, 25, 50, 75 and 95 percentiles, as well as the estimated expression levels.

### Usage

    gmhta(
          object
         ,background=FALSE
         ,gsnorm=c("median", "none", "mean", "meanlog")
         ,savepar=FALSE
         ,eps=1.0e-6
         ,addConstant = 0
         ,cl=NULL
         ,BatchFold=10
    )

### Arguments

`object` | an object of `FeatureSet`
`background` | Logical value. If `TRUE`, perform background correction before applying gmhta.
`gsnorm` | character. specifying the algorithm of global scaling normalisation.
`savepar` | Logical value. If `TRUE` the estimated parameters of the model are saved in file par\_gmhta.txt
`eps` | Optimisation termination criteria.
`addConstant` | numeric. This is an experimental feature and should not generally be changed from the default value.
`cl` | This function can be parallelised by setting parameter cl. For more details, please refer to the vignette.
`BatchFold` | we divide tasks into BatchFold*n jobs where n is the number of cluster nodes. The first n jobs are placed on the n nodes. When the first job is completed,the next job is placed on the available node. This continues until all jobs are completed. The default value is ten. The user also can change the value according to the number of cluster nodes n. We suggest that for bigger n BatchFold should be smaller.

### Details

The obtained expression measures are in log base 2 scale. Using the known relationships between genes, transcripts and probes, we propose a gamma model for hta2.0 data to calculate transcript and gene expression levels. The algorithms of global scaling normalisation can be one of "median", "none", "mean", "meanlog". "mean" and "meanlog" are mean-centered normalisation on raw scale and log scale respectively, and "median" is median-centered normalisation. "none" will result in no global scaling normalisation being applied. This function can be parallelised by setting parameter cl. For more details, please refer to the vignette.

### Value

A list of two object of class `exprReslt`.

### Author(s)

Xuejun Liu,WuJun Zhang,Zhenzhu gao, Magnus Rattray

### References

XueJun Liu. (2013) puma 3.0: improved uncertainty propagation methods for gene and transcript expression analysis, BMC Bioinformatics, 14:39.

Manhong Dai, Pinglang Wang,Andrew D. Boyd. (2005) Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data,Nucleic Acid Research 33(20):e175.

### See Also

Related class `exprReslt-class`

### Examples

    if(FALSE){
    ## The following scripts show the use of the method.
    #library(puma)
    ## load CEL files
    # object<-read.celfiles("celnames")
    #eset<-gmhta(object,gsnorm="none",cl=cl)

    }

[Package puma version 3.29.1 Index]