我正在尝试帮助我的女朋友完成 R 课程的任务。我对代码一无所知,我知道 C++,但 RI 不知道。
我正在阅读文档,但如果你能帮助我,那就太好了!
以下代码来自该网站。
代码第一部分:
# rm(list=ls(all=TRUE))
# It's discourteous to wipe out your responder's entire workspace
install.packages("abind")
library(abind)
simHPP<-function(lambda,T,N){
# lambda: scalar, intensity of the Poisson process
# T: scalar, time horizon
# N: scalar, number of trajectories
EN <- rpois(N,lambda*T)
y <- matrix(T,nrow=2*max(EN)+2,ncol=N)*matrix(1,nrow=2*max(EN)+2,ncol=N)
yy <- abind(y,matrix(1,nrow=2*max(EN)+2,ncol=N)*EN,along=3)
i=1
while(i<=N){
if(EN[i]>0){
yy[1:(2*EN[i]+1),i,1] <- c(0,rep(sort(T*runif(EN[i])),each=2))
} else {
yy[1,i,1]=0
}
yy[1:(2*EN[i]+2),i,2] <- c(0,floor((1:(2*EN[i]))/2),EN[i])
i=i+1
}
return(yy)
}
功能abind
,它是如何工作的?R中的赋值运算符(=)看起来像(<-)?
例如,这是在做什么?
yy[1:(2*EN[i]+1),i,1] <- c(0,rep(sort(T*runif(EN[i])),each=2))