R语言 - 逻辑回归

> library( ggplot2)

#设置随机种子为1

> set.seed( 1)

#将整数1,2,3分别定义为变量b0,b1,b2

> b0 <- 1 ; b1 <- 2 ; b2 <- 3

#产生1000个服从正态分布的随机数分别定义为变量x1,x2

> x1 <- rnorm( 1000) ; x2 <- rnorm( 1000)

#将以上产生的b0,b1,b2,x1,x2计算定义为浮点数变量z

> z <- b0 + b1* x1 + b2* x2

#将z带入逻辑回归计算函数，定义为变量pr

> pr <- 1/( 1+ exp(- z))

#使用模拟产生二项分布数据函数rbinom，定义样本容量为1000，点分布的试验次数为1，点分布成功概率为pr “0<pr<1”

> y <- rbinom( 1000, 1, pr)