similar to: Need help with time series

acf wants a time series, so tries to make one: as.ts(myts) Time Series: Start = 19357 End = 20027 Frequency = 1 [1] 24957 NA NA NA NA NA NA NA NA NA [11] NA NA NA NA NA NA NA NA NA NA [21] NA NA NA NA NA NA NA NA NA NA [31] NA 10577 NA NA NA NA NA

Dear all, I'm indirectly faced with the fact that setting the 'tsp' attribute of an object modifies its class definition: > class( structure(1:2, tsp=c(1,2,1), class=c("myts","ts")) ) [1] "ts" "myts" In general, this is of really little (ok, I admit: totally no) interest for me because 'myts' class is added just after assigning the

I have defined a daily timeseries for the 365 days of 2003 issuing: myts = ts(dati[,2:10],frequency=365,) > myts Time Series: Start = c(1, 1) End = c(1, 365) Frequency = 365 and mytime = as.POSIXct(strptime(as.character(dati[,1]),format="%Y-%m-%d")) contains the dates from "2003-01-01" to "2003-12-31" How can I combine mytime and myts in order to list

hello, I got a problem with plotting large time series, since I want to store the results in a .PDF file (I want to store several pages of plots). The PDF files get too large to be handled (> 10MB, one was even 200MB big). So I wonder, if there would be a possibilty to either - reduce the file size of the PDF - change the way the plot is generated to reduce the plot size? I use:

Hi dear R-users, I have a question about a function I'm trying to improve. How can I stop the function calculation at the last numeric value of my data? The problem is that the end of my data contains missing values (NAs). And the aim of my function is to compare the first numeric value with the next one (till the end). For the moment, It works well when my data doesn't contains any NAs

Sorry to keep posting but I want to do this right and I'm hoping for some pointers I now have two time series objects which I need to subtract. Unfortunatly the two series dont have the same sample rates. When I try to subtract them avgSub<-avg1-avg2 The time series object is clever enough to object. So I guess I need to write a function for subtraction of the time series objects which

I tried to assign values to specific elements of a time series and got in trouble. The code below should be almost self-explanatory. I wanted to assign 0 to the first element of x, but instead I assigned zero to the second element of x, which is not what I wanted. Is there a function that will allow me to do this without going into index arithmetic (which would be prone to errors)? FS >

I was looking for suggestions as to how to smooth a timeseries and, having accomplished that, how to find the fitted curve values for intermediate points. ? I've tried numerous examples of possible approaches in R that I've found on the web, but when applied to my simple data, R returns an error message in numerous cases. ? The main problem seems to be that I have monthly data, starting

Dear R-help list members, I am new to R, and having problems with plotting part of a time series. I have read in my data using read.table, and my 'time' column is recognised as a numeric variable. When I convert this to a time format, I am no longer able to use its values on my x-axis. The problem is that I only want part of the time series to appear on the plot, ignoring the first 23

help . I have a file which contains a table that has 2 columns of which the first one is a time sereis. I need to extract the first column alone and do an acf plot on it. I dont know how to extract the first column and pass it to the acf function in the ts package. Can anybody pls help me. end -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list

Hi you guys: I have been wondering if there is any way to change the labeling in plot.ts( ), for example , if I plot two sequences, i always got y labels as "series1", "series2", I tried to use ylab=c((expression(mu_1)),(expression(mu_2))), but it doesn't work. thanks in advance for any help best.

The following code works fine as long as I don''t try to run it through the distributed server. It doesn''t get there... It runs fine out of delayed_job, runs fine if called directly. But if ''distrib'' is true (the default) it runs right up to the call to the server and right past it without getting to the server or raising any errors. Am using Ruby 1.8.6.26 ,

Hi, Given test <- matrix(c(1, 1,2,2), 2,2) t <- apply(test, 1, t.test) How can I obtain a table of p-values, confidence interval etc, instead of [[1]] One Sample t-test data: newX[, i] t = 3, df = 1, p-value = 0.2048 alternative hypothesis: true mean is not equal to 0 95 percent confidence interval: -4.853102 7.853102 sample estimates: mean of x 1.5 [[2]]

Dear All, Need some help as I was trying to follow the guide below. https://wiki.samba.org/index.php/Samba4/HOWTO/Join_a_domain_as_a_DC Until the steps of ldbsearch -H /usr/local/samba/private/sam.ldb '(invocationid=*)' --cross-ncs objectguid and my ldbsearch reply with such a result. ldb: unable to dlopen /usr/lib64/samba/ldb/acl.so : /usr/lib64/ldb/libreplace.so: version

1. This is an imaginary data on monthly outcomes of 2 years and I want to forecast the outcome for next 12 months of next year. data Data1; input Yr Jan Feb Mar Apr May June July Aug Sept Oct Nov Dec; datalines; 2008 12 13 12 14 13 12 11 15 10 12 12 12 2009 12 13 12 14 13 12 11 15 10 12 12 12 ; run; I converted the above data into the below format to use it in R as it was giving error: asking

Hello to everyone,... I am a new R ambitious user. I would like to be the first at my department using R, but I have encountered a difficulty during the last days that I cannot overcome reading help() and searching over the net. Problem: I have multiple files with financial data like the following (header included): E.g.: filename: AOL.txt aol.txt 4 3 5 3... filename: IBM.txt ibm.txt 6 2 5 2...

Hello everybody When one is working with complex matrices, "transpose" very nearly always means *Hermitian* transpose, that is, A[i,j] <- Conj(A[j,i]). One often writes A^* for the Hermitian transpose. I have only once seen a "real-life" case where transposition does not occur simultaneously with complex conjugation. And I'm not 100% sure that that wasn't a

Hi everybody, I'm a new R french user. Sorry if my english is not perfect. Hope you'll understand my problem ;) I have to work on temperature data (35000 lines in one file) containing some missing data (N/A). Sometimes I have only 2 or 3 N/A following each other, but I have also sometimes 100 or 200 N/A following each other. Here's an example of my data, when I have only small gaps

I have a dataframe ("mydf") that contains "differences of means". I wish to test whether these differences are significantly different from zero. Below, I calculate the t-statistic for each column. What is a "good" method to calculate/look-up the p-value for each column? mydf=data.frame(a=c(1,-22,3,-4),b=c(5,-6,-7,9)) mymean=mean(mydf) mysd=sd(mydf)

Hi all, I am using JRI with the XTS library. Rengine rEngine = new Rengine(new String [] {/*"--vanilla"*/}, false, null); rEngine.eval("library(xts)"); When I create an xts object, using previously assigned double/date arrays with roughly 5000 elements, the eval method needs very long (up to a minute) to return. I would like to understand what is going on, but I find no