R 生成模型数据

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

#生成Logisitc模型数据，其中x~N(0,1), a=2, b=-0.8;
#每次生成n个，生成s组
generate_data<-function(n,s,a=2,b=-0.8){
for(k in 1:s){
set.seed(k)
x=array(0,dim=c(n,s));p=array(0,dim=c(n,s));y=array(0,dim=c(n,s))
x[,k]=rnorm(n);  
p[,k]=exp(a+b*x[,k])/(1+exp(a+b*x[,k]))
y[,k]=rbinom(n,1,prob=p[,k])
k=k+1
}
list(k=k,x=x,y=y)
}
generate_data(10,5)
结果如下：
$k
[1] 6
$x
      [,1] [,2] [,3] [,4]        [,5]
[1,]    0    0    0    0 -0.84085548
[2,]    0    0    0    0  1.38435934
[3,]    0    0    0    0 -1.25549186
[4,]    0    0    0    0  0.07014277
[5,]    0    0    0    0  1.71144087
[6,]    0    0    0    0 -0.60290798
[7,]    0    0    0    0 -0.47216639
[8,]    0    0    0    0 -0.63537131
[9,]    0    0    0    0 -0.28577363
[10,]    0    0    0    0  0.13810822
$y
      [,1] [,2] [,3] [,4] [,5]
[1,]    0    0    0    0    1
[2,]    0    0    0    0    0
[3,]    0    0    0    0    1
[4,]    0    0    0    0    1
[5,]    0    0    0    0    1
[6,]    0    0    0    0    1
[7,]    0    0    0    0    1
[8,]    0    0    0    0    0
[9,]    0    0    0    0    1
[10,]    0    0    0    0    0

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：generate logisitc logisit logis Data 模型

只能生成最后一组数据，请问如何让修改可以得到全部的数据？

哦哦，终于搞明白了：
把x=array(0,dim=c(n,s));p=array(0,dim=c(n,s));y=array(0,dim=c(n,s))放在循环体外，就对了。
放在循环里面，每次都生成n*s个0，所以只有最后一次的值。

看看

楼主我不太明白

是个好东西，我看看

感谢分享

点赞支持

thanks for sharing

thanks for 你的分享