[R] Re: [BioC] mva functions
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
[R] Postscript PBs
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
Re: [R] Integer manipulation
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] trellis plot
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] barplot legend size
-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] Trying to make a nested lme analysis
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] 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
[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
Re: [R] rsync
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] rsync
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] question
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] trellis plot
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] rsync
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]
-- = 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] Problem with read.table
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
RE: [R] Problem with read.table
-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
How about replacing /t by \t ? Jerome __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Re:problem with read.table
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
[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
RE: [R] manipulating ... inside a function
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
