R 这个简单的随机梯度下降代码有什么问题?
我正在尝试使用R SGD软件包学习SGD。然而,我无法正确理解这个简单的示例。怎么了?(注意,对于新加坡元X*beta,所有数字都是正数,这意味着我们将始终预测1)R 这个简单的随机梯度下降代码有什么问题?,r,logistic-regression,gradient-descent,R,Logistic Regression,Gradient Descent,我正在尝试使用R SGD软件包学习SGD。然而,我无法正确理解这个简单的示例。怎么了?(注意,对于新加坡元X*beta,所有数字都是正数,这意味着我们将始终预测1) 库(sgd) 种子集(0) d=ggplot2::钻石 d$price=ifelse(d$price>2500,1,0) 噪声_idx=样本(nrow(d),nrow(d)/5) d$price[noise_idx]=ifelse(d$price[noise_idx]==1,0,1) glm.fit=glm(价格~.-克拉,d,系
库(sgd)
种子集(0)
d=ggplot2::钻石
d$price=ifelse(d$price>2500,1,0)
噪声_idx=样本(nrow(d),nrow(d)/5)
d$price[noise_idx]=ifelse(d$price[noise_idx]==1,0,1)
glm.fit=glm(价格~.-克拉,d,系列=二项式)
菲特军士长
library(sgd)
set.seed(0)
d=ggplot2::diamonds
d$price=ifelse(d$price>2500,1,0)
noise_idx=sample(nrow(d),nrow(d)/5)
d$price[noise_idx]= ifelse(d$price[noise_idx]==1,0,1)
glm.fit=glm(price~.-carat,d,family=binomial)
sgd.fit <- sgd(price ~ .-carat, data=d,
model="glm", model.control=binomial(link="logit"))
glm.p=model.matrix(price~.-carat,d) %*% glm.fit$coefficients
sgd.p=model.matrix(price~.-carat,d) %*% sgd.theta$coefficients
par(mfrow=c(1,2))
hist(glm.p,50)
hist(sgd.p,50)