1. 矩阵乘法
用两个百分号%包住星号*(如果不用百分号,就是按元素的乘法)
p<-exp(X%*%beta_true)/(1+exp(X%*%beta_true))
2. 矩阵求逆
NR_step<--solve(t(X)%*%W%*%X)%*%t(X)%*%(y-p)
1. 矩阵乘法
用两个百分号%包住星号*(如果不用百分号,就是按元素的乘法)
p<-exp(X%*%beta_true)/(1+exp(X%*%beta_true))
2. 矩阵求逆
NR_step<--solve(t(X)%*%W%*%X)%*%t(X)%*%(y-p)