>
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)