similar to: Using ksmooth to produce the resulting plots

Hi, Assume that we may model the Nottingham temperature data (nottem) or Sunspot data (sunspot) set by a nonparametric autoregressive model of the form Yt = m(Yt-1) + et. Using the kernel estimation method, produce the resulting plots. We may use the fucntion tsmooth(x,y,"notmal",bandwidth=0.01). How can i define x and y using data nottem and sunspot?

Dear R Users, Before running ksmooth( ), a suitable bandwidth selection is needed. I use some functions for this task and receive these results for my data: width.SJ(y,nb=100,method="ste") : 40.25 bcv(y,nb=100) : 40.53 ucv(y) : 41.26 bandwidth.nrd(y) : 45.43 After implementing the function ksmooth(x,y, bandwidth= each of abovementioned bandwidths), I have some NAs

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

Apologies in advance if my comments don't help, in which case, no need to respond, but I noted in ?ksmooth: "bandwidth the bandwidth. The kernels are scaled so that their quartiles (viewed as probability densities) are at ? 0.25*bandwidth." So, could this be a source of the discrepancies you cited? Given that ?ksmooth explicitly says: "Note: This function was implemented for

Dear Sir, Madam, or to whom this may concern, my name is Jan Failenschmid and I am a Ph.D. student at Tilburg University. For my project I have been looking into different types of kernel regression estimators and corresponding R functions. While comparing different functions I noticed that stats::ksmooth returned different estimates for the same bandwidth as other kernel regression estimators

I am working with a model that I have to estimate a nonparametric function. The model is partial linear i.e. Y=X$\beta$ + f(z) + $\epsilon$ I am using the ' double residual methods' Robinson (1988) Speckman (1988) where I estimate a nonparametric function for each of the parametric variables in terms of the nonparametric one i.e. X[,i]=g(Z)+ u this is done because I need the $E(

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'd like to superimpose a ksmooth() onto a barplot(). My data is: > d 2009-06-20 2009-06-21 2009-06-22 2009-06-23 2009-06-24 2009-06-25 2009-06-26 2009-06-27 2009-06-28 2009-06-29 2009-06-30 2009-07-01 2009-07-02 Same Breed (B) 12.64 21.08 13.52 12.51 13.71 9.91 14.24 7.18 11.81 5.92 12.04 17.96

Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari --------------------------------- [[alternative HTML version deleted]]

Dear R user, My input data, positive and negative, is complete without missing data. After running ksmooth( ) , I receive for $y , many NAs. What could be the reasons and how can I receive complete output? So many thanks for any help. Amir Safari --------------------------------- [[alternative HTML version deleted]]

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

Le 31/05/2017 ? 17:30, Serguei Sokol a ?crit : > > More thorough reading revealed that I have overlooked this phrase in the > line's doc: "left and right /thirds/ of the data" (emphasis is mine). Oops. I have read the first ref returned by google and it happened to be tibco's doc, not the R's one. The layout is very similar hence my mistake. The latter does not

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