Published
Edited
May 18, 2022
1 fork
1 star
Attempt to add regression line:
longReadsB = Scatterplot(myFileTable, {
x: d => d[" A_fr"],
y: d => d[" qA"],
title: d => d.FileName,
xLabel: "A Frequency →",
yLabel: "↑ A Quality",
yDomain: [0,15],
stroke: "steelblue",
width: 450,
height: 450
}),
linearRegression(myFileTable, {x: " A_fr", y: " qA"})
import {Scatterplot} from "@d3/scatterplot"
import {Histogram} from "@d3/histogram"
import {LineChart} from "@d3/multi-line-chart"
d3Fetch = require('d3-fetch')
//data = FileAttachment(myFileTable).csv({typed: true})
function linearRegression(data, options = {}) {
let {stroke, x, y, z} = options;
let [vstroke, cstroke] = maybeColor(stroke, "currentColor");
if (z === undefined && vstroke != null) z = vstroke;
const X1 = [];
const Y1 = [];
const X2 = [];
const Y2 = [];
const S = vstroke ? [] : undefined;
return Plot.link(data, {
...Plot.transform(options, (data, facets) => {
const X = Plot.valueof(data, x);
const Y = Plot.valueof(data, y);
const Z = z !== undefined ? Plot.valueof(data, z) : undefined;
const [x1, x2] = d3.extent(X);
const regressionFacets = [];
for (const facet of facets) {
let F = facet.filter(i => defined(X[i]) && defined(Y[i]));
const regressionFacet = [];
for (const I of Z ? d3.group(F, i => Z[i]).values() : [F]) {
const f = linearRegressionLine(I, X, Y);
const i = I[0];
X1[i] = x1;
X2[i] = x2;
Y1[i] = f(x1);
Y2[i] = f(x2);
regressionFacet.push(i);
}
regressionFacets.push(regressionFacet);
}
return {data, facets: regressionFacets};
}),
x1: X1,
y1: Y1,
x2: X2,
y2: Y2
});
}
function linearRegressionLine(I, X, Y) {
const n = I.length;
if (n === 1) return () => Y[I[0]];
let sx = 0, sy = 0, sxx = 0, sxy = 0;
for (const i of I) {
const x = X[i];
const y = Y[i];
sx += x;
sy += y;
sxx += x * x;
sxy += x * y;
}
const m = (n * sxy - sx * sy) / (n * sxx - sx * sx);
const b = (sy - m * sx) / n;
return x => b + m * x;
}
function maybeColor(value, defaultValue) {
if (value === undefined) value = defaultValue;
return value === null ? [undefined, "none"]
: isColor(value) ? [undefined, value]
: [value, undefined];
}
function isColor(value) {
if (typeof value !== "string") return false;
value = value.toLowerCase().trim();
return value === "none"
|| value === "currentcolor"
|| (value.startsWith("url(") && value.endsWith(")")) // <funciri>, e.g. pattern or gradient
|| (value.startsWith("var(") && value.endsWith(")")) // CSS variable
|| d3.color(value) !== null;
}
function defined(x) {
return x != null && !Number.isNaN(x);
}
Purpose-built for displays of data
Observable is your go-to platform for exploring data and creating expressive data visualizations. Use reactive JavaScript notebooks for prototyping and a collaborative canvas for visual data exploration and dashboard creation.
