Performs leave-some-out measures for a regimix model. This includes a measure of how much effect leaving out an observation has on the probability of each site's RCP label. Also, this function can be used as a cross-validation workhorse.

` model` |
a regmix object whose fit you want to assess |

` ...` |
ignored |

` oosSize` |
the size of the withheld partitions (out-of-sample size). Use 1 (default) for leave-one-out statistics, such as Cook's distance and leave-one-out validation. |

` times` |
the number of times to perform the re-estimation (the number of leave out groups). For each 1:times a random partition of the data, of size oosSize, is taken and the model is fitted to one of the partitions. It is predicted to the other partition. The exception is when oosSize=1 and times=model$n (leave-one-out). In such cases (the default too), the observations are left out one-by-one and not randomly. |

` mc.cores` |
the number of cores to spread the workload over. Default is 1. Argument is useless on Windows machines – see |

` quiet` |
should printing be suppressed? Default is no, it should not. Note that in either case, printing of the iteration trace etc is suppressed for each regimix fit. |

An object of class `regiCooksD`

. It is a list of 4 elements:

` Y` |
the species data, |

` CV` |
the model$n by model$S by times array of out-of-sample predictions (this array contains a lot of NAs for where predictions would in-sample), |

` cooksD` |
a model$n by model$nRCP matrix of statistics that resemble Cook's distance. The statistic is the change in the prediction of RCP probability from the model with all the data to the model with only the in-sample data, and |

` predLogL` |
the predictive log-likelihood of each point in each withheld sample (log-likelihood contributions of withheld observations, again there will be many NAs). |

cooks.distance(model, ..., oosSize = 1, times = model$n, mc.cores = 1, quiet = FALSE, type = "cooksD")

`regimix`

`stability.regimix`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ```
## Not run:
#not run as R CMD check complains about the time taken.
#This code will take a little while to run (about 0.5 minutes on my computer)
#For leave-one-out cooks distance, use oosSize=1
#for serious use times will need to be larger.
system.time({
example( regimix);
cooksD <- cooks.distance( fm, oosSize=10, times=25)
})
## End(Not run)
example( regimix) #will fit a model and store in fm
#for serious use times will need to be larger.
#For leave-one-out cooks distance, use oosSize=1
cooksD <- cooks.distance( fm, oosSize=10, times=5)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.