similar to: smooth curves

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 all, I'm trying to find a way to get xyplot to produce a second set of axes outside the right hand side of the graph. This is my progress so far: EE <- equal.count(ethanol$E, number=9, overlap=1/4) xyplot(NOx ~ C | EE, data = ethanol, prepanel = function(x, y) prepanel.loess(x, y, span = 1), xlab = "Compression Ratio", ylab = "NOx (micrograms/J)",

Hello there, I am presently using the locfit package in "R" and would appreciate some help here. Could anyone tell me tell me how I can obtain the various components (x and y values)of the density estimation after using the "R" command "locfit"? For example, with the command "fhat<-density(somename)" I can obtain the x and y values by simply calling

Greetings and happy new year! I am in the process of converting some of the old S-PLUS scripts from Visualizing Data (Cleveland, 1993) into lattice. In fact, I did most of it several years ago, and at the time, all of the scripts that contained loess() worked fine. Tonight, I ran most of the scripts again, but every one that I tried with a loess() call crashed R. I tried it in two sessions, one

Dear all, In many cases, I need a plotting region much bigger than the screen (e.g. for maps or for graphs with many labels). A. MS-Windows if I try windows(width=25, height=25, rescale="fixed") it seems to be OK (a screen with scrollbars, exactly what I need) but if I try then plot(faithful$eruptions, faithful$waiting) I receive Error in plot.new() : Outer margins too large

I ran the following script from xyplot Examples using Tin-R on Windows and saw no plot produced. EE <- equal.count(ethanol$E, number=9, overlap=1/4) xyplot(NOx ~ C | EE, data=ethanol, prepanel = function(x,y) prepanel.loess(x, y, span=1), xlab="Compression Ratio", ylab="NOx (micrograms/J)", panel = function(x,y) { panel.grid()(h = -1, v=2)

Hello, I would be very grateful if anyone could help with what seems like a simple lattice task. I want to use xyplot, where the symbols for the plotted points are taken from another column in the data frame. So if the data frame looked like: a <- as.data.frame(matrix(data=c(1,1,10,2,2,20,3,3,30), nrow=3, ncol=3, byrow=TRUE)) a V1 V2 V3 1 1 1 10 2 2 2 20 3 3 3 30 you would get

Dear, I am doing some statistics analysis using R software. I am currently trying to implement your library locfit. I am using Win NT and R1.1.1 version obtained on the following Webmaster : http://cm.bell-labs.com/cm/ms/departments/sia/project/locfit/index.html When trying to run it, the following message appears: ********** Error in .C("guessnv", nvm = integer(1), ncm =

Hi, I wonder if anyone knows how to control color and darkness of the background grid generated by ' type = c("g", "p") ' in a lattice plot (e.g., in xyplot). The documentation does not seem to offer a way to change them. Thanks a lot! Wen

I have a question regarding the package LOCFIT of C Loader. After fitting a smooth surface (y modeled by two regressors), I can't seem to find a nice way of making a perspective (or other 3d plot). The contours don't look quite so good. It is meant to be used as illustration of smoothing to my students in regression analysis. Could anybody perhaps assist me in producing a

Hi: (Tried to find a bug report about this issue, but was unable to find it, let me know if this is a known issue) I have been working on an interface to highlight xyplot panels on mouse overs in JavaGD but I have stumbled with what seems to be a bug in trellis.focus. I am using R 2.8 with lattice 0.17-15 *** To replicate the bug: 1.- display an xyplot. For example, from the xyplot help page:

Dear R People: Could someone recommend a small example for spline regression, please? I would like to show this to my students, step by step. Thanks in advance, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess at gmail.com

Hi, I'd like to plot 4 groups of data using xyplot and panel.superpose so that the points are overlayed on a single plot. For each group of data I'd also like a loess smoothed function (using panel.loess). I have tried the following: xyplot(series ~ time | gr, data=etable, panel = function(x,y, ...) { panel.superpose(x,y, ...) panel.loess(x,y,span=.15)

Dear R listers, When using smooth.spline to interpolate data, results are generally good. However, some cases produce totally unreasonable results. The data are values of pressure, temperature, and salinity from a probe that is lowered into the ocean, and the objective is to interpolate temperature and salinity to specified pressures. While smooth.spline provides excellent values at the

Dear All, I've been?searching for?appropriate codes to compute the rate of change and the curvature?of ?nonparametric regression model whish was denoted by a smooth function?but?unfortunately?don't manage to?do?it. I presume that such characteristics from a smooth curve can be determined by the first and second derivative operators. The following are the example of fitting a

Dear list members, I am trying to understand this output from the smoothing package locfit (1.5-4, running on R 2.8.1 on Windows Vista 64 bit). # sample code x<-1:100 y<-rnorm(100) fit<-locfit(y~x,family="gaussian") #default parameters are fine plot(fit,band="global") #plot seems "reasonable", confidence bands use a global estimate of variance

Dear R users I want to compute Kendall's Tau between two vectors x and y. But x and y may have zeros in the same position(s) and I wrote the following function to be sure to drop out those "double zeros" "cor.kendall" <- function(x,y) { nox <- c() noy <- c() # for (i in 1:length(x)) if (x[i]!= 0 | y[i] != 0) nox[length(nox)+1]<- x[i] for (i in

Hi, I cannot find a 'vectorized' solution to this 'for loop' kind of problem. Do you see a vectorized, fast-running solution? Objective: Take the value of X at each timepoint and calculate the corresponding value of Y. Leading 0's and all 1's for X are assigned to Y; otherwise Y is incremented by the number of 0's adjacent to the last 1. The frequency and

Dear R-devel, I'm working on Prof. Loader's new version of locfit to try to get it pass R CMD check. I'm almost there, but I have a problem with some Rd files that I hope some one can help me resolve. Here's an example: In the package there's a function called locfit.censor(). This function can be used in a few different ways: locfit.censor(x, y, cens, ...)

DeaR list, I'm trying to represent some information via 3D plots. My data and session info are at the end of this message. So far, I have tried scatterplot3d (scatterplot3d), persp3d (rgl), persp (graphics) and scatter3d (Rmcdr) but any of them gave me what I'd like to have as final result (please see [1] for a similar 3D plot changing PF by ypred, pdn by h4 and pup by h11). In general