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

PHPMotion implementation - URL based videos (Hosted on separate location) -

-
i have implement video collection website users of specific isp. site hosted on public web-server movies hosted on local isp server. user of isp + general visitors can see available movies; isp's users able play them. player on site point url of local machine ( http://192.x.x.x/movies/mymovie.mov ) private class (192.x.x.x) accessible isp user. question: possible phpmotion script? or need custom work. if possible, can please suggest available mods? phpmotion require modifications able handle this. uploading, conversion, thumbnail generation, , playback pretty static in phpmotion. really though, rather modding entire upload , conversion process, if created quick custom form add entries phpmotion database (maybe include thumbnail upload), tweaked playback paths, should go on phpmotion side of things. assume had videos in format player handle (such flv or mp4).
Read more

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

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