Analyze data with the Statistics software R - Revision history

A WASSERMANN at 12:55, 18 October 2024

2024-10-18T12:55:52Z

← Older revision		Revision as of 12:55, 18 October 2024
Line 145:		Line 145:
	* [http://en.wikipedia.org/wiki/Robust_statistics Wikipedia on Robust Statistics]		* [http://en.wikipedia.org/wiki/Robust_statistics Wikipedia on Robust Statistics]

	[[Category:~~Examples]]~~		[[Category:Old]]
	~~[[Category:Statistics~~]]

Michael at 10:35, 21 February 2013

2013-02-21T10:35:29Z

← Older revision		Revision as of 10:35, 21 February 2013
Line 92:		Line 92:
	var est2 = a[5]*1.0;		var est2 = a[5]*1.0;

	~~console.log(t, res, m, sd, med, mad, est1, est2);~~
	if (!graph2) {		if (!graph2) {
	graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});		graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});

Michael at 10:35, 21 February 2013

2013-02-21T10:35:06Z

← Older revision		Revision as of 10:35, 21 February 2013
Line 115:		Line 115:
	}		}

	$('output').innerHTML = '<b><font size="+1">Normal location and scale:</font></b><br /><br />'+		document.getElementById('output').innerHTML = '<b><font size="+1">Normal location and scale:</font></b><br /><br />'+
	'<b>Estimates for location (true value = 10):</b><br />' +		'<b>Estimates for location (true value = 10):</b><br />' +
	'Mean = ' + ~~brd~~.round(m,2) + '<br />' +		'Mean = ' + Math.round(m,2) + '<br />' +
	'Median = ' + ~~brd~~.round(med,2) + '<br />' +		'Median = ' + Math.round(med,2) + '<br />' +
	'Radius-minimax estimator = ' + ~~brd~~.round(est1,2) + '<br /><br />' +		'Radius-minimax estimator = ' + Math.round(est1,2) + '<br /><br />' +
	'<b>Estimates for scale (true value = 3):</b><br />' +		'<b>Estimates for scale (true value = 3):</b><br />' +
	'Standard deviation = ' + ~~brd~~.round(sd,2) + '<br />' +		'Standard deviation = ' + Math.round(sd,2) + '<br />' +
	'MAD = ' + ~~brd~~.round(mad,2) + '<br />' +		'MAD = ' + Math.round(mad,2) + '<br />' +
	'Radius-minimax estimator = ' + ~~brd~~.round(est2,2) + '<br />';		'Radius-minimax estimator = ' + Math.round(est2,2) + '<br />';
	brd.update();		brd.update();
	};		};

Michael at 10:33, 21 February 2013

2013-02-21T10:33:57Z

← Older revision		Revision as of 10:33, 21 February 2013
Line 88:		Line 88:
	var sd = a[1]*1.0;		var sd = a[1]*1.0;
	var med = a[2]*1.0;		var med = a[2]*1.0;
	~~var mad = a[3]*1.0;~~
	var mad = a[3]*1.0;		var mad = a[3]*1.0;
	var est1 = a[4]*1.0;		var est1 = a[4]*1.0;
	var est2 = a[5]*1.0;		var est2 = a[5]*1.0;

	console.log(t);		console.log(t, res, m, sd, med, mad, est1, est2);
	if (!graph2) {		if (!graph2) {
	graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});		graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});

Michael at 10:32, 21 February 2013

2013-02-21T10:32:56Z

← Older revision		Revision as of 10:32, 21 February 2013
Line 92:		Line 92:
	var est1 = a[4]*1.0;		var est1 = a[4]*1.0;
	var est2 = a[5]*1.0;		var est2 = a[5]*1.0;

			console.log(t);
	if (!graph2) {		if (!graph2) {
	graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});		graph2 = brd.createElement('curve', [[x[0],x[x.length-1]],[m,m]], {strokecolor:'red'});

Michael at 18:03, 19 July 2012

2012-07-19T18:03:28Z

