similar to: panel argument in coplot

I am updating my lab manual in Spanish. How should I cite your supplements to Venables and Ripley 1997 and the scripts. Thank you very much Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232

I am using rw0641 with Windows 98. In the help for aov it states that the formula can sepecify multiple responses for a "maov". The help doesn't give an example, Venables and Ripley 1997:381 doesn't either. It isn't mentioned in the scripts nor in the section R Complements nor in Rnotes. I tried aov(c(y1,y2,y3)~x) which seemed reasonable to me but doesn't work. How

Is there any place where the Chambers and Hastie 1993 Statistical Models in S data sets can be found? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San

I am trying to assign values to a vector (pvalue). Similar code works in C but not in R. What am I doing wrong? r <- pvalue <- 0 for(i in (1:(k-1))){ for(j in (i+1):k){ r <- r+1 tstat <- (means[i]-means[j])/rms pvalue[r] <- 2*(1-pt(abs(tstat),df)) } } Thank you. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del

NCSS has a nice plot of the histogram and density plot superimposed. The S-PLUS 4 Guide to Statistics on page 45 suggests the following code for a general purpose density plot for a vector named data iqd <- summary(data)[5]-summary(data)[2] plot(density(data,width=2*iqd),xlab="x",ylab="",type="l") Is it possible to superimpose the density plot on a histogram

I have a vector Pes with 600 elements some of which are NA's. How can I form a vector of the indices of the NA's. for(i in 1:600) if(is.na(Pes[i])) print(i) prints the indices of the NA's but I can't figure out how to put the results in a vector. Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de

I am using rw0641 with Windows 98. To list just the influential repetitiones that result from "influence.measures", I am using the input result <- lm(y~x) and the code from the example in the help for "influence.measures" INFLM <- function(result){ inflm <- influence.measures(result) which(apply(inflm$is.inf,1,any)) } It works fine up to now with the

I am using rw0641. In my continuing quest to understand repeated measures analysis, I again return to lme. I exported the Potthoff and Roy data Orthodont.dat from S-PLUS 4.5 to avoid capture errors and ran the examples in the R help. I imported the data.frame with data <- read.table("Orthodont.dat",header=T) attach(data) and created the objects Orthodont.fit1 <-

I am using R 1.1 with Redhat 6.2 and RW 1.001 with Win98 (the upkey doesn't work on my IBM either as has been previously reported by others). The function aov doesn't return either the residuals or the residual degrees of freedom for split-plot designs. If you use the following code from Baron and Li's "Notes on the use of R for psycology experiments and questionnaires"

I am using RW 1.2.3. on an IBM PC 300GL. Using the data bp.dat which accompanies Helen Brown and Robin Prescott 1999 Applied Mixed Models in Medicine. Statistics in Practice. John Wiley & Sons, Inc., New York, NY, USA which is also found at www.med.ed.ac.uk/phs/mixed. The data file was opened and initialized with > dat <- read.table("bp.dat") >

Among other computers, I am using rw0642 on an IBM 300GL with 32MB RAM and Windows98. 1. If after opening rw0642 the first command is to help such as "? par", when the help window is closed by clicking X in the upper right corner of the help window the following message is shown This program has performed an illegal operation and will be shut down. If the problem persists contact the

I am using RedHat 6.1, the R1.1.0 binary download, and an HP Deskjet 692C. For some time I have been trying unsuccessfully to integrate R postscript graphics in LaTex. I consulted the Bug Tracking System. In the preamble of the LaTex file (test.tex) I have placed the line \usepackage{graphicx,color} and I use \begin{figure}[htbp] \begin{center} \includegraphics[height=4in]{graphic1.ps}

Hi, Does anyone know if the D'agostino test is available with R ? Alex

Hi! I would like to know, how do you plot more than one data series when using "coplot"? I think I know the answer if "plot" is being used: x <- 1:10 y <- 1:10 y2 <- 1:10* 1.1 plot ( y ~ x, type="n" ) points( y ~ x, type = "b", col = "red" ) # plot the points and lines for the first series points( y2 ~ x, type = "b", col =

R 1.9.1 on Win2000 or Win98SE. I am using coplot as follows: coplot(AVG~LRPI| REGION) the output seems normal but I get: "Warning message: calling par(new=) with no plot" This is the only explanation that I have for being unable to use par() with coplot for changing the way the xlab and ylab appears. From within coplot I can change the text itself but not the font, fontsize, etc,

Dear R-friends, For the following data: > xy x y i 1 731 0.313 2 2 739 0.340 2 3 790 0.373 2 4 855 0.451 2 5 980 0.608 2 6 575 0.156 1 7 608 0.207 1 8 630 0.249 1 9 670 0.332 1 10 838 0.377 1 11 964 0.466 1 > coplot(y ~ x|i, data=xy) coplot gives 3 panels, rather than 2, namely one for i=1 and one for i=2. Futhermore, when I extand data fram xy to have i=3 as follows:

Dear useRs, When I use coplot() and output to png/jpeg/bmp, the grid lines from the scatter plots disappear. If I output to pdf() the grid lines are there, however I can't use it - I have many points, and the resulting PDF file is large and very slow to open and scroll through. (By the way, if I click File-Save As-png/jpeg/bmp from Rgui.exe, the grid lines are preserved - but I need to use

Greetings: I am unsuccessful in suppressing "Given : myvariable" from a coplot. There was such a question in the past but the thread breaks down. I am sure this is a "for dummies" question :-). I tried: coplot(myvar~myvar | myvar, show.given=FALSE, xlab="....", ylab="...", main=" ") and some other variations (including without main=" ")

Try this: split(delta,cut(delta, breaks=c(0,2,4,7,10,max(delta)),include.lowest=T)) $"[0,2]" [1] 0.0 1.4 2.0 2.0 $"(2,4]" [1] 2.3 2.3 $"(4,7]" [1] 4.5 6.4 $"(7,10]" [1] 8.9 7.4 $"(10,12.4]" [1] 12.4 11.5 "Peter B. Mandeville" <mandevip at uaslp.mx>@stat.math.ethz.ch on 06/24/2002 13:11:41 Sent by: owner-r-help at

Hi all, I could not find a way to add the dashed horizontal separator line to the top panel of a coplot (the one indicating which panel belongs to which range of the given variable a in y~x | a). This line should separate the ranges plotted in a row of panels in the panel matrix below (hmmm, hope this is reasonably clear). I like it because it enhances readability. I'll try a sketch: >