linux - Fitting multiple curves to one data set -


i have data set receive outside source, , have no real control over. can edit ad infinitum, cannot choose parameters filter set. data, when plotted, shows 2 clumps of points several sparse, irrelevant points, well. here sample plot:

raw plot

(sorry, don't have adequate privileges embed picture)

there clump of points on left, clustered around (1, 16). clump part of set of points lies on (or near) line stretching (1, 17.5) (2.4, 13). also, there apparent curve (1.75, 18) (2.75, 12.5). finally, there sparse points above second curve, around (2.5, 17). visually, it's not difficult separate these groups of points. however, need separate these points within data file 3 groups, i'll call line, curve, , other (the curve group 1 need). i'd write program can reasonably without needing visually see plot.

now, i'm going add couple items make worse. sample set of data. while shapes of curve , line relatively constant 1 data set next, positions not. these regions can (and do) shift, both horizontally , vertically. real constant there's negative-slope line top-left bottom-right of plot, ~$\frac{1}{x^4}$ curve top-center bottom-right, , of sparse points in top-right corner, above curve.

i'm on linux, , i'm out of ideas. can tell approaches i've tried, though have not done well.

  • first, cleaned data set , sorted in ascending order x-coordinate. thought maybe points sorted in sort of logical way allow me 'head' or 'tail' data achieve desired result, not case.
  • i can write code in (python, fortran, c, etc.) removes point if it's not within x distance of previous point. fine, except scattering of points such 2 points near each other in x, separated appreciable distance in y. doesn't line , curve draw near 1 larger x-values.
  • i can fit curve partial data set. when sort data x-coordinate, example, can choose plot first 30 points, or last 200, or set of 40 in middle somewhere. that's not problem. line points tuck underneath curve points, causes problem.
  • if line points constant (which they're not), rotate plot angle line vertical , can @ points right of line, rotate back. may best way go doing this, in order that, need able isolate linear points, more or less essence of problem.

the other idea seems plausible me, try identify point density , split data separate files parameters. think best candidate problem, since independent of point location. however, i'm not sure how go doing this, because line , curve come quite close larger x-values (in sample plot, it's x-values greater 2).

i know not fall in request of mwe, don't know how i'd go providing more classical mwe. if there's else can provide help, please ask. thank in advance.


Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more