Defines functions rhoMax1dep

Documented in rhoMax1dep

#' To calculate the maximal allowed correlations
#' max for using cBern1dep to generate binary data
#' with 1-dependent structure
#' @param p the vector of marginal probabilities with dimension m
#' @return an (m-1)-dimensional vector rho, which is the maximum
#' the correlation between the adjacent variables
#' @export
rhoMax1dep <- function(p){
  m <- length(p)
  rholimit <- rep(1,m-1)

  rholimit[1] <- sqrt(p[2]*(1-p[1])/p[1]/(1-p[2]))
  for (i in 2:(m-1)){
    rholimit[i] <- rhoMaxMdep1(p[i-1],p[i],p[i+1])
  return( rholimit )}


Try the CorBin package in your browser

Any scripts or data that you put into this service are public.

CorBin documentation built on Nov. 15, 2020, 1:06 a.m.