Identification of differential genes using t-tests by R

Created: 2008-04-15 15:37:42 Last updated: 2008-07-01 17:21:18

Download Workflow

This workflow starts by retrieving the names of microarray datasets from the Maxd database. The user has to select sets of control and test data for analysis using t-tests by R. A list of significant differentially expressed genes is then analysed using the Go Term Finder tool which generates a list of GO terms associated with the genes. A CSV file containing the list of significant genes is also generated.

Preview

Download as scalable diagram (SVG)

Run

Run this Workflow in the Taverna Workbench...

Option 1:

Copy and paste this link into File > 'Open workflow location...'
http://www.myexperiment.org/workflows/181/download?version=2
[ More Info Expand ]

Workflow Components

Inputs (0)

Processors (20)

Name	Type	Description
Ttest_pvalue	stringconstant
GoTermFinder_pvalue	stringconstant
getNamesUsingXPath	local	Parses the measurement names of microarray datasets from the XML output from maxdBrowse
flatten	local	Removes a layer of lists
performT_tests	rshell	This R script performs a series of t-tests for each gene expression level between its 2 input data sets.
mergeOutput1	beanshell	Merges the selected input data sets into a single data set
createFinalCSV	beanshell	Creates a text file in CSV format contain a list of the all the differentially expressed genes identified by the R script.
cleanCSV	beanshell	Removes unneccesary prefix text in the affy identifiers from the CSV file output generated by the R script
createGeneList	beanshell	Create a list of genes for analysis by the GOTermFinder service. Spurious genes are removed including those which have been annotated as questioable and hypothetical and therefore will not have a GO term asscoiated with them.
selectTestData	beanshell	Uses Java Swing classes to enable the user to select data sets for t-test analysis
getAffyIds	beanshell	Parses the affy identifiers from the CSV file
mergeOutput2	beanshell	Combines the selected data sets into a single data set
selectControlData	beanshell	Uses Java Swing classes to enable the user to select data sets for t-test analysis
getMeasurementNames	arbitrarywsdl	Querys the maxdload database using maxdBrowse to obtain a list of measurement names from the big experiment data stored there
queryMaxd1	arbitrarywsdl	Will take a list of Gene names, a list of Measurement names and a list of Column Types and produce a table. All values are coloured in relation to the a reference. In the list of Measurement names, the first Measurement listed is taken to be the reference. Formatted as text.
getGeneName	arbitrarywsdl	Gets the gene name for a given yeast open reading frame number
analyseGenesPDFOutput	arbitrarywsdl	Analyses a list of genenames using the GOTermFinder service
getTargetDescription	arbitrarywsdl	Gets a description of a given affy identifer using the YeastAffymapper service
getTranscriptId	arbitrarywsdl	Gets the open reading frame number for a given affy identifier
queryMaxd2	arbitrarywsdl	Will take a list of Gene names, a list of Measurement names and a list of Column Types and produce a table. All values are coloured in relation to the a reference. In the list of Measurement names, the first Measurement listed is taken to be the reference. Formatted as text.

Beanshells (8)

Name	Description	Inputs	Outputs
mergeOutput1	Merges the selected input data sets into a single data set	in	out
createFinalCSV	Creates a text file in CSV format contain a list of the all the differentially expressed genes identified by the R script.	affyIds orfs genenames descriptions	csv
cleanCSV	Removes unneccesary prefix text in the affy identifiers from the CSV file output generated by the R script	in	out
createGeneList	Create a list of genes for analysis by the GOTermFinder service. Spurious genes are removed including those which have been annotated as questioable and hypothetical and therefore will not have a GO term asscoiated with them.	listOfOrfs listOfGenenames listOfDescriptions	orfs spurious_genes
selectTestData	Uses Java Swing classes to enable the user to select data sets for t-test analysis	input	output
getAffyIds	Parses the affy identifiers from the CSV file	inList	outList
mergeOutput2	Combines the selected data sets into a single data set	in	out
selectControlData	Uses Java Swing classes to enable the user to select data sets for t-test analysis	input	output

Outputs (3)

Name	Description
pdf	A PDF document generated by the GoTermFinder service
csv	A text file in CSV format showing a list of differentially expressed genes and their annotations
spurious_genes	A text file in CSV format shoing a list of genes annotated as hypothetical, questionable, genomic sequences, etc. These genes do not have a GO term associated with them and so are removed from the list sent to the GoTermFinder service.

Links (28)

