similar to: Function tsmooth

Hi, Consider the Nottingham temperature and the Sunspot data: >data(nottem) >data(sunspot) Assume that we may model each data set by a nonparametic autoregressice model of the form Yt = m(Yt-1) +et, What are the value of x and y for using function ksmooth(x,y,"normal",bandwidth=0.01)? Thanks a lot! Owen

There is now a preliminary version of a time series package in the R-devel snapshots, and we would welcome feedback on it. It is based in part on the packages bats (Martyn Plummer) and tseries (Adrian Trapletti) and in part on code I had or have written. (Thanks for the contributions, Martyn and Adrian!) Some of the existing ts code has been changed, for example to plot multiple time series, so

can you recommend a good manual for R that starts with a data set and gives demonstrations on what can be done using R? I downloadedR Langauage definition and An introduction to R but haven't found them overly useful. I'd really like to be able to follow some tutorials using a dataset or many datasets. The datasets I have available on R are Data sets in package 'datasets':

I noticed this in the NEWS: o Datasets 'sunspot.month' and 'sunspot.year' are available separately but not via data(sunspot) (which is used by package lattice to retrieve a dataset 'sunspot'). Now, the 'sunspot' dataset included in lattice was actually a mistake, I realized later that I instead wanted the 'sunspots' data from

Dear helpers, I am trying to write a function to return prediction values using package ts. I have written three different versions since I am not sure what's wrong with my func2. func and func1 return the same results.But func1 and func2 don't. In particular, the only difference between "func1" and "func2" is the function variable name being y and data, respectively.

Hi all, Does anyone know how ts objects ts(base) can be used within the 'seas' package? I would like to obtain seasonal statistics of regular time-series and for example look at the result of the plot.seas.var() function or use the change function() to look at change between periods or time-series. The nottem time-series are similar to the time-series we are analyzing (but with

Hi all, I need to (quickly!grrr) update IP addresses for which a given machine will answer. I did this: /sbin/ifconfig eth0:0 xx.xx.xx.10 netmask 255.255.255.0 /sbin/ifconfig eth0:1 xx.xx.xx.11 netmask 255.255.255.0 /sbin/ifconfig eth0:2 xx.xx.xx.30 netmask 255.255.255.0 /sbin/ifconfig eth0:3 xx.xx.xx.31 netmask 255.255.255.0 /sbin/ifconfig eth0:4 xx.xx.xx.13 netmask 255.255.255.0

Hi all, On my newly up-and-running nameserver (CentOS 5), I noticed the following alerts in /var/log/messages after restarting BIND. (lines inserted to aid in reading). As I'm new to SELinux, I'm hoping for some pointers on 1) if this is an issue which simply *must* be addressed, or if it's something I should live with, and 2) how to eliminate the warming messages without sacrificing

Hi am using rails 2.x and outoftime-sunspot gems I want to use morethan one solr cores and using sunspot in my application How to configure morethan one cores in sunspot.yml let''s assume Model1 using core1 Model2 using core2 How to configure them in my application -- Posted via http://www.ruby-forum.com/. -- You received this message because you are subscribed to the Google Groups

Problem: How can I format tickmark labels of log scaled axes of lattice graphics in the usual `xxx'-Format (and not in the scientific format). Example: (according to the help-page of xyplot): In the first plot I get the xxx-Format, in the second plot I get the scientific format (10^xxx): data(sunspot) plot( 1:37, sunspot, log='y',type='l') xyplot( sunspot ~ 1:37,

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 folks, I've put my new DNS server in place, told the primary that this IP would be pulling zones and restarted BIND on the primary. I configured my named.conf file to pull zones from the primary, started BIND on the new box. Here's where I get confused ++++++++++++++++++++++++ [root at sunspot etc]# /etc/init.d/named status number of zones: 169 debug level: 0 xfers running: 0 xfers

Here are another three problems with logarithmic scales: 1) segments() does not work with logarithmic scales. I suggest to change lines 962-973 in "plot.c": for (i = 0; i < n; i++) { if (FINITE(xt(x0[i%nx0])) && FINITE(yt(y0[i%ny0])) && FINITE(xt(x1[i%nx1])) && FINITE(yt(y1[i%ny1]))) { GP->col = INTEGER(col)[i % ncol];

"contour" can't return a list of contours (it would be nice) "arrows" has different options from S-PLUS (no "open" or "rel" options, no "lwd" parameter or graphics parameter pass-through; "size" in S corresponds to "length" in R). [I know we needn't follow S-PLUS syntax slavishly, but it would be nice at least to

Fritz just put the first version of a new time series package to the contrib section at CRAN. The package is called "tseries.tgz" and provides a library for time series analysis. It contains acf Autocorrelation Function adf.test Augmented Dickey-Fuller Test amif Auto Mutual Information Function bds.test BDS Test

Fritz just put the first version of a new time series package to the contrib section at CRAN. The package is called "tseries.tgz" and provides a library for time series analysis. It contains acf Autocorrelation Function adf.test Augmented Dickey-Fuller Test amif Auto Mutual Information Function bds.test BDS Test

hi all, i would like insert a fucntion call at the end of each basic bloc the fucntion i have defined and declared at the begining of the module ie i have written a pass derived from a module pass virtual bool runOnModule(Module & M){ // i declared the fucntion and defined it for (Module::iterator a = M.begin(), b = M.end(); a != b; for (Function::iterator i = a->begin(), e =

Hi all, Can I run a cointegrating regression, for example delta Xt=a1(Yt-1-cXt-1)+E1t and delta Yt=-b1(Yt-1-cXt-1)+E2t with R were Xt and Yt are non stationary time series at t a,b,c are parameters and E1t and E2t are error terms at t. Yt-Xt is stationary Any suggestions are welcome. Best regards, /fb -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Hi all, I want to estimate parameters from a VARMA(p,q)-Modell. The equations of the model or the model structures is given by: Xt=beta1+beta2*Xt-1+beta3*Yt-1+epsilon1 Yt=beta4+beta5*Yt-1+espilon2 epsilon1 and espilon2 are white noise. Xt is given by a vector of n elements e.g. (2, 4, 7, 9, …,n)’ and Yt is given by a vector of n elements e.g. (4,9,12,17,…,n)’. The lineVar from

Hi, This is actually a theory question. I'm a bit confused by the s.window parameter in the stl() function (which is in the ts package). For example, in the stl documentation it uses the nottem data, and then: plot(stl(nottem, s.win = 4, t.win = 50, t.jump = 1)) What does it mean by s.win = 4? Is it because a year has 4 seasons (namely Spring, Summer, Autumn and Winter)? If so will it