[BioC] Looking for R script to categorize genes using GO terms

Morten Mattingsdal morten.mattingsdal at student.uib.no
Thu Nov 3 10:52:49 CET 2005


Hi...
I have a script (which you can see is written by a biologist :) here 
:http://folk.uio.no/mortema/master.R

or code (depends limma object "MA" to work)

md<-get("GO:0007517", GOBPCHILDREN)
myGO<-append(md,"GO:0007517")
for (i in myGO){
 probe<-unique(lookUp(i, "mgug4121a", "GO2ALLPROBES"))
}
res<-NULL
for (i in probe){
 res<-append(res,grep(i,MA$gene$ID))
}
myM<-matrix(nrow=(length(res)),ncol=12) 
jens = 1
for(i in res){
 myM[jens,] = MA$M[i,]
 jens <- jens+1
}
name=NULL
for (i in res){
name<-append(name,MA$gene$GeneName[i])
jens<-jens+1
}
row.names(myM)<-name
write.table(myM,file="muscle development.txt",sep="\t",row.names=T) #expressions values for genes involved in mucle development


>Hi Dear list,
>
>I am wondering whether someone has written some R script to categorize
>genes into different Go categories (e.g., replication, cell wall,...)
>according to gene AffyID? I am working on ATH1 array. Or any hint/ quick
>way to do it?
>
>Thanks in advance
>Fangxin
>
>
>--------------------
>Fangxin Hong  Ph.D.
>Plant Biology Laboratory
>The Salk Institute
>10010 N. Torrey Pines Rd.
>La Jolla, CA 92037
>E-mail: fhong at salk.edu
>(Phone): 858-453-4100 ext 1105
>
>_______________________________________________
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch
>https://stat.ethz.ch/mailman/listinfo/bioconductor
>  
>



More information about the Bioconductor mailing list