RE: [R] manipulating "..." inside a function

2003-03-21 Thread Liaw, Andy
Thanks to those who responded.  I finally got by with some pretty ugly hack
(ceorcing ... to list, delete the unwanted part, add other arguments, and
use do.call).

Cheers,
Andy

> -Original Message-
> From: Ben Bolker [mailto:[EMAIL PROTECTED]
> Sent: Friday, March 21, 2003 3:59 PM
> To: Liaw, Andy
> Cc: '[EMAIL PROTECTED]'
> Subject: Re: [R] manipulating "..." inside a function
> 
> 
> 
>   I have a couple of slightly ugly functions to do this in my bbmisc 
> package:
> 
> clean.args <- function(argstr,fn,extrabad=NULL,dots.ok=TRUE) {
>   fnargs <- names(formals(fn))
>   if (length(argstr)>0 && !("..." %in% fnargs && dots.ok))  {
> badargs <- !sapply(names(argstr),"%in%",c(fnargs,""))
> argstr <- argstr[!badargs]
>   }
>   for (i in extrabad)
> argstr[[i]] <- NULL
>   argstr
> }
> 
> remove.args <- function(argstr,fn) {
>   fnargs <- names(formals(fn))
>   argstr[!(names(argstr) %in% fnargs)]
> }
> 
> 
> On Fri, 21 Mar 2003, Liaw, Andy wrote:
> 
> > Dear R-help,
> > 
> > Can some one tell me how to do the following (if it's possible)?
> > 
> > Suppose I have a function like this:
> > 
> > f <- function(x, y, ...) {
> > ## some code
> > g(x, y, ...)
> >## some more code
> > }
> > 
> > The problem is that g() may not understand everything that 
> comes through in
> > "...".  Is there a way to delete some component of "..." 
> and then pass it to
> > g()?
> > 
> > Here's the description of the real problem:  f() is a 
> panel.something
> > function, and g() is a model fitting function.  Lattice passes
> > "panel.number" as part of "..." to f(), and g() complains 
> about unused
> > argument "panel.number".
> > 
> > I'd be very grateful for any help!
> > 
> > Cheers,
> > Andy
> > 
> > 
> > 
> --
> 
> > 
> > __
> > [EMAIL PROTECTED] mailing list
> > https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> > 
> 
> -- 
> 318 Carr Hall[EMAIL PROTECTED]
> Zoology Department, University of Florida
> http://www.zoo.ufl.edu/bolker
> Box 118525   (ph)  352-392-5697
> Gainesville, FL 32611-8525   (fax) 352-392-3704
> 
> 

--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] manipulating "..." inside a function

2003-03-21 Thread Ben Bolker

  I have a couple of slightly ugly functions to do this in my bbmisc 
package:

clean.args <- function(argstr,fn,extrabad=NULL,dots.ok=TRUE) {
  fnargs <- names(formals(fn))
  if (length(argstr)>0 && !("..." %in% fnargs && dots.ok))  {
badargs <- !sapply(names(argstr),"%in%",c(fnargs,""))
argstr <- argstr[!badargs]
  }
  for (i in extrabad)
argstr[[i]] <- NULL
  argstr
}

remove.args <- function(argstr,fn) {
  fnargs <- names(formals(fn))
  argstr[!(names(argstr) %in% fnargs)]
}


On Fri, 21 Mar 2003, Liaw, Andy wrote:

> Dear R-help,
> 
> Can some one tell me how to do the following (if it's possible)?
> 
> Suppose I have a function like this:
> 
> f <- function(x, y, ...) {
> ## some code
> g(x, y, ...)
>## some more code
> }
> 
> The problem is that g() may not understand everything that comes through in
> "...".  Is there a way to delete some component of "..." and then pass it to
> g()?
> 
> Here's the description of the real problem:  f() is a panel.something
> function, and g() is a model fitting function.  Lattice passes
> "panel.number" as part of "..." to f(), and g() complains about unused
> argument "panel.number".
> 
> I'd be very grateful for any help!
> 
> Cheers,
> Andy
> 
> 
> --
> 
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> 

-- 
318 Carr Hall[EMAIL PROTECTED]
Zoology Department, University of Floridahttp://www.zoo.ufl.edu/bolker
Box 118525   (ph)  352-392-5697
Gainesville, FL 32611-8525   (fax) 352-392-3704

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] manipulating "..." inside a function

2003-03-21 Thread Wiener, Matthew
As Brian Ripley pointed out in a recent post, you can just give g() its own
"..." argument.

Regards, 

Matt Wiener

-Original Message-
From: Liaw, Andy [mailto:[EMAIL PROTECTED] 
Sent: Friday, March 21, 2003 3:37 PM
To: '[EMAIL PROTECTED]'
Subject: [R] manipulating "..." inside a function


Dear R-help,

Can some one tell me how to do the following (if it's possible)?

Suppose I have a function like this:

f <- function(x, y, ...) {
## some code
g(x, y, ...)
   ## some more code
}

The problem is that g() may not understand everything that comes through in
"...".  Is there a way to delete some component of "..." and then pass it to
g()?

Here's the description of the real problem:  f() is a panel.something
function, and g() is a model fitting function.  Lattice passes
"panel.number" as part of "..." to f(), and g() complains about unused
argument "panel.number".

I'd be very grateful for any help!

Cheers,
Andy



--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


--
Notice: This e-mail message, together with any attachments, contains
information of Merck & Co., Inc. (Whitehouse Station, New Jersey, USA) that
may be confidential, proprietary copyrighted and/or legally privileged, and
is intended solely for the use of the individual or entity named on this
message.  If you are not the intended recipient, and have received this
message in error, please immediately return this by e-mail and then delete
it.


==


--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] manipulating "..." inside a function

2003-03-21 Thread Liaw, Andy
Dear R-help,

Can some one tell me how to do the following (if it's possible)?

Suppose I have a function like this:

f <- function(x, y, ...) {
## some code
g(x, y, ...)
   ## some more code
}

The problem is that g() may not understand everything that comes through in
"...".  Is there a way to delete some component of "..." and then pass it to
g()?

Here's the description of the real problem:  f() is a panel.something
function, and g() is a model fitting function.  Lattice passes
"panel.number" as part of "..." to f(), and g() complains about unused
argument "panel.number".

I'd be very grateful for any help!

Cheers,
Andy


--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Re:problem with read.table

2003-03-21 Thread sabrina servanty
You were right!Thank a lot, I think that I'm really dumm!I've lost 3 hours trying to solve the problem!!Regards,Sabrina Servanty__Gagne une PS2 ! Envoie un SMS avec le code PS au 61166(0,35 Euro Hors coût du SMS)

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Problem with read.table

2003-03-21 Thread Spencer Graves
Have you tried "count.fields" to confirm that R thinks all records have 
the same number of fields?

Spencer Graves

sabrina servanty wrote:
Dear all,

I was used to work on R1.6 and I have now passed on R1.6.2 but I can't 
read my
file (and that is a big problem!!).
I made a data sheet with some
spreadsheet in Excell, and save it as separeted by tab .txt.
I write in R
read.table ("file.txt",h=T,sep="/t",dec=",")
But R consider that I have only one column (eG one variable)!!!

I have tried a lot of thing (I don't wrote the spreadsheet,I have 
verified in word
that my column was really separated by tabulation...) but I really don't 
find.
It may be really simple but I'm not really good to speak with R!

Regards

Sabrina Servanty.

__
Boîte aux lettres - Caramail - http://www.caramail.com


__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Problem with read.table

2003-03-21 Thread Jerome Asselin

How about replacing "/t" by "\t" ?

Jerome

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] Problem with read.table

2003-03-21 Thread Marc Schwartz
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of sabrina
servanty
Sent: Friday, March 21, 2003 1:22 PM
To: [EMAIL PROTECTED]
Subject: [R] Problem with read.table


Dear all, 

I was used to work on R1.6 and I have now passed on R1.6.2 but I can't
read my 
file (and that is a big problem!!).
I made a data sheet with some 
spreadsheet in Excell, and save it as separeted by tab .txt.
I write in R
read.table ("file.txt",h=T,sep="/t",dec=",")
But R consider that I have only one column (eG one variable)!!!

I have tried a lot of thing (I don't wrote the spreadsheet,I have
verified in word 
that my column was really separated by tabulation) but I really
don't find.
It may be really simple but I'm not really good to speak with R!

Regards

Sabrina Servanty.
-

Sabrina,

Try reversing the "/" in your 'sep = "/t"' argument.  You have it
reveresed at the moment, assuming that is is copied here as you are
using it.

The tab character is:  "\t"

So the command should be:

read.table ("file.txt",h=T,sep="\t",dec=",")

HTH,

Marc Schwartz

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Problem with read.table

2003-03-21 Thread Ben Bolker

  Should that be sep="\t" ?

On Fri, 21 Mar 2003, sabrina servanty wrote:

> Dear all,
> 
> I was used to work on R1.6 and I have now passed on R1.6.2 but I can't read my
> file (and that is a big problem!!).
> I made a data sheet with some
> spreadsheet in Excell, and save it as separeted by tab .txt.
> I write in R
> read.table ("file.txt",h=T,sep="/t",dec=",")
> But R consider that I have only one column (eG one variable)!!!
> 
> I have tried a lot of thing (I don't wrote the spreadsheet,I have verified in word
> that my column was really separated by tabulation...) but I really don't find.
> It may be really simple but I'm not really good to speak with R!
> 
> Regards
> 
> Sabrina Servanty.
> 
> __
> Boîte aux lettres - Caramail - http://www.caramail.com
> 
> 
> 

-- 
318 Carr Hall[EMAIL PROTECTED]
Zoology Department, University of Floridahttp://www.zoo.ufl.edu/bolker
Box 118525   (ph)  352-392-5697
Gainesville, FL 32611-8525   (fax) 352-392-3704

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R]

2003-03-21 Thread Wagner Silva


-- 
=
  Wagner Silva
[EMAIL PROTECTED]
  Bacharel em Ciência da Computação
Universidade Federal de Lavras

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] rsync