Source	Sink
GoTermFinder_pvalue:value	analyseGenesPDFOutput:pvalue
Ttest_pvalue:value	performT_tests:pvalue
cleanCSV:out	flatten:inputlist
createGeneList:orfs	analyseGenesPDFOutput:genes
createGeneList:spurious_genes	spurious_genes
flatten:outputlist	createFinalCSV:affyIds
flatten:outputlist	getAffyIds:inList
getAffyIds:outList	getGeneName:affyId
getAffyIds:outList	getTargetDescription:affyId
getAffyIds:outList	getTranscriptId:affyId
getGeneName:getGeneSymbolReturn	createFinalCSV:genenames
getGeneName:getGeneSymbolReturn	createGeneList:listOfGenenames
getMeasurementNames:return	getNamesUsingXPath:xml-text
getNamesUsingXPath:nodelist	selectControlData:input
getNamesUsingXPath:nodelist	selectTestData:input
getTargetDescription:getTargetDescriptionReturn	createFinalCSV:descriptions
getTargetDescription:getTargetDescriptionReturn	createGeneList:listOfDescriptions
getTranscriptId:getTranscriptIdReturn	createFinalCSV:orfs
getTranscriptId:getTranscriptIdReturn	createGeneList:listOfOrfs
mergeOutput1:out	performT_tests:control_csv
mergeOutput2:out	performT_tests:test_csv
performT_tests:file	cleanCSV:in
queryMaxd1:return	mergeOutput1:in
queryMaxd2:return	mergeOutput2:in
selectControlData:output	queryMaxd1:MeasurementNames
selectTestData:output	queryMaxd2:MeasurementNames
analyseGenesPDFOutput:analyseGenesPDFOutputReturn	pdf
createFinalCSV:csv	csv

Coordinations (1)

Controller	Target
mergeOutput1	selectTestData

Information Workflow Type

Taverna 1

Information Uploader

Peter Li

Information License

All versions of this Workflow are licensed under:

Information Version 2 (latest) (of 2)

Information Credits (1)

(People/Groups)

Peter Li

Information Attributions (0)

(Workflows/Files)

None

Information Tags (6)

Uploader tags

Log in to add Tags

Information Shared with Groups (0)

None

Information Featured In Packs (2)

Log in to add to one of your Packs

Information Attributed By (0)

(Workflows/Files)

None

Information Favourited By (2)

Information Statistics

4633 viewings

3252 downloads

[ see breakdown ]

Citations (1)

1. Li, P., Castrillo, J. I., Velarde, G., Wassink, I., Soiland-Reyes, S., Owen, S., Withers, D., Oinn, T., Pocock, M. R., Goble, C. A., Oliver, S. G. & Kell, D. B., Performing statistical analyses on quantitative data in Taverna workflows: an example using R and maxdBrowse to identify differentially expressed genes from microarray data, BMC Bioinformatics 2008, 9:334, 07 August 2008, http://www.biomedcentral.com/1471-2105/9/334

Version History

In chronological order:

Identification of differentially-expressed genes using R

Created by Peter Li on Tuesday 15 April 2008 15:37:42 (UTC)

Last edited by Peter Li on Tuesday 15 April 2008 15:41:25 (UTC)
Identification of differential genes using t-tests by R

Created by Peter Li on Tuesday 15 April 2008 15:37:42 (UTC)

Last edited by Peter Li on Tuesday 01 July 2008 17:21:18 (UTC)

Revision comment:

Replaced the cleanGenenames processor with a new Beanshell script called createGeneNames. This new processor generates a list of genes using their yeast gene name or open reading frame number if it does not have a standard yeast name.

Reviews (0)

No reviews yet

Be the first to review!

Comments (0)

View Timeline

No comments yet

Log in to make a comment

Other workflows that use similar services (3)

Only the first 2 workflows that use similar services are shown. View all workflows that use these services.

Taverna 1

Uploader

Peter Li

Identification of differential genes using... (1)

Download

This workflow starts by retrieving the names of microarray datasets from the Maxd database. The user has to select sets of control and test data which are then analysed by the LIMMA Bioconductor package in an R script. This produces a list of significant differentially expressed genes which is then analysed using the Go Term Finder tool to generate a PDF report of the common GO terms associated with the genes. A CSV file containing the list of sign...

Created: 2008-07-02 | Last updated: 2008-07-02

Credits: Peter Li

Taverna 1

Uploader

Peter Li

Mapping microarray data onto metabolic pat... (1)

Download

This workflow maps microarray data onto metabolic pathway diagrams represented as SBML models drawn using Cell Designer. To run this workflow requires libsbml to be installed into taverna - see http://www.mcisb.org/software/taverna/libsbml/index.html

Created: 2007-11-14 | Last updated: 2007-11-22

Credits: Peter Li