similar to: skip missing values in plots

I'm working on a graphic but have run into a road block about two issues. I don't have a color printer so I want to produce the graphic in black and white. I'm currently using this statement trellis.device(bg="white") but the body of the graphic contains color. What is the code to create the whole thing in black and white. The other issue might be a bit more tricky.

I'd like to do a simple scatter plot but instead of using the variable values on the X axis I would like to plot the percentiles. I searched in the manual for percentiles but did not find what I was looking for. I've been using SAS for several years but I new to R. -- Dean Sonneborn Programmer Analyst Department of Public Health Sciences University of California, Davis (916) 734-6656

I would like to add a footnote to this graph but do not see a "footnote" command in the package:lattice documentation. I would like to note the "span=.8" as the footnote. postscript(file= ?C:/Documents and Settings/dsonneborn/My Documents/Slovak/output/pcb_tables/smooth_PCB_lines_four.ps?, bg=?transparent?, onefile=FALSE, pointsize=20,paper=?letter?, horizontal=TRUE,

On the first issue, there's a recent post you can find in the archives. It's from Deepayan Sarkar on January 2 this year. It would probably pop up on a search for "black white lattice" or something similar. The key part of his answer: I'd do something like this as part of the initialization: <<...>> library(lattice) ltheme <- canonical.theme(color = FALSE)

This is a solution I am posting for a problem that others may have. If you want to: 1. Place lattice graphics from R into an Adobe InDesign document, and 2. Use the export as eps function in R to maximize resolution (it is much better than exporting as a metafile or bitmap), and 3. Use long strings of text in your titles or captions or to label your axes. Then you will have problems because:

I want to produce some boxplots and plot the logged values but have the axis scale in the original, not-logged scale. It seeming like I have the first few steps but I'm having trouble with the last. Here's what I'm doing (which I got for the documentation for boxplot and axis). How do I get the ticks to be labeled 2,5, and 9 of the original scale? a<-c(1,2,3,4,5,6,7,8,9,10)

I'm using this code to plot a smoothed line. These two columns of data really represent 4 groups and I'd like to plot a separate line for each group but have them all in the same plot. The R-Docs for lowess do not seem to indicate some type of "GROUPS=var_name" option. What would be the syntax for this? plot(AWGT ~ lipid ) lines(lowess(lipid , AWGT, f=.8)) -- Dean

I am trying to create a high low close style chart but I keep getting the following error statement; "x is not a open/high/low/close time series". I've used the function is.ts and R responds TRUE but the ohlcPlot function gives the above error statement. I'm actually planning to plot some odds ratios with their confidence intervals and the hi-low-close chart should do the

I am running this code to produce some boxplots. I have every thing that I need except that I would like the whisker line to be solid line, not dashes. I have reviewed the lattice docs but have not seemed to come across this point. print( bwplot( group ~ lpcb_tot, data= data7, xlab="Log PCB", + ylab= + " G B + S M S M

I have been trying to run a logistic regression using a number of studies. Below is the syntax, error message & data. Any advice regarding what I am doing wrong or solutions are appreciated, regards Bob Green > logreg <- read.csv("c:\\logregtest.csv",header=T) > attach(logreg) > names(logreg) [1] "medyear" "where" "who"

Hi there, I've looked through the very helpful advice about adding fitted lines to plots in the r-help archive, and can't find a post where someone has offered a solution for my specific problem. I need to plot logistic regression fits from three differently-sized data subsets on a plot of the entire dataset. A description and code are below: I have an unbalanced dataset

I currently have 3 separate boxplots but would like to put them all in the graphic so they would have the same scale. Below are the three statements and as you can see the Y axis is weight: bwplot(AWGT~ male2 .... bwplot(AWGT ~ bin_pcb2 ..... bwplot(AWGT ~ bin_pcb2 | male2 ..... Does anyone have some sample code where they have done something like this? Thanks, Dean Sonneborn M.S. Public

I am using the R code listed below to create 4 smoothed lines on a xyplot. I'm having trouble fine tuning it. First I think I may need a black and white plot so how do I get it to plot the lines with different characters, preferable the same characters used in the key (plus, X circle and triangle). I might also be interest in a version that draws four solid lines of different colors but

Hi, I was trying to plot the logistic regression from a regression "logreg" I just ran. I downloaded the "car" package from the R website and went to Packages -> install package from local zip file I checked in my library file and the package is there. I restarted R. I then ran the command: reg.line(logreg,col=palette()[2], lwd=2, lty=1) And I get the error: Error: could

Hi there: recently i try to use LogicReg package for a tree model(logistics fit ) . i list my code and error below: > dim(model.dat) [1] 48000 745 > fit1 <- logreg(resp = model.dat[,745], bin=model.dat[, 9:700], type = 3, select = 3, ntrees = c(1,2), nleaves=c(1,7), ) Insufficient declaration LGCn1MAX in logreg() is 20000 LGCn1MAX should be at least 48000 Please fix and

Dear R-users: How does one calculate in R the odds ratios for a CATEGORICAL predictor variable that has 4 levels. I see r-help inquiries regarding odds ratios for what looked like a continuous predictor variable. I was wondering how to get the pairwise odds ratios for comparisons of levels of a categorical predictor variable. I can't seem to get the correct output using: >

I want to use the grplasso package on a data set where I want to fit a linear model.? My interest is in identifying significant?beta coefficients.? The documentation is a bit cryptic so I'd appreciate some help. ? I know this is a strategy for large numbers of variables but consider a simple case for pedagogical puposes.? Say I have?two 3 category predictors (2 dummies each), a binary

Hello. I have run the command imp<-mice(mydata, im=c("","pmm","logreg","logreg"),m=5)  for a variable with no missing data, a numeric one and two variables with binary data. I got the following message: There were 37 warnings (use warnings() to see them) > warnings() Warning messages: 1: In any(predictorMatrix[j, ]) ... : coercing argument of

I want to get some "simple" logic regression examples to work before exploring a hard problem. I can get results, but I'm having some problems using "cat" to save the logic regression equation to a disk file. Consider this: # Simple Logic Regression Example # efg, 23 Feb 2005 library(LogicReg) # Create simulated data with known logic equation: # "noise"

I'm running some GEE models but when I request the summary(pcb.gee) all I get are rows and rows of intercorelations and they fill up the screen buffer so I can not even scroll back to see what else might be in the summary. How do I get the summary function to NOT print the intercorrelations? Thanks, -- Dean Sonneborn Programmer Analyst Department of Public Health Sciences University of