2003-03-21 Thread Erich Neuwirth


Seth Falcon wrote:

I am now downloading r-devel to compile it on windows XP. The CRAN
source code page says "you will prefere to use rsync". I am googling 
around, and cannot find anything about rsync on windows. 
   

You can use rsync on Windows via the Cygwin toolset (see
www.cygwin.com).  
 

If you do that and also want to use Brian Ripley's toolkit for compiling,
you might get into trouble because of incompatible versions of
cygwin1.dll
I managed to get a working solution by using rsync from

http://optics.ph.unimelb.edu.au/help/rsync/rsync_pc1.html

and putting rsync.exe into the directory with brian's binaries (sed, 
grep, diff ...)
this rsync vversion can work with brian's version of cygwin1.dll.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Problem with read.table

2003-03-21 Thread sabrina servanty
Dear all, I was used to work on R1.6 and I have now passed on R1.6.2 but I can't read my file (and that is a big problem!!).I made a data sheet with some spreadsheet in Excell, and save it as separeted by tab .txt.I write in Rread.table ("file.txt",h=T,sep="/t",dec=",")But R consider that I have only one column (eG one variable)!!!I have tried a lot of thing (I don't wrote the spreadsheet,I have verified in word that my column was really separated by tabulation...) but I really don't find.It may be really simple but I'm not really good to speak with R!RegardsSabrina Servanty.__Boîte aux lettres - Caramail - http://www.caramail.com

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] Install R on unix

