Description Usage Arguments Details Value See Also Examples

`simcam`

is used to combine several one-class SIMCA models for multiclass classification.

1 | ```
simcam(models, info = "")
``` |

`models` |
list with SIMCA models ( |

`info` |
text with information about the the object. |

Besides the possibility for multiclass classification, SIMCAM also provides tools for investigation of relationship among individual models (classes), such as discrimination power of variables, Cooman's plot, model distance, etc.

When create `simcam`

object, the calibration data from all individual SIMCA models is
extracted and combined for making predictions and calculate performance of the multi-class model.
The results are stored in `$calres`

field of the model object.

Returns an object of `simcam`

class with following fields:

`models ` |
a list with provided SIMCA models. |

`dispower ` |
an array with discrimination power of variables for each pair of individual models. |

`moddist ` |
a matrix with distance between each each pair of individual models. |

`classnames ` |
vector with names of individual classes. |

`nclasses ` |
number of classes in the object. |

`info ` |
information provided by user when create the object. |

`calres ` |
an object of class |

Methods for `simca`

objects:

`print.simcam` | shows information about the object. |

`summary.simcam` | shows summary statistics for the models. |

`plot.simcam` | makes an overview of SIMCAM model with two plots. |

`predict.simcam` | applies SIMCAM model to a new data. |

`plotModelDistance.simcam` | shows plot with distance between individual models. |

`plotDiscriminationPower.simcam` | shows plot with discrimination power. |

`plotModellingPower.simcam` | shows plot with modelling power for individual model. |

`plotCooman.simcam` | shows Cooman's plot for calibration data. |

`plotResiduals.simcam` | shows plot with Q vs. T2 residuals for calibration data. |

Methods, inherited from `classmodel`

class:

`plotPredictions.classmodel` | shows plot with predicted values. |

`plotSensitivity.classmodel` | shows sensitivity plot. |

`plotSpecificity.classmodel` | shows specificity plot. |

`plotMisclassified.classmodel` | shows misclassified ratio plot. |

Since SIMCAM objects and results are calculated only for optimal number of components, there is no sense to show such plots like sensitivity or specificity vs. number of components. However they are available as for any other classification model.

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 | ```
## make a multiclass SIMCA model for Iris data
library(mdatools)
# split data
caldata = iris[seq(1, nrow(iris), 2), 1:4]
se = caldata[1:25, ]
ve = caldata[26:50, ]
vi = caldata[51:75, ]
testdata = iris[seq(2, nrow(iris), 2), 1:4]
testdata.cref = iris[seq(2, nrow(iris), 2), 5]
# create individual models
semodel = simca(se, classname = 'setosa')
semodel = selectCompNum(semodel, 1)
vimodel = simca(vi, classname = 'virginica')
vimodel = selectCompNum(vimodel, 1)
vemodel = simca(ve, classname = 'versicolor')
vemodel = selectCompNum(vemodel, 1)
# combine models into SIMCAM objects, show statistics and plots
model = simcam(list(semodel, vimodel, vemodel), info = 'Iris data')
summary(model)
plot(model)
# show predictions and residuals for calibration data
par(mfrow = c(2, 2))
plotPredictions(model)
plotCooman(model, nc = c(1, 2))
plotResiduals(model, nc = 1)
plotResiduals(model, nc = 2)
par(mfrow = c(1, 1))
# show different plots for the model
par(mfrow = c(2, 2))
plotModelDistance(model, nc = 1)
plotDiscriminationPower(model, nc = c(1, 2))
plotModellingPower(model, nc = 1)
plotModellingPower(model, nc = 2)
par(mfrow = c(1, 1))
# apply the SIMCAM model to test set and show statistics and plots
res = predict(model, testdata, testdata.cref)
summary(res)
plotPredictions(res)
``` |

mdatools documentation built on July 7, 2018, 1 a.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.