plot.betaclust {betaclust} | R Documentation |

Visualise a `betaclust`

clustering solution by plotting the fitted and kernel density estimates, the uncertainty and the information criterion.

```
## S3 method for class 'betaclust'
plot(
x,
what = "fitted density",
plot_type = "ggplot",
data = NULL,
sample_name = NULL,
title = NULL,
patient_number = 1,
threshold = FALSE,
scale_param = "free_y",
...
)
```

`x` |
A |

`what` |
The different plots that can be obtained are either "fitted density","kernel density","uncertainty" or "information criterion" (default = "fitted density"). |

`plot_type` |
The plot type to be displayed are either "ggplot" or "plotly" (default = "ggplot"). |

`data` |
A dataframe of dimension |

`sample_name` |
The names of DNA samples in the dataset analysed by the K.R model. If no value is passed then default values of sample names, e.g. Sample 1, Sample 2, etc are used as legend text (default = NULL). |

`title` |
The title that the user wants to display. If no title is to be displayed the default is "NULL". |

`patient_number` |
The column number representing the patient in the patient-wise ordered dataset selected for visualizing the clustering solution of the K.. or KN. model (default = 1). |

`threshold` |
The "TRUE" option displays the threshold points in the graph for the K.. and the KN. model (default = "FALSE"). |

`scale_param` |
The position scales can be fixed or allowed to vary between different panels generated for the density estimate plots for visualizing the K.R clustering solution. Options are "fixed", "free_y","free_x" or "free" (default = "free_y"). The option "fixed" results in the x and y scales being fixed across all panels, "free" varies the x and y scales across the panels, "free_x" fixes the y scale and lets the x scale vary across all panels and "free_y" fixes the x scale and lets the y scale vary across all panels. |

`...` |
Other graphics parameters. |

The fitted density estimates can be visualized under the optimal clustering solution by specifying what = "fitted density" and kernel density estimates under the optimal clustering solution by specifying what = "kernel density".

The threshold inferred can be visualized by specifying threshold = TRUE. The KN. model calculates different pairs of threshold points for each patient as the shape parameters are allowed to vary for each patient. So the patient for whom the threshold needs to be displayed can be specified by inputting the column number representing the patient in the patient-wise ordered dataset in the parameter patient_number.

Interactive plots can also be produced using plot_type = "plotly". The uncertainty in the clustering solution can be plotted using what = "uncertainty". The information criterion values for all fitted models can be plotted using what = "information criterion".

This function displays the following plots as requested by the user:

fitted density estimates - Plot showing the fitted density estimates of the clustering solution under the optimal model selected.

kernel density estimates - Plot showing the kernel density estimates of the clustering solution under the optimal model selected.

uncertainty - A boxplot showing the uncertainties in the optimal clustering solution.

information criterion - Plot showing the information criterion values for all models fitted to support the selection of the optimal model.

```
my.seed <- 190
M <- 3
N <- 4
R <- 2
data_output <- betaclust(pca.methylation.data[1:100,2:9], M, N, R,
model_names = c("K..","KN.","K.R"), model_selection = "BIC",
parallel_process = FALSE, seed = my.seed)
plot(data_output, what = "fitted density", plot_type = "ggplot",
data=pca.methylation.data[1:100,2:9])
plot(data_output, what = "uncertainty", plot_type = "ggplot")
plot(data_output, what = "information criterion", plot_type = "ggplot")
```

[Package *betaclust* version 1.0.0 Index]