2003-03-21 Thread Huntsinger, Reid
It looks like the Makefiles "make" is looking at aren't the ones configure
created. I did the install on solaris 2.6 with no problem just now. The line
about "Stata 6 installation" is certainly suspicious... maybe some makefiles
got overwritten somehow?

Can you start clean and try again?

Reid Huntsinger


-Original Message-
From: pingzhao [mailto:[EMAIL PROTECTED] 
Sent: Thursday, March 20, 2003 5:00 PM
To: [EMAIL PROTECTED]
Subject: [R] Install R on unix


Hello,

I just tried to install R (unix version), but I have 
not successed.

I download the file and uncompresssed it in a R folder.
then I used the following commands:

./configure
make

The attached are error message:
Could any one tell me how to fix it??

Thanks


R is now configured for sparc-sun-solaris2.6

  Source directory:  .
  Installation directory:/usr/local

  C compiler:gcc  -g -O2
  C++ compiler:  g++  -g -O2
  Fortran compiler:  f77  -g

  X11 support:   yes
  Gnome support: no
  Tcl/Tk support:yes
  Readline support:  no

  R profiling support:   yes
  R as a shared library: no

  Recommended packages:  yes

configure: WARNING: you cannot build DVI versions of the R manuals
configure: WARNING: you cannot build info versions of the R manuals
configure: WARNING: you cannot build PDF versions of the R manuals
slri05:/home/pingzhao/R/R-1.6.2> make
creating src/scripts/R.fe
mkdir -p -- ../../bin


Stata 6 installation


The installation files are missing

(no action taken)
*** Error code 1
make: Fatal error: Command failed for target `install-cmds'
Current working directory /home/pingzhao/R/R-1.6.2/src/scripts
*** Error code 1
make: Fatal error: Command failed for target `R'
Current working directory /home/pingzhao/R/R-1.6.2/src/scripts
*** Error code 1
make: Fatal error: Command failed for target `R'
Current working directory /home/pingzhao/R/R-1.6.2/src
*** Error code 1
make: Fatal error: Command failed for target `R'

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] trellis plot

