Displaying 1 result from an estimated 1 matches for "silversignif".
2003 Jul 17
0
Silverman modality test
...nN1)>e)
{
modes <- nbmodes (x, winN)
winN1 <- winN
if (modes > n)
{
minw <- winN
winN <- (winN + maxw)/2
}
else
{
maxw <- winN
winN <- (winN - minw)/2
}
}
return (winN)
}
silversignif <- function (x, m, nboot=1000)
# Silverman's significance for the null that x has at
# most modes.
{
h <- hcrit (x, m)
h0 <- h
n <- 0
theta <- function (y, h0=h, sigma2=var(y))
{
(y + rnorm (1, mean=0, sd=h0^2)) / sqrt ( 1+ (h0^2/sigma2))
}
for (i in 1:nbo...