similar to: scatter.smooth error

Dear list, In using the scatter.smooth() function (modreg) on a small data set (100 obs) the following error was produced: > scatter.smooth(Na, S) Warning message: k-d tree limited by memory. ncmax= 200 I haven't used scatter.smooth much but when I have, I haven't seen this message before. gc() returns > gc() used (Mb) gc trigger (Mb) Ncells 417693 11.2 667722

Hi, I have a data.frame in wide format which I would like to reshape into a long format: example (nonsense) data: > dput(perus2) structure(list(id = c(30L, 38L, 21L, 12L, 22L, 28L, 31L, 44L, 8L, 47L, 23L, 20L, 41L, 42L, 29L, 50L, 5L, 33L, 4L, 17L, 11L, 1L, 18L, 6L, 9L, 32L, 16L, 14L, 39L, 48L, 37L, 43L, 25L, 27L, 35L, 40L, 45L, 49L, 2L, 34L, 10L, 15L, 24L, 26L, 46L, 7L, 13L, 3L, 19L, 36L),

I rewrote scatter.smooth to handle missing values, but I have a question about a move I had to make. Here's the code: Mscatter.smooth<-function (x, y, span = 2/3, degree = 1, family = c("symmetric", "gaussian"), xlab = deparse(substitute(x)), ylab = deparse(substitute(y)), ylim = range(y, prediction$y), evaluation = 50, ...) { if (inherits(x,

Dear List, Having had my appetite sufficiently whetted by Prof. Ripley's email about the new graphics capabilities in Unixes, I wanted to try them out. I updated to svn r44608, configured with the following options: R is now configured for x86_64-unknown-linux-gnu Source directory: .. Installation directory: /usr/local C compiler: gcc -O3 -g -std=gnu99

Hello, I have been using scatter.smooth and been trying to format my plots and been unable to adjust the line type or color of the smooth curve. When I specify col in scatter.smooth, I can adjust the color of the scatter plot points, but the line always stays black and solid. lty has no effect. When I specify par prior to calling scatter.smooth, col will change the color of the points,

Hi there, I would like to draw a scatter plot and fit a smooth line by loess. Below is the data. However, the curve line started from 0, which my "resid" list doesn't consist of 0 value. It returned some warnings which I don't know if this is the reason affecting such problem. Here I also attached the warning messages. Please let me know if there is a solution to fix this. Thank

Hi im trying to add error bars to my barplots, there very basic, i have a few grapghs where the y variable is different but on all the X variable is Age (Adult and Juvenile) however this is split into two levels so i have males and females, so my graph basically has four bars on it. I know how to add eror bars for instance when there is only one level eg lookng at the diffrence between male and

I'm updating the loess routines to allow for, among other things, arbitrary local polynomial degree and number of predictors. For now, I've given the updated package its own namespace. The trouble is, panel.loess still calls the original code in package:stats instead of the new loess package, regardless of whether package:loess or package:lattice comes first in the search list. If I

Hi, I have a scatter plot of the variables GNI and Lifeexp (Gross National Income and Life Expectancy, both metric). So I plotted them and I want to add a regression line and a lowess line. I use lowess and not loess because I have missing values. My code: plot(GNI,Lifeexp) abline(lm(Lifeexp~GNI), col="red") y.loess<-loess(Lifeexp~GNI,na. action = na.exclude)

Hello, I'm currently using a tool that provides a Loess fitting, but I obtained results that are slightly different from those provided by R implementation of the Loess. That's why I would like to know if you could give me a source (bibliography or web) that explains in a clear way each step of the algorithm, with the possible options to choose, etc.. in order for me to understand those

Has the loess.smooth() function changed? It used to work, but now it causes R to abort with a segmentation fault. I stole the function points.lines() from V&R 1st ed. pp. 67--68, but now it only works if I remove the line with loess.smooth. Here's the function I'm using: points.lines <- function(x, y, ...) { cor1 <-round(cor(x, y, use="pairwise"), digits=2)

In the HistData package, I have a data frame, PearsonLee, containing observations on heights of parent and child, in weighted form: library(HistData) > str(PearsonLee) 'data.frame': 746 obs. of 6 variables: $ child : num 59.5 59.5 59.5 60.5 60.5 61.5 61.5 61.5 61.5 61.5 ... $ parent : num 62.5 63.5 64.5 62.5 66.5 59.5 60.5 62.5 63.5 64.5 ... $ frequency: num 0.5 0.5

Suppose that I'm working on Hadley's diamond dataset and I want to review the relationship between price, colour and carat. I might run the following: library(ggplot2) #plot scatter and add some hex binning q<-qplot(carat,price,data=diamonds, geom=c("hex"), main="Variability of Diamond Prices by Carat and Colour") #facet to get one scatter for

Hi folks. If I have the following in my "data" event pH1 pH2 1 4.0 6.0 2 4.3 5.9 3 4.1 6.1 4 4.0 5.9 and on and on..... for about 400 events Is there a way I can get R to plot event vs. pH1 and event vs. pH2 and then do a loess or lowess line for each?? Thanks in advance David [[alternative HTML version deleted]]

Hello, I am fairly new to R am stumped on how to get the xyplot function in the lattice package to produce a scatter plot of count data versus time, such that the count data represent 8 different groups, and the plot produced has 8 unique "smoothed" lines for the different groups. This is the closest I can get to the desired plot: xyplot(masmean ~ mas, data = Pre96CR.masmeans, groups

I need to allocate (using C nomenclature) a set of "global" variables, some integer scalars and some double vectors. I have placed the name of such variables in the file containing the main script and the called functions: # ---------------------------------------------------------------------------- # --------------------------- GLOBAL DATA & CONSTANTS ----------------------

I am doing very regular stuff like the following: attach(wtana) fm<- lm(Body.Wt.on.SD1~Heart.Wt, data=wtana) #fm<- lm(wtana$Body.Wt.on.SD1~wtana$Heart.Wt) lrf<- loess(Body.Wt.on.SD1~Heart.Wt, wtana) #lrf<- loess(wtana$Body.Wt.on.SD1~wtana$Heart.Wt) plot(Body.Wt.on.SD1,Heart.Wt) #plot(wtana$Body.Wt.on.SD1,wtana$Heart.Wt) #lines(spline(Heart.Wt,fitted(lrf)), col=2) abline(fm, col=4)

Hello! I have the following problem. I would like to plot the hclust object "hcd" (bellow, at the end of the mail) with ticks at seq(0.05,0.25,by=0.05). I tried using the code plot(hcd) and plot(hcd,axes=FALSE) axis(2,seq(0.05,0.25,by=0.05)) In both cases, the resoult is the same, ticks at 0.05 and 0.25 are not printed. I tried changing the ylim argumet in plot, however I got a

Hi, is it possible to reverse the order in which panel.lmline() or panel.smooth() operation in xyplot() ? This type of situation might occur when plotting some variable with depth, but the relation we want to describe is variable ~ depth, and not depth ~ variable, as the plotting formula would suggest. # an example: d <- 1:100 v <- d * rnorm(100) xyplot(d ~ v, ylim=c(100,0),

Hi, The following call: xyplot(incidence ~ year, melanoma, panel=panel.smooth) produces a blank plot region with an error message: Error using packet 1 plot.new has not been called yet > sessionInfo() R version 2.11.1 (2010-05-31) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 LC_MONETARY=C