[R] change plotting symbol for groups in trellis graph

2007-01-30 Thread Tom Boonen
Hi,

how can I change the plotting symbol for the groups in a trellis panel dotplot.

My graph is similar to:

library(trellis)
dotplot(variety ~ yield | site, data = barley, groups = year,
 key = simpleKey(levels(barley$year), space = right),
 xlab = Barley Yield (bushels/acre) ,
 aspect=0.5, layout = c(1,6), ylab=NULL)

I'd like to plot the different years with different symbols (eg. one
with pch=19 the other with pch=21 or so). Thank you very much for your
help.

Tom

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] as.data.frame(cbind()) transforming numeric to factor?

2006-08-18 Thread Tom Boonen
Dear List,

why does as.data.frame(cbind()) transform numeric variables to
factors, once one of the other variablesused is a character vector?

#
x.1 - rnorm(10)
x.2 - c(rep(Test,10))
Foo - as.data.frame(cbind(x.1))
is.factor(Foo$x.1)

Foo - as.data.frame(cbind(x.1,x.2))
is.factor(Foo$x.1)
#

I assume there is a good reason for this, can somebody explain? Thanks.

Best,
Tom

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] as.data.frame(cbind()) transforming numeric to factor?

2006-08-18 Thread Tom Boonen
Thanks everybody. I recognize my mistake now.

I think as.data.frame(cbind(x.1,x.2),stringsAsFactors = FALSE) would
be a good idea.

Tom

On 8/18/06, Prof Brian Ripley [EMAIL PROTECTED] wrote:
 On Fri, 18 Aug 2006, Tom Boonen wrote:

  Dear List,
 
  why does as.data.frame(cbind()) transform numeric variables to
  factors, once one of the other variablesused is a character vector?
 
  #
  x.1 - rnorm(10)
  x.2 - c(rep(Test,10))
  Foo - as.data.frame(cbind(x.1))
  is.factor(Foo$x.1)
 
  Foo - as.data.frame(cbind(x.1,x.2))
  is.factor(Foo$x.1)
  #
 
  I assume there is a good reason for this, can somebody explain? Thanks.

 Only if you can explain the good reason why you did not just use
 data.frame(x.1, x.2)!

 cbind() makes a matrix out of vectors, here a character matrix.  And then
 as.data.frame() converts character columns to factors.

 --
 Brian D. Ripley,  [EMAIL PROTECTED]
 Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
 University of Oxford, Tel:  +44 1865 272861 (self)
 1 South Parks Road, +44 1865 272866 (PA)
 Oxford OX1 3TG, UKFax:  +44 1865 272595


__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] ARMA(1,1) for panel data

2006-08-14 Thread Tom Boonen
Dear List,

I am new to TS-Modeling in R. I would like to fit an ARMA(1,1) model
for a balanced panel, running Y on a full set of unit and year dummies
using an arma(1,1) for the disturbance:

y_it=unit.dummies+yeardummies+e_it

where: e_it=d*e_it-1+u_it+q*u_it-1

How can I fit this model in R? arma() does not seem to take covariates
(or I don't understand how to specify the function so that it would).
Thank you very much.

Best, Tom

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.