search for: cryan

...1 other 2001 12 2 UHS 2001 17 3 Lourdes 2001 10 4 Endwell 2001 3 5 other 2002 13 6 UHS 2002 19 > What is the easiest way to make a line graph like this (made this one in Stata): http://bingweb.binghamton.edu/~cryan/Workload.png Thank you. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu PGP public keys available at http://home.stny.rr.com/ryancw/ "If you want to build a ship, don't drum...

...e fitting and assessing loglinear models for several nominal variables . . . something akin to step or drop1 or add1 for linear or logistic regression? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu "If you want to build a ship, don't drum up the men to gather wood, divide the work and give orders. Instead, teach them to yearn for the vast and endless sea." [Antoine de St. Exupery]

Running R 2.13.1 on Windows XP. I would like to get week of the year (1-52) for each date. library(chron) dts <- dates(c("02/27/92", "02/27/92", "01/14/92","02/28/92", "02/01/92")) dts dts.chron <- as.chron(dts) dts.chron class(dts.chron) # all of these component extractions work: months(dts.chron) weekdays(dts.chron) years(dts.chron)

...cal vector x <- as.logical(c(0,0,0,1,0,0,1,1,0)) x How would I make the words TRUE appear on the screen in a different color from the words FALSE? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu "If you want to build a ship, don't drum up the men to gather wood, divide the work and give orders. Instead, teach them to yearn for the vast and endless sea." [Antoine de St. Exupery]

...n use the default method: > > plot(as.vector(time(ttt)), as.vector(ttt), type = "p", col=ifelse(ttt<8, > "black", "red")) > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Tue, May 1, 2018 at 1:17 PM, Christopher W Ryan <cryan at binghamton.edu> > wrote: > >> How would I color points conditional on their value in a plot of a time >> series. Something like this: >> >> ## demonstration data >> ttt <- ts(rpois(12, lambda = 8), start = c(2000, 1), freq = 4) >> ttt >> plo...

The ts method for plot() is quirky. You can use the default method: plot(as.vector(time(ttt)), as.vector(ttt), type = "p", col=ifelse(ttt<8, "black", "red")) Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, May 1, 2018 at 1:17 PM, Christopher W Ryan <cryan at binghamton.edu> wrote: > How would I color points conditional on their value in a plot of a time > series. Something like this: > > ## demonstration data > ttt <- ts(rpois(12, lambda = 8), start = c(2000, 1), freq = 4) > ttt > plot(ttt, type = "p") > &gt...

...int on the graph, I get an error, "no point within 0.25 inches." But in some areas, I can click where there is no visible point anywhere close, and an identify() label will pop up. The troublesome code is shown below. I've posted an example graph here: http://bingweb.binghamton.edu/~cryan/junkidentify2.pdf I'd be grateful if anyone could take a look at it. The text() labels are in red, and the results of pointing-and-clicking are shown in black. Notice point 4, for example: near the origin in black, but upper center in red (which is where the real data point should be: point 4...

How would I color points conditional on their value in a plot of a time series. Something like this: ## demonstration data ttt <- ts(rpois(12, lambda = 8), start = c(2000, 1), freq = 4) ttt plot(ttt, type = "p") ## doesn't work--all points the same color plot(ttt, type = "p", col = ifelse(ttt < 8, "black", "red")) ## also doesn't work--all

Based on the feedback received, I did the following: a) moved my lib sub-directory from the existing installed R version to c:\myRLib b) installed the updated R version c) created .Renviron file in the home directory (C:\R-2.5.1) with the line R_LIBS=c:/myRLib d) used .libPaths() command to confirm that the new R installation was recognizing the myRLib sub-directory e) deleted my old R

...gt; plot(as.vector(time(ttt)), as.vector(ttt), type = "p", col=ifelse(ttt<8, >> "black", "red")) >> >> >> Bill Dunlap >> TIBCO Software >> wdunlap tibco.com >> >> On Tue, May 1, 2018 at 1:17 PM, Christopher W Ryan <cryan at binghamton.edu> >> wrote: >> >>> How would I color points conditional on their value in a plot of a time >>> series. Something like this: >>> >>> ## demonstration data >>> ttt <- ts(rpois(12, lambda = 8), start = c(2000, 1), freq =...

