第５回MIKU勉強会その２

第５回MIKU勉強会 - MIKU勉強会の続き

①について

x<- c(165,170,172,175,170,172,183,187,180,185)
y<- c(50,60,65,65,70,75,80,85,90,95)
X<- sum(x)
Y<- sum(y)
S<- rep(0,100)
a<- 1:10
b<- 1:10
A<- as.matrix(expand.grid(a,b))   #(a,b)は1~10までの整数の組み合わせ100通り
for(i in 1:100){　　　　　　　　　#以下の操作を繰り返す　
  a<- A[i,1];b<- A[i,2]
  S[i]<- Y-a*X-b
}

 image(matrix(S,10,10,byrow=TRUE))　　#byrow=TRUEは行ごとに要素を埋める

#次はa,bに-1000〜1000までの100個の数を入れる

X<- sum(x)　　　　　　　　　　　　
Y<- sum(y)
S<- rep(0,100^2)
a<- b<- seq(from=-1000,to=1000,length.out=100)
A<- as.matrix(expand.grid(a,b))
for(i in 1:(100^2)){
  a<- A[i,1];b<- A[i,2]
  csum<- 0
  for(k in 1:10){
    f<- a*x[k]+b
    c<- (y[k]-f)^2
    csum<- csum + c
  }
  S[i]<- csum
}
image(matrix(S,100,100,byrow=T))

persp(matrix(S,100,100,byrow=TRUE),phi=20,theta=60,col='pink')