using library pd.hg.u95av2
1
0
Entering edit mode
@pau-marc-munoz-torres-6731
Last seen 7.5 years ago
Barcelona

Hello everybody.

 

 I'm traying to use some cel files that i download from internet. In the paper that I found that the hibridation probe set in the experiment was Affymetrix Hum95Av2. Now, I would like to relate the experiment result with the gene annotaton. To this objective i looked in bioconductor fro the correct database. I found three that i thought were correct, hgu95av2 , u95av2 and pd.hg.u95av2.

 

 To relate de data with the database i used the following commands (I tested it with all the databases) :

library('affy')
library ("hgu95av2")
setwd("/home/paumarc/Kopriva/cel/")
data <- ReadAffy()
my_frame <- data.frame(exprs(data))

Annot <- data.frame(ACCNUM=sapply(contents(hgu95av2ACCNUM), paste, collapse=","),SYMBOL=sapply(contents(hgu95av2SYMBOL),paste,collapse=","),DESC=sapply(contents(hgu95av2.GENENAME), paste,collapse=","),DESC=sapply(contents(hgu95av2UNIPROT), paste, collapse=", "))
all <- merge(Annot, my_frame, by.x=0, by.y=0, all=T)
write.csv(all, file = "hgu95av2.csv")

 

Never the less it did not work. Can you tell me what do I do wrong or if I should use another database?

thanks.

 

affy databases • 1.1k views
ADD COMMENT
0
Entering edit mode
@pau-marc-munoz-torres-6731
Last seen 7.5 years ago
Barcelona

problem solved

 

 i forgot to normalitze data and i changed to the database hgu95av2.. Now it works

 

thanks

ADD COMMENT

Login before adding your answer.

Traffic: 492 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6