I participate peripherally on a listserve for middle- and high-school science teachers. Sometimes questions about graphing or data analysis come up. I never miss an opportunity to advocate for R. However, the teachers are often skeptical that their students would be able to issue commands or write a little code; they think it would be too difficult. Perhaps this stems from the Microsoft- and

I'm trying to learn Sweave. So far things are going well with the chunks of code identified by << >>= But I'm having trouble with the in-line text use of \Sexpr. Here is a short example .Rnw file: \documentclass[12pt]{article} \usepackage[margin=1.25in]{geometry} \usepackage{graphicx} \usepackage{Sweave} \begin{document} <<test>>= pnorm(1) @ \Sexpr{2+6}

...ologies, forgot to copy R-help on this response. Begin forwarded message: From: Roger Koenker <rkoenker at illinois.edu> Subject: Re: [R] effects() extractor for a quantile reqression object: error message Date: September 6, 2024 at 8:38:47?AM GMT+1 To: "Christopher W. Ryan" <cryan at binghamton.edu> Chris, This was intended to emulate the effects component of lm() fitting, but was never implemented. Frankly, I don?t quite see on first glance how this works for lm() ? it seems to be (mostly) about situations where X is not full rank (see lm.fit) and I also never bothere...

running R 2.5.1 on WinXP. plot(model formula here) produces 4 useful residual diagnostic plots. Two questions about them: 1. How can I isolate just one of the four, and have only that one be shown? str(plot(model)) didn't give me any insight. 2. How can I make the lines (I guess they are lowess lines) not appear on the plots? Thanks. --Chris Ryan

I've written a simple script that does some surveillance analysis on daily counts of walk-in clinic visits, for our county health department. (Actually, Michael Hohle's surveillance package does all the work; I just customized it a little to work with the way our data are recorded.) The output consists of a graph and a little table. My colleague will run it once in a while, on

...jaw$age^i } for (i in 2:6) { jaw$age."i" <- jaw$age^i } What is the proper syntax to create and name these variables all at once? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu PGP public keys available at http://home.stny.rr.com/ryancw/ "If you want to build a ship, don't drum up the men to gather wood, divide the work and give orders. Instead, teach them to yearn for the vast and endless sea." [Antoine de St. Exupery]

...t;<chunkname>>= and @. Several of them don't contain any R code. Can I compile the whole document with Sweave("Master.Rnw") ? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu PGP public keys available at http://home.stny.rr.com/ryancw/ "If you want to build a ship, don't drum up the men to gather wood, divide the work and give orders. Instead, teach them to yearn for the vast and endless sea." [Antoine de St. Exupery]

Great question. What do I want? I want my co-workers to stop using Excel spreadsheets for data entry, storage, and sharing! I want them to understand the value of data discipline. But alas . . . . I work in a county health department in the US. Between dplyr, stringr, grep, grepl, and the base R read() functions, I'm doing OK. I need to learn more about APIs, so I can see if I can make R

...16 16 5 12 10 6 13 11 7 18 NA 8 8 NA . . . . cut off for brevity . . . -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu PGP public keys available at http://home.stny.rr.com/ryancw/ "If you want to build a ship, don't drum up the men to gather wood, divide the work and give orders. Instead, teach them to yearn for the vast and endless sea." [Antoine de St. Exupery]

I have a dataframe that looks like this: > str(chr) 'data.frame': 84 obs. of 7 variables: $ county: Factor w/ 3 levels "Broome","Nassau",..: 3 3 3 3 3 3 3 3 3 3 ... $ item : Factor w/ 28 levels "Access to healthy foods",..: 21 19 20 18 16 3 2 6 17 8 ... $ value : num 8644 15 3.5 3.9 7.7 ... $ low : num 7897 9 2.5 2.6 7 ... $ high : num 9390