How to get gene expression value from Arrayexpress

ArrayExpress has a wonderful R package for data search, download and analysis but it doesn't always work in perfection. Therefore it always good to have an alternative to download the raw data/CEL and analyze it. This R script will simply take the accession number in command line argument and save the expression data in a file named as  Gene_expression.txt.
How to download expression data set from NCBI GEO HERE

Prerequisite 


We need the following R libraries to run the script
  • ArrayExpress 
  • aff

Uses 


Rscript script_name accession_number

Limitations 

  • This script is written specifically for Arabidopsis data sets, therefore, you have to modify it as per your requirement. 

Script 


How to perform parallel BLAST

Cheat Sheet to Install and work with R on Ubuntu