# NOT RUN {
#library(prozor)
data(protpepmetashort)
library(Matrix)
colnames(protpepmetashort)
head(protpepmetashort)
dim(protpepmetashort)
count = prepareMatrix( protpepmetashort, peptideID = "peptideSeq" )
dim(count)
inverse = prepareMatrix( protpepmetashort, peptideID = "peptideSeq" , weight = "inverse")
#aa = prepareMatrix(protpepmetashort, peptideID = "peptideSeq" , weight = "AA")
#xx = prepareMatrix(protpepmetashort, peptideID = "peptideSeq" , weight = "coverage")
image( as.matrix(count) )
corProt = cor( as.matrix(count) )
par(mfrow =c(1,2))
image(corProt)
#penalise peptides matching many proteins
corProtn = cor( as.matrix(inverse) )
image(corProtn)
# }
Run the code above in your browser using DataLab