# JointRegBC.default: Joint Modelling of Mixed Correlated Binary and Continuous... In JointRegBC: Joint Modelling of Mixed Correlated Binary and Continuous Responses : A Latent Variable Approach

## Description

A joint regression model for mixed correlated binary and continuous responses is presented. In this model binary response can be dependent on the continuous response. With this model, the dependence between responses can be taken into account by the correlation between errors in the models for binary and continuous responses.

## Usage

 ```1 2``` ```## Default S3 method: JointRegBC(ini = NA, X, y, z, p, q, ...) ```

## Arguments

 `ini` Initial values `X` Design matrix `z` Continuous responses `y` Binary responses `p` Order of dimension of Binary responses `q` Order of dimension of continuous responses `...` Other arguments

## Details

Models for JointRegBC are specified symbolically. A typical model has the form response1 ~ terms and response2 ~ terms where response1and response2 are the (numeric) binary and continuous responses vector and terms is a series of terms which specifies a linear predictor for responses. A terms specification of the form first + second indicates all the terms in first together with all the terms in second with duplicates removed. A specification of the form first:second indicates the set of terms obtained by taking the interactions of all terms in first with all terms in second. The specification first*second indicates the cross of first and second. This is the same as first + second + first:second.

## Value

 `Binary response` Coefficient of ordinal response `Continuous Response` Coefficient of continuous response `Variance of Countinuous Response` Variance of continuous response `Correlation` Coefficient of continuous response `Hessian` Hessian matrix `convergence` An integer code. 0 indicates successful convergence. `objective` -loglikelihood.

## Note

Supportted by Shahid Beheshti University

## Author(s)

Ehsan Bahrami Samani and Zhale Tahmasebinejad

## References

Bahrami Samani, E. and Tahmasebinejad. Zh.(2011). Joint Modelling of Mixed Correlated Nominal, Ordinal and Continuous Responses. Journal of Statistical Research. 45(1):37-47.

`nlminb`,`fdHess,clogit`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88``` ```function (ini = NA, X, y, z, p, q, ...) { options(warn = -1) f <- function(ini, X, y, z, p, q) { X = cbind(1, X) y <- as.vector(y) z <- as.vector(z) ini <- as.vector(ini) X <- as.matrix(X) n = nrow(X) muz = muy = muygivenzx = q2 = q1 = l1 = l2 = l3 = muygivenzx = as.vector(0) sez <- ini[p + q + 2] seygivenzx <- (1 - (ini[p + q + 1])^2) mz=matrix(0,n,p) my=matrix(0,n,q) for(i in 1:n){ for(j in 1:p){ mz[i,j]=ini[1:p][[j]]*X[i, ][[j]] }} for(i in 1:n){ for(k in 1:q){ my[i,k]=ini[(p + 1):(p + q)][[k]]*X[i, -1][[k]] }} for (i in 1:n) { muz[i] <- sum(mz[i,]) muy[i] <- sum(my[i,]) muygivenzx[i] <- muy[i] + (ini[p + q + 1] * (z[i] - muz[i]))/sez q1[i] <- ( - muygivenzx[i])/sqrt(seygivenzx) l1[i] <- log(pnorm(q1[i])) + log(dnorm(z[i], muz[i], sez)) l2[i] <- log(1 - pnorm(q1[i])) + log(dnorm(z[i], muz[i], sez)) } data0 <- cbind(y, l1) data1 <- cbind(y, l2) data0[data0[, 1] == 1, 2] <- 0 data1[data1[, 1] == 0, 2] <- 0 t0 <- sum(data0[, 2]) t1 <- sum(data1[, 2]) t <- c(t0, t1) Tfinal <- sum(t) return(-Tfinal) } n = nlminb(ini, f, X = X, y = y, z = z, p = p, q = q, lower = c(rep(-Inf, p+q), -0.999, 0), upper = c(rep(Inf, p+q), 0.999, Inf), hessian = T) h = fdHess(n\$par, f, z = z, y = y, X, p, q) h1 = h\$Hessian ih = ginv(h1) se = sqrt(abs(diag(ih))) n\$Hessian <- h1 n\$p <- p n\$q <- q n\$se <- as.vector(se) n\$call <- match.call() class(n) <- "JointRegBC" object = n Co.Re <- data.frame(Parameter = object\$par[1:p], S.E = object\$se[1:p], `Confidence Interval` = paste("(", round(object\$par[1:p] - 2 * object\$se[1:p], 3), ",", round(object\$par[1:p] + 2 * object\$se[1:p], 3), ")", sep = "")) Binary.Re <- data.frame(Parameter = object\$par[(p + 1):(p + q)], S.E = object\$se[(p + 1):(p + q)], `Confidence Interval` = paste("(", round(object\$par[(p + 1):(p + q)] - 2 * object\$se[(p + 1):(p + q)], 3), ",", round(object\$par[(p + 1):(p + q)] + 2 * object\$se[(p + 1):(p + q)], 3), ")", sep = "")) Cor <- data.frame(Parameter = object\$par[p + q + 1], S.E = object\$se[p + q + 1], `Confidence Interval` = paste("(", round(object\$par[p + q + 1] - 2 * object\$se[p + q + 1], 3), ",", round(object\$par[p + q + 1] + 2 * object\$se[p + q + 1], 3), ")", sep = "")) Var <- data.frame(Parameter = object\$par[p + q + 2], S.E = object\$se[p + q + 2], `Confidence Interval` = paste("(", round(object\$par[p + q + 2] - 2 * object\$se[p + q + 2], 3), ",", round(object\$par[p + q + 2] + 2 * object\$se[p + q + 2], 3), ")", sep = "")) res <- list(call = object\$call, `Continuos Response` = Co.Re, `Variance Of Countinous Response` = Var, `Binary Response` = Binary.Re, Correlation = Cor) res\$Hessian <- h1 res\$convergence <- n\$convergence res\$objective<- n\$objective res\$call <- match.call() class(res) <- "JointRegBC" res } ```