r - ScatterPlot and ONLY one Histogram plot together -

-

i want visualize time series data 'scatter plot' , histogram on right side, haven't been able figure out how turn off histogram on upper side.

code example:

install.packages("psych") library(psych) data = matrix(rnorm(n=100000,mean=2,sd=1.5), nrow = 100, ncol=1000) fs = list() fs$p_z = 1*(data>2) n_p = 1;   for(i in floor(seq(1,dim(data)[2],length.out=n_p))) {   scatter.hist(x = rep(1:length(data[,i])), y = data[,i],                xlab = 'observations',                ylab = 'log(tpm)',                title = 'mixture plot',                col = c("red","blue")[fs$p_z[,i]+1],                 correl = false, ellipse = false, smooth = false) }

result:

scatter plot 2 histograms, 1 upper side , 1 on right side

expected result: same 1 have no histogram on upper side. i.e., histogram on right side log(tpm).

note: using psych package, scatter.hist function seemed easy , nice use, couldn't find how turn off 1 histogram.

where flexibility ends, hacking begins. if @ scatter.hist function, see pretty basic usage of r base graphics. following modified code create plot want:

scat.hist <- function(x, y, xlab = null, ylab = null, title = "", ...) {  ## create layout layout(matrix(c(1,2),1,2,byrow=true), c(3,1), c(1,3))  ## plot scatter par(mar=c(5,5,3,1)) plot(x= x, y = y, xlab = xlab, ylab = ylab, main = title, ...)  ## right histogram yhist <- hist(y, plot = false, breaks = 11) par(mar=c(5,2,3,1)) mp <- barplot(yhist$density, space=0, horiz=true, axes = false) ## density d <- density(y, na.rm = true, bw = "nrd", adjust = 1.2) temp <- d$y         d$y <- (mp[length(mp)] - mp[1] + 1) * (d$x - min(yhist$breaks))/(max(yhist$breaks) - min(yhist$breaks))         d$x <- temp lines(d) }

let's try first row:

i = 1 scat.hist(x = seq_along(data[,i]), y = data[,i], col = c("red", "blue")[fs$p_z[,i]+1], xlab = 'observations', ylab = 'log(tpm)', title = 'mixture plot')

enter image description here


Comments

Post a Comment

Popular posts from this blog

c# - Unity IoC Lifetime per HttpRequest for UserStore -

-
i'm trying clean default implementation of accountcontroller.cs comes out of box in new mvc5/owin security implementation. have modified constructor this: private usermanager<applicationuser> usermanager; public accountcontroller(usermanager<applicationuser> usermanager) { this.usermanager = usermanager; } also, have created lifetime manager unity looks this: public class httpcontextlifetimemanager<t> : lifetimemanager, idisposable { private httpcontextbase _context = null; public httpcontextlifetimemanager() { _context = new httpcontextwrapper(httpcontext.current); } public httpcontextlifetimemanager(httpcontextbase context) { if (context == null) throw new argumentnullexception("context"); _context = context; } public void dispose() { this.removevalue(); } public override object get...
Read more

Change the color of an oval at click in Java AWT -

-
i have draw ovals in java, , @ click change color. beginning tried change color after 20 ms, doesn't work. my code is: public class mycomponentnew extends frame { public graphics2d g2d; public mycomponentnew(string title) { super(title); setsize(400, 550); } @override public void paint(graphics g) { this.g2d = (graphics2d) g; this.g2d.setcolor(color.red); this.g2d.filloval(10, 55, 50, 100); } public void changecolor () { this.g2d.setcolor(color.blue); this.g2d.filloval(10, 55, 50, 100); } } and in class main method have: mycomponentnew m; m = new mycomponentnew("fereastra cu baloane"); m.setvisible(true); m.addwindowlistener(new windowadapter() { @override public void windowclosing(windowevent we) { system.exit(0); } }); try { thread.sleep(20); } catch(interruptedexception e) {} m.changecolor(); the color of oval remains red. ...
Read more

I am trying to solve the error message 'incompatible ranks 0 and 1 in assignment' in a fortran 95 program. -

-
i writing program in fortran find velocity of parachuting person in relation time. keep getting error can't fix. new programming , appreciated. the error is v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) 1 error: incompatible ranks 0 , 1 in assignment @ (1) and program is program para integer :: real :: v(11) !velocity real :: q !initial velocity real :: h !time step real :: c !drag coefficient real :: m !mass ! gravity equal 32 ft/s^2 write (*,*)'enter time step' read(*,*)h write(*,*)'enter initial velocity' read(*,*)q write(*,*)'enter drag coefficient' read(*,*)c write(*,*)'enter mass' read(*,*)m i=1,10 ! 1 10, 1 being interval. end v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) q=v(1) end program you cannot use [] normal parenthesis in expressions. array constructor, [ items ] means array items elements. end do should after line.
Read more