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.

1 | ```
JointRegBC(ini = NA, X, y, z, p, q, ...)
``` |

`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 |

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.

`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. |

Supportted by Shahid Beheshti University

Bahrami Samani and Zhale Tahmasebinejad

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 | ```
data("Bahrami1")
gender<-Bahrami1$ GENDER
age<-Bahrami1$AGE
duration <-Bahrami1$ DURATION
y<-Bahrami1$ STEATOS
z<-Bahrami1$ BMI
sbp<-Bahrami1$ SBP
X=cbind(gender,age,duration ,sbp)
P<-lm(z~X)[[1]]
names(P)<-paste("Con_",names(P),sep="")
Q<-clogit(y~X)[[1]]
names(Q)<-paste("Binary",names(Q),sep="")
W=c(cor(y,z),var(z))
names(W)=c("Corr","Variance of Continous Response")
ini=c(P,Q,W)
p=5;
q=4;
JointRegBC(ini,X=X,y=y,z=z,p=p,q=q)
## The function is currently defined as
structure(function (x, ...)
UseMethod("JointRegBC"), class = "JointRegBC")
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.