RE: [R] Yet another beginner question

2003-09-09 Thread Heywood, Giles
The Irregular Time-Series ('its') class could be helpful to you. The following does the basics of what you have asked. require(its) its.format("%m/%d/%y") myits <- its(cbind(reelect,new),dates=as.POSIXct(x=strptime(dateofpoll,format=its.for mat( plot(myits,type="b") You may want to customise

Re: [R] Yet another beginner question

2003-09-09 Thread Uwe Ligges
Peter Flom wrote: Thanks for all the help on my earlier questions. How do you plot a simple time series with unequal intervals? I have the following dateofpoll <- as.ts(c("6/1/02", "7/1/02", "10/1/02", "1/4/03", "1/25/03", "6/7/03", "7/16/03", "8/17/03", "9/4/03")) reelect <- c(51, 47, 49, 5

RE: [R] memory problem in exporting data frame

2003-09-09 Thread Henrik Bengtsson
Hi, I replied to a related question yesterday (Mon Sept 8, 2003) with subject "RE: [R] cannot allocate vector of size...". That was as also about running low of memory, but about *reading* data from file and not writing. However, the problem is likely to be due to the same thing. You pass a large

Re: [R] R video

2003-09-09 Thread Thomas Petzoldt
David Khabie-Zeitoune schrieb: Hi Does anybody know of any R packages under Windows to produce video files from a sequence of R graphs -- e.g. in .wmv or avi format? Hello, some times ago I've testet several methods and found the following to be useful and effective: 1.) Create a numbered sequ

[R] ANOVA

2003-09-09 Thread Britta Lintfert
I have Data like this: Tone OQtil4 OQ0 GOtil4GO0SKrhsJ SK0RCrhsJ RC0 1 HsLG -9.1347877 -2.97 -7.409590 -6.40 9.389357 20.60 10.688854 24.40 2 HsLG -7.9270569 -2.03 -7.861541 -6.90 10.165324 26.40 10.640183 23.10 3 HsLG -7.0394106 -1.26 -7.509566 -6.53 12.033194

Re: [R] R video

2003-09-09 Thread Erich Neuwirth
Quicktime also can convert sequences of images into quicktime and avi movies. It accepts quite a lot of different image formats. Thomas Petzoldt wrote: David Khabie-Zeitoune schrieb: Hi Does anybody know of any R packages under Windows to produce video files from a sequence of R graphs -- e.

Re: [R] ANOVA

2003-09-09 Thread Winfried Theis
Hi Britta! On Tue, 2003-09-09 at 10:57, Britta Lintfert wrote: [snip] > > > now I want to know, if there are significant differents in the variables > OQ, GO, SK, RC depending on Tone. This I can do with an ANOVA , or?? Yes, this is correct. > > But when I start aov in R I get the following me

Re: [R] Matrix resampling (bootstraps)

2003-09-09 Thread Ernesto Jardim
On Tue, 2003-09-09 at 05:11, Hector L. Ayala-del-Rio wrote: > Dear all, > I am trying to generate bootstrap replicate matrixes (rows=samples, > column=species, sampling with replacement) from a matrix dataset, but I do > not know how to do it in R. I have tried boot() and bootstrap(), but th

[R] lattice.xyplot: adding grid lines

2003-09-09 Thread Marc Mamin
Hallo, I'd like to add grid lines to a lattice graph having 2 series of Y data. See these 2 examples: data(iris) [1] xyplot(Sepal.Length + Sepal.Width ~ Petal.Length , data = iris, allow.multiple = TRUE, scales = "same",type="l", ) [2] xyplot(Sepal.Length + Sepal.Width ~ Petal.

[R] lattice plot - portrait / landscape

2003-09-09 Thread Mahbub Latif
Hi, How can I use portrait/landscape option in lattice bwplot? Is there any option in trellis.device where I can define this? Thanks in advance, Mahbub. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] ANOVA

2003-09-09 Thread Spencer Graves
To expand on Winfried's comments, you may wish to try the following before considering multivariate: anova(lm(OQtil4 ~ Tone, data=gesamt)) and similarly for your other response variables. There are ways of putting this kind of thing in a loop, storing only the p values or whatever you want

Re: [R] No joy installing R with shared libs

2003-09-09 Thread Laurent Faisnel
>> Can some kind soul please give me a fool proof recipe for building R >> and RSPython so that it actually works? > I don't have a recipe, but one thought to help debug the process: Try > installing RPy [1]. RPy also provides access to R via Python and uses > the libR.so library. If you can i

[R] problem trying to use a superclass on R-1.7.1

2003-09-09 Thread Laurent Faisnel
Hi dear helpers, I have been using for quite a long time a script where a class extends another one (without trouble). But I now have problems of that kind while trying to run the script : Error in insertMethod(methods, sig, args, def, TRUE) : inserting method corresponding to empty si

Re: [R] lattice plot - portrait / landscape

2003-09-09 Thread Andrew C. Ward
Dear Manbub, I assume that you mean you want the boxes to be vertical rather than horizontal (which is the default). Compare the following two uses of bwplot: data(singer) bwplot(voice.part ~ height, data=singer, xlab="Height (inches)") bwplot(height ~ voice.part, data=singer, ylab="Heig

Re: [R]: RODBC column length>255

2003-09-09 Thread simon_gibbon
Dear Wayne, I had similar problems and as Andrew Ward suggests you can overcome them by recompiling the RODBC code, which is what I did after Professor Ripley's kind advice for Windows. I am now able to read in 64K character strings. However, as of version 1.0.2 the source code sets a 64K lim

[R] How to convert a Rd file into multi html files?

2003-09-09 Thread Yiming Zhou
Hi, I wrote all documents for each objects in my package in one .Rd file. Now I can't convert the .Rd file into multi html files such that a html file corresponds a object in package. The command I used was ' R CMD Rdconv -t html foo.Rd'. It just converted the document of first object in .Rd in

Re: [R] How to convert a Rd file into multi html files?

2003-09-09 Thread Uwe Ligges
Yiming Zhou wrote: Hi, I wrote all documents for each objects in my package in one .Rd file. Now I can't convert the .Rd file into multi html files such that a html file corresponds a object in package. The command I used was ' R CMD Rdconv -t html foo.Rd'. It just converted the document of first

[R] Re: Hierarchical clustering

2003-09-09 Thread ucgamdo
I think are looking for the function 'cutree' from package mva checkout its documentation: > require(mva) > ?cutree pleanty of examples to do what you want. # Hi R lovers! I am using the agnes function of the package

Re: [R] How to convert a Rd file into multi html files?

2003-09-09 Thread Philipp Pagel
Hi! On Tue, Sep 09, 2003 at 04:17:25AM -0700, Yiming Zhou wrote: > The command I used was ' R > CMD Rdconv -t html foo.Rd'. It just converted the document of first > object in .Rd into html format. I'm not exactly an expert in writing R documentation but while playing arround with packag

RE: [R] No joy installing R with shared libs

2003-09-09 Thread Arne.Muller
Hi, I've experienced similar failures with the RSperl installation. So I'd be interested if someone sorts out the library misery ... ;-) Arne > -Original Message- > From: Laurent Faisnel [mailto:[EMAIL PROTECTED] > Sent: 09 September 2003 12:48 > To: [EMAIL PROTECTED] > Subject:

[R] Making R packages

2003-09-09 Thread Gattuso, Jean-Pierre
Hi: I am posting this message for a colleague who has a lot of trouble to build an R package on Windows. He did not find a solution to his problems on the R-help archives and hopes that one of the R gurus will be able to help. He has a directory "CO2" which should contain all the required fi

[R] Building XML package for MacOS X

2003-09-09 Thread Sean Davis
I am working to build the XML package for R on MacOS X. I have installed libxml2-2.5.9 into /usr/local. I set the LIBXML_INCDIR=/usr/local/include/libxml2. I use R INSTALL, I get the following: R INSTALL -c -l /usr/local/R/library XML_0.94-1.tar.gz : : : gcc -bundle -flat_namespace -undef

Re: [R] Making R packages

2003-09-09 Thread Thomas Petzoldt
Gattuso, Jean-Pierre wrote: The following command is run in a windows console: C:\Program Files\R\rw1071\bin\rcmd INSTALL D:\CO2 Then he gets the following error (approximate translation of a French message) : "Perl is not recognized as an internal or external command, an executable fil

Re: [R] Making R packages

2003-09-09 Thread Duncan Murdoch
On Tue, 9 Sep 2003 14:39:53 +0200, you wrote: >Then he gets the following error (approximate >translation of a French message) : > "Perl is not recognized as an internal or >external command, an executable file or a command >file" > >Perl has of course been installed on the PC and >my co

[R] Hardware oddity - linux/windows

2003-09-09 Thread michaell taylor
Hi, I am trying to figure out why a new machine is running my R script so slowly. The script was developed on a Linux , dual 2.8 xeon processor machine with 4GB ram. On this machine, the script runs in about an hour (it creates lots of multilevel simulations). While running, 100% of a single

Re: [R] Making R packages

2003-09-09 Thread Uwe Ligges
Gattuso, Jean-Pierre wrote: Hi: I am posting this message for a colleague who has a lot of trouble to build an R package on Windows. He did not find a solution to his problems on the R-help archives and hopes that one of the R gurus will be able to help. He has a directory "CO2" which should

Re: [R] Making R packages

2003-09-09 Thread Ernesto Jardim
On Tue, 2003-09-09 at 13:39, Gattuso, Jean-Pierre wrote: > Hi: > > I am posting this message for a colleague who has > a lot of trouble to build an R package on > Windows. He did not find a solution to his > problems on the R-help archives and hopes that > one of the R gurus will be able to he

Re: [R] lattice.xyplot: adding grid lines

2003-09-09 Thread Deepayan Sarkar
On Tuesday 09 September 2003 04:49, Marc Mamin wrote: > Hallo, > > I'd like to add grid lines to a lattice graph having 2 series of Y data. > > See these 2 examples: > > > > data(iris) > > [1] > xyplot(Sepal.Length + Sepal.Width ~ Petal.Length , > data = iris, allow.multiple = TRUE, sca

RE: [R] Making R packages

2003-09-09 Thread Wayne Jones
Try adding C:\Perl\bin\; to your path variable. -Original Message- From: Gattuso, Jean-Pierre [mailto:[EMAIL PROTECTED] Sent: 09 September 2003 13:40 To: [EMAIL PROTECTED] Subject: [R] Making R packages Hi: I am posting this message for a colleague who has a lot of trouble to build an

Re: [R] Matrix resampling (bootstraps)

2003-09-09 Thread James MacDonald
Putting aside the issue of whether you should be using boot() or not, you can resample your matrix by doing something like this. > a <- matrix(1:10, nrow=10, ncol=10, byrow=TRUE) > a [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,]12345678910

Re: [R] lattice plot - portrait / landscape

2003-09-09 Thread Deepayan Sarkar
On Tuesday 09 September 2003 05:48, Andrew C. Ward wrote: > Dear Manbub, > > I assume that you mean you want the boxes to be vertical > rather than horizontal (which is the default). Compare the > following two uses of bwplot: >data(singer) >bwplot(voice.part ~ height, data=singer, xlab="He

RE: [R] Making R packages

2003-09-09 Thread Wayne Jones
The document http://www.stat.auckland.ac.nz/~kwan022/rinfo.php gives an excellent step by step description of how to make an R package for Windows. -Original Message- From: Gattuso, Jean-Pierre [mailto:[EMAIL PROTECTED] Sent: 09 September 2003 13:40 To: [EMAIL PROTECTED] Subject: [R] Maki

RE: [R] lattice.xyplot: adding grid lines

2003-09-09 Thread Marc Mamin
Hallo, thank you for your answer, there is still one problem in your solution: Sepal.Length + Sepal.Width are considerd as a single serie of data... Marc Mamin -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 09, 2003 3:03 PM To: Marc Mamin;

[R] Changing Tick Marks for Date Plots

2003-09-09 Thread Sean Carmody
I have been experimenting with various approaches to plotting (irregular) time-series with reasonable success. One thing I have been able to do is change the number of tick marks on the axis. Consider the following simple example: > x <- as.data.frame(matrix(ncol=2,nrow=500)) > names(x) <- c("date

Re: [R] Building XML package for MacOS X

2003-09-09 Thread Prof Brian Ripley
The simplest way is probably to edit the libxml2 library and remove the module containing xmlParserError. That's what the Windows maintainer has to do. On Tue, 9 Sep 2003, Sean Davis wrote: > I am working to build the XML package for R on MacOS X. I have installed > libxml2-2.5.9 into /usr/loc

Re: [R] lattice.xyplot: adding grid lines

2003-09-09 Thread Deepayan Sarkar
Sorry, that's what you get when you don't run your code before sending it. Should be panel.superpose instead of panel.xyplot, since you want a grouped display. e.g., xyplot(Sepal.Length + Sepal.Width ~ Petal.Length , data = iris, allow.multiple = TRUE, scales = "same",type="l", pa

Re: [R] Hardware oddity - linux/windows

2003-09-09 Thread Prof Brian Ripley
Do you have an anti-virus software running? That's about what I see when Sophos is switched on. On 9 Sep 2003, michaell taylor wrote: > > Hi, > > I am trying to figure out why a new machine is running my R script so > slowly. > > The script was developed on a Linux , dual 2.8 xeon processo

Re: [R] How to convert a Rd file into multi html files?

2003-09-09 Thread Paul Gilbert
The R packaging programs expect one .Rd file for each function, all in the man directory. Installing the package converts each .Rd file into an html file as well as other documentation formats. There was once some discussion of allowing multiple functions to be documented in a single file, but

[R] tree mortality risk model using glm()

2003-09-09 Thread Tord Snall
Dear all, I've used glm(family=binomial(link="logit")) several times, but now I think that a log link is more appropriate. I want to fit a model for probability of tree fall (TF)), with tree diameter (dbh) and soil moisure (soil) as predictors. A large number of trees have been checked every sec

[R] ploting 2 variables in a lattice xyplot

2003-09-09 Thread Ernesto Jardim
Hi, I want to plot 2 variables (confidence intervals) in a single xyplot. I have a dataframe with variables Yup, Ylo, X and Z and I want to have a xyplot with both variables ploted. I'm trying with xyplot(Yup~X|Z, panel=function(x){...}) but this way I'm not able to pass variable Ylo into the f

Re: [R] How do I coerce numeric factor columns of data frame to vector?

2003-09-09 Thread Michael A. Miller
> "Thomas" == Thomas W Blackwell <[EMAIL PROTECTED]> writes: > Michael - Because these columns are factors to begin with, > using as.numeric() alone will have unexpected results. See > the section "Warning:" in help("factor"). Ah, thanks for pointing that out. When I've used as.

RE: [R] Changing Tick Marks for Date Plots

2003-09-09 Thread Heywood, Giles
The way to control number of tickmarks in plot is via par(lab) par(lab=c(5,5,7)) #the default plot(rnorm(20),rnorm(20)) par(lab=3*c(5,5,7)) plot(rnorm(20),rnorm(20)) However this does not work for axis.POSIXct, which the function called by the plot method for 'its'. I am not sure why this is. T

[R] charge a vector with variables and to use as variable in a checkbutton?

2003-09-09 Thread solares
hello, how i cant to charge in form dynamic a checkbutton, try to do it with a vector be charged automaticamente but not works, for example library(tcltk) tt<-tktoplevel() f<-tkframe(tt) tkpack(f) i<-2 if (i==1) {b1<-tkcheckbutton (f,text="b1",variable="b1",relief="raised");tkpack(b1);print(tclva

Re: [R] How do I coerce numeric factor columns of data frame to vector?

2003-09-09 Thread Na Li
On 9 Sep 2003, Michael A. Miller outgrape: > > > > > > "Thomas" == Thomas W Blackwell <[EMAIL PROTECTED]> writes: > > > Michael - Because these columns are factors to begin with, > > using as.numeric() alone will have unexpected results. See > > the section "Warning:" in help("factor"). >

Re: [R] tree mortality risk model using glm()

2003-09-09 Thread Thomas Lumley
On Tue, 9 Sep 2003, Tord Snall wrote: > Dear all, > > I've used glm(family=binomial(link="logit")) several times, but now I think > that a log link is more appropriate. > > I want to fit a model for probability of tree fall (TF)), with tree > diameter (dbh) and soil moisure (soil) as predictors. A

Re: [R] ploting 2 variables in a lattice xyplot

2003-09-09 Thread Douglas Bates
Ernesto Jardim <[EMAIL PROTECTED]> writes: > Hi, > > I want to plot 2 variables (confidence intervals) in a single xyplot. I > have a dataframe with variables Yup, Ylo, X and Z and I want to have a > xyplot with both variables ploted. I'm trying with > > xyplot(Yup~X|Z, panel=function(x){...})

Re: [R] How do I coerce numeric factor columns of data frame to vector?

2003-09-09 Thread Martin Maechler
> "NaLi" == Na Li <[EMAIL PROTECTED]> > on Tue, 09 Sep 2003 10:15:22 -0500 writes: NaLi> On 9 Sep 2003, Michael A. Miller outgrape: >> > > > > > "Thomas" == Thomas W Blackwell <[EMAIL PROTECTED]> writes: >> >> > Michael - Because these columns are factors to begin wi

Re: [R] Stifling REprintf() output

2003-09-09 Thread Thomas Lumley
On Mon, 8 Sep 2003, Jeff Gentry wrote: > > In some code that I have written, use of url() is generating the output > line: > "cannot open: HTTP status was `404 Not Found`" > My fault. This used to use error(), but as you will recall, the result was that the connection wasn't freed and getBioC()

Re: [R] charge a vector with variables and to use as variable in a checkbutton?

2003-09-09 Thread Dirk Eddelbuettel
On Tue, Sep 09, 2003 at 12:10:06PM -0300, [EMAIL PROTECTED] wrote: > hello, how i cant to charge in form dynamic a checkbutton, try to do it > with a vector be charged automaticamente but not > works, for example [...] > while (i<=j){ > b[i]<-tkcheckbutton(f,text="b[i]",variable="b[i]",relief="r

Re: [R] Building XML package for MacOS X

2003-09-09 Thread Sean Davis
I know this is not entirely on-topic, but assuming that several users have had to do this, could someone give me some specifics? If it is off-topic, feel free to reply directly rather than to the group. Thanks > The simplest way is probably to edit the libxml2 library and remove > the module con

[R] Cross-platform compatibility with R?

2003-09-09 Thread Joseph Ortiz
Dear All, I'm a newbie to R having just switched from S-Plus. I was able to port my files from S-Plus for Linux to R for OSX using the S-Plus "data.dump", and the R "data.restore" commands. I've decided that I like R enough to install it on my Linux box, but I can't figure out an easy way to g

Re: [R] How do I coerce numeric factor columns of data frame to vector?

2003-09-09 Thread Peter Dalgaard BSA
Martin Maechler <[EMAIL PROTECTED]> writes: > >> as.numeric (levels (y)[as.numeric (y)]) > NaLi> [1] 1 3 1 > > NaLi> which is a bit awkward. > as.numeric(as.character(y)) ! > > {in some cases you might consider using as.integer() instead of as.numeric()} Better to convert first and

[R] Making R packages (Unix)

2003-09-09 Thread Gattuso, Jean-Pierre
Hi: I have have taken over from a colleague who prepared an R package and failed to build it on Windows. I am doing this with unix as I am a mac user. Below is the output I get when I use the build command: [gattuso:unix/R/CO2.Rcheck] gattuso% R CMD build CO2 * checking for file 'CO2/DESCRIPT

Re: [R] Cross-platform compatibility with R?

2003-09-09 Thread Uwe Ligges
Joseph Ortiz wrote: Dear All, I'm a newbie to R having just switched from S-Plus. I was able to port my files from S-Plus for Linux to R for OSX using the S-Plus "data.dump", and the R "data.restore" commands. I've decided that I like R enough to install it on my Linux box, but I can't figure

Re: [R] Cross-platform compatibility with R?

2003-09-09 Thread Peter Dalgaard BSA
Joseph Ortiz <[EMAIL PROTECTED]> writes: > I tried to use tar to archive my directories and binary .Rdata files > and simply extract them to a new machine. This works fine going from > one OSX machine to another, but does not work going from OSX to Linux, > presumably due to binary file or directo

Re: [R] Cross-platform compatibility with R?

2003-09-09 Thread Douglas Bates
Joseph Ortiz <[EMAIL PROTECTED]> writes: > Dear All, > > I'm a newbie to R having just switched from S-Plus. > > I was able to port my files from S-Plus for Linux to R for OSX using > the S-Plus "data.dump", and the R "data.restore" commands. I've > decided that I like R enough to install it on

Re: [R] ploting 2 variables in a lattice xyplot

2003-09-09 Thread Ernesto Jardim
On Tue, 2003-09-09 at 16:35, Douglas Bates wrote: > Ernesto Jardim <[EMAIL PROTECTED]> writes: > > > Hi, > > > > I want to plot 2 variables (confidence intervals) in a single xyplot. I > > have a dataframe with variables Yup, Ylo, X and Z and I want to have a > > xyplot with both variables ploted

RE: [R] memory problem in exporting data frame

2003-09-09 Thread array chip
Hi all, Thanks for all the suggestions. I was able to get the data frame out by first deleting some other large objects in the directory, and then changing the data frame into matrix by as.matrix(), splitting the matrix into 4 blocks and finally using write.table() to write the matrix into 4 files

Re: [R] memory problem in exporting data frame

2003-09-09 Thread Peter Dalgaard BSA
array chip <[EMAIL PROTECTED]> writes: > Hi all, > > Thanks for all the suggestions. I was able to get the > data frame out by first deleting some other large > objects in the directory, and then changing the data > frame into matrix by as.matrix(), splitting the matrix > into 4 blocks and finall

Re: [R] How to convert a Rd file into multi html files?

2003-09-09 Thread Paul Gilbert
Paul Gilbert wrote: The R packaging programs expect one .Rd file for each function, ... Sorry, this was not exactly correct. It expects one .Rd file for each documentation item (man page) indicated by \name{}, but that can document more than one function. Paul Gilbert _

Re: [R] How do I coerce numeric factor columns of data frame to vector?

2003-09-09 Thread Michael A. Miller
> "Na" == Na Li <[EMAIL PROTECTED]> writes: > x <- c(1, 3, 1) > y <- factor (x) > y > [1] 1 3 1 > Levels: 1 3 > as.numeric (y) > [1] 1 2 1 Ah, now I get it - thanks! Mike __ [EMAIL PROTECTED] mailing list https://www

Re: [R] ploting 2 variables in a lattice xyplot

2003-09-09 Thread Deepayan Sarkar
On Tuesday 09 September 2003 10:35 am, Douglas Bates wrote: > Ernesto Jardim <[EMAIL PROTECTED]> writes: > > Hi, > > > > I want to plot 2 variables (confidence intervals) in a single xyplot. I > > have a dataframe with variables Yup, Ylo, X and Z and I want to have a > > xyplot with both variables

RE: [R] memory problem in exporting data frame

2003-09-09 Thread Henrik Bengtsson
> -Original Message- > From: array chip [mailto:[EMAIL PROTECTED] > Sent: den 9 september 2003 19:04 > To: Henrik Bengtsson; 'Thomas W Blackwell'; Patrick Burns > Cc: [EMAIL PROTECTED] > Subject: RE: [R] memory problem in exporting data frame > > > Hi all, > > Thanks for all the suggest

[R] Computing a CDF or many quantiles

2003-09-09 Thread Kevin S. Van Horn
Given f, a pdf over a finite interval, is there any existing R function that can efficiently tabulate the cumulative distribution function for f, or produce all N+1 quantiles of the form i/N? "Efficiently" here means better than doing repeated integrations for each point. _

Re: [R] Computing a CDF or many quantiles

2003-09-09 Thread Jerome Asselin
Perhaps approx() or approxfun() can help you create an efficient CDF. HTH, Jerome On September 9, 2003 01:38 pm, Kevin S. Van Horn wrote: > Content-Length: 417 > Status: R > X-Status: N > > Given f, a pdf over a finite interval, is there any existing R function > that can efficiently tabulate th

[R] Newbie question about plotting density objects

2003-09-09 Thread cmprobst
Dear List, I have an array of 6400 x 56 elements. I want to calculate the density function for each column and plot all 56 density functions in one plot. I have tried several procedures, but they all failed. What can I do? TIA Christian [[alternative HTML version deleted]] __

RE: [R] Computing a CDF or many quantiles

2003-09-09 Thread Ted Harding
On 09-Sep-03 Kevin S. Van Horn wrote: > Given f, a pdf over a finite interval, is there any existing R function > that > can efficiently tabulate the cumulative distribution function for f, or > produce all N+1 quantiles of the form i/N? "Efficiently" here means > better than doing repeated integ

[R] recording and taking mean of a set of matrices

2003-09-09 Thread Ross Boylan
I'm looking for a good form in which to store matrix results of a simulation. I am doing a simulation study. Each simulation generates some data and then analyzes it. I want to record the results of many simulations and analyze them. Say r has the results of one simulation, and I care about r$c

Re: [R] recording and taking mean of a set of matrices

2003-09-09 Thread Spencer Graves
In simulation like you describe, it is best to avoid using rbind in a loop, as that has more overhead than creating objects of the size required to store the results before you start the loop. Also, if all your results are numbers, it may be better to avoid data.frames as they require more

RE: [R] recording and taking mean of a set of matrices

2003-09-09 Thread Adaikalavan RAMASAMY
" mean(list(m1, m2)) " will not work. mylist <- list(m1, m2) sapply( mylist, FUN=mean ) gives will give you the mean of m1 and m2 sapply( mylist, FUN= function(x) apply(x, 2, mean) ) will give you the column means of m1 and m2 in a matrix format. Double check the resulting dimension. Here are

Re: [R] recording and taking mean of a set of matrices

2003-09-09 Thread Gabor Grothendieck
Just one small comment. R stores arrays in reverse odometer order, i.e. the leftmost subscript varies fastest. Thus, you might want to use the last subscript to represent the simulation number rather than the first. That way, the entire vector or matrix is stored together (i.e. contiguously in

Re: [R] Newbie question about plotting density objects

2003-09-09 Thread Jason Turner
cmprobst wrote: Dear List, I have an array of 6400 x 56 elements. I want to calculate the density function for each column and plot all 56 density functions in one plot. I have tried several procedures, but they all failed. What can I do? 56 lines could even confuse a spider, so you might want t

Re: [R] recording and taking mean of a set of matrices

2003-09-09 Thread Spencer Graves
If all you want are the means of the statistics, then you may not need to store all the individual simulation results. What about the following: N <- 9 # number of simulations k <- 3 # number of coefficients CoefSum <- rep(0, k) names(CoefSum) <- letters[1:k] VarSum <- array(0, dim=c(k, k))

RE: [R] recording and taking mean of a set of matrices

2003-09-09 Thread Richard A. O'Keefe
If there are a _lot_ of results, and they are not needed until after the simulations have all run, it's always possible to write them out to a file or files and then use scan() or read.table() or something like that to read them back in. __ [EMAIL PROTEC

Re: [R] Making R packages (Unix)

2003-09-09 Thread Michael Mader
"Gattuso, Jean-Pierre" wrote: > > Hi: > > I have have taken over from a colleague who > prepared an R package and failed to build it on > Windows. I am doing this with unix as I am a mac > user. Below is the output I get when I use the > build command: > > [gattuso:unix/R/CO2.Rcheck] gattuso% R