计算R中的克矩阵

Question

我需要一些帮助来计算给定内核的nXn克矩阵K.这是我生成模拟数据的R代码。我可以是任何正定矩阵。以对角矩阵为简单。

     set.seed(3)
     n=20
     x=runif(n)
     y=rnorm(n)
    df<-cbind(x,y)

  I=diag(2)

kernel<-function(x,y) {
    t(x)%*%I%*%y
     }
   # for example 
#K[1,1]
    t(df[1,])%*%I%*%df[1,]
          [,1]
    [1,] 0.5829376
#K[1,2]
t(df[1,])%*%I%*%df[2,]
          [,1]
[1,] 0.978207