similar to: Plots spanning columns

Hi! Maybe this is a silly question, but I need the column rank (http://en.wikipedia.org/wiki/Rank_matrix) of a matrix and R function 'rank()' only gives me the ordering of the elements of my matrix. How can I compute the column rank of a matrix? Is there not an R equivalent to Matlab's 'rank()'? I've been browsing for a time now and I can't find anything, so any

Hi there! Still struggling to translate Matlab code into R's tsDyn package. Here is my question: Is there in R an equivalent function to Matlab's prctile()? To the moment I thought it was quantile(), but I just realized I was wrong. The definition of the Matlab function: prctile Percentiles of a sample SyntaxY = prctile(X,p) Description Y = prctile(X,p) calculates a value that is

Hi there! I need some assistance here with vector orderings. I have a set of q vectors of length p, grouped by rows in a matrix A, q?p, that I need to order lexicographically (http://en.wikipedia.org/wiki/Lexicographical_order). I also have another matrix B, p?r, and a vector c, that should be ordered according to the order of A. So far, I was doing ordering <-

Hi there! I can't find any information about creating combinations with repetitions in R. The function combn() does create combinations, but _without_ repetitions. Here is what I need to do: svIter <- 1000 xx <- matrix(rnorm(m*n), c(m, n)) sequence <- seq(range(xx)[1], range(xx)[2], length.out = svIter^(1/q)) expand.grid(secuence, secuence, .../q times/..., secuence) That is,

Hi! I've been browsing through the last months' archive and I can't find an answer to my question, so here it is (let's hope it's not too obvious): I'm working on extensions of an R library, and I would be very surprised if everyone developing R packages is doing the following, as I do: 1.- Write down a modification of an R file 2.- Exit the current R

Hi! I'm translating some code from Matlab to R and I found a problem. I need to translate Matlab's function 'lsqnonlin' (http://www-ccs.ucsd.edu/matlab/toolbox/optim/lsqnonlin.html) into R, and at the beginning I thought it would be the same as R's 'optim'. But then I looked at the definition of 'lsqnonlin' and I don't quite see how to make

Dear all, I'm having trouble using the 'apply' function to multiply some arrays. Let A be a 500x2 array and B a 500x3 array. I need a vector C which has 6 columns being the first three the result of A[,1] * B and the second three the result of A[,2] * B. What is the most efficient way to express that? I'm trying to use 'apply' with no success... Thanks for your

Hi again! I'm using the 'optim' method to fix the parameters of a model. I have written the function to be minimised and another function which returns the gradient of the error. Now my problem is that, in order to compute that gradient, it would be extremely convenient to be able to pass some parameters to the gradient function. I don't see how to do it given the fixed syntax

Estimados. Hace un tiempo que tengo una duda, estaba pensando en los problemas como optimización de costos, donde hay varias alternativas y librerías, pasando por soluciones inspiradas en energía, genética o algo matemático como matrices y álgebra. Luego aparecen tensorflow, cntk, y otros tantos, de los cuáles https://keras.rstudio.com/ ofrece alternativas para mezclar los mundos por decirlo de

Hi there! Happy new year to everyone! There's a piece of code that I must write that's driving me crazy. Maybe any of you has previous experiences in something similar. Any help will be greatly appreciated! The problem is as follows. I have a matrix Z with dimensions c(m, p). I have two vectors Gamma and Th with length p. What I want to do is, for j in 1:p, for each pair

Estimado Jesús Para Fernández En teoría es ese material, lo vi muy rápido y en la parte genética tiene cosas que biológicamente no son así, hay un libro de Falconer, Introducción a la genética cuantitativa, que tiene escrita la parte matemática, hay un abismo entre la biología y los ingenieros que se inspiran en la biología. Yo pensaba en la resolución de un problema real, aunque relativamente

pacf in devel seems by default to return a different number of lags than 1.7.1 for $pacf. I don't see any mention of this in the NEWS file, or any change in the documentation, so I suspect it is and error, though it may be an undocumented improvement. (Newbie question: How is the simplest way to display a function like pacf.default that is not exported from a namespace?) Paul

hi, in 1.9.1, the return value from ARMAacf(pacf=TRUE) is not named by lags, contrary to ?ARMAacf. the simple fix is to move names(Acf) <- down after if(pacf), with an appropriate starting lag as pacf=TRUE appears to start at lag 1 (whereas pacf=FALSE starts at lag 0). for consistency, one could argue to append 1 for lag 0 for pacf=TRUE (or start pacf=F at lag 1). however, given the

How can I make sure the residual signal, after subtracting the trend extracted through some technique, is actually trend-free ? I would greatly appreciate any suggestion about some Stationarity tests. I'd like to make sure I have got the difference between ACF and PACF right. In the following I am citing some definitions. I would appreciate your thoughts. ACF(k) estimates the correlation

I am getting the following: *** Weave Errors *** Error in driver$runcode(drobj, chunk, chunkopts) : Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" *** Source Errors *** Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" make[1]: *** [checkVignettes] Error 1 I don't really understand the new namespace mechanism,

Dear list, according to the documentation of acf{ts} "the partial correlation coefficient is estimated by fitting autoregressive models of successively higher orders up to lag.max. " However, R seems to return the Yule-Walker estimates of the PACF by default. You can check this using c(1:10) as the series: the YW estimates are 0.7000000 and -0.1527035 for lags 1 and 2 . If the PACF

When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]]

R: Is there an easy way to get the acf and pacf for an irregular times series? That is, the acf and pacf with lag lengths that are in units of time, not observation number. Thanks, Jason Higbee Research Associate Federal Reserve Bank of St. Louis The views expressed in this email are the author's and not necessarily those of the Federal Reserve Bank of St. Louis or the Federal Reserve

Estimados amigos Estoy dibujando las funciones acf y pacf de una variable de una serie "zoo": > ls.str(pat="T0.5") T0.5 : 'zoo' series from 2017-11-08 23:00:00 to 2017-11-15 06:59:00   Data: num [1:9120, 1:3] 55 49.8 51 50.1 36.5 ...   Index:  POSIXct[1:9120], format: "2017-11-08 23:00:00" "2017-11-08 23:01:00" "2017-11-08

Hi i am doing an environmental research The equation is as follow: gam(y1 ~ x1 + s(x2) + s(x3) + s(x4), family = gaussian, fit = true) I would like to obtain the beta coefficient and 95CI of x4 (or s(x4)), what should I do? Thanks, Lung -- View this message in context: http://r.789695.n4.nabble.com/mgcv-beta-coefficient-and-95-CI-tp3320491p3320491.html Sent from the R help mailing list