Skip to content

Instantly share code, notes, and snippets.

@jacobw56
Last active October 13, 2017 19:09
Show Gist options
  • Save jacobw56/2fd529120462c8ee044bccc3b0836547 to your computer and use it in GitHub Desktop.
Save jacobw56/2fd529120462c8ee044bccc3b0836547 to your computer and use it in GitHub Desktop.
Gaussian curve in d3
<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8">
<meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1">
<title>Normal Plot</title>
<meta name="description" content="">
<script src="http://d3js.org/d3.v3.min.js" charset="utf-8"></script>
<style type="text/css">
body {
font: 10px sans-serif;
}
.axis path,
.axis line {
fill: none;
stroke: #000;
shape-rendering: crispEdges;
}
/*.x.axis path {
display: none;
}*/
.line {
fill: none;
stroke: lightsteelblue;
stroke-width: 12px;
}
.area {
fill: steelblue;
stroke-width: 0;
}
</style>
</head>
<body>
</body>
<script type="text/javascript">
//setting up empty data array
var data = [];
getData(); // popuate data
// line chart based on http://bl.ocks.org/mbostock/3883245
var margin = {
top: 20,
right: 20,
bottom: 30,
left: 50
},
width = 960 - margin.left - margin.right,
height = 960 - margin.top - margin.bottom;
var x = d3.scale.linear()
.range([0, width]);
var y = d3.scale.linear()
.range([(height*(2/3)), 6]);
var line = d3.svg.line()
.x(function(d) { return x(d.q); })
.y(function(d) { return y(d.p); });
var area = d3.svg.area()
.x(function(d) { return x(d.q); })
.y0(0)
.y1(function(d) { return y(d.p); });
var svg = d3.select("body").append("svg")
.attr("width", width + margin.left + margin.right)
.attr("height", height + margin.top + margin.bottom)
.append("g")
.attr("transform", "translate(" + margin.left + "," + margin.top + ")");
x.domain(d3.extent(data, function(d) {
return d.q;
}));
y.domain(d3.extent(data, function(d) {
return d.p;
}));
/*
svg.append("g")
.attr("class", "x axis")
.attr("transform", "translate(0," + height + ")")
.call(xAxis);
svg.append("g")
.attr("class", "y axis")
.call(yAxis);
*/
svg.append("path")
.datum(data)
.attr("class", "area")
.attr("d", area);
svg.append("path")
.datum(data)
.attr("class", "line")
.attr("d", line);
function getData() {
// loop to populate data array with
// probabily - quantile pairs
for (var i = 0; i < 1000; i++) {
q = normal() // calc random draw from normal dist
p = gaussian(q) // calc prob of rand draw
el = {
"q": q,
"p": p
}
data.push(el)
};
// need to sort for plotting
//https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Array/sort
data.sort(function(x, y) {
return x.q - y.q;
});
}
// from http://bl.ocks.org/mbostock/4349187
// Sample from a normal distribution with mean 0, stddev 1.
function normal() {
var x = 0;
var y = 0;
var rds, c;
do {
x = Math.random() * 2 - 1;
y = Math.random() * 2 - 1;
rds = x * x + y * y;
} while (rds == 0 || rds > 1);
c = Math.sqrt(-2 * Math.log(rds) / rds); // Box-Muller transform
return (Math.abs(x * c) > 3) ? normal() : x * c; // throw away extra sample y * c
}
//taken from Jason Davies science library
// https://github.com/jasondavies/science.js/
function gaussian(x) {
var gaussianConstant = 1 / Math.sqrt(2 * Math.PI);
var mean = 0;
var sigma = 1;
x = (x - mean) / sigma;
return gaussianConstant * Math.exp(-.5 * x * x) / sigma;
};
</script>
</html>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment