pt88 - 経済統計練習帳

var sol_ols=new Array();
var x_ols = new Array();
var y_ols = new Array();
var number_ols,datanumber_ols;
 
var number;
var loandi=new Array();
var topix=new Array();
var irl=new Array();
var irloan=new Array();
var topix=new Array();
var id; 
var ex;
 
function data_b(){
 
var sheet = ex.getSheetByName("b");
 
var v=sheet.getRange(1, 1,67,7).getValues();
 
 for(s=1;s<number+1;s++){
 
 irl[s]=v[s][0]; 
 irloan[s]=v[s][1];  
   loandi[s]=v[s][2]; 
 topix[s]=v[s][3];
 }
 
 
}
 
function myFunction() {
 
 
  id="1PL9vcWZgA8b7LAjMV5j0FeHRCJF-xW1eFrvI2Fgi9E4";
  ex = SpreadsheetApp.openById(id);
 
  number=66;
 
 data_b();
 
  ex_b2();
 
  ex_b3();
 
 
 
}
 
function ex_b2(){
 
number_ols=2;
 datanumber_ols=number;
  for(s=1;s<datanumber_ols+1;s++){
  x_ols[s]=new Array();
  }
 
 
  for(s=1;s<datanumber_ols+1;s++){
  x_ols[s][1]=irl[s];
  x_ols[s][2]=loandi[s];
  y_ols[s]=irloan[s];
  }
 
 
  ols();
 
  Logger.log(sol_ols[1]);
  Logger.log(sol_ols[2]);
 
 
}
 
 
function ex_b3(){
 
number_ols=3;
 datanumber_ols=number;
  for(s=1;s<datanumber_ols+1;s++){
  x_ols[s]=new Array();
  }
 
 
  for(s=1;s<datanumber_ols+1;s++){
  x_ols[s][1]=1;
  x_ols[s][2]=loandi[1];
  z=Math.log(topix[2])-Math.log(topix[1]); 
  if(s>1)z=Math.log(topix[s])-Math.log(topix[s]);  
  x_ols[s][3]=z;
  if(s>1)x_ols[s][2]=loandi[s-1];
   y_ols[s]=loandi[s];
  }
 
 
  ols();
 
  Logger.log(sol_ols[1]);
  Logger.log(sol_ols[2]);
  Logger.log(sol_ols[3]);
 
 
}
 
 
function ols(){
 
var a_ols=new Array();
var b_ols=new Array();
 
for(s=1;s<number_ols+1;s++){
a_ols[s]=new Array();
}
 
 
for(s=1;s<number_ols+1;s++){    
y1 = 0;
for (sx = 1; sx <datanumber_ols;sx++ ){
y1=y1+x_ols[sx][s]*y_ols[sx];    
}
b_ols[s]=y1;
}
 
 
 
for(s1=1;s1<number_ols+1;s1++){    
for(s2=1;s2<number_ols+1;s2++){ 
y1 = 0;
for (sx = 1; sx <datanumber_ols;sx++ ){
y1=y1+x_ols[sx][s1]*x_ols[sx][s2];    
}
a_ols[s1][s2]=y1;
}} 
 
 
for(s=1;s<number_ols+1;s++){
for(sx=1;sx<number_ols+1;sx++){
z=a_ols[sx][s]/a_ols[s][s];
if(s==sx)z=0;
for(tx=1;tx<number_ols+1;tx++){
a_ols[sx][tx]=a_ols[sx][tx]-z*a_ols[s][tx];
}
b_ols[sx]=b_ols[sx]-z*b_ols[s];
}
}
 
 
for(s=1;s<number_ols+1;s++){
sol_ols[s]=b_ols[s]/a_ols[s][s];
}
 
 
 
}