小编Nic*_*ick的帖子

关于绘制置信区间有很多答案.

我正在阅读Lourme A.等人(2016年)的论文,我想从论文中得出90%的置信区间和10%的特殊点,如图2所示:.

我不能使用LaTeX并插入带有置信区域定义的图片:

library("MASS")
library(copula)
set.seed(612)

n <- 1000 # length of sample
d <- 2    # dimension

# random vector with uniform margins on (0,1)
u1 <- runif(n, min = 0, max = 1)
u2 <- runif(n, min = 0, max = 1)

u = matrix(c(u1, u2), ncol=d)

Rg  <- cor(u)   # d-by-d correlation matrix
Rg1 <- ginv(Rg) # inv. matrix 

# round(Rg %*% Rg1, 8) # check

# the multivariate c.d.f of …

我有一个矩阵 M：

n = 3    
x=c(0.85, 0.1, 0.05)
M <- matrix(NA, n, n); 

for(i in 1:n){
for(j in 1:n){
M[i,j] = x[i] * x[j]
}}

#       [,1]  [,2]   [,3]
# [1,] 0.7225 0.085 0.0425
# [2,] 0.0850 0.010 0.0050
# [3,] 0.0425 0.005 0.0025

我需要找到所有反对角线的总和，包括 M[1,1] 和 M[n, n]。我的尝试是

d <-matrix(c(0, 1, 2, 1, 2, 3, 2, 3, 4), n)
tapply(M, d, sum)

     0      1      2      3      4 
0.7225 0.1700 0.0950 0.0100 0.0025 

结果对我来说是正确的。

问题。如何定义矩阵d的条目？可以作为 col(M) 和 row(M) 上的函数。

我有一个矩阵：

A<-t(matrix(
c(0, 0, 1,
  0, 0, 0,
  0, 0, 1,
  0, 0, 1,
  0, 0, 0,
  1, 1, 0), 3, 6))

我需要保留只出现一次的列。因此，预期结果只是第三列：(1, 0, 1, 1, 0, 0)。

我已经找到了唯一和重复的函数，但我需要更强大的东西来删除出现多次的所有列（在我的示例中是第一个和第二个）。

我有两个4乘4阵列:

a1 = [ 1 2 3 4; 5 6 7 8; 9 10 11 12; 13 14 15 16 ]
a2 = [ 17 18 19 20; 21 22 23 24; 25 26 27 28; 29 30 31 32 ]

我需要创建16×8阵列C:

1   2   3   4   17  18  19  20
1   2   3   4   21  22  23  24
1   2   3   4   25  26  27  28
1   2   3   4   29  30  31  32
5   6   7   8 …

我想做10次迭代,但有时标志大于1.在这种情况下,我必须将变量i减少1并再次计算标志.怎么做这个减量？

    for(i in (1:n)){ 
    flag <- ... # some code
        if (flag > 1) {
         # some code
        } 
        #else decrement i <- i - 1 ??
    }

谢谢.