Description Usage Arguments Details Value Author(s) See Also Examples
Sumamrize and provide p-value after association test of 2 or more variables. Percentages are reported for discrete variables. Users can choose the type of test they want. For continuous variable median and range are reported if the test specifid is non parametric, mean and standard deviation are reported for a parametric test. For discrete categories we perform Fisher's Exact Test, for continous following options are available - parametric - Students' T-Test or one way ANNOVA, for non-paramteric Wilcoxon Rank sum test or Kruskal Wallis test.
1 | get.summary2(fac,var,type,test.type="np",skip.test=FALSE, var.n=NULL)
|
fac |
vector, a dicrete vector of interest for which association is to be tested against other categorical or continuous variables |
var |
vector, can be discrete or continuous |
type |
|
test.type |
possible values - |
skip.test |
logical, Default is |
var.n |
character, If a variable name is provided, it will add this to output. |
Summary statistics are rounded to nearest 2 decimal points.
Note that length of fac
and var
should be the same. For Wilcoxon Rank sum test ties are handled through wilcox.test
function in coin
package.
workspace parameter in Fisher's Exact test is set to workspace=2e8
. All p-values < 0.0001 are replaced with a character string P<0.0001 or otherwise rounded to three decimal points.
fmat |
returns a character matrix summarixing |
Arshi Arora
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 | #two categorical variables
set.seed(100)
fac<-sample(1:2, 50, replace=TRUE)
var<-sample(1:3, 50, replace=TRUE)
get.summary2(fac,var, type=1,var.n="Example, Type=1" )
#Performing Fisher's Exact test
#| |1(n=23;46%) |2(n=27;54%) |RowTotal |pval |
#|:---------------|:-----------|:-----------|:--------|:------|
#|Example, Type=1 | | | |0.5971 |
#|1 |9(39%) |8(30%) |17(34%) | |
#|2 |7(30%) |12(44%) |19(38%) | |
#|3 |7(30%) |7(26%) |14(28%) | |
#continuous with two categories of fac
var<-rnorm(50)
get.summary2(fac, var, type=2, test.type="p", var.n="2-Way, test.type=p")
#Performing T-Test
#| |1 |2 |RowTotal |pval |
#|:------------------|:-----------|:-----------|:--------|:-----|
#|2-Way, test.type=p |-0.11(1.42) |-0.05(0.99) | |0.855 |
#|NA |0 |0 | | |
get.summary2(fac, var, type=2, test.type="np", var.n="2-Way, test.type=np")
#Performing Wilcoxon Rank Sum test
# 1 2 RowTotal pval
#2-Way, test.type=np "-0.1[-2.27-2.58]" "-0.07[-1.93-1.73]" "" "0.676"
#NA "2.58" "1.73" "" ""
fac<-sample(1:3, 50, replace=TRUE)
get.summary2(fac, var, type=2, test.type="p",var.n="Example, test.type=p")
#Performing one way ANOVA
#| |1 |2 |3 |RowTotal |pval |
#|:--------------------|:-----------|:----------|:---------|:--------|:-----|
#|Example, test.type=p |-0.46(1.31) |-0.16(0.9) |0.38(1.3) | |0.132 |
#|NA |0 |0 |0 | | |
get.summary2(fac, var, type=2, test.type="np",var.n="Example, test.type=np")
#Performing Kruskal-Wallis test
#| |1 |2 |3 |RowTotal |pval |
#|:---------------------|:-----------------|:-----------------|:----------------|:--------|:-----|
#|Example, test.type=np |-0.53[-2.27-2.58] |-0.07[-1.93-1.82] |0.42[-1.74-2.45] | |0.131 |
#|NA |2.58 |1.82 |2.45 | | |
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.