R 如何在总面积中找到LGCP随机场λ值?
RandomField包中有一个rLGCP模型示例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 -
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