-1

我希望我正在开发的应用程序自动转到 Excel 表并计算转移概率矩阵。

 action<-actions
    state<-states
    p<-array(0,c(length(state),length(state),length(action)))
    colnames(p)=state
    rownames(p)=state
    # transition probability matrix based on the action that we have choosen
    empiricala1<-read.xlsx("www/dynamic model.xlsx",1)
    empiricala2<-read.xlsx("www/dynamic model.xlsx",2)


    #show(empirical) calculate transition probability from historical data
    em1<-as.data.frame(empiricala1)
    em2<-as.data.frame(empiricala2)
    tab1 <- table(em1)
    tab2 <- table(em2)

    tab1<-addmargins(prop.table(table(em1$current,em1$nextstate),1),2)
    tab2<-addmargins(prop.table(table(em2$current,em2$nextstate),1),2)
    transitionprob1<-p[,,1]<-prop.table(table(em1$current,em1$nextstate),1)
    transitionprob2<-p[,,2]<-prop.table(table(em2$current,em2$nextstate),2)
    print(transitionprob1)
    print(transitionprob2)


    for(i in 1:length(action)){

      p[,,i]<-prop.table(table(em[i]$current,em[i]$nextstate),i)

    }

我得到的错误如下：

Error in em[i] : object of type 'closure' is not subsettable

我该如何解决这个问题。

4

1 回答 1

1

好的，所以要扩展评论...

您有两个数据框em1和em2. 您似乎希望将相同的操作应用于两个数据帧（例如table(em1)和table(em2)。编写起来真的很乏味，尤其是当您获得更多变量时。

你试图做的是：

for (i in 1:2) em[i]

这样做的问题是你没有得到em1and em2。取而代之的是，您得到em[1]并em[2]因此引用了该em对象（一开始就不存在）。

有几种方法可以解决这个问题。

1. 将数据框移动到列表中

lis <- list(em1, em2)

*apply这样，您可以使用例如家庭或for-loop来遍历列表：

sapply(lis, nrow)
for (i in 1:length(lis)) nrow(lis[[i]])

2.使用get

另一种选择是使用get，它允许您提供一个字符串，随后它将返回该字符串中描述的变量。

vec <- c("em1","em2")
sapply(vec, function(x) nrow(get(x)))

请注意，get通常不鼓励使用。我也会选择第一个选项。

于 2016-02-27T01:33:42.770 回答