search for: ethanol

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)", panel = function(x, y) { panel.grid(h=-1, v= 2) panel.xy...

I would like to draw smooth curves instead of polygons. I could not find any spline function to do that : given an x and a y vectors, they all take the x in increasing order. Is there a function to draw a smooth curve through a set of points in any order ? with many thanks in advance Alain Guerreau CNRS Paris [[alternate HTML version deleted]]

...cument I found on http://www.locfit.info However, first, it is a local regression method. I am not sure how to specify the degree of the regression model. second, i dont know what are "scale" and "alpha". Are they associated with the bandwidth? Thanks very much! data(ethanol) # a bivariate local regression with smaller smoothing parameter fit <- locfit(NOx~E+C, data=ethanol, scale=0, alpha=0.5) plot(fit) --------------------------------- [[alternative HTML version deleted]]

I am trying to read and write database tables that have a date field in them. I am constructing a data.frame, then using dbWriteTable to create the table and dbReadTable to read it. >datatbl<-data.frame( dates, trialnum, fooddel, ethdel, trialtime, trialtype, deliveries, food, ethanol, fcumrec, dcumrec, rrf, rrd) >ifelse(startdate==filelist, filetbl<-datatbl,filetbl<-rbind(filetbl, datatbl)) >dbWriteTable(con, subj, as.data.frame(filetbl), overwrite=T) This creates the table named for the contents of subj...the problem is in the class of the first column, date...

...on. I have made two attempts. The first, below, almost works, but only takes the first character from the V3 column: library(lattice) xyplot(V1~V2, data=a, pch=as.character(a$V3)) I've also tried this example, which is given in an online user guide for trellis. It uses the built-in "ethanol" data set (which is also in lattice), and subscripts... but when I try the same code in lattice, I get an "invalid graphics state" error. library(lattice) xyplot(NOx ~ E | C, data = ethanol, aspect = 1/2, panel = function(x, y, subscripts) text(x, y, subscripts,...

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) panel.xyplot(x,y) pan...

...t seems to be a bug in trellis.focus. I am using JGR R 2.8 with lattice 0.17-15 Note: It's important to you use JGR to use the Java Graphics Device. *** To Replicate the Bug ***: 1.- display an xyplot. For example, from the xyplot help page: library(lattice) require(stats) EE <- equal.count(ethanol$E, number=9, overlap=1/4) ## Constructing panel functions on the fly; prepanel 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) {...

...at 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 with an existing .Rdata and in another directory without one. The call below, isolated from its script, is one of the culprits: with(ethanol, loess(NOx ~ C * E, span = 1/3, degree = 2, parametric = "C", drop.square = "C", family="s")) In response to the Posting Guide, when I say 'crash' I mean that the system abnormally terminated R when this line or any other one that invokes loess() is ru...

...it(0.8,"npc"), height=unit(0.2,"npc"), just=c("left","bottom"), x=unit(0.8,"npc"), y=unit(0.8,"npc"), angle=-90, name="v3" ) pushViewport(v1) grid.rect() print( xyplot( NOx ~ E, data=ethanol, ), newpage=FALSE ) upViewport(); pushViewport(v2) grid.rect() print(histogram( ethanol$E ),newpage=FALSE) upViewport() pushViewport(v3) grid.rect() print( histogram( ethanol$NOx ) , newpage=FALSE ) [[alternative HTML version deleted]]

...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: library(lattice) require(stats) EE <- equal.count(ethanol$E, number=9, overlap=1/4) ## Constructing panel functions on the fly; prepanel 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) {...

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)

...: C/Fortran function name not in load table ********** Do you know what this message means and do u know what procedure shall I implement in order to be able to use best fit. I look forward to receiving your advices. Thank you very much for your time in advance locfit. N.B : i use the example of Ethanol in pckg (locfit). Nour-Eddine. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subjec...

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

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

...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 "good" graph in 3d? It probably works in S-Plus, but I am using R. (I have tried plot (locfit object, type="persp") on the ethanol dataset - but this gives a poor result, with a few warnings.) THANK YOU!. Jacob Jacob L van Wyk Department of Statistics University of Johannesburg APK P O Box 524 Auckland Park 2006 South Africa Tel: +27-11-489-3080 Fax: +27-11-489-2832

Hi all: There was a concern raised by reviewers of a manuscript of mine over the proper execution of a Pearson's correlation. In brief, this was undertaken in order to determine the relationship between the extent of wheel running (y axis) and ethanol intake (x axis) across three, separate 10 day periods in 7 animals. In the paper, the correlational plots for each 10 day-period had 70 data points: One point for each day and each animal across 10 days of experimentation. The reviewers, however, appropriately pointed out that this is a violation...

Is there any place where the Chambers and Hastie 1993 Statistical Models in S data sets can be found? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San

I've looked every way I can think of for help on trellis.datasets, but nothing comes up for me. Please help me find information on what is included, and how to get at those data. Thanks. Dave Parkhurst

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

...using the leaps function, and I run across the following error: > leaps(data[,3:7], data[,1], nbest = 10) Error in leaps.setup(x, y, wt = wt, nbest = nbest, nvmax = NCOL(x) + int, : character variables must be duplicated in .C/.Fortran Here is a sample of what 'data' looks like: ethanol flask batch delabso delgluc delglyc ph 1 0.00 1 0 1.41 0.0 0.7 1 2 0.00 2 0 1.33 0.0 0.6 9 3 0.00 2 0 1.18 0.0 1.1 1 4 0.00 3 1 1.58 0.0 3.5 1 5 0.00 4 0 1.25 0.0 5.0 1 6 0.00...