similar to: Newbie Time Series Questions

I have a dataframe that I'm trying to reshape, and need advice. My data: > klam.merge[200:225,] stream lulc x sumlength pct.lgth 200 1223030419685 92 0.25000000 9.89 2.52780586 201 1223030419686 23 0.00274154 4.73 0.05796068 202 1223030419686 41 0.75009917 4.73 15.85833341 203 1223030419686 42 2.65000000 4.73 56.02536998 204

Hello all! I have a problem to group my data (years) in 10 years classes. For example for year year decade 1598 1590-1600 1599 1590-1600 1600 1590-1600 1601 1600-1610 --- my is like this> [1] 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 [16] 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 [31] 1628 1629 1630 1631 1632 1633

Hi, I have a monthly time series with missing values and I would use stl function to identify seasonality. I tried all settings of na.action but the result is the same: stl(tm245,s.window=11, na.action=na.pass) Error in stl(tm245, s.window = 11, na.action = na.pass) : NA/NaN/Inf in foreign function call (arg 1) Can you help me? Thanks Andrea Toreti [[alternative HTML version

Hi, I'm dealing with time series. I usually use stl() to estimate trend, stagionality and residuals. I test for normality of residuals using shapiro.test(), but I can't test for autocorrelation and heteroskedasticity. Is there a way to perform Durbin-Watson test and Breusch-Pagan test (or other simalar tests) for time series? I find dwtest() and bptest() in the package lmtest, but it

Hello, I'm a neurophysiologist, and actually quite new to R. In fact, to analyse the responses of my neurons, I would like to have R doing an non-orthogonal ANOVA with a type III sum of squares. Could anybody help me? Thank you a lot in advance. Francois -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'm using R 1.6.2 on a Windows 2000 machine. I've plotted the results of an MDS run labeled by a numerical ID, and color coded by a group code: plot(cv.mds.spr$points, type="n", main="Non-Metric Multidimensional Scaling of SprRun CV Watersheds") text(cv.mds.spr$points, labels = as.character(cv.wshed.id.spr), col = codes(cv.wshed.grp), cex=.75) Question is, how do I

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

I am running R 1.2.2 under Linux. When using the function stl in the ts library, how can I save the seasonal component? What I would like was something like: library(ts) data(nottem) data.stl <- stl(nottem, "per") x <- data.stl$sea This what I get: > x NULL I would, however, like to store in x the seasonal component. Thanks in advance. Francisco. -- Francisco

Hi, again i am stuck in my presentation, and i have never learn R before in my life but need this to be done, so please help me out for a favour: http://www.nabble.com/file/p20333155/kew.dat kew.dat run this in R and these comes up: Month Year Rain 1 Jan 1900 74.400000 2 Feb 1900 80.500000 3 Mar 1900 23.600000 4 Apr 1900 23.600000 5 May 1900 25.100000 6

Easy question that I can't find an answer for. I'm trying to subset a data frame and want to exclude the positive values, i.e. I want the NA values. My data: > summary(temp$tuna) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 1 2 3 3 4 5 1211 Querying for subset(temp, tuna %in% "NA", select.... subset(temp, tuna == NA,

Hello all! I have a problem to extract values greater that for example 1820. I try this code: x[x[,1]>1820,]->x1 Please help me! Thank you! The data structure is: structure(c(2.576, 1.728, 3.434, 2.187, 1.928, 1.886, 1.2425, 1.23, 1.075, 1.1785, 1.186, 1.165, 1.732, 1.517, 1.4095, 1.074, 1.618, 1.677, 1.845, 1.594, 1.6655, 1.1605, 1.425, 1.099, 1.007, 1.1795, 1.3855, 1.4065, 1.138, 1.514,

This isn't an R issue per se, but given the recent issues with Mozilla and help.search(), I figured I'd float it. I use mozilla on a WIN2000 box. I was using 1.0 and the help.search() wasn't working as has been discussed over the past few days. Yesterday I upgraded to 1.1, and the javascript help still does not work. (I checked and it is enabled.) Another problem I noticed this morning

Hi all, How can the label of the x-axes in the plot() of a stl.object be adapted? e.g., When plotting: plot(stl(nottem, "per")) In the labels of the x-axes is “time”. How can this be changed to e.g., “Time (dekade) “? It does not work with xlab or others anymore… Thanks, Jan _______________________________________________________________________ Ir. Jan Verbesselt Research

Hi, I have the panel data of income statement of several banks. The date 9803 means 1998-Q1, 9806 means 1998-Q2, etc. I transform the date code to 1 (for 1998Q1), 2 (for 1998Q2), ...., 16 (for 2011Q4) where 1, 2, .... are placed in Col1. Now the income statement of a specific quarter is actually the accumulation from the beginning of the year. For example, the cost data of 1999Q3 is the

I'm using 1.6.2 on a win 2k box. (I know I'm due for an upgrade.) I've used brute force to label a series of vectors, and I'm wondering if there's a better way to do this. I wanted to create a biplot of the 2nd & 3rd PCs, and did this: test <- edit(loadings(cv.prc.spr)) test.1 <- test[,1] test.2 <- test[,2] test.3 <- test[,3] x0 <- 0 y0 <- 0 i <-

Hello all! I have a problem with R. I try to merge data like this: structure(c(2.1785, 1.868, 2.1855, 2.5175, 2.025, 2.435, 1.809, 1.628, 1.327, 1.3485, 1.4335, 2.052, 2.2465, 2.151, 1.7945, 1.79, 1.6055, 1.616, 1.633, 1.665, 2.002, 2.152, 1.736, 1.7985, 1.9155, 1.7135, 1.548, 1.568, 1.713, 2.079, 1.875, 2.12, 2.072, 1.906, 1.4645, 1.3025, 1.407, 1.5445, 1.437, 1.463, 1.5235, 1.609, 1.738, 1.478,

Dear R-users, I have a complex line by xy-values (ordered by z). And I would like to get interpolated y-values on the positions of x = 0:600. How do I get the correct points? x=c(790,790,790,790,790,786,783,778,778,766,763,761,761,761,715,628,521,350,160,134,134,129,108,101,93,111,161,249,288,243,139,45,7)

Does the polymars library exist anymore? I see reference to it at http://rweb.stat.umn.edu/R/doc/html/packages.html, but not on CRAN (looked in contrib and contrib/Old). -- Rob Schick Ecologist NOAA Fisheries Santa Cruz Lab 110 Shaffer Road Santa Cruz, CA 95060 Phone: 831.420.3960 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'm using R 1.6.2 on Windows 2000. I have two similar sets of files in two different directories. One contains species presence data and environmental measurements; the other contains species absence data and environmental measurements for the same variables as the presence data. The absence datasets contain many more observations than the presence datasets, and I'd like to draw a

Hello, >From the date time stamps I see that not only do I have the results of the build in 'gnuwin32/front-ends' but I can also see the time stamps have changed in the main bin directory so I assume that the result of 'make all recommended' also copies the results there. Now on to the packages. I was interested in following the flow of 'stl' so I added the following