我正在尝试从中创建条件概率矩阵：
我已经重写了没有循环的代码
a # Signal Vector
b # Price Change Vector
Signalt<- seq(0, 1, 0.05) # Produce the 5% tiles
abst <- c(seq(1, 1.02, by = 0.0025), 2) #Produce the 0% to 2% tiles with 0.25% increments. Added 1 to include price change of 0 in `temp`
xbool = ((Signal >= Signalt[1] & a < Signalt[1 + 1]) *1) # 1 for True 0 for False
temp = (PercChange + 1) * xbool
temp2 <- temp[which(temp > 0)]
CondProb <- cut(temp2, abst, include.lowest = T)
table(CondProb)
这将输出带有带有出现次数的绝对列的表。我当然需要它占每行总数的百分比，但我希望首先能够运行循环并获得矩阵输出。
Loop Original - 应该在很大程度上被忽略，因为我已经更改了大部分编码设置
Signal <- runif(100)
PercChange <- abs((rnorm(100)/100))
signalt <- seq(0, 1, 0.05)
abst <- seq(0, c(0.02:1), 0.0025)
CondDistMat <- matrix(0, nrow = length(signalt), ncol = length(abst))
for(j in 1:length(signalt - 1)){
xbool = (is.na((Signal >= signalt[j] & Signal < signalt[j + 1]) ) * 1)
ysubset = (PercChange * xbool[j] )
CondProb = hist(ysubset, breaks = abst, freq = TRUE)
CondDistMat[signalt, abst] <- CondProb$density
}
列将是由 定义的百分位数，abst而行将是由 定义的 5% 平铺signalt。这个想法是通过布尔向量产生1，其中绝对回报PercChange应该在列中，然后绘制每个概率signalt。
但是，我无法产生输出 - 任何人都可以发现错误吗？提前致谢
所需的输出应该类似于附加的图像
。