mux009 - 経済統計練習帳

var sig,mu;
var d=new Array();
var p=new Array();
var ps=new Array();
var mux=new Array();
var sigx=new Array();
 
 
function sample(){
 
 
for(s=1;s<101;s++){
p[s]=0.01;    
}    
 
tx=0;  
for(s=1;s<11;s++){
for(sx=1;sx<11;sx++){
tx=tx+1; 
mux[tx]=4.5+0.1*s;
sigx[tx]=0.5+0.1*sx;
}        
}
 
 
d[1]=5;
d[2]=4.9;
d[3]=4.7;
d[4]=5.3;
d[5]=4.8;
 
 
var timer;
var delay = 1000;
 
var loop = function () {
 
 for(s=1;s<101;s++){
p1=p[s];
for(sx=1;sx<6;sx++){
p1=p1*f(d[sx],mux[s],sigx[s]);
}   
ps[s]=p1;
}
 
p1=0;   
for(s=1;s<101;s++){
p1=p1+ps[s];    
}  
 
for(s=1;s<101;s++){
p[s]=ps[s]/p1;    
}      
 
 
    paint();
 
    clearTimeout(timer);
    timer = setTimeout(loop, delay);
}
 
loop();
 
 
 
 
 
}
 
 
  function f(x,mu,sig){
var pi,f1,f2;
 
        pi=Math.PI;
 
        f2=Math.sqrt(2*pi*sig*sig);
        f1=-(x-mu)*(x-mu)/(2*sig*sig);
 
        return Math.exp(f1)/f2;
}
 
function paint(){
 
var sub=new svgpaint();
 
str="<svg width=\"500\" height=\"500\">";
for(s=1;s<101;s++){
str=str+sub.rect(5*s,500-500*p[s],5,500*p[s],"red");    
} 
str=str+sub.line(0,0,0,500,3,"blue");    
str=str+sub.line(0,500,500,500,3,"blue");    
 
 
str=str+"</svg>";
 
$("#memo").html(str);   
 
}