library(Rfit)
x<-c(-4,-2,5,9)
x
walsh(x)
mywalsh<-function (x)
{
n <- length(x)
w <- vector(n * (n + 1)/2, mode = "numeric")
ind <- 0
for (i in 1:n) {
for (j in i:n) {
ind <- ind + 1
w[ind] <- 0.5 * (x[i] + x[j])
}
}
return(w)
}
mywalsh(x)
