Last data update: 2014.03.03

R: Select Ranked Lists
selectrankedlistsR Documentation

Select Ranked Lists

Description

Given a list with two elements, one containing ranklists for a set of regression models and the second containing the associated p-values. This function is used to extract a subset of the models.

Usage

selectrankedlists(ranklist, colsinc)

Arguments

ranklist

List with two elements (as output from generate profiles): Ranklist: Matrix containing the ranks of gene expression. Rows containing the genes, columns the different profiles.

Pvalues: Matrix containing the associated p-values to the differential expression profiles in Ranklist.

colsinc

Vector of integers containing the column references of the profiles to select.

Details

Format of list provided to selectrankedlists is the same as is output from generateprofiles generateprofiles. The output from selectrankedlists can be used as input to the classify profile function classifyprofile.

Value

Ranklist

Matrix containing the ranks of gene expression. Rows containing the genes, columns the selected profiles

Pvalues

Matrix containing the associated p-values to the differential expression profiles in Ranklist

Author(s)

C. Pacini

Source

http://www.ebi.ac.uk/arrayexpress/experiments/E-GEOD-22528

See Also

generateprofiles, classifyprofile.

Examples

data(profiles)
selectprofile<-selectrankedlists(profiles,1)
classification<-classifyprofile(data=selectprofile$ranklist,noperm=10,signif.fdr=1)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(DrugVsDisease)
Loading required package: affy
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colnames, do.call, duplicated, eval, evalq,
    get, grep, grepl, intersect, is.unsorted, lapply, lengths, mapply,
    match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, rank,
    rbind, rownames, sapply, setdiff, sort, table, tapply, union,
    unique, unsplit

Loading required package: Biobase
Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.

Loading required package: limma

Attaching package: 'limma'

The following object is masked from 'package:BiocGenerics':

    plotMA

Loading required package: biomaRt
Loading required package: ArrayExpress
Loading required package: GEOquery
Setting options('download.file.method.GEOquery'='auto')
Setting options('GEOquery.inmemory.gpl'=FALSE)
Loading required package: DrugVsDiseasedata
Loading required package: cMap2data

Attaching package: 'cMap2data'

The following object is masked from 'package:DrugVsDiseasedata':

    genelist

Loading required package: qvalue
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/DrugVsDisease/selectrankedlists.Rd_%03d_medium.png", width=480, height=480)
> ### Name: selectrankedlists
> ### Title: Select Ranked Lists
> ### Aliases: selectrankedlists
> ### Keywords: ~kwd1 ~kwd2
> 
> ### ** Examples
> 
> data(profiles)
> selectprofile<-selectrankedlists(profiles,1)
> classification<-classifyprofile(data=selectprofile$ranklist,noperm=10,signif.fdr=1)
Number of Significant results greater than 10 Using top 10 hits - consider using average linkage instead
> 
> 
> 
> 
> dev.off()
null device 
          1 
>