Sub-grouping in a plot in R -


i want represent data different individual cells in xyplot , give them color in base of different category. however, when have been able represent dots:

xyplot( signal ~ time | as.factor(treatment), data=data,groups=cell,         fill.color = as.character(data$color),         panel = function(x, y,fill.color,...,subscripts){           fill = fill.color [subscripts]           panel.xyplot(x, y,pch=19, col=fill, type ="p")} )

but in way imposible visually track cells. therefore, want same lines , polygons error area happens each time try, lattice override second group should assign color (or not able tell lattice. here way that: 

cell<-rep(x = c("a","b","c","d"),50) signal<-rep(sample(seq(from = 0, = 50, = 1), size = 50, replace = true),4) time<-sort(rep(seq(1,50),4),decreasing = f) treatment<-rep(c("hard","soft"),50*2) color<-rep(c("red","orange"),50*2) data<-data.frame(cell,signal,time,treatment,color)  my.panel2 <- function(x, y, subscripts, col, pch,cex,sd,fill.color,...) {   low95 <- y-sd[subscripts]   up95 <- y+sd[subscripts]   fill=fill.color [subscripts]   panel.xyplot(x, y, col=fill.color, pch=pch,cex=cex, ...)   panel.arrows(x, low95, x, up95, angle=90, code=3,lwd=3,                 length=0.05, alpha=0.2,col=col) }   xyplot(signal~time|as.factor(treatment), groups=as.factor((data$cell)),         data=data, type='l',        color.line=as.character((data$color)))

thanks; santi

thanks @vince, gave me interesting idea. ggplot2 based on layers decided try similar in lattice as.layer(). first spitted data in 3 groups, 1 per level 

fm1<-filter(fmeans, group=="ave-int") fm2<-filter(fmeans, group=="mini-int") fm3<-filter(fmeans, group=="high-int") line1<-fm2[43,] fm3<-rbind(fm3,line1)

then decided use polygon error band using panel.polygon():

my.panel.pol <- function(x, y, subscripts, col,sd,...) {   plot.line <- trellis.par.get("plot.line")   xs <- if(is.factor(x)) {     factor(c(levels(x) , rev(levels(x))), levels=levels(x))   } else {     xx <- sort(unique(x))     c(xx, rev(xx))   }   low95 <- y-sd[subscripts]   up95 <- y+sd[subscripts]   panel.xyplot(x, y, col=col,...)   panel.polygon(xs, c(up95, rev(low95)), col=col, alpha=0.2, border=f) }

then represent data , merge in 1 graphic

   a<-       xyplot(mean ~ slice*12 |treatment , fm1,               layout=c(2,2),col="red",              grid=t,              ylim = c(0,max(fmeans$mean)),               group = stemcell, type = "l",              sd=fm1$sd,              panel.groups= "my.panel.pol",              panel="panel.superpose")     b<-xyplot(mean ~ slice*12 |treatment , fm2,               layout=c(2,2),col="blue",               group = stemcell, type = "l",               sd=fm2$sd,               panel.groups= "my.panel.pol",               panel="panel.superpose")     c<-xyplot(mean ~ slice*12 |treatment , fm3,               layout=c(2,2),col="green",               group = stemcell, type = "l",               sd=fm3$sd,               panel.groups= "my.panel.pol",               panel="panel.superpose")      a+as.layer(b)+as.layer(c)

so here result:

enter image description here


-

Comments

Post a Comment

Popular posts from this blog

c# - Binding a comma separated list to a List<int> in asp.net web api -

the default binder in web api expecst http://url.com/webapi/report/?pageids=3243&pageids=2365 to bind public ihttpactionresult report(list<int> pageids){ // exciting webapi code} i wish bind http://url.com/webapi/report/?pageids=3243,2365 as running out of space in url get . i have created class public class commaseparatedmodelbinder : system.web.http.modelbinding.imodelbinder { public bool bindmodel(httpactioncontext actioncontext, modelbindingcontext bindingcontext) { //binding in here } } and registered in webapiconfig.cs var provider = new simplemodelbinderprovider( typeof(list<int>), new commaseparatedmodelbinder()); config.services.insert(typeof(modelbinderprovider), 0, provider); i have altered method signature use model binder so public ihttpactionresult report( [modelbinder] list<int> pageids){ // exciting webapi code} however break point in binder not being
Read more

Delphi 7 and decode UTF-8 base64 -

in delphi 7, have widestring encoded base64(that received web service widestring result) : pd94bwwgdmvyc2lvbj0ims4wij8+dqo8c3ryaw5nptiq2lpyqjwvc3ryaw5npg== when decoded it, result not utf-8: <?xml version="1.0"?> <string>طھط³طھ</string> but when decoded base64decode.org, result true : <?xml version="1.0"?> <string>تست</string> i have use encddecd unit decodestring function. the problem have using decodestring . function, in delphi 7, treats decoded binary data being ansi encoded. , problem text utf-8 encoded. to continue encddecd unit have couple of options. can switch decodestream . instance, code produce utf-8 encoded text file data: {$apptype console} uses classes, encddecd; const data = 'pd94bwwgdmvyc2lvbj0ims4wij8+dqo8c3ryaw5nptiq2lpyqjwvc3ryaw5npg=='; var input: tstringstream; output: tfilestream; begin input := tstringstream.create(data); try output := tfilestream.
Read more

html - Is there any way to exclude a single element from the style? (Bootstrap) -

so, have img logo, overlaid on div styled background-image: url("...") , , both elements inside 2 div classes .to-color , .to-overlay respectively, black transparent overlay. is there way exclude logo style given these 2 classes? here's jsfiddle example: link html , css example: - html - <div class="to-color"> <div class="to-overlay"> <div class="respo"> <div class="row"> <div class="col-xs-4 col-xs-offset-4 text-center"> <img src="http://ansonalex.com/wp-content/uploads/2011/06/google-logo-768x260.png" alt="..." class="img-responsive"> </div> </div> </div> </div> </div> - css - .to-color { background-color: #000; } .to-overlay { opacity : 0.5; } .respo { background-image: url("https://source.unsplash.com/category/objects/1600x900"); background-position: center center;
Read more