How to draw a polygon around NA values in R? -

-

i'm trying draw error region around data using base graphics. i've figured out how polygon, starts act badly if there na values in data.

dat <- rnorm(10, mean = 1:10) depth <- 11:20 sd <- rnorm(10, mean = 1.5, sd = 0.5) col <- "blue" alpha <- .2  col <- adjustcolor(col, alpha.f = alpha) par(mfrow = c(1,2)) plot(dat, depth, type = "o", main = "no nas in dat or depth") x <- c(dat - sd, rev(dat + sd)) y <- c(depth, rev(depth)) polygon(x = x, y = y, col = col, border = na)  dat[7] <- na plot(dat, depth, type = "o", main = "nas in dat or depth") x <- c(dat - sd, rev(dat + sd)) y <- c(depth, rev(depth))  polygon(x = x, y = y, col = col, border = na)

this gives me following image: the desired outcome if there no nas (left panel) , error when there nas (right panel)

it seems na value divides lower polygon 2 polygons. i'd keep 1 polygon.

here's possible solution using rle function:

set.seed(123) # added reproducibility dat <- rnorm(10, mean = 1:10) depth <- 11:20 sd <- rnorm(10, mean = 1.5, sd = 0.5) col <- "blue" alpha <- .2  col <- adjustcolor(col, alpha.f = alpha) par(mfrow = c(1,2)) plot(dat, depth, type = "o", main = "no nas in dat or depth") x <- c(dat - sd, rev(dat + sd)) y <- c(depth, rev(depth)) polygon(x = x, y = y, col = col, border = na)  dat[7] <- na plot(dat, depth, type = "o", main = "nas in dat or depth") x <- c(dat - sd, rev(dat + sd)) y <- c(depth, rev(depth))  ############################################ ## code print error range starts here: ## ############################################ enc <- rle(!is.na(dat)) endidxs <- cumsum(enc$lengths) for(i in 1:length(enc$lengths)){   if(enc$values[i]){     endidx <- endidxs[i]     startidx <- endidx - enc$lengths[i] + 1      subdat <- dat[startidx:endidx]     subsd <- sd[startidx:endidx]     subdepth <- depth[startidx:endidx]      x <- c(subdat - subsd, rev(subdat + subsd))     y <- c(subdepth, rev(subdepth))      polygon(x = x, y = y, col = col, border = na)   } }

enter image description here

the idea plot 1 polygon each consecutive non-na block.

since rle, given vector, returns lengths , values of consecutive blocks having same value, use identify blocks not na , plot polygon subsetting original dat, depth , sd vectors.


Comments

Post a Comment

Popular posts from this blog

javascript - Chart.js (Radar Chart) different scaleLineColor for each scaleLine -

-
Image
i try , create radar chart using chart.js has various colours each scaleline, or coloured between scalelines. wondering if possible? from: to: i have working graph, though there doesn't seem method change individual scale lines. kind regards leigh you can extend radar chart type this, so chart.types.radar.extend({ name: "radaralt", initialize: function (data) { chart.types.radar.prototype.initialize.apply(this, arguments); var originalscaledraw = this.scale.draw; var ctx = this.chart.ctx; this.scale.draw = function () { var linewidth = this.linewidth; // bypasses line drawing in originalscaledraw this.linewidth = linewidth; originalscaledraw.apply(this, arguments); ctx.linewidth = this.linewidth; var scale = this; // draw chart.helpers.each(scale.ylabels, function (label, index) { // color of ...
Read more

apache - Error with PHP mail(): Multiple or malformed newlines found in additional_header -

-
suddenly have started receiving above error without changes having been made script. host 1and1 (i know...) the script still works fine on different server, , suspicion there must have been server config change has lead this, although hosts plead ignorance. there's no information on above error @ in google can find - have ideas? server running apache if helps. had similar problem. came out of blue. no php code changed. what changed: php upgraded 5.5.25-1 5.5.26. a security risk in php mail() function has been fixed , newlines in additional_headers allowed no more. because newlines mean: starts email message (and surely don't want inject newlines through headers followed evil message). what have worked fine, e.g. having newlines after headers or passing whole message additional_headers , function no more. solution : sanitize headers. no multiple newlines in additional_headers argument. these count "multiple or malformed newlines":...
Read more

java - Android – MapFragment overlay button shadow, just like MyLocation button -

-
Image
i'm struggling add shadow custom overlay button on google map. goal make custom button google's mylocation button. here's how app looks now: layout.xml <relativelayout xmlns:android="http://schemas.android.com/apk/res/android" xmlns:tools="http://schemas.android.com/tools" xmlns:map="http://schemas.android.com/apk/res-auto" android:layout_height="match_parent" android:layout_width="match_parent" style="@style/apptheme"> <fragment android:layout_width="match_parent" android:layout_height="match_parent" android:id="@+id/map" tools:context=".mapsactivity" android:name="com.google.android.gms.maps.supportmapfragment" android:layout_alignparentend="true" android:layout_alignparentstart="true" android:layout_alignparentbottom="true" ...
Read more