我的变量筛选标准是估计系数不等于零,所以我用的连乘,若这十个变量的乘积不等于零,则认为筛选出了所有重要变量。假设我有四种方法,每种方法模拟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