similar to: Binding a matrix to a matrix

Hello. I have the following object which is a list of length NumSim with each entry being a matrix of dimensions Ncurvas x 3: dW = replicate(NumSim,cbind(rnorm(Ncurvas),rnorm(Ncurvas),rnorm(Ncurvas)),simplify=F) I would like to transform it into an array of dimension Ncurvas x 3 x NumSim. Does anybody does how to do this? or how to generate directly and array composed of independent random

Hello I am trying to import the attached file Curva LIBOR to R. I am trying to use the following commands and obtaining the following errors > res <- read.xlsx("C:\\Users\\FELIPE PARRA\\Documents\\Quantil\\Federacion\\Curva LIBOR.xlsx", 4) Error en .jcall(rowCells[[ic]], "I", "getColumnIndex") : RcallMethod: invalid object parameter > res <-

Hello all, I have been working with b-splines and noted that the splines package is not available in CRAN. Does any body know what happened with it? Or, is there any package that replaces it? Thank you Felipe Parra [[alternative HTML version deleted]]

Hello, I have stings which have all sort of characters (numbers, letters, punctuation marks, etc) I would like to stay only with the numbers in them, does somebody know how to do this? Thank you Felipe Parra [[alternative HTML version deleted]]

Hello I want to unlist the attached element getting only the first element in each element of the list. The last element of the list looks as this: [[5065]] [[5065]]$Pluv3Meses [1] 274.4 [[5065]]$PluvMesesMedio [1] 378.2667 [[5065]]$Pluv2UltimosMeses [1] 23.33333 So I would like to get for each element of the list the element called Pluv3Meses. The whole list has 5065 elements but when I try to

Hello I am new to RExcel and I would like to run a source code form the excel worksheet. I would like to run the following code source("C:\\Quantil Aplicativos\\Genercauca\\BackwardSelectionNC.r") from the excel wroksheet. Does anybody know how to do this? Thank you Felipe Parra [[alternative HTML version deleted]]

Hello, does somebody know in a boxplot, what does each element in the boxplot represent? 1. lines at the extremes of the dotted lines? 2. Extremes of the boxes 3. Black line in the middle of the box? 4. notches? Thank you Felipe Parra [[alternative HTML version deleted]]

Hello I am tryin to use the apply functions with two data frames I've got and I am getting the following error message Error en HistRio$SecSte : $ operator is invalid for atomic vectors I don't understand why. when I use the apply I am doing: PromP <- function(HistRio,AnaQuim){ xx <- c(0,0,0) if(length(which(AnaQuim$SecSte==HistRio$SecSte))>0){ xx[1]<-1 }

Hello, i have the following list strsplit(as.character(Elecciones$Municipios),"\\.") [[1]] [1] "ANTIOQUIA" "ABEJORRAL" [[2]] [1] "META" "ACACIAS" [[3]] [1] "CASANARE" "AGUAZUL" [[4]] and I would like to make a vector of the first element of each of the list items, in this case ANTIOQUIA, META, CASANARE, etc . Do you know how

Hello, I am runnning a program on R with a "big" number of simulations and I am getting the following error: Error: no se puede ubicar un vector de tamaño 443.3 Mb I don't understand why because when I check the memory status in my pc I get the following: > memory.size() [1] 676.3 > memory.size(T) [1] 1124.69 > memory.limit() [1] 4000 which should in theory allow to

Hello. I would like to know if there exists in R a function like trunc but where i can choose how many decimal places can I stay with in the number I have (sort of the same as the trunc function in excel). I would like, for example if I have the number 0.974678 and I choose to stay with 3 decimal places to have as ouput 0.974. Thank you Felipe Parra [[alternative HTML version deleted]]

Hello to all, I have been using an optim with the following call: optim(param_ini,fun_errores2,Precio_mercado=Precio,anos_pagosE2=anos_pagos,control=list(maxit=10000,reltol=1e-16)) depending on the intial values I'm getting the same solution but once I get the convergence message=10 (no convergence) and for the others I get convergence message = 0 Solution1: $par beta1

Hello I am trying to use predict.lm, but I am having trouble getting out of sample predictions. I am getting the same output if I use the following three commands: predict(ModeloLineal,predictors[721:768,]) predict(ModeloLineal,predictors[1:768,]) predict(ModeloLineal) where ModeloLineal is the output from ModeloLineal<-lm(dataTS[,6] ~ predictors[1:720,]), so the first 720 observations of

Hello, I am trying to create a panel with the attached data frame. using the following code: > PanelRio = DataRiopaila[which(duplicated(DataRiopaila$SEC_STE)==T),] > PanelRio=plm.data(PanelRio,index=c("SEC_STE","FechaSiembra")) series Pluv3Meses,PluvMes4al10,Pluv2UltimosMeses,Rad3Meses,RadMes4al10,Rad2UltimosMeses are NA and have been removed It tells me it removes

Hello. I am trying to use the portfolio backtesting function in fPortfolio package, but I don't now why in my version of fPortfolio I don't have either the portfolioBactest nor the portfolioBacktesting functions. Does anybody knows what might be going on? thank you Felipe Parra [[alternative HTML version deleted]]

Hello. I am using fGarch to estimate the following model: Call: garchFit(formula = fmla, data = X[, i], trace = F) Mean and Variance Equation: data ~ arma(1, 1) + garch(1, 1) Conditional Distribution: norm Coefficient(s): mu ar1 ma1 omega alpha1 beta1 -0.94934 1.00000 -0.23211 54.06402 0.45709 0.61738 Std. Errors: based on Hessian Error Analysis:

Hello I am trying to compare a vector with a Matrix's rows.The vector has the same length as the number of columns of the matrix, and I would like to find the row numbers where the matrix's row us the same as the given vector. What I am doing at the moment is using apply as follows: apply(Matrix,1,function(x)all(x%in%LHS)) but this isn't too fast actually. I would like to know if

Hello. I have the following funtion: fechasEntrega = function(FechasEntrega,fecha){ if(length(which(FechasEntrega<fecha))>0){ tkmessageBox(title = "Error en Fecha de Valoracion",message="Hay una fecha de entrega anterior a la fecha de valoracion. Todas las fechas de entrega deben ser posteriores a la fecha de valoración para el correcto funcionamiento del

Hello , I have estimated the following model, a sarima: p=9 d=1 q=2 P=0 D=1 Q=1 S=12 In R 2.12.2 Call: arima(x = xdata, order = c(p, d, q), seasonal = list(order = c(P, D, Q), period = S), optim.control = list(reltol = tol)) Coefficients: ar1 ar2 ar3 ar4 ar5 ar6 ar7 ar8 ar9 0.3152 0.8762 -0.4413 0.0152 0.1500 0.0001 -0.0413 -0.1811

Hello I am trying to import an Excel file but I am loosing the headers, My headers are in the first to rows of the EXCEL file. In the following R output, the NA are supposed to be the second item in the Header. Is there any way to Import more than one row as headers?. Thank you Felipe Parra > Data<-odbcConnectExcel('Curva LIBOR.xlsx',readOnly=T,) >