Re: [R] pmvt problem in multcomp

2004-05-26 Thread Chihiro Kuroki
Dear Mr.Torsten: At Sun, 23 May 2004 11:40:51 +0200 (CEST), Torsten Hothorn wrote: Yes, `pmvt' returns NaN without indicating this error. We need to check. Thanks for the report (and *please* cc emails reporting problems with packages to the maintainer!), I am dreadfully sorry I did not cc

RE: [R] Metafiiles into Word R 1.9.0

2004-05-26 Thread Jesus Frias
Hi again, Thanks a lot to Jason Turner for the hint to the FAQ-windows. savePlot() seems to be in trouble with lattice objects as much as any other device, which I did not know. I can generate now the metafiles: library(lattice) example(xyplot) win.metafile() print(dotplot(variety ~

[R] Using R in C++

2004-05-26 Thread Jörgen Wallerman
Hello, Is it possible to use R functions (in my case: ks.test()) from C++ -applications? That is, I get the impression R can execute C/C++ code, but is there any possibility to do the opposite? Where can I find help? --- Ph. D. Jörgen Wallerman Swedish

[R] is.weekend() odd behaviour

2004-05-26 Thread Matthew Dowle
It seems is.weekend() is unsure ? # Start R 1.9.0 with --vanilla on windows xp # load package chron table(is.weekend(sapply(1:100, function(i){Sys.sleep(0.05);Sys.time()}))) FALSE TRUE 6832 date() [1] Wed May 26 11:18:56 2004 version _ platform

[R] Once again in search of memory

2004-05-26 Thread Eric Lecoutre
Hello any R user, Working on R 1.8.1 or R 1.9 on Windows, 256 Mb RAM. I am trying to debug/correct/speedup the execution of one of our student's program. Basically, we perform simulations. R is launched with high memory-limits. One execution of the program requires nearly 200 Mbs and is OK the

[R] Common principle components

2004-05-26 Thread J. Pedro Granadeiro
Dear all, Can anyone point me to a package that can perform Common principle components? Thank you. José P. Granadeiro __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

AW: [R] Question

2004-05-26 Thread TEMPL Matthias
A small comment: The code of Agnes is written in Fortran. Following book give more details: Kaufman, L. and Rousseeuw, P.J. (1990). _Finding Groups in Data: An Introduction to Cluster Analysis_. Wiley, New York. The 'hclust' function is based an Algorithm contributed to STATLIB by F.

Re: [R] is.weekend() odd behaviour

2004-05-26 Thread Uwe Ligges
Matthew Dowle wrote: It seems is.weekend() is unsure ? It is completely sure, if an object is given is.weekend() knows about as in: table(is.weekend(sapply(1:100, function(i){Sys.sleep(0.05);as.chron(Sys.time())}))) (note the as.chron() call!) Uwe Ligges # Start R 1.9.0 with --vanilla on

[R] multi read.table function and read.table function not accepting col.names

2004-05-26 Thread Laura Quinn
I am using R-1.8.0 on Debian. I'm trying to read in a large table (1441*16) which currently has no header line. I have set up a list of column names which is 16 names long. when i try the following: myfiledate.01-read.table(filenamedate.01,row.names=NULL,col.names=names,na.strings=-999.00) I

RE: [R] is.weekend() odd behaviour

2004-05-26 Thread Matthew Dowle
Thanks. When is.weekend() is given an object it doesn't know about, could a warning or error be added? At the moment, I think you're agreeing, that it silently returns a random TRUE/FALSE. -Original Message- From: Uwe Ligges [mailto:[EMAIL PROTECTED] Sent: 26 May 2004 11:55 To: Matthew

[R] Page's trend test?

2004-05-26 Thread Vaclav Petricek
Hello Is there an R implementation of the Page's trend test as described in http://en.wikipedia.org/wiki/Page%27s_trend_test or something equivalent? Thanks for your help, Vaclav __ [EMAIL PROTECTED] mailing list

Re: [R] multi read.table function and read.table function not accepting col.names

2004-05-26 Thread Wolski
Hi! 1. assuning that names is the vector with the names. try col.names=names instead of names 2. Several ways a) use lapply(listwithpathstofiles,read.tabel,remaining, options) //u will get it stored in the list. b) or use for loop and append to list. mylist-list() c) or look ?assing to

Re: [R] is.weekend() odd behaviour

2004-05-26 Thread Gabor Grothendieck
Uwe Ligges ligges at statistik.uni-dortmund.de writes: Matthew Dowle wrote: table(is.weekend(sapply(1:100, function(i){Sys.sleep(0.05);Sys.time()}))) FALSE TRUE 6832 It seems is.weekend() is unsure ? It is completely sure, if an object is given is.weekend() knows about

[R] Venn diagrams

2004-05-26 Thread Wayne Jones
Hello there, Does anyone know how to plot venn diagrams in R? Ive searched the mail archive lists but to no avail. Regards, Wayne Jones KSS Ltd Seventh Floor St James's Buildings 79 Oxford Street Manchester M1 6SS England Company Registration Number 2800886 Tel: +44 (0) 161 228 0040

Re: [R] Venn diagrams

2004-05-26 Thread Chuck Cleland
http://www.stats.uwo.ca/faculty/murdoch/software/ Also, a google search [Venn r-project] returns a number of promising hits. hope this helps, Chuck Cleland Wayne Jones wrote: Does anyone know how to plot venn diagrams in R? Ive searched the mail archive lists but to no avail. -- Chuck

Re: [R] multi read.table function and read.table function not accepting col.names

2004-05-26 Thread Laura Quinn
Thank you for your earlier help, I have a couple of related questions. With the lapply function, is it still possible to assign column names to each column of every data frame, again I have tried and it doesn't seem to like it.. and is it also possible to assign names to each data frame when I

Re: [R] thanks again

2004-05-26 Thread partha_bagchi
Uh-oh .. the Dalgaard's Box has been reopened Peter Dalgaard [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 05/25/2004 04:17 PM To: Philippe Grosjean [EMAIL PROTECTED] cc: Chad Shaw [EMAIL PROTECTED], [EMAIL PROTECTED] Subject:Re: [R] thanks again

Re: [R] Turning pass/fail results into a proportion

2004-05-26 Thread Gabor Grothendieck
Matthew Walker m.g.walker at massey.ac.nz writes: : I have a data frame that contains a number of pass/fails for certain : variable sizes. From that, I would like to form another data frame that : contains the proportions of pass/fails per variable. : : So, for example: : df - data.frame(

[R] I/O fortran instructions and dyn.load

2004-05-26 Thread Gilles GUILLOT
I have the following Fortran code subroutine sub(path) character*100 path open(10,file=path) end saved as test.f which I compile with g77 -c test.f then I make the shared libary in R (Version 1.9.0) with system(R CMD SHLIB test.o) so far, everything OK. But dyn.load(test.so) returns the

[R] power analysis

2004-05-26 Thread Ruben Real
Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 7bit Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted sender hosts) X-Virus-Scanned: by amavisd-new X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on hypatia.math.ethz.ch

Re: [R] cor and missing values. Bug?

2004-05-26 Thread Robert W. Baer, Ph.D.
Not to put too fine a point on it, but did you consider checking the NEWS file for the most recent version (1.9.0, http://cran.r-project.org/src/base/NEWS)? o The cor() function did not remove missing values in the non-Pearson case. There is still something a little strange

Re: [R] I/O fortran instructions and dyn.load

2004-05-26 Thread Prof Brian Ripley
Try R CMD SHLIB test.f If SHLIB know it has to deal with Fortran code, it adds the appropriate libraries. On Wed, 26 May 2004, Gilles GUILLOT wrote: I have the following Fortran code subroutine sub(path) character*100 path open(10,file=path) end saved as test.f which I compile

FW: [R] is.weekend() odd behaviour

2004-05-26 Thread Matthew Dowle
Kurt, Uwe suggested I write to you as maintainer of chron ... at the start of is.weekend, in the check on the argument type, maybe just change chron() to as.chron()? This would ensure as.chron.POSIXt gets called on POSIXt arguments, and (I think) fixes the problem. I tested and it seems ok. No

[R] (BRMA) NTI - Internet Controlada

2004-05-26 Thread BRmultiaccess
NTI - Internet Controlada E-mail bloqueado devido a restrições no servidor de e-mail, o e-mail possuia um virus, qualquer duvida contactar o Departamento de Informática. Gaia Importação e Exportação Ltda Para:

Re: [R] I/O fortran instructions and dyn.load

2004-05-26 Thread Gilles GUILLOT
It works with R CMD SHLIB test.f Thanks. Gilles Try R CMD SHLIB test.f If SHLIB know it has to deal with Fortran code, it adds the appropriate libraries. On Wed, 26 May 2004, Gilles GUILLOT wrote: I have the following Fortran code subroutine sub(path) character*100 path

[R] Opening help pages in new tab of existing Mozilla Firebird

2004-05-26 Thread Kevin Wright
Subject pretty much says it all. I currently have options()$browser set to open help pages in Mozilla Firebird, but it starts a new window each time and I would like a new 'tab' in an existing window. Sorry if this is obvious, but I can't find anything. Kevin Wright

Re: [R] Need help on Mclus output

2004-05-26 Thread Bhaskar S. Manda
On Wed, 26 May 2004 03:34:14 +0200, Max Marinucci wrote I have fitted a mixture with 4 normal components on a univariate distribution using the Mclust package. Now, I would like to get a variable with the cluster membership of each class, or in If you've loaded library(mclust) you can do

Re: [R] No direct or inherited method for function update

2004-05-26 Thread Dick Beyer
Thanks very much for your help. Following your suggestion, I found that the offending bioconductor package is affy. I will report that to the bioconductor folks. *** Richard P. Beyer, Ph.D. University of Washington

RE: [R] Using R in C++

2004-05-26 Thread Vadim Ogranovich
Look at Writing R Extensions guide. It covers both R-from-C and C-from-R. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jörgen Wallerman Sent: Wednesday, May 26, 2004 3:09 AM To: '[EMAIL PROTECTED]' Subject: [R] Using R in C++ Hello, Is

Re: [R] Opening help pages in new tab of existing Mozilla Firebird

2004-05-26 Thread Michael Morley
Try this Mozilla Firefox extension, http://www.intraplanar.net/projects/tabprefs/ If you more control over tabs there is this extension as well http://white.sakura.ne.jp/~piro/xul/_tabextensions.html.en -Mike Kevin Wright wrote: Subject pretty much says it all. I currently have options()$browser

Re: [R] Opening help pages in new tab of existing Mozilla Firebird

2004-05-26 Thread Marc Schwartz
On Wed, 2004-05-26 at 10:32, Kevin Wright wrote: Subject pretty much says it all. I currently have options()$browser set to open help pages in Mozilla Firebird, but it starts a new window each time and I would like a new 'tab' in an existing window. Sorry if this is obvious, but I can't

Re: [R] Common principle components

2004-05-26 Thread Martin Maechler
JosePG == J Pedro Granadeiro [EMAIL PROTECTED] on Wed, 26 May 2004 16:01:22 +0100 writes: JosePG I am sorry for not being clear. I meant the methods detailed in: JosePG Flury, B. (1988). Common Principle Components Analysis and Related JosePG Multivariate Models, John Wiley and

Re: [R] Common principle components

2004-05-26 Thread Spencer Graves
Google just returned about 1.8M hits for principle components and 3.2M hits for principal components. Evidently, eigenvectors are more your friend (pal) than a structure for the space by a factor of roughly 2 to 1. spencer graves Martin Maechler wrote: JosePG == J Pedro Granadeiro [EMAIL

Re: [R] thanks again

2004-05-26 Thread Achim Zeileis
Chad: Thanks again for the conference. We had good fun. Thanks, we too ;-) If only I could get some sleep now... I want a useR t-shirt. After lots of drinking and discussion with Fritz Leisch at dinner on saturday, I feel there will be no chance for an official T-shirt. My question is:

Re: FW: [R] is.weekend() odd behaviour

2004-05-26 Thread Gabor Grothendieck
I think the user should be forced to convert explicitly from POSIXct to chron before is.weekend is used. That way the user will know for sure what time zone he is converting with respect to. Having it done silently will just cause subtle errors. IMHO the correct functionality is to give an

Re: [R] thanks again

2004-05-26 Thread Dirk Eddelbuettel
On Wed, May 26, 2004 at 06:30:40PM +0200, Achim Zeileis wrote: As this useR! is over, I guess it's a bit too late for a useR! 2004 shirt. In general, I would agree with you that it would be nice (and not only for fun) to have shirts (and coffee mugs and basecaps and ...) with R logos or

Re: [R] thanks again

2004-05-26 Thread Andrew Robinson
http://www.cafepress.com/ is one such option. I have no personal experience with it. Andrew On Wed, 26 May 2004, Dirk Eddelbuettel wrote: On Wed, May 26, 2004 at 06:30:40PM +0200, Achim Zeileis wrote: As this useR! is over, I guess it's a bit too late for a useR! 2004 shirt. In

Re: [R] Subtracting number of days from a date

2004-05-26 Thread Gabor Grothendieck
Hoeven, Maarten van der Maarten.van.der.Hoeven at knmi.nl writes: : suppose I have a date, say May 15 2004, and I want to now what date it is 23 days before that date. The way to : calculate the new date should (...) take account of leap years :) : : In pseudocode: : : olddate - May 15 2004 :

Re: [R] Subtracting number of days from a date

2004-05-26 Thread Prof Brian Ripley
On Wed, 26 May 2004, Hoeven, Maarten van der wrote: suppose I have a date, say May 15 2004, and I want to now what date it is 23 days before that date. The way to calculate the new date should (...) take account of leap years :) In pseudocode: olddate - May 15 2004 newdate -

Re: [R] thanks again

2004-05-26 Thread Wolski
Hi! The magic cauldron http://www.catb.org/~esr/writings/magic-cauldron/magic-cauldron.html and http://www.catb.org/~esr/writings/magic-cauldron/magic-cauldron-9.html#ss9.4 Sincerely Eryk *** REPLY SEPARATOR *** On 5/26/2004 at 9:40 AM Andrew Robinson wrote:

[R] (no subject)

2004-05-26 Thread Moustafa ElHousinie
I have data on disappearance of virus C from blood of acutely infected humans. Blood samples were taken approximately 2, 3 and 6 months after the expected date of exposure. It is noted that : 1- not all timepoints are available 2- in few cases the person turned negative early (2-3

[R] identify() - image()

2004-05-26 Thread javier garcia - CEBAS
Hi all; Just to ask if you know about any available function in R to identify points in a image plotted in X11. Something like the function identify(), but able to extract (x,y,value) groups from the image. Thanks and best regards, Javier __ [EMAIL

[R] 2 way repeated measures ANOVA using R: syntax and reporting question

2004-05-26 Thread GIDEON WASSERBERG
Dear Friends I have a technical question about conducting 2 way repeated measures ANOVA analysis using R. 1. Data set: repeated measurement of activity over night (2 hr. intervals) repeated (within subject)factor: Hours Between subject: Species, Sex Dependent variables: specimens Here is how

[R] Outlier identification according to Hardin Rocke (1999)

2004-05-26 Thread Kevin Wright
I'm trying to use a paper by Hardin Rocke: http://handel.cipic.ucdavis.edu/~dmrocke/Robdist5.pdf as a guide for a function to identify outliers in multivariate data. Attached below is a function that is my attempt to reproduce their method and also a test to see what fraction of the data are

[R] A data selection problem, suggestions highly appreciated

2004-05-26 Thread Yong Wang
Hi, All I get following question: A data format like following: [Day time x y] Jan1 18:56:24 x1 y1 Jan1 18:56:25 x2 y2 Jan1 18:56:27 x3 y3 Jan1 18:56:28 x4 y4 Jan1 18:56:31 x5 y5 . . what I wanna do is to partion the time interval by unit of 5 seconds.

[R] Aid on two-way ANOVA with contrasts

2004-05-26 Thread joseclaudio.faria
Dears members of R list, It would like that a more experienced statician in R helped me to complete the analysis to follow: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4,

Re: [R] cor and missing values. Bug?

2004-05-26 Thread Peter Dalgaard
Robert W. Baer, Ph.D. [EMAIL PROTECTED] writes: Not to put too fine a point on it, but did you consider checking the NEWS file for the most recent version (1.9.0, http://cran.r-project.org/src/base/NEWS)? o The cor() function did not remove missing values in the

Re: [R] Common principle components

2004-05-26 Thread Peter Dalgaard
Spencer Graves [EMAIL PROTECTED] writes: Google just returned about 1.8M hits for principle components and 3.2M hits for principal components. Evidently, eigenvectors are more your friend (pal) than a structure for the space by a factor of roughly 2 to 1. spencer graves Of course,

Re: [R] cor and missing values. Bug?

2004-05-26 Thread Frank E Harrell Jr
On 27 May 2004 00:20:17 +0200 Peter Dalgaard [EMAIL PROTECTED] wrote: Robert W. Baer, Ph.D. [EMAIL PROTECTED] writes: Not to put too fine a point on it, but did you consider checking the NEWS file for the most recent version (1.9.0, http://cran.r-project.org/src/base/NEWS)?