关于程序模拟的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

    在模拟研究中，设定变量个数为100个，其中有10个为重要变量，其回归系数不等于零，其他90个变量的回归系数设为零。当模拟次数为100次时，我想统计某种方法在这100次模拟中正确筛选出重要变量的次数，应该用什么函数实现呢？例如这个方法在第一次模拟中筛选出20个变量，包括了10个重要变量，则记为1次，在第二次模拟中筛选出15个变量，但是不全包括10个重要变量，则记为零次，依次类推，想统计100次模拟中有多少次包含重要变量。跪求答案，请高手不吝赐教啊。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：回归系数 变量个数 筛选出 求答案 程序

先把每一次筛选出的变量名都存下来，然后看是否包含那10个重要变量。

帮顶

我的变量筛选标准是估计系数不等于零，所以我用的连乘，若这十个变量的乘积不等于零，则认为筛选出了所有重要变量。假设我有四种方法，每种方法模拟100次，变量个数为1000，重要变量为前10个，我是这样写的程序：c<-array(1:24,c(2,3,4))###4表示四种方法，3表示变量个数，2表示模拟次数。现在我要将前两个变量相乘，看他的乘积是否为零。
P1<-array(0,c(2,4))
P2<-array(0,c(2,4))
for(i in 1:2)
        {for (k in 1:4)
     d<-c[i,1:2,k]
         P1[i,k]<-prod(d)
         P2[i,k]<-as.numeric(P1[i,k]!=0)
}
P<-apply(P2,2,sum)/2