mlegp-nugget-related | R Documentation |
Functions for detecting replicates and for calculating sample variance at specific design points
varPerReps(X, Y) estimateNugget(X, Y) anyReps(X)
X |
the design matrix |
Y |
a vector (or 1 column matrix) of observations |
varPerReps
returns a 1-column matrix where element i corresponds to the sample variance in observations corresponding to design point X[i]
estimateNugget
returns a double calculated by taking the mean of the matrix returned by varPerReps
anyReps
returns TRUE
if two or more rows of X
are identical
These functions are used by mlegp
to set an initial value of the nugget when a constant nugget is being estimated. The function varPerReps
may also be useful for specifying the form of the nugget matrix for use with mlegp
.
Garrett M. Dancik dancikg@easternct.edu
https://github.com/gdancik/mlegp/
x = matrix(c(1,1,2,3,3)) # the design matrix y = matrix(c(5,6,7,0,10)) # output anyReps(x) varPerReps(x,y) estimateNugget(x,y)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.