similar to: scatter plot

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 really like these Trellis graphics but how do I get this code to skip the missing? logreg<-read.csv("logreg.csv", header=TRUE, sep=",", na.string=" ") attach(logreg) bwplot(yesno~bc_pcb_tot |varlist, data=logreg, main="Box Cox PCB transformation", auto.key=TRUE, fontfamily = "HersheySans" ) Dean Sonneborn M.S. Public Health Sciences *

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,

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)

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'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 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 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 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

I am using lattice to make figures as pdfs: trellis.device(device = "pdf",file = "Figure6.pdf",color = FALSE) I need to specify some blank space on the left-hand margins (the pages will be bound so we need about 0.5 inch)). I have tried a number of solutions but none seems to work (e.g. par.set). Can this be done when initiating the plotting device? Or is the some other way

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)

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

Hello! The question is simple: What is the escape sequence for a new line when using Hershey fonts? I obviously tried '\n' but it didn't work (see the sample below). I looked at 'demo(Hershey)' but all it only shows escape sequences for printable characters. The sample I've been using to try to find the escape sequence is below. You can comment or un-comment the

Hi all, I am trying to make tiff files of lattice plots at a resolution greater than 300 dpi required by a journal (PLoS ONE). I have tried both the tiff and bitmap functions. tiff keeps panel colors but reduces axes and tick labels so they are nearly invisible. bitmap maintains correct label size but only produces greyscale. Regular plots work fine with tiff; the problem is only with lattice

Would anyone with experience with the map functions know how to divide Czechoslovakia into the Czech Republic and Slovakia. They have been two separate countries for some time now. I'm thinking about the worldhires map database in particular. -- Dean Sonneborn, MS Programmer Analyst Department of Public Health Sciences University of California, Davis (530) 754-9516

I'm using this syntax to create data subsets for plots: subset=source=="Both". Now I would like to create a subset defined by two different variables like this: subset=source=="Both" subset=site=="home" but this syntax is not correct. The documents in the manual for subset seem to be creating whole new data files not just selecting rows based on the

Greetings Everybody: I generated a 1.2MB dta file based on the general social survey with Stata8 for linux. The file can be re-opened with Stata, but when I bring it into R, it says all the values are missing for most of the variables. This dataset is called "morgen.dta" and I dropped a copy online in case you are interested http://www.ku.edu/~pauljohn/R/morgen.dta looks like this

Dear List, I have some discrete data and want to calculate the percentiles and the percentile ranks for each of the unique scores. I can calculate the percentiles with quantile(). I know that "ecdf" can be used to calculate the empirical cumulative distribution. However, I don't know how to exact the cumulative probabilities for each unique element. The requirement is similar

I am attempting to read in a SAS 9.1 data file. After starting R I change to the directory containing the sas data file and use the "dir" command to confirm that it is there. Then I run the following R-code: library(foreign) sashome <- "/Program Files/SAS/SAS 9.1" test<-read.ssd(file.path(sashome), "pcb", sascmd = file.path(sashome,