Description Usage Arguments Details Value Author(s) References See Also Examples

Main initialization functions.

1 2 3 |

`PARAM` |
a list of uninitialized parameters, as usual, the returned
values of |

`initial.em.gbd()`

takes in a template of `PARAM`

(uninitialized),
and usually is available by calling `set.global()`

, then
return an initialized `PARAM`

which is ready for EM runs.

Internally, there are six different initializations implemented for
the function `initial.em.gbd()`

including `prob.extend`

,
`prob.simple`

, `qnorm.extend`

, `qnorm.simple`

, `extend`

,
and `simple`

. These methods are mainly based on transformation of
original space of data (p-values and voxel locations) into more linear
space such that the Euclidean distance more makes sense (fairly) to
classify data in groups.

`initial.RndEM.gbd()`

implement RndEM initialization algorithm
based on repeated calling `initial.em.gbd()`

.

Note that all configurations are included in `PARAM`

set by
`set.global()`

.

These functions return an initialized `PARAM`

for EM runs based on
pre-stored configuration within the input uninitialized `PARAM`

.

Wei-Chen Chen and Ranjan Maitra.

http://maitra.public.iastate.edu/

`set.global()`

, `fclust()`

, `PARAM`

.

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 | ```
library(MixfMRI, quietly = TRUE)
library(EMCluster, quietly = TRUE)
# .FC.CT$algorithm <- "em"
# .FC.CT$model.X <- "V"
# .FC.CT$ignore.X <- TRUE
.FC.CT$check.X.unit <- FALSE
### Test toy1.
set.seed(1234)
X.gbd <- toy1$X.gbd
PV.gbd <- toy1$PV.gbd
PARAM <- set.global(X.gbd, PV.gbd, K = 2)
PARAM.new <- initial.em.gbd(PARAM)
PARAM.toy1 <- em.step.gbd(PARAM.new)
id.toy1 <- .MixfMRIEnv$CLASS.gbd
print(PARAM.toy1$ETA)
RRand(toy1$CLASS.gbd, id.toy1)
### Test toy2.
set.seed(1234)
X.gbd <- toy2$X.gbd
PV.gbd <- toy2$PV.gbd
PARAM <- set.global(X.gbd, PV.gbd, K = 3)
PARAM.new <- initial.em.gbd(PARAM)
PARAM.toy2 <- em.step.gbd(PARAM.new)
id.toy2 <- .MixfMRIEnv$CLASS.gbd
print(PARAM.toy2$ETA)
RRand(toy2$CLASS.gbd, id.toy2)
``` |

MixfMRI documentation built on April 26, 2018, 5:03 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.