get_top_pattern function

get top patterns from a slice of bronze-standard measurement

get_top_pattern(BrS_dat, Y, case_status, n_pat, exclude_missing = TRUE)

Arguments

• BrS_dat: bronze-standard data, which is usually data_nplcm$Mobs$MBS[[1]]
• Y: A vector of case/control status: 1 for case; 0 for control
• case_status: 1 for case; 0 for controls
• n_pat: the number of top patterns one wants to show
• exclude_missing: DEFAULT is TRUE for excluding any individual with missing measurements.

Returns

a list of results: obs_pat - observed rates; pattern_names; exist_other - if actual no. of patterns is larger than n_pat; N- No. of individuals with Y = case_status.

Examples

data(data_nplcm_noreg)
get_top_pattern(data_nplcm_noreg$Mobs$MBS[[1]],data_nplcm_noreg$Y,1,5,FALSE)

data(data_nplcm_noreg)
get_top_pattern(data_nplcm_noreg$Mobs$MBS$MBS1,data_nplcm_noreg$Y,case_status=1,n_pat=5)

See Also

Other exploratory data analysis functions: plot_logORmat(), show_individual(), summarize_BrS(), summarize_SS(), visualize_season()