2003-03-21 Thread Jerome Asselin

Hi Latif,

You need to build your own panel function that will fit your purpose.
This will do what you want, but it's not very elegant. A better solution 
would have the panel function depend on the value of z. Any other suggestions 
on own to do this?

Jerome

library(lattice)

x <- c(rep(LETTERS[1:4],13), rep(LETTERS[4:1],12))

y <- rnorm(100)
z <- rep(1:2,50)
x <- c("A","B",x)
y <- c(-1.5,-2.5,y)
z <- c(1:2,z)
bwplot(y~factor(x)|z,layout=c(2,1), panel=function(x,y) 
 {
   panel.bwplot(x[-1],y[-1],horizontal=F)
   panel.xyplot(x[1],y[1],pch=20,cex=2,col="red")
 })


On March 21, 2003 06:00 am, you wrote:
> 
> Hi there,
> 
> I need some help about trellis plot. I have the following plot.
> 
> x <- c(rep(LETTERS[1:4],13), rep(LETTERS[4:1],12))
> 
> y <- rnorm(100)
> z <- rep(1:2,50)
> bwplot(y~factor(x)|z,layout=c(2,1),  panel=function(x,y) 
panel.bwplot(x,y,horizontal=F))
> 
> 
> Now I want to place "*" on the positions (1,-1.5) in the first panel and 
(2,-2.5) in the second panel. I need help on this.
> 
>  
> 
> Thanks.
> 
> 
> 
> Mahbub.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] question

