search for: pzs6

...d guy and new to this list... But have been using R for years now. I want to make a map of counties in the US with shaded colors that depend on the level of variable "Y" that I want to map. I have the US county and state fips codes and the Y variable. How do I do this? Please reply to pzs6@cdc.gov Thanks, Phil Smith pzs6@cdc.gov [[alternative HTML version deleted]]

...e map package, and then used: library(maps) I can see that by using map( 'state' ) you get the state boundaries, also. How do I fill in different colors for the different states? I see there is a col parameter for map, but can't see how to get it to work! Please reply directly to: pzs6@cdc.gov Many Thanks! Phil Smith Centers for Disease Control and Prevention Atlanta [[alternative HTML version deleted]]

Hi R Community: I want to take the difference in two dates: dt2 - dt1. But, I want the answer in months between those 2 dates. Can you advise me? Please respond to: pzs6 at cdc.gov Thank you! Phil Smith Centers for Disease Control and Prevention

...lowed ... and I have trouble!! Finally, I check to see if ests$stfips has the same length as the number of rows in "ests"... > nrow( ests ) == length( as.character( ests$stfips ) ) [1] TRUE ... and it does!! Does anybody have a suggestion to solve this problem? Thanks! Phil Smith pzs6 at cdc.gov National Immunization Program Centers for Disease Control and Prevention Atlanta, GA

...strings of those labels (by using srt or crt). There is a "names" argument to boxplot, but I haven't had much luck controlling what it prints, the size of the font, and the character rotation. Can somebody enlighten me on how to do this? Please respond to my work email address: pzs6@cdc.gov Many thanks! Phil Smith Centers for Disease COntrol and Prevention [[alternative HTML version deleted]]

...only the long dash shows up in the box. The consequence of this is that the race/ethnic group that corresponds to the line type that is only a long dash cannot be distinguished from the legend. How do I stretch that legend box out so as to allow lty to draw longer line segments? Please reply to: pzs6 at cdc.gov Many thanks! Phil Smith Centers for Disease Control and Prevention

Hi R-people: I need a function to compute the number of months between 2 dates, in the same way a mother would do it. For example, if a kid is born on February 6, the number of months between that date and March 7 is exactly 1 month, although it is only 29 days. Thank you! Phil Smith CDC [[alternative HTML version deleted]]

...wrong in R, and 2) what the distinction between R and Splus is in this case. Thank you, Philip J. Smith, PhD Centers for Disease Control and Prevention National Immunization Program Immunization Services Division MS E-32 1600 Clifton Road, NE Atlanta, GA 30333 ph: 404 639 8729 fax: 404 639 3266 pzs6 at cdc.gov

Dear R-people: François Michonneau's method to obtain the special character plus/minus works on Windows 2000 professional. Many Thanks to François for his work! Phil Smith Centers for Disease Control and Prevention Atlanta, GA _____ From: François MICHONNEAU [mailto:francois.michonneau@gmail.com] Sent: Monday, March 20, 2006 1:58 PM To: r-help@stat.math.ethz.ch; Smith, Phil

Folks: I'm having trouble doing a forward variable selection using step() First, I fit an initial model: fit0 <- glm ( est~1 , data=all, subset=c(n>=25) ) then I invoke step(): fit1 <- step( fit0 , scope=list(upper=est~ pcped + pchosp + pfarm ,lower=est~1)) I get the error message: Error in eval(expr, envir, enclos) : invalid 'envir' argument I looked at the

Hi R-people: When I use the command to fit a model with an intercept, only: glm ( formula=haspdata ~ 1, data=dat, family=binomial, weights= dat$hy.wgt.s, subset=(dat$haspdat0!=3) ) I get the message: Warning message: non-integer #successes in a binomial glm! in: eval(expr, envir, enclos) Does anyone know what this means?? The data for this command is listed below. Thanks, Phil Smith CDC

Hi R people! :-) I am printing a table of percentages (P) and their 95% confidence interval half-widths (CI). I would like to (i) create a character string of P and CI that has the plus/minus character to the left of CI, then (ii) use write.csv to save that file with the P, CI and the special plus/minus character. I poured over the r-help archives and did not find advice on this. That

Hi All: I'm using win.metafile() to produce windows metafiles. When I use Word to insert the wmf picture, Word gives an error! I did my homework in posting this question, and couldn't find a fix. Any suggestions? Phil Smith CDC

Folks: With the help of David L. Reiner, I've developed a function that computes the number of months between 2 dates, x and y. num.months <- function ( x , y ) { x <- as.Date( x ) y <- as.Date( y ) seeq <- seq(from=x , to=y , by="months") ans <- length( seeq ) - 1 if ( max( seeq) > y ) ans <- ans -1 ans } To ease your reading this function, I've