← Older revision		Revision as of 18:03, 19 July 2012
Line 76:		Line 76:
	t += p[i].Y() + ';';		t += p[i].Y() + ';';
	}		}
	new Ajax.Request('/~alfred/jsxgraph~~/trunk~~/examples/rserv.php', {		new Ajax.Request('/~alfred/jsxgraph/examples/rserv.php', {
	method:'post',		method:'post',
	parameters:'input='+escape(t),		parameters:'input='+escape(t),

Michael at 10:39, 17 January 2012

2012-01-17T10:39:52Z

← Older revision		Revision as of 10:39, 17 January 2012
Line 76:		Line 76:
	t += p[i].Y() + ';';		t += p[i].Y() + ';';
	}		}
	new Ajax.Request('/~~~mkohl~~/rserv.php', {		new Ajax.Request('/~alfred/jsxgraph/trunk/examples/rserv.php', {
	method:'post',		method:'post',
	parameters:'input='+escape(t),		parameters:'input='+escape(t),

Michael at 09:16, 7 June 2011

2011-06-07T09:16:13Z

← Older revision		Revision as of 09:16, 7 June 2011
Line 33:		Line 33:
	y[6+i] = 10+z2*3;		y[6+i] = 10+z2*3;
	x[6+i] = 6+i;		x[6+i] = 6+i;
	p[i] = brd.createElement('point', [x[i],y[i]],{name:' ',fixed:true,~~style~~:7});		p[i] = brd.createElement('point', [x[i],y[i]],{name:' ',fixed:true,size:2,face:'[]'});
	p[6+i] = brd.createElement('point', [x[6+i],y[6+i]],{name:' ',fixed:true,~~style~~:7});		p[6+i] = brd.createElement('point', [x[6+i],y[6+i]],{name:' ',fixed:true,size:2,face:'[]'});
	}else{		}else{
	y[i] = 10+z1*3;		y[i] = 10+z1*3;
	p[i] = brd.createElement('point', [x[i],y[i]],{name:' ',fixed:true,~~style~~:7});		p[i] = brd.createElement('point', [x[i],y[i]],{name:' ',fixed:true,size:2,face:'[]'});
	}		}
	}		}
Line 61:		Line 61:
	brd.removeObject(p[9]);		brd.removeObject(p[9]);

	p[0] = brd.createElement('glider', [x[0],y[0],l[0]],{name:' ',~~style~~:6});		p[0] = brd.createElement('glider', [x[0],y[0],l[0]],{name:' ',size:4,face:'o'});
	p[10] = brd.createElement('glider', [x[10],y[10],l[1]],{name:' ',~~style~~:6});		p[10] = brd.createElement('glider', [x[10],y[10],l[1]],{name:' ',size:4,face:'o'});
	p[8] = brd.createElement('glider', [x[8],y[8],l[2]],{name:' ',~~style~~:6});		p[8] = brd.createElement('glider', [x[8],y[8],l[2]],{name:' ',size:4,face:'o'});
	p[9] = brd.createElement('glider', [x[9],y[9],l[3]],{name:' ',~~style~~:6});		p[9] = brd.createElement('glider', [x[9],y[9],l[3]],{name:' ',size:4,face:'o'});
	brd.unsuspendUpdate();		brd.unsuspendUpdate();

Line 127:		Line 127:
	}});		}});
	}		}
			brd.addHook(doIt, 'mouseup');
	</jsxgraph>		</jsxgraph>

Michael at 09:11, 7 June 2011

2011-06-07T09:11:32Z

← Older revision		Revision as of 09:11, 7 June 2011
Line 8:		Line 8:
	* radius-minimax estimator: green (optimally robust; cf. Rieder et al. (2008))<br /><br />		* radius-minimax estimator: green (optimally robust; cf. Rieder et al. (2008))<br /><br />
	By changing the y-position of the four movable points you should recognize the instability (non-robustness) of mean and standard deviation in contrast to the robust estimates; e.g., move one of the four movable points to the top of the plot.<br /><br />		By changing the y-position of the four movable points you should recognize the instability (non-robustness) of mean and standard deviation in contrast to the robust estimates; e.g., move one of the four movable points to the top of the plot.<br /><br />
	~~<html>~~
	~~<link rel="stylesheet" type="text/css" href="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraph.css" />~~
	~~<script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/prototype.js"></script>~~
	~~<script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraphcore.js"></script>~~
	~~<div id="jxgbox" class="jxgbox" style="width:700px; height:400px;" onmouseup="doIt();"></div>~~
	~~</html>~~
	===Online results:===		===Online results:===
			<html><script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/prototype.js"></script></html>
	<pre id='output'>Statistics:<br></pre>		<pre id='output'>Statistics:<br></pre>
	<~~html>~~		<jsxgraph width="700" height="400">
	~~<script language~~="~~JavaScript~~">		brd = JXG.JSXGraph.initBoard('jxgbox', {boundingbox: [-0.15, 60, 11.15, -20],axis:true});
	brd = JXG.JSXGraph.initBoard('jxgbox', {~~originX~~: 20, ~~originY: 300~~, ~~unitX: 60~~, ~~unitY: 5~~,axis:true});
	brd.suspendUpdate();		brd.suspendUpdate();
	var graph1,graph2,graph3,graph4,graph5,graph6,graph7,graph8,graph9;		var graph1,graph2,graph3,graph4,graph5,graph6,graph7,graph8,graph9;
Line 133:		Line 127:
	}});		}});
	}		}
	</~~script>~~		</jsxgraph>
	~~</html~~>

	=== The underlying source code ===		=== The underlying source code ===

A WASSERMANN at 10:07, 26 November 2009

2009-11-26T10:07:03Z

← Older revision		Revision as of 10:07, 26 November 2009
Line 10:		Line 10:
	<html>		<html>
	<link rel="stylesheet" type="text/css" href="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraph.css" />		<link rel="stylesheet" type="text/css" href="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraph.css" />
			<script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/prototype.js"></script>
	<script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraphcore.js"></script>		<script type="text/javascript" src="http://jsxgraph.uni-bayreuth.de/distrib/jsxgraphcore.js"></script>
	<div id="jxgbox" class="jxgbox" style="width:700px; height:400px;" onmouseup="doIt();"></div>		<div id="jxgbox" class="jxgbox" style="width:700px; height:400px;" onmouseup="doIt();"></div>