Nothing
KEnvelope <-
function(X, r = NULL, NumberOfSimulations = 100, Alpha = 0.05,
ReferenceType = "", NeighborType = ReferenceType, SimulationType = "RandomPosition",
Precision = 0, Global = FALSE, verbose = interactive()) {
CheckdbmssArguments()
# Choose the null hypothesis
SimulatedPP <- switch (SimulationType,
RandomPosition = expression(rRandomPositionK(X, Precision=Precision, CheckArguments = FALSE)),
RandomLabeling = expression(rRandomLabeling(X, CheckArguments = FALSE)),
PopulationIndependence = expression(rPopulationIndependenceK(X, ReferenceType, NeighborType, CheckArguments = FALSE))
)
if (is.null(SimulatedPP))
stop(paste("The null hypothesis", sQuote(SimulationType), "has not been recognized."))
# local envelope, keep extreme values for lo and hi (nrank=1)
Envelope <- envelope(X, fun=Khat, nsim=NumberOfSimulations, nrank=1,
r=r, ReferenceType=ReferenceType, NeighborType=NeighborType,
CheckArguments = FALSE,
simulate=SimulatedPP, verbose=verbose, savefuns=TRUE
)
attr(Envelope, "einfo")$H0 <- switch (SimulationType,
RandomPosition = "Random Position",
RandomLabeling = "Random Labeling",
PopulationIndependence = "Population Independence"
)
# Calculate confidence intervals
Envelope <- FillEnvelope(Envelope, Alpha, Global)
# Return the envelope
return (Envelope)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.