我不知道如何從這個回圈中提取 estquant 的每個值。我應該在最后添加什么代碼來提供所有值而不是一個值！
p <- 0.5
m <- 2
d1 <- as.matrix(d);d1
for (i in 1:m){
Xj <- d1[,i]
nj <- length(Xj)
Fj <- pbeta(Fx,i,nj 1-i)
a <- pbeta(p,i,nj 1-i)
estFj <- knots(ecdf(Xj))
estquant <- min(estFj[estFj >= a])
}
uj5u.com熱心網友回復：
你想estquant成為一個 length 的向量m。
所以：
p <- 0.5
m <- 2
d1 <- as.matrix(d);d1
estquant <- numeric(m)
for (i in 1:m){
Xj <- d1[,i]
nj <- length(Xj)
Fj <- pbeta(Fx,i,nj 1-i)
a <- pbeta(p,i,nj 1-i)
estFj <- knots(ecdf(Xj))
estquant[i] <- min(estFj[estFj >= a])
}
estquant
（在回圈中為物件 1 對 1 賦值時預定義物件很重要，否則 R 必須為每次迭代重新定義物件，這很耗時。）
轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/386702.html