如何使用 R 提取 estquant 的所有值？
How do I extract all the values of estquant using R?
我不知道如何从此循环中提取 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])
}
您希望 estquant 成为长度为 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 必须为每次迭代重新定义对象，这很耗时。）
我不知道如何从此循环中提取 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])
}
您希望 estquant 成为长度为 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 必须为每次迭代重新定义对象，这很耗时。）