similar to: automatic exploration of all possible loglinear models?

I can't seem to find just what I'm looking for in R help, Everitt and Hothorn HSAUR, Murrell's book, or the R graphics gallery at http://addictedtor.free.fr/graphiques/. Probably not looking efficiently, but anyway, If my data look like this: > head(data) cat startyear studentid 1 other 2001 12 2 UHS 2001 17 3 Lourdes 2001 10 4

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

suppose I have some logical 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

R 2.5.1 on WinXP I'm trying to create new variables in a dataframe called jaw, as powers of jaw$age up to the sixth power, and name them jaw$age.(--the digit corresponding to the power--) Obviously none of these work (silly for me to try, I suppose): for (i in 2:6) { jaw$age.'i' <- jaw$age^i } for (i in 2:6) { jaw$age.i <- jaw$age^i } for (i in 2:6) {

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)

Suppose I have a Master.Rnw file that looks something like this: \documentclass[12pt]{mypaper} \usepackage[margin=1in]{geometry} \usepackage{setspace} \usepackage{url} \usepackage{indentfirst} \usepackage{fancyhdr} \usepackage{Sweave} \pagestyle{fancy} \lhead{sonographic rectal diameter and ADHD} \rhead{ } \usepackage{abbrevs} %\usepackage{natbib} %\usepackage{apacite}

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

I start with a dataframe called xrays. It contains scores on films from each of two radiologists. It is in "long" format. I used the reshape package to melt a data frame and then cast it into "wide" format, one line for each patient (identified by redlognumb) with scores from both radiologists for a given patient on the same line. I named the result of the casting xrays.data.

I can't seem to use the point-and-click identify() function properly. I'm running R 2.5.1 (I know, I need to get around to upgrading) under Win XP. The problem is, when I click on a point 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

Excellent! Worked like a charm. Thanks. --Chris Ryan On Tue, May 1, 2018 at 4:33 PM, William Dunlap <wdunlap at tibco.com> wrote: > 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

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

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

Hi all, I've been running loglinear models for three-way tables: one of the variables having three levels, and the other two having two levels each. An example looks like below: > yes.no <- c("Yes","No") > switch <- c("On","Off") > att <- c("BB","AA","CC") > L <- gl(2,1,12,yes.no) > T <-

You may also want to check this out: plot(ttt, type = "p") points(ttt, col = ifelse(ttt < 8, "black", "red")) Eivind K. Dovik Bergen, NO On Tue, 1 May 2018, Christopher W Ryan wrote: > Excellent! Worked like a charm. Thanks. > > --Chris Ryan > > On Tue, May 1, 2018 at 4:33 PM, William Dunlap <wdunlap at tibco.com> wrote: > >> The

Dear netters, Sorry for cross-posting this question. I am sure R-Help is not a research methods discussion list, but we have many statisticians in the list and I would like to hear from them. Any function/package in R would be able to deal with the problem from this researcher? ---------- Forwarded message ---------- From: Heidi Bertels Date: Tue, Jun 5, 2012 at 4:31 PM Subject: How to best

Hi all: Here's a question about result of loglinear analysis. There're 2 factors:area and nation.The raw data is in the attachment. I fit the saturated model of loglinear with the command: glm_sat<-glm(fre~area*nation, family=poisson, data=data_Analysis) After that,I extract the coefficients: result_sat<-summary(glm_sat) result_coe<-result_sat$coefficients I find that all the

Dear R Users, A would like to fit a loglinear analysis to a three dimensional contingency table. But I Don't want to run a full saturated modell. Is there any package in R that could handle somekind of stepwise search to choose out the best soultion? And how can I fit a non fully saturated modell, which only use the important interactions? Best Regards Zoltan Kmetty [[alternative HTML

I am using the loglin function of the base package to fit log-linear models. I am interested in obtaining the parameter values and their standard errors. Parameters are easily obtained, but I haven't found the way of obtaining their standad errors. Is this possible with the loglin function? If not, is there any other function to get them? Many thanks, -- Vicente Piorno Departamento de Ecolox?a

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}

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