simulateSpatialData: simulate normal distributed data

View source: R/simulateSpatialData.R

simulateSpatialDataR Documentation

simulate normal distributed data

Description

Simulate normal distributed data with spatial correlation structure

'theta' ($\theta$) describes how rapidly the correlation declines with respect to the distance between two voxels. The three-dimensional coordinates of the voxels are defined as all combinations of vector $c = 1, ..., m1/3$, then $\Sigma_\theta = \exp(-\theta K)$ where $K$ is the matrix containing the euclidean distances between the three-dimensional coordinates' voxels. So, $m^1/3$ must be an integer value.

Usage

simulateData(pi0,m,n, theta, seed = NULL, power = 0.8, alpha = 0.05)

Arguments

pi0

numeric value in '[0,1]'. Proportion of true null hypothesis.

m

numeric value. Number of variables.

n

numeric value. Number of observations.

theta

numeric value in '[0,1]'. Level of correlation between pairs of variables. See details

seed

integer value. If you want to specify the seed. Default to @NULL

power

numeric value in '[0,1]'. Level of power. Default 0.8.

alpha

numeric value in '[0,1]'. It expresses the alpha level to control the family-wise error rate. Default 0.05.

Value

Returns a matrix with dimensions m \times n.

Author(s)

Angela Andreella


angeella/ARIpermutation documentation built on Aug. 24, 2023, 3:36 p.m.