R 如何在总面积中找到LGCP随机场λ值?

R 如何在总面积中找到LGCP随机场λ值?,r,random,attributes,R,Random,Attributes,RandomField包中有一个rLGCP模型示例 if(require(RandomFields)) { # homogeneous LGCP with exponential covariance function X <- rLGCP("exp", 3, var=0.2, scale=.1) # inhomogeneous LGCP with Gaussian covariance function m <- as.im(function(x, y){5 -

RandomField包中有一个rLGCP模型示例

 if(require(RandomFields)) {
  # homogeneous LGCP with exponential covariance function
  X <- rLGCP("exp", 3, var=0.2, scale=.1)

  # inhomogeneous LGCP with Gaussian covariance function
  m <- as.im(function(x, y){5 - 1.5 * (x - 0.5)^2 + 2 * (y - 0.5)^2}, W=owin())
  X <- rLGCP("gauss", m, var=0.15, scale =0.5)
  plot(attr(X, "Lambda"))
  points(X)

  }
if(需要(随机字段)){
#具有指数协方差函数的齐次LGCP

X我不完全确定这是否是您要查找的,但绘图中每个点的Lambda值矩阵存储在由
spatstat::rLGCP
创建的模型的
Lambda
属性中

您可以这样访问它们:

m <- as.im(function(x, y){5 - 1.5 * (x - 0.5)^2 + 2 * (y - 0.5)^2}, W=owin())
X <- rLGCP("gauss", m, var=0.15, scale = 0.5)
lambda_matrix <- attr(X, "Lambda")$v
m