ｐｐｐ０９９９ｘ - 経済統計練習帳

var x=new Array();
var y=new Array();
var c=new Array();
var mu,sig 
var d=new Array();
var mus,sigs;
var mm=new Array();
var ss=new Array();
var number;
 
function sample(){
 
number=1;    
 
 
 
 
var timer;
var delay = 1000;
 
d[1]=50;
d[2]=54;
d[3]=47;
d[4]=60;
d[5]=44;    
 
mus=40;
sigs=10;
maxz=0;
 
mm[number]=mus;
ss[number]=sigs;  
 
var loop = function () {
 
dmu=1;
if(Math.random>0.5)dmu=-1;
dsig=0.1;
if(Math.random>0.5)dsig=-0.1;   
 
mu=mus+dmu;
sig=sigs+dsig;   
 
number=number+1;
 
mm[number]=mu;
ss[number]=sig;      
 
for(s=1;s<10001;s++){
x[s]=randx(mu,sig);    
}    
 
for(s=1;s<10001;s++){
y[s]=Math.floor(x[s]);  
}    
 
for(s=0;s<101;s++){
c[s]=0;    
}    
 
for(s=1;s<10001;s++){
sx=y[s];
if(sx>0)c[sx]=c[sx]+1;  
}     
 
z=1;    
for(s=1;s<6;s++){
sx=d[s];
z=z*c[sx];    
}    
 
if(z>maxz)mus=mu;
if(z>maxz)sigs=sig;
if(z>maxz)maxz=z;    
 
paint();    
 
 
clearTimeout(timer);
    timer = setTimeout(loop, delay);
}
 
loop();
 
 
}
 
 
function randx(mu,sig) {
 
 
  pi = Math.PI;
 
  p1 = 0;
  p2 = 0;
 
  x1=Math.random();
  x2=Math.random();
 
  f1=Math.log(x1);
  r=Math.sqrt(-2*f1)*Math.cos(2 * pi * x2);
 
 
 
p3=mu+sig*r;
 
return p3;
 
 
}
 
 
function paint(){
 
      var sub=new svgpaint();
 
str="<svg width=\"700\" height=\"700\">";
for(s=1;s<number+1;s++){
str=str+sub.rect(5*mm[s],500-20*ss[s],5,5,"blue");      
}
str=str+"</svg>";
$("#memo").html(str);
 
 
 
  }