で８８８８６ - 経済統計練習帳

var x=new Array();
var y=new Array();
var number;
var sol=new Array();
var datanumber;
 
function myFunction() {
 
number=2;  
datanumber=5;  
 
  for(s=1;s<datanumber+1;s++){
  x[s]=new Array();
  }
 
 
    id="1dX9slF7uA1HOco_x0SGHmBjKlJxd2JjzSR8Hfj7k8qI";
 
 var ex1 = SpreadsheetApp.openById(id);
 var sheet = ex1.getSheetByName("data");
 
 
    for(s=1;s<datanumber+1;s++){
   x[s][2]=sheet.getRange(s,1).getValue();
   x[s][1]=1;
   y[s]=sheet.getRange(s,2).getValue(); 
    }
 
  ols();
 
 
  sheet.getRange(7,1).setValue(sol[1]);
  sheet.getRange(7,2).setValue(sol[2]);
 
 
}
 
 
 
 
 
function ols(){
 
var a=new Array();
var b=new Array();
 
for(s=1;s<number+1;s++){
a[s]=new Array();    
}
 
 
for(s1=1;s1<number+1;s1++){ 
b1=0;    
for(s=1;s<datanumber+1;s++){
b1=b1+y[s]*x[s][s1];   
}    
 
b[s1]=b1;
}
 
 
for(s1=1;s1<number+1;s1++){ 
for(s2=1;s2<number+1;s2++){ 
a1=0;    
for(s=1;s<datanumber+1;s++){
a1=a1+x[s][s2]*x[s][s1];    
}    
a[s1][s2]=a1;
}
}
 
 
 
 
 
for(s1=1;s1<number+1;s1++){
for(s2=1;s2<number+1;s2++){
z=a[s2][s1]/a[s1][s1];
if(s1==s2)z=0;
for(s=1;s<number+1;s++){
a[s2][s]=a[s2][s]-z*a[s1][s];
}
b[s2]=b[s2]-z*b[s1]
 
}
}
 
for(s=1;s<number+1;s++){
sol[s]=b[s]/a[s][s];
}
 
 
 
 
}