R: Display the criteria associated with a collection of...

plot.PLNmixturefamily {PLNmodels} | R Documentation

## Display the criteria associated with a collection of PLNmixture fits (a PLNmixturefamily)

### Description

Display the criteria associated with a collection of PLNmixture fits (a PLNmixturefamily)

### Usage

    ## S3 method for class 'PLNmixturefamily'
    plot(
      x,
      type = c("criteria", "diagnostic"),
      criteria = c("loglik", "BIC", "ICL"),
      reverse = FALSE,
      ...
    )

### Arguments

`x` | an R6 object with class `PLNmixturefamily`
`type` | a character, either `"criteria"` or `"diagnostic"` for the type of plot.
`criteria` | vector of characters. The criteria to plot in c("loglik", "BIC", "ICL"). Default is c("loglik", "BIC", "ICL").
`reverse` | A logical indicating whether to plot the value of the criteria in the "natural" direction (loglik - 0.5 penalty) or in the "reverse" direction (-2 loglik + penalty). Default to FALSE, i.e use the natural direction, on the same scale as the log-likelihood.
`...` | additional parameters for S3 compatibility. Not used

### Details

The BIC and ICL criteria have the form 'loglik - 1/2 * penalty' so that they are on the same scale as the model log-likelihood. You can change this direction and use the alternate form '-2*loglik + penalty', as some authors do, by setting `reverse = TRUE`.

### Value

Produces either a diagnostic plot (with `type = 'diagnostic'`) or the evolution of the criteria of the different models considered (with `type = 'criteria'`, the default).

### Examples

    data(trichoptera)
    trichoptera <- prepare_data(trichoptera$Abundance, trichoptera$Covariate)
    myMixtures <- PLNmixture(Abundance ~ 1 + offset(log(Offset)),
               data = trichoptera, control_main = list(smoothing = "forward", iterates = 1))
    plot(myMixtures, reverse = TRUE)

[Package PLNmodels version 0.11.6 Index]