1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
#' @importFrom graphics mtext par plot.default rect segments text title
#' @rdname plotMat
#'
#' @param main.title Main title in \code{plotMatNm} version.
#' @param title.row Title for row normalized version in \code{plotMatNm}.
#' @param title.col Title for column normalized version in \code{plotMatNm}.
#' @param title.col Title for column normalized version in \code{plotMatNm}.
#' @param main.title.line Used instead of \code{title.line} in \code{plotMatNm}.
#' @param par.set Used instead of \code{title.line} in \code{plotMatNm}.
#'
#' @export
plotMatNm <-
function(M=x,x=M,...,main.title=NULL,title.row="Row normalized",title.col="Column normalized",main.title.line=-2,par.set=list(mfrow=c(1,2))){
if(is.null(main.title)){
objName<-deparse(substitute(M))
if(objName=="x")objName<-deparse(substitute(x))
main.title <- paste("Matrix",objName)
}
if(!is.null(par)){
par.def<-par(no.readonly = TRUE)
par(par.set)
}
row.normalized<-sweep(M, 1, apply(M, 1, sum),FUN="/")
row.normalized[is.nan(row.normalized)]<-0
plot.mat(M=row.normalized,main=title.row,outer.title=FALSE,...)
column.normalized<-sweep(M, 2, apply(M, 2, sum),FUN="/")
column.normalized[is.nan(column.normalized)]<-0
plot.mat(M=column.normalized,main=title.col,outer.title=FALSE,...)
title(main=main.title,outer=TRUE,line=main.title.line)
if(!is.null(par.set))par(par.def)
}
|
