similar to: problem with optim: (list) object cannot be coerced to 'double'

Hello, I am currently testing a Vector AR of dim 3 over not a lot of data (135 * 3 observations) . To test the adequation of my vecot ar, I use the Schwarz Bayesian Criterion and the classic modified Portmanteau test on the residuals (it can be found for instance in http://www.iue.it/PUB/ECO2004-8.pdf , page 15) -> the null hypothesis is "the residuals process are a vectorila white

Hey all, I''ve had these functions for some time now, and would like to offer them as two new dollar-sign functions - elements to be extended by Prototype geniuses. :-) makeText(string) as $T() - return text node element Does just what it says... I''m sure someone could extend it nicely when via Prototype. (example) var x = $T(''hello world'');

When I do t.test on two distributions (see example below), it outputs numerous data about the t.test. What I'd like to do is individually capture some of this data and assign it to other variables. However, I am unable to find anything in the help section. In the example below, the t value is -4.0441 and the p-value is 0.006771 How can I assign these values to two variables, let's

---------- Forwarded message ---------- Date: Wed, 07 Jul 2004 00:31:21 -0400 From: Declan McCullagh <declan@well.com> To: politech@politechbot.com Subject: [Politech] House bill exports analog phone regs to VoIP http://www.politechbot.com/docs/boucher.voip.bill.070604.pdf There's a new bill in the House of Representatives to regulate phone calls made over the Internet. It was

Good day, all, I'm trying to transfer a 174M file to a system with 1.2G free. Other files in the tree come over just fine, but this transfer dies partway through: rsync -avvz -e ssh --partial --progress server.with.the.file:/server/directory /local/directory opening connection using ssh server.with.the.file rsync --server --sender -vvlogDtprz --partial . /local/directory receiving file

Dear All, 1)Can I use lm() to fit more than one response in single expression. e.g data is a matrix of these variables R1 R2 R3 X Y Z 1 2 1 1 2 3 .... Now i wnat to fit R1:R3 ~ X+Y+Z. 2) How can i use Singular Value decomposition (SVD) as an alternate to lsq. Regards,

When I import the library timeSeries I get (at least) the variable USDCHF imported too. I would like to delete it, but I cannot. As you can see below. Clearly I am doing something wrong. What is it? > library(timeSeries) Loading required package: timeDate > class(USDCHF) [1] "timeSeries" attr(,"package") [1] "timeSeries" > rm(list=c("USDCHF"))

Record plot which stores a plot to an internal R data structure and Replay Plot which replays the plot enables one to keep plots around and use them accross sessions at least until version 3. If you cannot restore old plots and treat them as data and save accross sessions, there is little reason for recordPlot. Either R needs to dispense with internal formats entirelly, a bad move in my opinion

I?m wondering which is the most efficient (time, than memory usage) way to obtain a multivariate time series object from a data frame (the easiest data structure to get data from a database trough RODBC). I have a starting point using timeSeries or xts library (these libraries can handle time zones), below you can find code to test. Merging parallelization (cbind) is something I?m thinking at

I retrieve for a few hundred times a group of time series (10-15 ts with 10000 values each), on every group I do some calculation, graphs etc. I wonder if there is a faster method than what presented below to get an appropriate timeseries object. Making a query with RODBC for every group I get a data frame like this: > X ID DATE VALUE 14 3 2000-01-01 00:00:03 0.5726334

Greetings -- I've got some sensor data of the form t1_1, t1_2 t2_1, t2_2 ... tN_1,tN_2 -- time intervals measuring starts and stops of sensor activity. I'd like to see whether there's any regularity in it. Seems natural to consider these data timeseries -- except most of the timeseries packages and models assume regular ones, with a fixed frequency. I wonder what's a

Given S4 methods [ and [<-, how do I write the Rd-file usage clause for the latter one? What I have now is: \S4method{[}{TimeSeries,TimeDate,missing}(x, i, j, ..., drop) \S4method{[<-}{TimeSeries,TimeDate,missing,ANY}(x, i, j, ..., value) which results in the following output: ## S4 method for signature 'TimeSeries, TimeDate, missing': x[i, j, ..., drop]

i have time series of momentum signal. I want to get the date of each of the "-1" signal period. for example , the first period of -1 signal begins on 2005-9-21 and ends on 2005-9-28. 2nd period of -1 signal begins on 2005-09-30 and ends on 2005-10-28. Thx Cameron date Px 200MA Signals 2005-09-15 26.27 25.83865 1 2005-09-16 26.07 25.83275 1

Dear R Users and maintainers of packages zoo and timeSeries, I believe there is a recently introduced inconsistency between timeSeries and zoo which is causing a problem with rbind. I had previously reported that I was having problems with rbind in the following code: library(zoo) foo<-zoo(1,order.by=as.Date("2007-10-09")) bar<-zoo(2,order.by=as.Date("2007-10-10"))

Hi guys, I am logging data about my body (weight, body fat, blood pressure, ..) in a .csv file and would like to plot this as a time series. I uploaded the (noisified) .csv, you can see the link in the code I have so far (you can run the code directly as-is): df.raw <- read.csv("http://www.chaotic-neutral.de/temp/stats-noised.csv", sep=";", dec=".")

Hi there I have some timeseries data which I plot in a OHLC Plot. In the same plot I'd like to have the EMA of this timeseries. I tried to add the EMA point to OHLC with lines(), but this doesn't work. Has anyone an idea how to handle it? Regards, Michael Zak

Hi R users: When I try to merge or bind (cbind or rbind) two series, both with a "FinCenter" different that GMT, the result is "GMT" not the original financial center? What am I doing wrong? ###################################################### require(timeSeries) getRmetricsOptions("myFinCenter") setRmetricsOptions(myFinCenter = "America/Bogota")

Hi All, Does someone have an idea of how to cleverly convert a categorical timeseries into a transition matrix? Ie, I have something like: x<- c(1,1,2,1,1,2,2,2,1,2), And I want a matrix with counts and/or probabilities: > tr <- matrix(c(2,3,2,2),2,2) > tr [,1] [,2] [1,] 2 2 [2,] 3 2 Meaning that there are two transitions from 1 to 1, two from 1 to 2, three from 2 to 1

Hi, I am trying to get the following two timeseries (these are small subsets of the whole thing) into R so I can merge them using zoo. Timeseries 1=[ Date Count 9/28/2003 1505 10/5/2003 1535

Hi, I didn't find any post on this subject so I ll ask you some advices. Let's say that I have two library trees. Number 1 is the default R library tree on path1 Number 2 is another library tree on a server with all packages on path2. When I set library(aaMI,lib.loc=paths2) it loads the package even if its not on default R library When I set library(fOptions,lib.loc=paths2) it