similar to: rw0642

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

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

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

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

How can I add a least squares line to each panel in coplot? 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 Luis Potos? Fax: 48 28-23-52 Av.

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") >

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}

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"

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

1) rwinst Brian Ripley has provided me a new version of 'rwinst.exe' which should solve the signaled problem (installing in the root of the hard-disk). It should be available tomorrow from Wien (and in few days from the other mirrors). 2) gif There was some problems in rw0641. It used to work before and it should work again in rw0642 (at least, it does on my system). So, please, if you

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

thanks to all the help so far I have code now that works on windows (albeit very clunky) that can read the pertinent info from the shape header and return a list with all the necessary parameters to read info out of the shape file. Now come the big design questions 1) I have considered building a map object, basically a list that would hold the info from the shape file and a data frame as part of

I just downloaded the NT version of R and cannot find anything concrete on loading packages. There is alot on how to do it under Unix. The windows FAQ briefly mentions it but does not explain it. I have donwloaded windows binaries of the package "clust.zip" and unzipped it into my c:\rw0632\library directory... So I have c:\rw0632\library\cluster Now what do I do ?

Using rw0642 under win95, on executing from a directory on the c: drive: system("command.com /c dir", intern=T) there is a whirring from drive a: (empty), then the message appears in a dos box: General failure reading device aLP!!? Abort, Retry, Fail? Repeating the command a few times finally results in success with no whirring, after which there are no failures (even

Greetings, Cohen, Cohen, West, and Aiken 2003 (Applied Multiple Regression-Correlation Analysis for the Behavioral Sciences, Third Edition) on page 273 state the covariance matrix of the regression coefficients is provided by standard programs for multiple regression, including SAS, SPSS, and SYSTAT. How can it be calculated with R. Thank you very much. pbm Peter B. Mandeville cel:

I have been trying to get Xgobi to work with the function xgobi( ), but I keep on getting the message ""Xgobi requires a filename or some data from stdin" in the DOS-shell. My Xgobi.bat file is : @echo off PATH=c:\windows\ set DISPLAY=195.23.18.11:0.0 set XGOBIDIR=c:/Xgobi3 c:\Xgobi3\xgobi.exe In R I get: > xgobi(D1) C:\RW0642/library/xgobi/scripts/xgobi.bat -vtitle

Dear friends. I have trobpe with path. Below I have tried to paste a path and file name but it will not work > path <- "D:/rw0641/own/own/procedures/GFR_TIME" > path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good