R: Ошибка: неверное количество измерений в функции - PullRequest
Новая
       15

R: Ошибка: неверное количество измерений в функции

0 голосов
/ 25 марта 2020

Я пытаюсь создать функцию в R для метода ядра, который является многомерным расширением метода, представленного в https://shareok.org/bitstream/handle/11244/11006/Chen_okstate_0664D_12928.pdf?sequence=1

library(mvtnorm)
library(MASS)
Kij<- function(x,y,n,H)
{
  temp1<-temp1b<-array(,dim = c(n-1,p))
  temp1=(1/det(H^(1/2)))*dmvnorm((x-y)%*%(solve(H^(1/2))))
  temp2=(x^2)*mean(temp1 %*% t(temp1))  ##KKij
  temp1b=((1/2)*(x+y))*temp1   

  return (matrix(c(x*mean(temp1),temp2, c(sum(temp1b[1,]^2), sum(temp1b[2,]^2))),ncol=p))
}

Qomega<- function(x,n,H){
temp1<- array(,dim=c(p,3))
temp2=temp3=temp4=matrix(0,nrow=p)
   for (j in 1:n){
     temp1=Kij(x[j,],x[-j,],n,H)
     temp2=temp2+temp1[1,]        #Kij to \hat{Q}i
     temp3=temp3+temp1[2,]        # int_{x^2 f^3(x)}
     temp4=temp4+temp1[3,]        #Aij^2 to Ai^2
   }
return(matrix(c(temp2*(1/n),4*(temp3/n-(temp2/n)^2),temp4),nrow=p))
}

Затем, когда я попытался проверить используя функцию 

H <- n[1]^(-1/(p+4))*solve(sigma^(1/2))
p <- 2
n<-c(20,20,20) 
mu1=rep(1,p)
mu2=rep(1,p)
mu3=rep(1,p)
sigma <- matrix(c(1.0, 0,
                  0, 1.0), nrow = 2)
groups<-rep(1:3, n)
x1=mvrnorm(n[1], mu1,sigma)
x2=mvrnorm(n[2], mu2,sigma)
x3=mvrnorm(n[3], mu3,sigma)
x=rbind(x1,x2,x3)

я получил следующую ошибку

enter image description here

Но когда я попытался проверить ошибку с фактическим номером, работал очень хорошо, как показано ниже

enter image description here

Я не могу понять, в чем именно проблема.

...