我想知道逻辑回归如果做不到这一点,有没有办法将它们转换为近似正态分布的 rv?
我的数据生成过程模型是最简单的模型,例如
在哪里
\ mathcal是伯努利分布。
所以在R代码中:
n<-1000
p<-5
x<-matrix(rnorm(n*(p-1)),nc=p-1);
y<-cbind(1,x)%*%rep(1,p)
y0<-1/(exp(-y)+1);
y0<-rbinom(n,1,y0)
fit<-glm(y0~x,family="binomial")
y<-residuals(fit,type="deviance")
qqnorm(y)
qqline(y)