Displaying 1 result from an estimated 1 matches for "tp3477339p3477339".
2011 Apr 27
1
matrix of higher order differences
...diff(x, dif=i);
for(j in 1:length(differences)){
M[j,i+j] <- differences[j]
}
}
M[lower.tri(M)] <- t(M)[lower.tri(M)];
return(M);
}
x <- c(1,2,3,5,7,11,13,17,19);
diffmatrix(x);
--
View this message in context: http://r.789695.n4.nabble.com/matrix-of-higher-order-differences-tp3477339p3477339.html
Sent from the R help mailing list archive at Nabble.com.