extract data from txt files

[sunnycqcn]

I want to extract gene expression values from a txt file based on gene names.

for example:

A 1500

B 2500

C 3500

D 3400

I want to extract A and C gene expression values using gene name file extract it.

A 1500

C 3500

I have wroten a perl code, but it is not work well. Could you help me modify it or give me a new code?

Thanks,

#!/usr/bin/perl -w
use strict;
open(IN,"file1");
open(NAME,"file2");
open(OUT,">result");
while(<IN>){
chomp;
@data=split / /,$_;
$hash{$data[0]}=$data[1];
}
while(<NAME>){
chomp;
print OUT "$_ $hash{$_}\n";
}
close IN;
close NAME;
close OUT;

[Michael.Ante]

Why not simply grep:

grep -f name.txt gene_expression.txt > out.txt

P.S.: Please be aware of the fact, that Cross-Postings are forbidden/not welcome..

[sunnycqcn]

Thanks, you are right.

Why not simply grep:

grep -f name.txt gene_expression.txt > out.txt

P.S.: Please be aware of the fact, that Cross-Postings are forbidden/not welcome..

[ashishbansal]

After RMA processing, you should end up with an ExpressionSet object. Let's assume it's called eset, and your probe set IDs are in a character vector called yourProbeSetIDs. Then:

exprs(eset) # gives a dataframe with expression levels, probe set IDs are the row names
exprs(eset)[yourProbeSetIDs, ] # gives you of subset dataframe with your probe sets of interest
To write a tab-delimited file:

write.table(yourDataFrame, file="filename.txt", sep="\t", quote=F, row.names=T)

[Urmila785]

I'm looking for something like that, Thank you.
Online Clinical Research Training