similar to: aov with Error terms

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --0-488534781-901753860=:6003 Content-Type: TEXT/PLAIN; CHARSET=US-ASCII Content-ID: <Pine.NEB.4.00.9807291716501.6003@mathfs.math.montana.edu> Attached is R code for an aov

R developers, 2 things: a bug in model.frame and a question about setup of generic functions. I don't understand the following behavior for generic functions: Suppose I'm working with the cats data in the MASS library and I want to create a formula object to model Hwt on Sex: # This works: > formula(Hwt ~ Sex) Hwt ~ Sex # But the following does not: > formula(Hwt ~ Sex,

R developers, Last week I asked several questions which I think I can now answer myself in large part. > > Q 1) Is aov supposed to avoid qr decomposition and matrix inversion? > In S it is claimed to be faster than lm() for large datasets. Is > that due to avoidance of qr()? > (My main goal is to get the Error strata working.) My use of matrix

Dear R-help-ers (Currently running Version 0.65.1 Release (October 07, 1999) on Solaris 2.6) My problem has to do with updates of the R language. When I install a new version, I would like to just change one directive, or a soft link to the current doc/html/ directory, to update a bookmark I have pointing to the current R doumentation directory. I've tried building a soft link:

Dear Rhelp folks, I think I read this in recent Rhelp, but I can't find it in the archive. I need to know how to get lg2c in order to build the nlme package. I'm trying to install nlme_3.1-2.tar.gz on a Linux machine (6.1 Redhat) under R Version 0.99.0 (February 7, 2000). Thanks, Jim Jim Robison-Cox ____________ Department of Math Sciences | |

Just after R-0.63.0 came out I reported a problem to the list. I was unable to get ./configure to complete it's job, because it couldn't find a working fortran. It stopped after trying to test for underscore in subroutine names. This was especially perplexing, since I could still configure the older version (I had 62.2) with no error messages. I tried suggestions of B.D. Ripley and

Dear R-helpers, I just upgraded from version 62.1 to 62.2 and now am faced with the same problem I've run into before when upgrading: Is there an easy way to take the libraries I have been using in the older version and pull then into the new version. My organization looks like this on Solaris 2.5.2: /opt -- /R -- /R-0.62.1 --/library -- subdirectories for lots of

When I run configure for the new release, R-0.63, I get an error: loading cache ./config.cache checking for a BSD compatible install... aux/install-sh -c checking whether ln -s works... yes checking for ranlib... ranlib checking for bison... bison -y checking for ar... ar checking for ratfor... no checking for latex... /psoft/TeX/bin/latex checking for dvips... /psoft/TeX/bin/dvips checking for

When I run configure for the new release, R-0.63, I get an error: loading cache ./config.cache checking for a BSD compatible install... aux/install-sh -c checking whether ln -s works... yes checking for ranlib... ranlib checking for bison... bison -y checking for ar... ar checking for ratfor... no checking for latex... /psoft/TeX/bin/latex checking for dvips... /psoft/TeX/bin/dvips checking for

R-crew, I'm trying to install the R-Snapshot which I downloaded a few hours ago and I'm having a problem I can't trace. My system is Sun SPARC, SunOS 5.3 : 1996 >From running ./configure:*************************************** loading cache ./config.cache checking for a BSD compatible install... /opt/local/R/R-snapshot/etc/install-sh -c checking whether ln -s works... yes

Dear R-help folks, I am seeing unexpected behaviour from the function mean with option na.rm =TRUE (which is removing a whole column of a data frame or matrix. example: testcase <- data.frame( x = 1:3, y = rep(NA,3)) mean(testcase[,1], na.rm=TRUE) [1] 2 mean(testcase[,2], na.rm = TRUE) [1] NaN OK, so far that seems sensible. Now I'd like to compute both means at once:

Hello, Perhaps this is a bug, but I thought I'd start with R-Help. I'm trying to compile R-1.4.1 on a new Sun-Blade running Solaris 2.8 ( I got it running with no problem on older SUNs) I get the following error after doing the usual ./configure and make: gcc -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -g -O2 -c arithmetic.c -o arithmetic.o

Hello List, Knowing I would eventually need to ask for help on the list, I've held off my thanks and congratulations on the 1.+ releases of R. I am amazed at the dedication of the R-core team, and do add my applause and thanks. I've had a compile problem with the latest version which has not occurred before. I'm compiling under: platform sparc-sun-solaris2.6 arch sparc

R developers, I have a first attempt to make an aov function. Eventually I want to build in Error() structure, but first I am trying to get this presentable for balanced data with only a single stratum, just using residual error. I am following R. M. Heiberger's Computation for the Analysis of Designed Experiments, Wiley (1989) I a using a wrapper (aov.bal) to call the

Hi, I'm trying to install rjags and R2jags with R version 2.14.1 (2011-12-22) on Linux Mint 12 machine using software repository: deb http://cran.fhcrc.org/bin/linux/ubuntu oneiric/ jags seems to install OK, but rjags gives a weird error about not finding fprec symbol. Here is a history of the problem including several attempts to work around it. (Because of Ben Bolker's recent

To those who might be interested -- following is the solution to my previous post regarding reproducing output from SAS Proc Mixed for a two-factor crossed random effects ANOVA model. I am graciously endebted to the kind replys from two statisticians for this solution whose names I will refrain from mentioning for the sake of privacy. I hope this helps someone?! -- Phil Turk > hw7 <-

I have data that looks like lake,loglength,logweight 1,2.369215857,1.929418926 1,2.426511261,2.230448921 1,2.434568904,2.298853076 1,2.437750563,2.298853076 1,2.442479769,2.230448921 1,2.445604203,2.356025857 ... 102,2.722633923,3.310268367 102,2.781755375,3.502153893 102,2.836324116,3.683407299 102,2.802773725,3.583312152 102,2.790285164,3.546419267 102,2.806179974,3.599118565

I am trying to produce frequencies in defined intervals however I can't seem to figure out how to get R to bin my data the way I want it to. I have several thousand lengths of fish that I want to be binned as follows: Ex. Length Bin 209 200 219 210 431 430 727 720 That is, bins with any length equal to or greater than the lower

What code can i use to convert a table like this: Tag# Date 1 1 2 1 3 1 4 1 2 2 4 2 1 3 2 3 4 4 Into one like this: Tag 1 2 3 4 #Date header 1 1 0 0 1 2 1 1 1 0 3 1 0 0 0 4 1 1 0 1 Thanks, Ben Cox Research Assistant (M.S.) Montana Cooperative

All - Does anyone know if there is a method to calculate a goodness-of-fit statistic for quantile regressions with package quantreg? Specifically, I'm wondering if anyone has implemented the goodness-of-fit process developed by Koenker and Machado (1999) for R? Though I have used package quantreg in the past, I may have overlooked this function, if it is included. Citation: Koenker, R. and