[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