Re: [R] cross-classified random factors in lme without blocking - success

2003-11-22 Thread Gordon Smyth
Sundar's advice seems to do the trick. Here is a small simulation of a cross-classified random model with extraction of the fitted variance components from lme: > library(nlme) > set.seed(18112003) > na <- 20 > nb <- 20 > sigma.a <- 2 > sigma.b <- 3 > sigma.res <- 1 > mu <- 0.5 > n <- na*nb > a

Re: [R] plotting and overlay

2003-11-22 Thread Prof Brian Ripley
On Sat, 22 Nov 2003, Duncan Murdoch wrote: > On Sat, 22 Nov 2003 13:54:34 -0800, Spencer Graves > <[EMAIL PROTECTED]> wrote: > > >Start the plot using "plot", add lines and points using "lines" and > >"points", as described in the examples with "?plot". > > That's how you do plots that you con

Re: [R] options(width= ) in rw1081

2003-11-22 Thread Prof Brian Ripley
On Sat, 22 Nov 2003 [EMAIL PROTECTED] wrote: > I just downloaded and installed the factory-fresh > rw1081. Windows XP. > > In my Rprofile file in the \etc subdirectory I put > (among others) > > > options(continue=" ") > options(width=65) > options(scipen=5) > options(show.signif.stars=FALSE)

RE: [R] where to get the "leaps" package

2003-11-22 Thread Liaw, Andy
Once you start up Rgui, click on the menu "Packages" -> "Install packages from CRAN", and then select "leaps" from the list and click OK. Andy > From: [EMAIL PROTECTED] > > Dear all > I am supposed to use Mallow's Cp creterion to select a model > which require a "leaps" package. the version ri

[R] options(width= ) in rw1081

2003-11-22 Thread kjetil
Hola! I just downloaded and installed the factory-fresh rw1081. Windows XP. In my Rprofile file in the \etc subdirectory I put (among others) options(continue=" ") options(width=65) options(scipen=5) options(show.signif.stars=FALSE) R is honoring all of this on startup, exept the width=65 op

Re: [R] where to get the "leaps" package

2003-11-22 Thread kjetil
On 22 Nov 2003 at 18:12, Yong Wang wrote: You should find it on CRAN, the Comprehensive R Archival Network. Search google for CRAN. or if you are on windows, go to menu Packages --- install packages from CRAN. Kjetil Halvorsen > Dear all > I am supposed to use Mallow's Cp creterion to select a

[R] where to get the "leaps" package

2003-11-22 Thread Yong Wang
Dear all I am supposed to use Mallow's Cp creterion to select a model which require a "leaps" package. the version right now I am using is R 1.7.1(os:windows), the "leaps" package is not included. please let me know how and where can I get this package and add it to the current version, any caution

Re: [R] plotting and overlay

2003-11-22 Thread kjetil
On 22 Nov 2003 at 13:37, christoff pale wrote: In R many graphics functions take an argument add, and there are some functions (points, lines) which add data to an existing graph. Some examples: plot( dnorm, from=-3, to=3 ) plot( function(x) dt(x, 2), from=-3, to=3, col="red2", add=TRUE) x <-

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Al Piszcz
That is what the last suggestion was and it works, see previous message from me. Thx again. I'm changing my code to use this method now. On Sat, 22 Nov 2003, Duncan Murdoch wrote: > Date: Sat, 22 Nov 2003 17:10:39 -0500 > From: Duncan Murdoch <[EMAIL PROTECTED]> > To: Al Piszcz <[EMAIL PROTECT

Re: [R] plotting and overlay

2003-11-22 Thread Duncan Murdoch
On Sat, 22 Nov 2003 13:54:34 -0800, Spencer Graves <[EMAIL PROTECTED]> wrote: >Start the plot using "plot", add lines and points using "lines" and >"points", as described in the examples with "?plot". That's how you do plots that you construct yourself; you can get a lot of the higher level fun

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Duncan Murdoch
On Sat, 22 Nov 2003 15:14:31 -0500 (EST), you wrote: > >I would like to have one handle or reference to >'n' matrices. The matrices vary in size. >All data is floating point. >The input files have 21 columns and a varying >number of rows. > >I am open to any data structure that will >support this.

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Al Piszcz
Thank you very much! This method works with R 1.8.0 and R 1.8.1 R : Copyright 2003, The R Foundation for Statistical Computing Version 1.8.1 (2003-11-21), ISBN 3-900051-00-3 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions.

Re: [R] plotting and overlay

2003-11-22 Thread Spencer Graves
Start the plot using "plot", add lines and points using "lines" and "points", as described in the examples with "?plot". hope this helps. spencer graves christoff pale wrote: Hi, I am use to matlab and when I want to overlay multiple plots on the same graphics devise I type in "hold on" and

[R] plotting and overlay

2003-11-22 Thread christoff pale
Hi, I am use to matlab and when I want to overlay multiple plots on the same graphics devise I type in "hold on" and when I am done I type in "hold off" Does R have an equivalent feature? can someone give me an example where they plotted 2 functions on top of each other? thanks _

Re: [R] 1.8.1 behavior change? W EXAMPLE

