0

这是我的问题的一个例子:

    #data#
    g<-c(1,1,1,2,2,2)
    A<-runif(6,min=1,max=5)
    B<-runif(6,min=100,max=1000)
    C<-runif(6,min=30,max=31)
    D<-runif(6,min=67,max=98765)
    var<-cbind(A,B,C,D)
    label<-colnames(var)

    store<-data.frame(matrix(ncol=2))
    colnames(store)=c("usedVar","prediction")

    library(MASS)#get lda

    for (i in c(1:4))
    {
    for (k in c(1:4))
    {
    if(i==k)break()#using the same variable will not work in lda

    dis<-lda(g~var[,i]+var[,k],CV=TRUE,fold=6)#linear discriminant function
    ct<-table(g,dis$class)
    pre<-sum(diag(prop.table(ct)))# Prediction of cv-reclassification

    #store the results
    a<-i+4*(k-1)
    store[a,1]<-paste(label[i],label[k])
    store[a,2]<-round(pre,4)
    }
    }
    store

我有两个 for 循环被 if-break-case 所“扰乱”(扰乱的意思是,我不需要这些情况的信息)。通过使用 a<-i+4*(k-1) 将结果存储到我的“结果”data.frame,中断案例也将作为 NA 包含在内。有没有办法在没有 NA 案例的情况下存储我的循环结果以获得连续数据集?

4

1 回答 1

0

你主要需要改变你的计数器变量,a. 而是随着循环的进行增加一个计数器。我还更改了 if 语句以封装内部循环的整个内容,并预先分配store到适当的维度(我认为):

store<-data.frame(matrix(ncol=2,nrow=c(i*k-min(i,k))))
colnames(store)=c("usedVar","prediction")

library(MASS)#get lda
row.i <- 0 # start counter
for (i in c(1:4)) {
  for (k in c(1:4)) {
    if(i != k) {
      row.i <- row.i + 1 # only increment if condition met
      dis<-lda(g~var[,i]+var[,k],CV=TRUE,fold=6)#linear discriminant function
      ct<-table(g,dis$class)
      pre<-sum(diag(prop.table(ct)))# Prediction of cv-reclassification

      #store the results
      store[row.i,1]<-paste(label[i],label[k])
      store[row.i,2]<-round(pre,4)
    }
  }
}
store
于 2012-08-15T18:22:01.233 回答