similar to: Projects

G''day I''ve been performing an analysis of some data by a follow the leader method (Bad I know but it seems to work and I think I understand all the steps now) and I''ve run into a problem. The analysis is for all practical purposes identical to the non-linear mixed effects analysis performed on Ludbrooks Rabbit data in the second edition of Venables & Ripley (page

A non-text attachment was scrubbed... Name: not available Type: text Size: 2378 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980901/d2289901/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 2378 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980901/d2289901/attachment.pl

I am using R0.63.1 under Windows 95. Settings of <number> and <interval> when coplot is called are ignored, except for generating a warning. For example: > data(quakes) > coplot(long ~ lat | depth * mag, data = quakes, number=3) Warning: parameter "number" couldn't be set in high-level plot() function The warning is repeated ad nauseum. The following changes

A non-text attachment was scrubbed... Name: not available Type: text Size: 3126 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19991130/5cb00c0f/attachment.pl

Hello everybody, I'm writting some functions for experimental designs. The one I'm working on is similar to "fac.design" on Splus. The problem I have is with the form of the output : When the argument fraction is submitted, Splus gives something like : A B C 1 - - - 2 + + - 3 + - + 4 - + + Fraction : ~A:B:C The first part of this output is the

Hi, this is not a new doubt, but is a doubt that I cant find a good response. Look this output: > m.lme <- lme(Yvar~Xvar,random=~1|Plot1/Plot2/Plot3) > anova(m.lme) numDF denDF F-value p-value (Intercept) 1 860 210.2457 <.0001 Xvar 1 2 1.2352 0.3821 > summary(m.lme) Linear mixed-effects model fit by REML Data: NULL AIC BIC

I am currently plotting a graph with the function coplot() conditional on two factors which gives me 36 sub-panels. I have a few questions in general about changing the size of the graph. My question follows : a) I have a fair bit of text to put on the margin sides using mtext(). What is the best way to shrink the whole graph to say 80 - 90% so that I can include a bit more words. b) Since I am

Is there a quick way to make boxplots groups by two variables? By that I mean, that if x axes have values ("A","B","C"), than at each value there would be a few boxplots each for a value of second variable (say ("1","2","3")). Thank You. -- Vadim Kutsyy http://www.kutsyy.com vadim at kutsyy.com The University of

Is there a quick way to make boxplots groups by two variables? By that I mean, that if x axes have values ("A","B","C"), than at each value there would be a few boxplots each for a value of second variable (say ("1","2","3")). Thank You. -- Vadim Kutsyy http://www.kutsyy.com vadim at kutsyy.com The University of

Consider the following analysis of a class experiment done as a Latin Square: > spinner <- gl(4,4,16,label=c("Murray","Angela","Shasha","Stephen")) > order <- gl(4,1,16) > treat <- scan() 1: 1 2 4 3 5: 4 3 1 2 9: 3 4 2 1 13: 2 1 3 4 17: Read 16 items > coin <-

I am using R for the first time in one of my classes. My students have alerted me to a problem for which we have not found an answer. We find that some means returned by the model.tables function are not correct when missing data is present in analysis of variance problems. We have duplicated the problem using R 1.2.0, 1.2.1, and 1.2.2 under Windows 98 and several distributions of Linux (Redhat

There seems no default setting for tck, causing problems for plot.times() > library(chron) > dts <- dates(c("02/27/92", "02/27/92", "01/14/92", + "02/28/92", "02/01/92")) > y<-chron(dts,out.format=c("dmy","h:m:s")) > y [1] 270292 270292 140192 280292 010292 > plot(y,1:5) Error in axis(n, at =

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:

Following p.206 of "Statistical Models in S", I wish to change the code for summary.glm() so that it estimates the dispersion for binomial & poisson models when the parameter dispersion is set to zero. The following changes [insertion of ||dispersion==0 at one point; and !is.null(dispersion) at another] will do the trick: "summary.glm" <- function(object, dispersion =

I wonder whether there is any comment on the following. It is a summary of a 10 minute contribution that I made to a biostatistical workshop help yesterday in Sydney. My aim was to highlight the benefits of the Linux/R development model, in a session with the title: "How can statistical packages be improved?" =================================================================== The R

Depending on the setting of par()$usr, the 'at' setting in mtext(.., adj=0, outer=T) may cause the text to appear in an anomalous position (e. g. in the first instance below, at the left of the plot region rather than at 'at=0' in the figure region), or the text may not appear at all. If one does not set the 'at' parameter the text appears (with 'adj=0') on the

In another thread on this list, various wild allegations have been made, relating to the New York Times article on R. I object both to the subject line and to the content of several of the messages, and will not repeat or quote any of that content. It smacks to me of mischief making. Discussion has centered around the following quote from the NY Times article: ?According to them, the

Hello ! I'm struggling again against lattice graprhics. ;) I'm trying to produce a conditionnal xyplot with two conditionning factors (let's say A and B). I want the levels of those factors (A1, A2, etc) to show in the margins of the lattice plot, not in the strips between the panels. A1 A2 A3 plot11 plot12 plot13 B1 plot21 plot22 plot23 B2 I managed to remove the

COPLOT: how do i define the given.values for two variables? Thomas Pesl -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch