silver-standard data summary
Arguments
- SS_dat
a data frame of silver-standard data. It can usually be obtained by
data_nplcm$Mobs$MSS[[1]], meaning the first SS measurement slice.- Y
a vector of case control status: 1 for case; 0 for control.
See also
Other exploratory data analysis functions:
get_top_pattern(),
plot_logORmat(),
show_individual(),
summarize_BrS(),
visualize_season()
Examples
data(data_nplcm_noreg)
summarize_BrS(data_nplcm_noreg$Mobs$MBS[[1]], data_nplcm_noreg$Y)
#> $Nd
#> [1] 300
#>
#> $Nu
#> [1] 300
#>
#> $MBS_mean
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> case 0.60 0.31333 0.37000 0.16333 0.11667 0.13667
#> control 0.21 0.21667 0.23667 0.12000 0.13333 0.13000
#>
#> $MBS_q1
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> case 0.5436223 0.2634318 0.3173072 0.12559292 0.08477613 0.10209823
#> control 0.1675776 0.1736517 0.1919746 0.08764294 0.09919185 0.09629289
#>
#> $MBS_q2
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> case 0.6538491 0.367955 0.4259800 0.2095867 0.1582502 0.1804223
#> control 0.2597553 0.266846 0.2880173 0.1619657 0.1767463 0.1730629
#>
#> $MBS_nm
#> [1] "A" "B" "C" "D" "E" "F"
#>
summarize_SS(data_nplcm_noreg$Mobs$MSS[[1]], data_nplcm_noreg$Y)
#> $count
#> [1] 28 7
#>
#> $NA_count
#> A B
#> 0 0
#>
#> $Nd
#> [1] 300
#>
#> $MSS_mean
#> [,1] [,2]
#> [1,] 0.09333 0.02333
#>
#> $MSS_q1
#> [,1] [,2]
#> [1,] 0.0649724 0.01037832
#>
#> $MSS_q2
#> [,1] [,2]
#> [1,] 0.1319772 0.04834131
#>
#> $MSS_nm
#> [1] "A" "B"
#>