2003-03-21 Thread Roger Peng
Try using the `pixmap' package from CRAN.

-roger
___
UCLA Department of Statistics
[EMAIL PROTECTED]
http://www.stat.ucla.edu/~rpeng

On Fri, 21 Mar 2003, vincent deschodt wrote:

> How can i open a  bmp image and transform it in a matrix of pixels?
> how can i save it after transform on matrix?
> (Sorry for my english i don't speak it very well)
> 
>   [[alternate HTML version deleted]]
> 
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] rsync

2003-03-21 Thread Seth Falcon
> I am now downloading r-devel to compile it on windows XP. The CRAN
> source code page says "you will prefere to use rsync". I am googling 
> around, and cannot find anything about rsync on windows. 

You can use rsync on Windows via the Cygwin toolset (see
www.cygwin.com).  

I've also seen an rsync for windows written in Python, but haven't
tried it.  

+ seth

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] rsync

2003-03-21 Thread Ko-Kang Kevin Wang
Hi,

Perhaps have a look at 
http://optics.ph.unimelb.edu.au/help/rsync/rsync_pc1.html "Installing 
rsync on a Windows machine"?

FYI, I did a google search on "Using rsync on Windows" and found it.

On Fri, 21 Mar 2003, kjetil brinchmann halvorsen wrote:

> Date: Fri, 21 Mar 2003 13:35:20 -0400
> From: kjetil brinchmann halvorsen <[EMAIL PROTECTED]>
> To: [EMAIL PROTECTED]
> Subject: [R] rsync
> 
> Hola!
> 
> I am now downloading r-devel to compile it on windows XP. The CRAN
> source code page says "you will prefere to use rsync". I am googling 
> around, and cannot find anything about rsync on windows. 
> 
> Anybody has any experience with rsync on windows?
> 
> Kjetil
> 
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> 

-- 
Cheers,

Kevin

--
/* Time is the greatest teacher, unfortunately it kills its students */

--
Ko-Kang Kevin Wang
Master of Science (MSc) Student
SLC Tutor and Lab Demonstrator
Department of Statistics
University of Auckland
New Zealand
Homepage: http://www.stat.auckland.ac.nz/~kwan022
Ph: 373-7599
x88475 (City)
x88480 (Tamaki)

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] Savitzky-Golay Derivative and Smoothing

2003-03-21 Thread Liaw, Andy
If I'm not mistaken, that's sort of local polynomial with even degree and
fixed bandwidth (based on my own interpretation of description in Numerical
Recipes).  You can do that with functions in the KernSmooth package.

HTH,
Andy

> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> Sent: Friday, March 21, 2003 1:43 PM
> To: [EMAIL PROTECTED]
> Subject: [R] Savitzky-Golay Derivative and Smoothing
> 
> 
> Hello,
> 
> Is there any libary with the algorithms of the Savitzky-Golay 
> Derivative
> and Smoothing. I found the calculation on the web site
> "www.galactic.com/algorithms/" but I'm to new in R so I cant 
> programm it in
> R.
> Can someone help me?
> 
> Thanx
> Andreas
> 
> micro-biolytics GmbH
> Andreas Wolf
> 
> Georges-Köhler-Allee 102
> D - 79110 Freiburg
> Fon: +49-761-2037524
> Fax: +49-761-2037522
> 
> eMail: [EMAIL PROTECTED]
> Web:   http://www.micro-biolytics.com
> 
> - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> CONFIDENTIALITY NOTICE
> 
> This e-mail is intended only for the addressee named above 
> and may contain
> confidential information. If you are not the intended 
> recipient, please let
> us know immediately by return e-mail and then delete this 
> e-mail, without
> disclosing or copying the contents.
> 
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> 

--

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Savitzky-Golay Derivative and Smoothing

2003-03-21 Thread wolf
Hello,

Is there any libary with the algorithms of the Savitzky-Golay Derivative
and Smoothing. I found the calculation on the web site
"www.galactic.com/algorithms/" but I'm to new in R so I cant programm it in
R.
Can someone help me?

Thanx
Andreas

micro-biolytics GmbH
Andreas Wolf

Georges-Köhler-Allee 102
D - 79110 Freiburg
Fon: +49-761-2037524
Fax: +49-761-2037522

eMail: [EMAIL PROTECTED]
Web:   http://www.micro-biolytics.com

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
CONFIDENTIALITY NOTICE

This e-mail is intended only for the addressee named above and may contain
confidential information. If you are not the intended recipient, please let
us know immediately by return e-mail and then delete this e-mail, without
disclosing or copying the contents.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] rsync

2003-03-21 Thread kjetil brinchmann halvorsen
Hola!

I am now downloading r-devel to compile it on windows XP. The CRAN
source code page says "you will prefere to use rsync". I am googling 
around, and cannot find anything about rsync on windows. 

Anybody has any experience with rsync on windows?

Kjetil

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Trying to make a nested lme analysis

2003-03-21 Thread Ronaldo Reis Jr.
Hi,

I'm trying to understand the lme output and procedure.
I'm using the Crawley's book.

I'm try to analyse the rats example take from Sokal and Rohlf (1995).
I make a nested analysis using aov following the book.

> summary(rats)
Glycogen   Treatment  Rat  Liver  
 Min.   :125.0   Min.   :1   Min.   :1.0   Min.   :1  
 1st Qu.:135.8   1st Qu.:1   1st Qu.:1.0   1st Qu.:1  
 Median :141.0   Median :2   Median :1.5   Median :2  
 Mean   :142.2   Mean   :2   Mean   :1.5   Mean   :2  
 3rd Qu.:150.0   3rd Qu.:3   3rd Qu.:2.0   3rd Qu.:3  
 Max.   :162.0   Max.   :3   Max.   :2.0   Max.   :3  

> attach(rats)
> Treatment <- factor(Treatment)
> Rat <- factor(Rat)
> Liver <- factor(Liver)

> model <- aov(Glycogen~Treatment/Rat/Liver+Error(Treatment/Rat/Liver))
> summary(model)

Error: Treatment
  Df  Sum Sq Mean Sq
Treatment  2 1557.56  778.78

Error: Treatment:Rat
  Df Sum Sq Mean Sq
Treatment:Rat  3 797.67  265.89

Error: Treatment:Rat:Liver
Df Sum Sq Mean Sq
Treatment:Rat:Liver 12  594.049.5

Error: Within
  Df Sum Sq Mean Sq F value Pr(>F)
Residuals 18 381.00   21.17   
> 

OK,

Then I try to make this analysis using lme.

> model <- lme(Glycogen~Treatment, random=~1|Treatment/Rat/Liver)
> summary(model)
Linear mixed-effects model fit by REML
 Data: NULL 
   AIC  BIClogLik
  233.6213 244.0968 -109.8106

Random effects:
 Formula: ~1 | Treatment
(Intercept)
StdDev:3.541272

 Formula: ~1 | Rat %in% Treatment
(Intercept)
StdDev: 6.00658

 Formula: ~1 | Liver %in% Rat %in% Treatment
(Intercept) Residual
StdDev:3.764883 4.600247

Fixed effects: Glycogen ~ Treatment 
Error in if (any(wchLv <- (as.double(levels(xtTab[, wchPval])) == 0))) { 
: 
missing value where logical needed
In addition: Warning message: 
NaNs produced in: pt(q, df, lower.tail, log.p) 
> 

The random effects are correct, the variance component is OK:

In nested aov | In nested lme
Residual
21.1666   | 21.16227
Liver in Rats
14.16667  | 14.17434
Rats in Treatment
36.0648   | 36.079

But I not understand why the Fixed effects error?

What is the problem in my formula to make this analysis using lme?

Thanks for all
Inte
Ronaldo


--
|   //|\\   [*]
|| ( õ õ )  [Ronaldo Reis Júnior  ]
| V [ESALQ/USP-Entomologia, CP-09 ]
||  / l \   [13418-900 Piracicaba - SP]
|  /(lin)\  [Fone: 19-429-4199 r.229  ]
||/(linux)\ [EMAIL PROTECTED]  ]
|/ (linux) \[ICQ#: 5692561]
||  ( x )   [*]
||| _/ \_ Powered by Gnu/Debian Woody
---
Insecta - Entomologia
Departamento de Biologia Animal
Universidade Federal de Viçosa
---

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] ArrayAnalyzer and Bioconductor

2003-03-21 Thread Robert Gentleman
Today, Insightful Corporation (www.insightful.com) is announcing the
availability of S+ArrayAnalyzer, a new, integrated module for S-PLUS
based on collaboration with the BioConductor Project
(www.bioconductor.org) - an open source and open development software
project for the analysis and comprehension of genomic data.

The collaboration between Insightful and BioConductor delivers
benefits to both commercial and academic researchers analyzing
microarray experiments:

- BioConductor's advanced analytics for genomic data will be available
  to a wider audience via a commercially supported product from
  Insightful, a long-established vendor with fully staffed tech
  support, training, and consulting.

- S+ArrayAnalyzer adds numerous features that can improve productivity
  for many users, such as installation, more data access and import
  options, a guided-workflow interface, interactive graphs with
  hyperlinked annotation, Web deployment of applications.

- Insightful will sponsor a graduate student position in the
  BioConductor project to help drive continued advancements and
  innovation in both the open source and commercial offerings.

- A new differential expression library (lpetest), written by
  Insightful and the University of Virginia, will be ported to R and
  made available on the Comprehensive R Archive Network (CRAN) and
  BioConductor in Summer 2003.

Spokespeople for BioConductor see the collaboration as helping ensure
the distribution and high-level end-user support of key research tools
to the broadest possible population of researchers. These people
include: S.Dudoit, Division of Biostatistics, University of
California, Berkeley; R.A. Irizarry, Department of Biostatistics,
Johns Hopkins University; V.J. Carey, Harvard Medical School;
R. Gentleman, Harvard School of Public Health.

Shawn Javid, Insightful's CEO had this to add, "The collaboration with
BioConductor is a blueprint for how Insightful can work with the open
source community to bring highly innovative data analysis applications
to the widest possible user base.  This cooperation improves the
analytic solutions available to our common base of S programmers and
the non-statisticians who benefit from using applications developed
with S-PLUS and R."

S+ArrayAnalyzer is available now with pricing for commercial and
academic organizations available by calling (800)569-0123 x479, or via
email at [EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] barplot legend size

2003-03-21 Thread Marc Schwartz
>-Original Message-
>From: [EMAIL PROTECTED] 
>[mailto:[EMAIL PROTECTED] On Behalf Of Ott Toomet
>Sent: Friday, March 21, 2003 8:44 AM
>To: [EMAIL PROTECTED]
>Subject: [R] barplot legend size
>
>
>Dear R-people,
>
>are there any way to change the size of legend in barplot?  I 
>have tried various versions of cex, both as par(cex.*= ) and 
>barplot(..., cex.*= ).  So long without success.
>
>Sincerely,
>
>Ott

Ott,

When you say the size of the legend I presume that you mean the text
within the legend given that you are trying to use par("cex").

Instead of creating the legend within the barplot() call (ie.
barplot(, legend = ) ), remove the 'legend' argument and call
legend() after the barplot call with your legend specifications.

See ?legend, which provides additional formatting flexibility,
including position, text size and spacing options among others.

HTH,

Marc Schwartz

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] barplot legend size

2003-03-21 Thread Ott Toomet
Dear R-people,

are there any way to change the size of legend in barplot?  I have
tried various versions of cex, both as par(cex.*= ) and barplot(...,
cex.*= ).  So long without success.

Sincerely,

Ott

> version
 _
platform i686-pc-linux-gnu
arch i686 
os   linux-gnu
system   i686, linux-gnu  
status
major1
minor5.1  
year 2002 
month06   
day  17   
language R
-- 
Ott Toomet
PhD Student

Dept. of Economics
Århus University
Building 322
Universitetsparken
8000 Århus C
Denmark

[EMAIL PROTECTED]
ph: (+45) 89 42 20 27
---

 (o_ (*_ (O_ (o< -!  
//\ //\ //\ //\  
V_/_V_/_V_/_V_/_ 
 
standarddrunken shocked noisy
penguin penguin penguin penguin

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Plot multi series on one plot

2003-03-21 Thread Michael A. Miller
> "mhoward" == mhoward  <[EMAIL PROTECTED]> writes:

> R help, How can I plot the below data table by Data ~ Site
> and group by Grinder and Equip Id so I get a chart like
> this Excel version? I have tried coplot with little success
> and lattice makes a pretty good chart like I want, but I am
> using the DCOM so it does not display correctly. I would
> like to make this using the base library if possible.

>  <<...OLE_Obj...>>

Your excel plot did not come along with your message, so I'm not
quite sure what you're looking for.  I won't let that stop me
from suggesting something though!

I know you've said lattice is not working for you, but I still
recommend it - something like

 df <- read.table('tmp.dat',header=T)
 dotplot(SITE~DATA|GRINDER*EQUIPID, data=df)

Other wise you can use dotchart, which is in the base library.
Try something like this:

 dotchart(tapply(df$DATA, list(df$SITE,paste(df$GRINDER,df$EQUIPID)), mean))

Mike

-- 
Michael A. Miller   [EMAIL PROTECTED]
  Imaging Sciences, Department of Radiology, IU School of Medicine

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Redundancy in Variables (and canonical correlation)

2003-03-21 Thread Barker, Chris

 Base on the earlier question: there is a "redundancy" measure for
Canonical Correlation "Canonical Correlation Redundancy Analysis". You
may be able to adapt it to your situation. 
 
See the 
paper by 
Stewart David W., Love W. (1968), A General Canonical Correlation Index,
Psychological Bulletin, 70, 160-163.

There are dozens of references to it and "canonical correlation
redundancy "
on a Google Search (many are pointers to the online SAS help). 





Chris Barker
   Director of Statistical Research
MEDTAP International, Inc.
   Redwood City, Ca
 
www.medtap.com
 
  650 632 4218

>

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] trellis plot

2003-03-21 Thread Mahbub Latif
Hi there,
I need some help about trellis plot. I have the following plot.
x <- c(rep(LETTERS[1:4],13), rep(LETTERS[4:1],12))
y <- rnorm(100)z <- rep(1:2,50)bwplot(y~factor(x)|z,layout=c(2,1),  panel=function(x,y) panel.bwplot(x,y,horizontal=F)) 
Now I want to place "*" on the positions (1,-1.5) in the first panel and (2,-2.5) in the second panel. I need help on this.
 
Thanks.
 
Mahbub.

 
 

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] Integer manipulation

2003-03-21 Thread Thomas Lumley
On Fri, 21 Mar 2003, Wayne Jones wrote:

> Hi there,
>
> When I enter a particularly long number into R it rounds it down into
> scientific notation.
> For example,
>
> > 251002679
>
> [1] 2.51e+12
>
> How can I preserve the precision of the original number?

It has been preserved:
> a<-251002679
> 25100-a
[1] -2679

For numbers above 2^31 (or whatever .Machine$integer.max is), you have to
use floating point, which gives you about 15 digits accuracy.

If you want to print a number to more digits than the default, use the
digits= option to print(), or change options(digits)

-thomas

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Integer manipulation

2003-03-21 Thread Wayne Jones
Hi there, 

When I enter a particularly long number into R it rounds it down into
scientific notation.
For example, 

> 251002679

[1] 2.51e+12

How can I preserve the precision of the original number?

Regards, 

Wayne



Dr Wayne R. Jones
Statistician / Research Analyst
KSS Group plc
St James's Buildings
79 Oxford Street
Manchester M1 6SS
Tel: +44(0)161 609 4084
Mob: +44(0)7810 523 713



KSS Ltd
A division of Knowledge Support Systems Group plc
Seventh Floor  St James's Buildings  79 Oxford Street  Manchester  M1 6SS  England
Company Registration Number 2800886 (Limited) 3449594 (plc)
Tel: +44 (0) 161 228 0040   Fax: +44 (0) 161 236 6305
mailto:[EMAIL PROTECTED]http://www.kssg.com


The information in this Internet email is confidential and may b... [[dropped]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] question

2003-03-21 Thread vincent deschodt
How can i open a  bmp image and transform it in a matrix of pixels?
how can i save it after transform on matrix?
(Sorry for my english i don't speak it very well)

[[alternate HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Postscript PB

2003-03-21 Thread Poizot Emmanuel
Hi,
I use R 1.6.2 under Mandrake9.0.
I've got a problem with the postscript files I try to creat.
When I look to the file with ghostview it's ok.
When I want to print it, I've got a blank page or a black page (fill of black 
encre)
I changed the printer (guessing it was my driver printer), it was the same.
Does anyone had the same problem and resolved it ?

-- 
Cordialement

Emmanuel POIZOT
Cnam/Intechmer
Digue de Collignon
50110 Tourlaville
Tél : (33)(0)2 33 88 73 42
Fax : (33)(0)2 33 88 73 39
-

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Postscript PBs

2003-03-21 Thread Poizot Emmanuel
Hi,
I use R 1.6.2 under Mandrake9.0.
I've got a problem with the postscript files I try to creat.
When I look to the file with ghostview it's ok.
When I want to print it, I've got a blank page or a black page (fill of black 
encre)
I changed the printer (guessing it was my driver printer), it was the same.
Does anyone had the same problem and resolved it ?

-- 
Cordialy

Emmanuel POIZOT
Cnam/Intechmer
Digue de Collignon
50110 Tourlaville
Tél : (33)(0)2 33 88 73 42
Fax : (33)(0)2 33 88 73 39
-

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] CCA (Curvilinear Component Analysis)

2003-03-21 Thread Anne-Laure Boulesteix
Hi,
Does anyone know if  Curvilinear Component Analysis (CCA) or Curvilinear 
Distance Analysis (CDA) has already been implemented in R ? I couldn't 
find it.

Many thanks
Anne-Laure
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Re: [BioC] mva functions

2003-03-21 Thread alessandro . semeria
Probably you have to set a new max value for the max 
amount of the "heap" memory, look at the R-help on 'memory.limit' 
function.

Good luck 
A.S.



Alessandro Semeria 
Models and Simulations Laboratory
The Environment Research Center - Montecatini (Edison Group), 
Via Ciro Menotti 48,
48023 Marina di Ravenna (RA), Italy
Tel. +39 544 536811
Fax. +39 544 538663
E-mail: [EMAIL PROTECTED]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help