0

我正在尝试帮助我的女朋友完成 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))
4

1 回答 1

1

有一个代码审查网站,这可能是合适的,因为您(部分)要求解释您没有给出上下文的代码。

(t 使用大写 T 作为向量名称虽然不是错误,但令人困惑,因为对该名称的默认分配是逻辑 TRUE。)

abind将数组作为参数并返回增广数组。如果您有两个 2 x 2 x 2 数组,abind那么您将得到一个 2 x 2 x 4 数组

abind( array(1:8, c(2,2,2) ),  array (9:16, c(2,2,2)))
#------------
, , 1
     [,1] [,2]
[1,]    1    3
[2,]    2    4

, , 2
     [,1] [,2]
[1,]    5    7
[2,]    6    8

, , 3
     [,1] [,2]
[1,]    9   11
[2,]   10   12

, , 4
     [,1] [,2]
[1,]   13   15
[2,]   14   16
于 2012-06-04T03:04:36.703 回答