Analysis of multiple traits can provide additional information beyond analysis of a single trait, allowing better understanding of the underlying genetic mechanism of a common disease. To accommodate multiple traits in familial correlation analysis adjusting for confounders, we develop a regression model for canonical correlation parameters and propose joint modeling along with mean and scale parameters. The proposed method is more powerful than the regression method modeling pairwise correlations because it captures familial aggregation manifested in multiple traits through maximum canonical correlation.