2b:null 2c:["$","title",null,{"children":"CAP 6737 - Assignment 06 Visual Cluster Analysis of Multidimensional data | Joshua Harkness | Observable"}] 31:T47f,{ const div = DOM.element('div'); const data = []; data.push({ x: irisSpeciesGrouped.get('setosa').map(d => d[selectDropDown.x]), y: irisSpeciesGrouped.get('setosa').map(d => d[selectDropDown.y]), mode: 'markers', marker: { size: 4, }, name: 'setosa' }); data.push({ x: irisSpeciesGrouped.get("versicolor").map(d => d[selectDropDown.x]), y: irisSpeciesGrouped.get("versicolor").map(d => d[selectDropDown.y]), mode: 'markers', marker: { size: 4, }, name: 'versicolor' }); data.push({ x: irisSpeciesGrouped.get("virginica").map(d => d[selectDropDown.x]), y: irisSpeciesGrouped.get("virginica").map(d => d[selectDropDown.y]), mode: 'markers', marker: { size: 4, }, name: 'virginica' }); const layout = { title: `Scatter Plot of Attribute ${selectDropDown.x} vs ${selectDropDown.y} `, xaxis: { autorange: true, title: `${selectDropDown.x}` }, yaxis: { autorange: true, title: `${selectDropDown.y}` } }; Plotly.newPlot(div, data, layout); return div; }32:Tae9,{ const myDiv = html`

`; const XY = iris.map(d => [d.sepal_length, d.petal_length]); const k = 3; //const { clusters, clusterCenters } = computeKmeanCluster(nClusters, xyData); let clusterCenters = initializeClusterCenters( k, XY, clusterIntitializationMethod ); let clusters; const colors = ["steelblue", "chocolate"]; // differnt color for each cluster for (let i = 0; i < nIterations; i++) { // Associate data points to cluster clusters = clusterCenters.map(_ => []); XY.forEach(d => { const distanceToClusterCenters = clusterCenters.map(c => squaredDistance(c, d) ); // Choose the closest center clusters[d3.minIndex(distanceToClusterCenters)].push(d); }); // Compute New Cluster center for (let j = 0; j < k; j++) clusterCenters[j] = updateClusterCenter(clusterCenters[j], clusters[j]); const data = clusters.map((cluster, i) => { // generate Trace data for the clusters return { x: cluster.map(d => d[0]), y: cluster.map(d => d[1]), mode: 'markers', marker: { size: 4, color: colors[i] } }; }); const layout = { showlegend: false, autosize: false, width: 800, height: 550, margin: { t: 50 }, hovermode: 'closest', xaxis: { title: 'sepal_length', showgrid: false, zeroline: false }, yaxis: { title: 'petal_length', showgrid: false, zeroline: false } }; layout.shapes = clusters.map((cluster, i) => { const convexHull = d3.polygonHull(cluster); if (clusterOutlineShape == "polygon") { const path = d3.line().curve(d3.curveCatmullRomClosed)(convexHull); return { type: 'path', path: path, fillcolor: 'rgba(128, 128, 128, 0.4)', line: { color: "black" } }; } else { const deltaX = 0.6 * d3.max(convexHull.map(p => Math.abs(p[0] - clusterCenters[i][0]))); const deltaY = 0.6 * d3.max(convexHull.map(p => Math.abs(p[1] - clusterCenters[i][1]))); return { type: 'circle', xref: 'x', yref: 'y', x0: clusterCenters[i][0] - deltaX, y0: clusterCenters[i][1] - deltaY, x1: clusterCenters[i][0] + deltaX, y1: clusterCenters[i][1] + deltaY, //opacity: 0.2, fillcolor: 'rgba(128, 128, 128, 0.4)', line: { color: 'black' } }; } }); Plotly.newPlot(myDiv, data, layout, { staticPlot: true }); await Promises.delay(1000); yield myDiv; } }33:Tb53,initializeClusterCenters = (k, XY, method) => { const N = XY.length; if (k == 0 || k > N) return null; const clusterCenters = XY.slice(0, k).slice(); switch (method) { case "Random": const dim = XY[0].length; for (let j = 0; j < k; j++) { for (let i = 0; i < dim; i++) clusterCenters[j][i] = d3.randomUniform( ...d3.extent(XY, d => d[i]) )(); } break; case "Forgy": // Choose k distinct items from XY array using Reservoir Sampler. // See algorithm: https://www.geeksforgeeks.org/reservoir-sampling1) // Step I: copy first k items of data to clusterCenters. for (let i = k; i < N; i++) { // a) Generate a random number from 0 to i where i is index of current item in data. // Let the generated random number is j. const j = d3.randomInt(i + 1)(); // b) If j is in range 0 to k-1, replace reservoir[j] with arr[i] if (j < k) clusterCenters[j] = XY[i].slice(); } break; case "Random Partition": // Randomly assign a partition to each element of XY. const clusterIndices = XY.map(_ => d3.randomInt(k)()); const lengths = [0, 0]; for (let j = 0; j < k; j++) { const cluster = clusterIndices .filter(index => index == j) .map(index => XY[index]); if (cluster.length > 0) clusterCenters[j] = updateClusterCenter(clusterCenters[j], cluster); } break; case "K-Mean++": //Randomly select the first cluster Center from the data points. clusterCenters[0] = XY[d3.randomInt(k)()].slice(); for (let j = 1; j < k; j++) { //For each data point compute its distance from the nearest, previously chosen centroid. const nearestDistances = XY.map(xy => d3.min( d3 .range(j) .map(index => Math.sqrt(squaredDistance(xy, clusterCenters[index])) ) ) ); const sumNearestDistances = d3.sum(nearestDistances); const probabilities = nearestDistances.map( d => d / sumNearestDistances ); const cumProbabilities = probabilities.slice(); for (let i = 1; i < N; i++) cumProbabilities[i] += cumProbabilities[i - 1]; // Select the next centers from the data points such that the probability of choosing a point as centroid is // directly proportional to its distance from the nearest, previously chosen centroid. // the point having maximum distance from the nearest centroid is most likely to be selected next as a centroid const x = Math.random(); const sampleIndex = d3.bisect(cumProbabilities, x); mutable feedback = sampleIndex; clusterCenters[j] = XY[sampleIndex].slice(); } break; } return clusterCenters; }34:[] 2f:{"id":"176c054c1c166f77","slug":"cap-6737-assignment-06-visual-cluster-analysis-of-multidim","title":"CAP 6737 - Assignment 06 Visual Cluster Analysis of Multidimensional data","thumbnail":"4f43e6b22498186b74904b75505ab12916b966f3ef8957dcb8a9d7b12a86ada6","access_level":"public","create_time":"$D2021-03-01T23:26:34.019Z","update_time":"$D2021-03-02T01:54:09.308Z","view_time":null,"first_public_time":"$D2021-03-02T00:38:46.936Z","has_unpublished_edits":null,"fork":false,"stars":0,"starred":false,"owner":{"type":"team","id":"73f33157ee13ba42","login":"darkhark","name":"Joshua Harkness","avatar_url":"https://avatars-next.observableusercontent.com/523d70e7f553a5ade6f8fb2ace128fcf1e933dcb4bef533e928a4452c0687855","tier":"starter_2024"},"version":88,"publish_version":88,"role":null,"authored":false,"authors":[{"type":"individual","id":"4514e759801885f0","login":"darkhark","name":"Joshua Harkness","avatar_url":"https://avatars-next.observableusercontent.com/523d70e7f553a5ade6f8fb2ace128fcf1e933dcb4bef533e928a4452c0687855","tier":null}],"members":[{"type":"individual","id":"4514e759801885f0","login":"darkhark","name":"Joshua Harkness","avatar_url":"https://avatars-next.observableusercontent.com/523d70e7f553a5ade6f8fb2ace128fcf1e933dcb4bef533e928a4452c0687855","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"titling":"auto","title":"CAP 6737 - Assignment 06 Visual Cluster Analysis of Multidimensional data","cells":[{"id":0,"value":"md`# CAP 6737 - Assignment 06 Visual Cluster Analysis of Multidimensional data`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":33,"value":"viewof selectDropDown = columns({\n x: select({\n options: irisDims,\n value: 'sepal_length'\n }),\n y: select({\n options: irisDims,\n value: 'sepal_width'\n }),\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":28,"value":"$31","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":69,"value":"viewof clusterIntitializationMethod = radio({\n options: [\"Random\",\"Forgy\", \"Random Partition\", \"K-Mean++\"],\n value: \"Forgy\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":67,"value":"viewof clusterOutlineShape = radio({\n options: [\"polygon\",\"ellipse\"],\n value: \"polygon\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":65,"value":"$32","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":12,"value":"md`## Imports`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":16,"value":"import { soFetch } from '@alecglassford/so-fetch'","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":19,"value":"d3 = require(\"d3@6\", \"d3-regression@1.3.4\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":30,"value":"Plotly = require(\"plotly.js-dist\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":35,"value":"import { radio, select, slider } from \"@jashkenas/inputs\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":37,"value":"import { columns } from \"@bcardiff/observable-columns\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":8,"value":"md`## Load the Data`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":10,"value":"iris = d3.csvParse(\n await (await soFetch(\n \"https://raw.githubusercontent.com/plotly/datasets/master/iris-id.csv\"\n )).text(),\n d3.autoType\n )","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":22,"value":"irisSpecies = [...new Set(iris.map(d => d.species))]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":25,"value":"irisDims = iris.columns.slice(0,4)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":42,"value":"irisSpeciesGrouped = d3.group(iris, d => d['species']);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":72,"value":"$33","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":74,"value":"updateClusterCenter = (center, cluster) => {\n for (let i = 0; i < center.length; i++)\n center[i] = d3.mean(cluster, d => d[i]);\n return center;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":76,"value":"squaredDistance = (p, q) =>\n d3.sum(p.map((_, i) => (p[i] - q[i]) * (p[i] - q[i])))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":80,"value":"mutable feedback = \"\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":83,"value":"nIterations = 5","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W34","stdlib":"1"}} 30:"$2f:body"