Roleswitch (version 1.10.0)

diagnosticPlot: Create diagnostic plot for understanding the Roleswitch model outputs.

Description

Create diagnostic plot for understanding the roleswitch outputs. Create a 2 by 4 panels of plots. From left to right, the top panel displays the observed N mRNA and M miRNA expression (x.o and z.o), the N by M seed-match matrix (c), and the inferred total mRNA expression; the bottom panel displays the inferred probability of the M miRNAs targeting the N mRNA (miRNA-mRNA; p.x), the probability of the N mRNA "targeting" the M miRNAs (mRNA-miRNA; p.z), the dot product of the above two matrices (Joint) and the convergence rate (delta.p.all).

Usage

diagnosticPlot(pred)

Arguments

pred

Results obtained from roleswitch.

References

Li, Y., ..., Zhang, Z., Inference of personalized miRNA-mRNA interactions toward redefining cancer signatures (in preparation).

See Also

roleswitch

Examples

x.o <- matrix(abs(rnorm(10, mean=3)), 
  dimnames=list(c(1:10),"mRNA")) # mRNA expression
  
z.o <- matrix(abs(rnorm(4, mean=3)),
  dimnames=list(c(1:4),"miRNA"))  # miRNA expression
  
c <- matrix(rpois(40, lambda=3),nrow=nrow(x.o), 
  dimnames=list(c(1:10),c(1:4))) # seed match matrix

rs.pred <- roleswitch(x.o, z.o, c)

diagnosticPlot(rs.pred)