Description Usage Arguments Value Examples
This function estimates the parameters used for visualization of model-based clustering performs with R package Rmixmod. To achieve the parameter infernece, it automatically samples probabilities of classification from the model parameters
1 2 3 | clusvisMixmod(mixmodResult, sample.size = 5000, maxit = 10^3,
nbrandomInit = 4 * mixmodResult@bestResult@nbCluster, nbcpu = 1,
loccont = NULL)
|
mixmodResult |
[ |
sample.size |
numeric. Number of probabilities of classification sampled for parameter inference. |
maxit |
numeric. It limits the number of iterations for the Quasi-Newton algorithm (default 1000). |
nbrandomInit |
numeric. It defines the number of random initialization of the Quasi-Newton algorithm. |
nbcpu |
numeric. It specifies the number of CPU (only for linux). |
loccont |
numeric. Index of the column containing continuous variables (only for mixed-type data). |
Returns a list
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 51 52 53 54 | ## Not run:
## First example: R package Rmixmod
# Package loading
require(Rmixmod)
# Data loading (categorical data)
data("congress")
# Model-based clustering with 4 components
set.seed(123)
res <- mixmodCluster(congress[,-1], 4, strategy = mixmodStrategy(nbTryInInit = 500, nbTry=25))
# Inference of the parameters used for results visualization
# (specific for Rmixmod results)
# It is better because probabilities of classification are generated
# by using the model parameters
resvisu <- clusvisMixmod(res)
# Component interpretation graph
plotDensityClusVisu(resvisu)
# Scatter-plot of the observation memberships
plotDensityClusVisu(resvisu, add.obs = TRUE)
## Second example: R package Rmixmod
# Package loading
require(Rmixmod)
# Data loading (categorical data)
data(birds)
# Model-based clustering with 3 components
resmixmod <- mixmodCluster(birds, 3)
# Inference of the parameters used for results visualization (general approach)
# Probabilities of classification are not sampled from the model parameter,
# but observed probabilities of classification are used for parameter estimation
resvisu <- clusvis(log(resmixmod@bestResult@proba),
resmixmod@bestResult@parameters@proportions)
# Inference of the parameters used for results visualization
# (specific for Rmixmod results)
# It is better because probabilities of classification are generated
# by using the model parameters
resvisu <- clusvisMixmod(resmixmod)
# Component interpretation graph
plotDensityClusVisu(resvisu)
# Scatter-plot of the observation memberships
plotDensityClusVisu(resvisu, add.obs = TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.