similar to: Computing the rank of a matrix.

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! I want to plot three graphs (residuals, ACF and PACF of a model). Ideally I would use a c(2,2) disposition where the residuals plot would start at position 1,1 and span to position 1,2. Then I would plot the ACF in position 2,1 and the PACF in position 2,2. Maybe is clearer like this: -------------------------- | | | residuals | |

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'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

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

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

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

Dear sir, we would like to know if there exist any R package containing the computational performance of the n-fold Stieljes' convolution of functions. We look forward to hearing from you. Thank you in advance. ____________________________________ M.Luz G?miz P?rez Dpt. Estad?stica e Investigaci?n Operativa Facultad de Ciencias Universidad de Granada Telf.: 958-243156 e-mail:

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

Hi. Sorry: my level of english is not enough good. I hope you understand what I mean. My question is: it's possible to share a subdirectory (and only one subdirectory for every user) in the section [homes]? The version of samba which i'm working is 2.2.7a, over Solaris 8. Thanks in advance. --

Hello! At the University of Granada (Spain) we use a Samba Server for aprox. 1000 users and runs ok. But (there's always a but) we need to know how encript password for the smbpasswd file. I know that we can use the smbpasswd program but we need get the encrypted password in the stdout. Has someone a program than does this?? TIA.

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

Buenas tardes a todos. Mi problema es el siguiente: Calculo una muestra de unos índices (s) y a partir de esa muestra elijo aquellos elementos de una tabla de datos que tengo (en este caso XP4) que cumplen la condición de que el índice de XYP sea el de la muestra. Ahora comienza mi problema, la única forma que he encontrado para que me R me tome todos los elementos que necesito es creando un

I tried to intall R in a iMac with Mac OS E1- 9.0.4 (Mac OS ROM VE-4.6.1). When I tried to run R for the first time I got a message about "Carbonlib-- CF RangeMake" not beeing found. Please, could you help me to tackle the problem. Thanks very much. Barahona -- Enrique Barahona Fernandez ESTACION EXPERIMENTAL DEL ZAIDIN, CSIC Profesor Albareda, 1 18008, GRANADA (Spain) Tel. + 34

Hi, I'm trying to use a similarity matrix (triangular) as input for pam() or fanny() clustering algorithms. The problem is that this algorithms can only accept a dissimilarity matrix, normally generated by daisy(). However, daisy only accept 'data matrix or dataframe. Dissimilarities will be computed between the rows of x'. Is there any way to say to that your data are already a

Dear all, I am using function glmer from package lme4 to fit a generalized linear mixed effect model. My model is as follows: model1 <- glmer(fruitset ~ Dist*wire + (1|Site), data, binomial) summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: fruitset ~ Dist * wire + (1 | Site) Data: data AIC BIC logLik deviance 68.23 70.65 -29.11 58.23 Random

I'm preparing some code to compute the optimal geometry of stressed solids. The core of the calculations is the optimization of elastic energy using the simulated annealing method implemented in the R optim() rutine. I've defined a function to compute this "energy" scalar (the fn parameter for optim) and prepared a list with the arrays defining the geometry and the elastic

Las dos Enviado desde mi iPhone > El 28/11/2013, a las 18:07, "Carlos J. Gil Bellosta " <cgb en datanalytics.com> escribió: > > ¿Numérica o simbólicamente? > > Un saludo, > > Carlos J. Gil Bellosta > http://www.datanalytics.com > > El día 28 de noviembre de 2013 18:03, jmcontreras <jmcontreras en ugr.es> escribió: >> Hola a todos

Hola: El 02/07/14 a las #4, Jose Luis Cañadas Reche escribió: > [...] > Ya de paso pregunto si alguien está interesado en > hacer un GIL (Grupo de Interés Local) en Andalucía.. ¿En que zona tienes interés? En Almería no hay muchos eRReros pero alguno que otro estamos. Salud y Revolución. Lobo. -- Libertad es poder elegir en cualquier momento. Ahora yo elijo GNU/Linux, para no