2003-11-22 Thread Al Piszcz
Here is a working example with R 1.8.0, followed by the test files, and finally the R 1.8.1 error messages. [1] R 1.8.0 WORKING EXAMPLE R R : Copyright 2003, The R Development Core Team Version 1.8.0 (2003-10-08) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Al Piszcz
I would like to have one handle or reference to 'n' matrices. The matrices vary in size. All data is floating point. The input files have 21 columns and a varying number of rows. I am open to any data structure that will support this. On Sat, 22 Nov 2003, Prof Brian Ripley wrote: > Date:

Re: [R] ISOdate() and strptime()

2003-11-22 Thread Prof Brian Ripley
Confirmation that this *is* an OS-specific problem: A professional implementation of the POSIX standard (Solaris) gets all of these correct. Your so-called OS lacks any implementation of strptime, so we borrowed one from glibc. Unfortunately, that is buggy, even to the extent that unclass(strp

[R] arima {ts}

2003-11-22 Thread Gopal Annasundaram
I'm trying to get more info about seasonal parameter of arima(ts). Can someone explain on its usage? I'm also interested in building a model for weekly sales for an item with covariates like consumer confidence index, promotion flag etc. Which package has functions to model and do predictions

[R] Warning: X11 protocol error: BadWindow (invalid Window parameter)

2003-11-22 Thread Korponai János
Hi all, I've got this using Rcmdr package only. R-cmdr> scatterplot.matrix(~PC1+PC2+PC3+PC4 | Species, reg.line=FALSE, smooth=FALSE, span=0.5, diagonal= 'density', by.groups=TRUE, data=iris) Warning: X11 protocol error: BadWindow (invalid Window parameter) Warning: X11 protocol error: BadWindow (

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Duncan Murdoch
On Sat, 22 Nov 2003 11:57:32 -0500 (EST), you wrote: > >In it creates the following warning/error: > >Any advice appreciated. > > > stt <- data.frame() > # load all datasets into a dataframe > for (ds in 1:n) { >stt[ds] <- as.matrix(read.table(fileList[ds])) > } I don't know what you are

Re: [R] 1.8.1 behavior change?

2003-11-22 Thread Prof Brian Ripley
On Sat, 22 Nov 2003, Al Piszcz wrote: > In stt [1] V1 <0 rows> (or 0-length row.names) so is that what `worked properly' means? (It is also what S+6.1 does, but I am pretty sure it is not what anyone wanted. Note the number of columns is wrong, too.) In 1.7.1 it gave a similar error to 1.8.1

RE: R crashes with package SJava; was [R] Memory leakage?

2003-11-22 Thread Takeshi Takama
Dear All, <> R-1.81 j2sdk1.4.0_03 SJava_0.66-1 Windows NT4.01 I try to run R from Java with SJava and I seem to have the same problem as discussed before in this mail list a few month ago: https://www.stat.math.ethz.ch/pipermail/r-help/2003-May/031960.html https://www.stat.math.e

[R] 1.8.1 behavior change?

2003-11-22 Thread Al Piszcz
Instt <- data.frame() > # load all datasets into a dataframe > for (ds in 1:n) { + stt[ds] <- as.matrix(read.table(fileList[ds])) + } Error in "[<-.data.frame"(`*tmp*`, ds, value = as.matrix(read.table(fileList[ds]))) : replacement has 358 rows, data has 0 >

Re: [R] R-1.8.1 is released

2003-11-22 Thread Duncan Murdoch
On 21 Nov 2003 13:45:53 +0100, you wrote: >I've rolled up R-1.8.1.tgz a short while ago. This is a patch version >mostly fixing a number of issues in 1.8.0, some of which were quite >serious. As usual, a few new features have crept in as well. (See below >for details.) > >You can get it from > >ht

[R]: how to plot smooth function estimate from gam (mgcv package) in other program

2003-11-22 Thread Stephanie Klot
Hi all, I would like to export the smooth function estimate I got from gam to plot it in another graphics software. In S-plus I use the function preplot() for that, but it seems not to work in R. Has somebody an idea how to solve that? Thanks Stephanie Stephanie

RE: [R] lm with ordered factors

2003-11-22 Thread Ted Harding
On 22-Nov-03 Ted Harding wrote: > [...] > I have a quantitative variable Y and a 4-level ordered factor A > (with very unequal numbers at the different levels, by the way). > > The command > > lm(Y ~ A) > > returns (amongst other stuff) an intercept, and coefficients > A.L, A.Q and A.C for the

Re: [R] lm with ordered factors

2003-11-22 Thread Prof Brian Ripley
On Sat, 22 Nov 2003 [EMAIL PROTECTED] wrote: > Hi Folks, > > No doubt a question with a well-known answer, but I'm unfortunately > not managing to find it readily ... ! > > I have a quantitative variable Y and a 4-level ordered factor A > (with very unequal numbers at the different levels, by th

[R] lm with ordered factors

2003-11-22 Thread Ted Harding
Hi Folks, No doubt a question with a well-known answer, but I'm unfortunately not managing to find it readily ... ! I have a quantitative variable Y and a 4-level ordered factor A (with very unequal numbers at the different levels, by the way). The command lm(Y ~ A) returns (amongst other st

[R] RAqua 1.8.1 binaries available on CRAN

2003-11-22 Thread Stefano Iacus
New binaries for RAqua 1.8.1 are now available on CRAN (main site) and will be propagated soon over the network. Stefano __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help