梯度下降算法-R语言

#读取自变量、因变量数据
x<-read.table("ex2x.txt")
y<-read.table("ex2y.txt")
Y=y[[1]]
m<-length(y[[1]])
#自变量增加一列构造矩阵
x1<-rep(1,m)
x2<-x[[1]]
X<-data.frame(x1,x2)
X<-as.matrix(X)
#初始化参数
theta<-matrix(c(0,0),ncol=1)
#设定迭代次数
maxiteration=2000
#设定学习速度
alpha=0.07
for ( n in c(1:maxiteration)){
  #计算梯度
  grad=(1/m)*t(X)%*%((X%*%theta)-Y)
  #参数更新或改进
  theta=theta-alpha*grad
  print(theta)
}