集研發功能錯誤

Question

我有下列R-代碼：

第一部分工作：

denSim <- function(n,nrow,ncol,v,ar) { 
XX <- matrix(NA,nrow=nrow,ncol=ncol) 
h1 <- h2 <- kknum <- kksum <- abc <- NULL 

for(i in 1:n){ 
XX[,i] <- arima.sim(list(order=v,ar=ar), n=nrow) 
h1[i]<-length(XX[,i])^(-1/5) 
h2[i]<-sd(XX[,i])*h1[i] 
kknum[[i]] <- exp(-0.5*((XX[nrow,i]-XX[,i])/bw.nrd0(XX[,i]))^2) 
kksum[i]<-sum(kknum[[i]][-nrow]) 
abc[i]<- 1/(h2[i]*sqrt(2*pi)) 
} 
Simls <- list(XX = XX,h1 = h1,h2=h2,kknum = kknum,kksum = kksum,abc = abc) 
return(Simls) 
} 

Sim <- denSim(n=200,nrow=500,ncol=200,v=c(1,0,0),ar=0.6) 

第二部分是不工作：

pdf <- function(z,i){ 
ABC <-Sim$abc[i] 
Z <- z - Sim$XX[-1,i] 
H2 <- Sim$h2[i] 
eta <- ((Z)/H2)^2 
KKN <- Sim$kknum[[i]][-500] 
KKS <- Sim$kksum[i] 
return(ABC*sum(KKN*exp(-0.5*eta))/KKS) 
} 
cdf <- function(i) integrate(pdf,0.6*mean(Sim$XX[,i]),20,i)$value 
cdf(1) 

是什麼問題？

Answer 1

integrate(...)採用「矢量化函數」，例如，一個將矢量作爲參數並返回相同長度的矢量的函數。你的函數pdf(...)需要一個標量值爲x並返回一個標量。

幸運的是，有一個簡單的解決方案，R中的Vectorize(...)函數是一個將標量函數轉換爲矢量函數的包裝函數。

pdf.vect <- Vectorize(pdf) 
cdf <- function(i) integrate(pdf.vect,0.6*mean(Sim$XX[,i]),20,i)$value 
cdf(1) 
# [1] 0.6820746 

編輯（應答到OP的評論）

是的，但integrate(...)也沒有矢量，所以你必須要做到這一點爲好。

cdf <- Vectorize(function(x,i) integrate(pdf.vect,0.6*mean(Sim$XX[,i]),x,i)$value) 

x <- seq(0,5,0.1) 
plot(x,pdf.vect(x,1),type="l", ylab="Y") 
par(new=T) 
plot(x,cdf(x,1),yaxt="n", xaxt="n", col="red",type="l",ylab="",)