similar to: locpoly

is there a way of estimating derivative curves, similar to the ones we get from 'locpoly', from 'loess' estimation. i am interested in estimation of 1st and 2nd derivatives... --------------------------------- [[alternate HTML version deleted]]

Hello All, Is there a local bandwidth selection routine for local polynomial regression (locpoly) ? Thanks Chinthaka Kuruwita

Dear list, I've got a question concerning difference between loess and locpoly. I have to use a plug-in method to chose a bandwith so I take locpoly method to fit a curve. My problem is:I know how to get predicted values in loess: m=loess(y~x) y_fitted=predict(m). But how to get the same in locpoly? I computed like this: bw=dpill(x, y, blockmax = 5, divisor = 20,trim = 0.01, proptrun =

Hi: I am using the package "KernSmooth" to do the local polynomial regression. However, it seems the function "locpoly" can only deal with univariate covaraite. I wonder is there any kernel smoothing package in R can deal with bivariate covariates? I also checked the package "lcofit" in which function "lcofit" can indeed deal with bivariate case. The

Hi people, I am new in this list and could not find a FAQ for it in particular, furthermore I could not find my question answered in the official R FAQ or docs. I have simply something like this: > f<-approxfun(data[,1],data[,2]) and f is: > f function (v) .C("R_approx", as.double(x), as.double(y), as.integer(n), xout = as.double(v), as.integer(length(v)),

Dear R users, I am currently using the locpoly function from the KernSmooth package to estimate densities. However, I have some trouble understanding how this estimation technique is implemented in R. My main concern comes from the fact that this function gives negative estimates when the bandwidth is sufficiently large (mainly in the tails of distributions). I have read some articles on this

Hi all, I have a data frame called "rst", see below: ------------------------------------------------------------------------------------------ # This is a paste able example # In case you don't have "KernSmooth" package installed, please uncomment below line. # install.packages("KernSmooth") library(KernSmooth) rst <- data.frame(hsp = rnorm(23), dal =

1. Please do make use of the subject line. 2. Please (re-)read the description of the `range.x' argument in ?locpoly: it's suppose to be a vector of min and max x values. 3. I hope you realize how grossly inefficient this computation is... Andy > From: klea lambrou > > > hello R-users.could you please help me on this one?i have > 2 vectors x1 > and y1 and

Hello!! At first thank you all for your suggestions and your help. But I would like to get another group of datas from the density estimation. I don't know if these information are provided by locfit or I should use an other function like plot(). Description: Each number between 1 and 100 (only integer) for example should have only one density value. But I also would like to have the density

Dear R Users, Suppose we are interested for generating a new vector ( x ) from a current vector (y) of length 1000 so that x includes the sum of every 5 values in y respectively from the first to the end of length y. The same length of y for x is desired, so that other 4 positions (indices) in x are filled out with NA. For generating such a new vector, I have no idea. I tried in some

Dear R Users, For running wavelet functions using dwt( ), modwt( ), and mra( ), a wavelet filter algorithm is applied. For all these functions, default is "la8" and other possibility is "haar". In related documents, another possibilities like as symlet and coiflet ... are not cited. Besides "la8" and "haar", which wavelet filters can be used? Thank

Hi R Users, I don't know how much is difficult my problem and even it is possible to solve in R or not. Given a vector with 2000 observations. I want to creat a new vector from that vector so that new vector be the sum of every 5 observations sequently. That is , each new observation is sum of every 5 sequent observations, from the first observation of previous vector to the end. So

Hi Does R have a function or has someone written a function to draw a 3d surface from a scatter plot of values using either ksmooth or locpoly. OR a transform function a that merges x relation z and y relation z to (x,y) relation z? I tried out scatterplot3d but it seems it would take a bit of work to get scatterplot3d to draw a curved surface. Lawrence

I would like to get the fortran code for the functions locpol and dpill in the kernsmooth library , the purpose is to incorporate other selection criterion for the bandwidth in which a varying bandwidths are replaced instead of the constant bandwidth. I will send you back the improved version as soon as I get the fortran code for the upper mention functions. ((LOCPOLY and DPILL )))

Updating my 1.8.0 R installation (>update.packages() ) I obtain the following (SORRY FOR THE LENGTH OF THE LOG BUT IT HELPS!!!): ................ downloaded 135Kb KernSmooth : Version 2.22-11 in /usr/lib/R/library Version 2.22-12 on CRAN Update (y/N)? y mgcv : Version 0.9-3.1 in /usr/lib/R/library Version 0.9-6 on CRAN Update (y/N)? y trying URL

Dear R Friends , I want to name my variables( more than 2 variables in a single plot) within a plot to distinct them from each other, but I cann't. How it is possible? I don't mean x and y axis using xlab or ylab. At the below , it follows some lines, only as an example that you could try please, if it is possible. I really thanks for your attention. Amir library(graphics) y<-

Hi dear R users, Is it basically correct that a problem is ( linearly on nonlinearly ) modeled so that the slack variable is bounded by an upper bound ? If so, how it can be handled and coded practically ? for example: x1+ x2 =< b so ----> x1 + x2 + s=b s=b- x1 - x2 b- x1 - x2 =< upper value But algorithms can not calculate b- x1 - x2 , because

Dear R Users, Suppose we want to creat a new vector ( x ) from a current vector (y) of length 1000. The current vector y includes negative, zero and positive values. We want our new vector x includes the negative values in y, otherwise NA with the same length as y. For this, we have x=y[y<0] . Now x includes a subset of y with shorter length than y. With x=match(y,x) we would

Dear R users, Sometimes it is needed to read only some columns from a table, in particulare for high frequency data. How it is possible to read just some certain columns using read.table ( ). The reason could be keeping space in R and in particular accelerating in reading data when the number of rows are huge and some of them are not needed. Thank you very much, Amir

> Background: > OS: Linux Mandrake 9.1 > release: R 1.9.0 > editor: Xemacs 21.4 > frontend: ESS 5.1.23 > --------------------------------- > > Colleagues > Is there a function in R that is similar to polyval in matlab? (y = polyval(p,x) returns the value of a polynomial of degree n evaluated at x. The input argument p is a vector of length n+1 whose elements are the