findintercorr_cont_nb

the method used to generate the <code>k_cont</code> continuous variables. "Fleishman" uses a third-order polynomial transformation
and "Polynomial" uses Headrick's fifth-order transformation.

method

a matrix with <code>k_cont</code> rows, each a vector of constants c0, c1, c2, c3 (if <code>method</code> = "Fleishman") or
c0, c1, c2, c3, c4, c5 (if <code>method</code> = "Polynomial"), like that returned by <code><a rd-options="SimMultiCorrData" href="/link/find_constants?package=SimMultiCorrData&version=0.2.2&to=SimMultiCorrData" data-mini-rdoc="SimMultiCorrData::find_constants">find_constants</a></code>

constants

a <code>k_cont x k_nb</code> matrix of target correlations among continuous and Negative Binomial variables

rho_cont_nb

a vector of size parameters for the Negative Binomial variables (see <code><a rd-options="stats" href="/link/NegBinomial?package=SimMultiCorrData&version=0.2.2&to=stats" data-mini-rdoc="stats::NegBinomial">NegBinomial</a></code>)

size

a vector of success probability parameters

prob

a vector of mean parameters (*Note: either <code>prob</code> or <code>mu</code> should be supplied for all Negative Binomial variables,
not a mixture; default = NULL)

the number of random numbers to generate in calculating the bound (default = 10000)

nrand

the seed used in random number generation (default = 1234)

seed

This function calculates a <code>k_cont x k_nb</code> intermediate matrix of correlations for the <code>k_cont</code> continuous and
 <code>k_nb</code> Negative Binomial variables. It extends the method of Amatya &amp; Demirtas (2015, 10.1080/00949655.2014.953534) to
 continuous variables generated using
 Headrick's fifth-order polynomial transformation and Negative Binomial variables. Here, the intermediate correlation
 between Z1 and Z2 (where Z1 is the standard normal variable transformed using Headrick's fifth-order or Fleishman's
 third-order method to produce a continuous variable Y1, and Z2 is the standard normal variable used to generate a
 Negative Binomial variable via the inverse cdf method) is calculated by dividing the target correlation by a correction factor.
 The correction factor is the product of the upper Frechet-Hoeffding bound on the correlation between a Negative Binomial
 variable and the normal variable used to generate it (see <code><a rd-options="SimMultiCorrData" href="/link/chat_nb?package=SimMultiCorrData&version=0.2.2&to=SimMultiCorrData" data-mini-rdoc="SimMultiCorrData::chat_nb">chat_nb</a></code>) and the power method
 correlation (described in Headrick &amp; Kowalchuk, 2007, 10.1080/10629360600605065) between Y1 and Z1. The function is used in
 <code><a rd-options="SimMultiCorrData" href="/link/findintercorr?package=SimMultiCorrData&version=0.2.2&to=SimMultiCorrData" data-mini-rdoc="SimMultiCorrData::findintercorr">findintercorr</a></code> and <code><a rd-options="SimMultiCorrData" href="/link/rcorrvar?package=SimMultiCorrData&version=0.2.2&to=SimMultiCorrData" data-mini-rdoc="SimMultiCorrData::rcorrvar">rcorrvar</a></code>.
 This function would not ordinarily be called by the user.

correlation

continuous

binomial

negative

Fleishman

Headrick

intermediate

Generate continuous (normal or non-normal), binary, ordinal, and count (Poisson or Negative
Binomial) variables with a specified correlation matrix. It can also produce a single continuous
variable. This package can be used to simulate data sets that mimic real-world situations (i.e.
clinical or genetic data sets, plasmodes). All variables are generated from standard normal
variables with an imposed intermediate correlation matrix. Continuous variables are simulated
by specifying mean, variance, skewness, standardized kurtosis, and fifth and sixth standardized
cumulants using either Fleishman's third-order (<DOI:10.1007/BF02293811>) or Headrick's
fifth-order (<DOI:10.1016/S0167-9473(02)00072-5>) polynomial transformation. Binary and
ordinal variables are simulated using a modification of the ordsample() function from 'GenOrd'.
Count variables are simulated using the inverse cdf method. There are two simulation pathways
which differ primarily according to the calculation of the intermediate correlation matrix. In
Correlation Method 1, the intercorrelations involving count variables are determined using a
simulation based, logarithmic correlation correction (adapting Yahav and Shmueli's 2012 method,
<DOI:10.1002/asmb.901>). In Correlation Method 2, the count variables are treated as ordinal
(adapting Barbiero and Ferrari's 2015 modification of GenOrd, <DOI:10.1002/asmb.2072>).
There is an optional error loop that corrects the final correlation matrix to be within a
user-specified precision value of the target matrix. The package also includes functions to
calculate standardized cumulants for theoretical distributions or from real data sets, check
if a target correlation matrix is within the possible correlation bounds (given the distributions
of the simulated variables), summarize results (numerically or graphically), to verify valid power
method pdfs, and to calculate lower standardized kurtosis bounds.

Allison Fialkowski

SimMultiCorrData

Simulation of Correlated Data with Multiple Variable Types

findintercorr_cont_nb function

a matrix with <code>k_cont</code> rows, each a vector of constants c0, c1, c2, c3 (if <code>method</code> = "Fleishman") or
c0, c1, c2, c3, c4, c5 (if <code>method</code> = "Polynomial"), like that returned by <code><a rd-options='SimMultiCorrData' href='find_constants'>find_constants</a></code>

a vector of size parameters for the Negative Binomial variables (see <code><a rd-options='stats' href='NegBinomial'>NegBinomial</a></code>)

This function calculates a <code>k_cont x k_nb</code> intermediate matrix of correlations for the <code>k_cont</code> continuous and
 <code>k_nb</code> Negative Binomial variables. It extends the method of Amatya &amp; Demirtas (2015, 10.1080/00949655.2014.953534) to
 continuous variables generated using
 Headrick's fifth-order polynomial transformation and Negative Binomial variables. Here, the intermediate correlation
 between Z1 and Z2 (where Z1 is the standard normal variable transformed using Headrick's fifth-order or Fleishman's
 third-order method to produce a continuous variable Y1, and Z2 is the standard normal variable used to generate a
 Negative Binomial variable via the inverse cdf method) is calculated by dividing the target correlation by a correction factor.
 The correction factor is the product of the upper Frechet-Hoeffding bound on the correlation between a Negative Binomial
 variable and the normal variable used to generate it (see <code><a rd-options='SimMultiCorrData' href='chat_nb'>chat_nb</a></code>) and the power method
 correlation (described in Headrick &amp; Kowalchuk, 2007, 10.1080/10629360600605065) between Y1 and Z1. The function is used in
 <code><a rd-options='SimMultiCorrData' href='findintercorr'>findintercorr</a></code> and <code><a rd-options='SimMultiCorrData' href='rcorrvar'>rcorrvar</a></code>.
 This function would not ordinarily be called by the user.

findintercorr_cont_nb: Calculate Intermediate MVN Correlation for Continuous - Negative Binomial Variables: Correlation Method 1

Description

Usage

Arguments

Value

References

See Also