Functions for detecting replicates and for calculating sample variance at specific design points
varPerReps(X, Y) estimateNugget(X, Y) anyReps(X)
the design matrix
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
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
Garrett M. Dancik firstname.lastname@example.org
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.