Sample cells with replacement to find bootstrapped distribution of coefficients
pbootVcov1(cl, zlmfit, R = 99)
bootVcov1(zlmfit, R = 99, boot_index = NULL)a cluster object created by makeCluster
class ZlmFit
number of bootstrap replicates
list of indices to resample. Only one of R or boot_index can be offered.
array of bootstrapped coefficients array of bootstrapped coefficients
pbootVcov1: parallel version of bootstrapping
data(vbetaFA)
zlmVbeta <- zlm(~ Stim.Condition, subset(vbetaFA, ncells==1)[1:5,])
#>
#> Done!
#Only run 3 boot straps, which you wouldn't ever want to do in practice...
bootVcov1(zlmVbeta, R=3)
#> , , (Intercept), C
#>
#> B3GAT1 BAX BCL2 CCL2 CCL3
#> [1,] 18.01107 17.51153 17.29382 20.91156 21.30163
#> [2,] 18.19441 17.44289 17.26292 20.23155 19.48139
#> [3,] 18.19441 17.44753 16.74632 19.63031 21.79991
#>
#> , , Stim.ConditionUnstim, C
#>
#> B3GAT1 BAX BCL2 CCL2 CCL3
#> [1,] -1.613563 -0.6570005 NA -4.300377 NA
#> [2,] NA -0.3006537 0.4797757 -3.620360 NA
#> [3,] -1.796906 -0.9377056 0.9963785 -3.019123 NA
#>
#> , , (Intercept), D
#>
#> B3GAT1 BAX BCL2 CCL2 CCL3
#> [1,] -4.643666 -0.9642210 -3.080827 -4.164635 -2.773660
#> [2,] -5.177913 -0.7923866 -2.773319 -3.383684 -2.632035
#> [3,] -5.068849 -0.9123317 -3.091916 -3.628102 -3.116119
#>
#> , , Stim.ConditionUnstim, D
#>
#> B3GAT1 BAX BCL2 CCL2 CCL3
#> [1,] 0.9194129 1.1834269 -2.27584386 0.4781265 -2.493911
#> [2,] -0.9609553 0.1372586 -1.34213135 0.1837434 -2.568610
#> [3,] 0.7573917 0.8149932 0.03162636 -0.4806393 -2.139224
#>