Re: [R] How to do aggregate operations with non-scalar functions
On Apr 7, 2005 1:18 AM, Itay Furman [EMAIL PROTECTED] wrote: On Tue, 5 Apr 2005, Gabor Grothendieck wrote: On Apr 5, 2005 6:59 PM, Itay Furman [EMAIL PROTECTED] wrote: Hi, I have a data set, the structure of which is something like this: a - rep(c(a, b), c(6,6)) x - rep(c(x, y, z), c(4,4,4)) df - data.frame(a=a, x=x, r=rnorm(12)) The true data set has 1 million rows. The factors a and x have about 70 levels each; combined together they subset 'df' into ~900 data frames. For each such subset I'd like to compute various statistics including quantiles, but I can't find an efficient way of [snip] I would like to end up with a data frame like this: a x 0%25% 1 a x -0.7727268 0.1693188 2 a y -0.3410671 0.1566322 3 b y -0.2914710 -0.2677410 4 b z -0.8502875 -0.6505710 [snip] One can use do.call(rbind, by(df, list(a = a, x = x), f)) where f is the appropriate function. In this case f can be described in terms of df.quantile which is like quantile except it returns a one row data frame: df.quantile - function(x,p) as.data.frame(t(data.matrix(quantile(x, p f - function(df, p = c(0.25, 0.5)) cbind(df[1,1:2], df.quantile(df[,r], p)) Thanks! Just what I wanted. A minor point is that for some reason the row numbers in the final data frame are not sequential (see below -- this is not a consequence of my changes). These are the original row numbers of the first row of each combo of a and x. If z is the result of do.call you can always do this: row.names(z) - 1:nrow(z) if this its needed. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] density estimation with weighted sample
Dear all I would like to perform density estimation with a weighted sample (output of an Importance Sampling procedure) in R. Could anybody give me an advice on what function to use (in which package)? Thanks a lot, Lorenzo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] density estimation with weighted sample
On Thu, 7 Apr 2005, Tomassini, Lorenzo wrote: I would like to perform density estimation with a weighted sample (output of an Importance Sampling procedure) in R. Could anybody give me an advice on what function to use (in which package)? This could mean 1) You have a sample with weights w, so `w=4' means `I have 4 of those'. 2) You have a sample from a density proportional to w(x)f(x) and want to estimate f. Your title suggests the first, your comment the second. If it is the second, use any package (even density() in R) to estimate the density g of the sampled distribution, for ghat/w and rescale to unit area. If you know a lot about w (e.g. in stereology) there are specialized methods which are better. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] build rpvm under cygwin
I tried ot build rpvm in my own makefile. But runs into some linker errors like e.g. undefined reference to `_R_alloc' My enviornment looks like this: CYGWIN pvm 3.4 compiled under cygwin myself R installed from the rw2001.exe setup file. I guess that the R under rw2001.exe was build with some other compiler? I then tried to compile R myself under CYGWIN but runs into the following problem: building from src typing make gcc -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CON FIG_H -D__NO_MATH_INLINES -g -O2 -c dynload.c -o dynload.o In file included from dynload.c:35: ../../src/include/Defn.h:60:22: psignal.h: No such file or directory I found the psignal.h header file under gnuwin32/fixed/h/psignal.h how do incoperate this udner cygwin? I also tried to type build under src/gnuwin32 but get another error: $ make make: ./Rpwd.exe: Command not found make[1]: ./Rpwd.exe: Command not found make --no-print-directory -C front-ends Rpwd make -C ../../include -f Makefile.win version make[3]: Nothing to be done for `version'. make Rpwd.exe gcc -O2 -Wall -pedantic -I../../include -c rpwd.c -o rpwd.o rpwd.c:22:20: direct.h: No such file or directory rpwd.c: In function `main': rpwd.c:38: warning: implicit declaration of function `chdir' rpwd.c:41: warning: implicit declaration of function `getcwd' make[3]: *** [rpwd.o] Error 1 make[2]: *** [Rpwd] Error 2 make[1]: *** [front-ends/Rpwd.exe] Error 2 make: *** [all] Error 2 Ideas would be appreciated. Regards Lars Schouw --- Lars Schouw [EMAIL PROTECTED] wrote: Dear Professor Ripley The good news is that I fot PVM up and running one two Windows nodes now. I had to connect them with each other manually . for now not using rsh or ssh. Now building RPVM for Windows might not be so easy as it sounds. Did anyone try this out before successfully? Also the SNOW package but that did not look so bad. Regards Lars --- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Thu, 24 Mar 2005, A.J. Rossini wrote: Looks like you are trying to install source tarball on Windows without the relevant toolset (compiler, etc)? To save further hassle, rpvm is not going to build on Windows unless you have PVM installed and working on Windows. If that is the case, this looks like the use of the wrong make, with the wrong shell (that message is coming from a Windows shell, not sh.exe). Do see the warnings in README.packages about the MinGW make. On Thu, 24 Mar 2005 00:11:34 -0800 (PST), Lars Schouw [EMAIL PROTECTED] wrote: I am trying to install the rpvm package doing this: C:\R\rw2000\binrcmd install rpvm_0.6-2.tar.gz '.' is not recognized as an internal or external command, operable program or batch file. '.' is not recognized as an internal or external command, operable program or batch file. make: *** /rpvm: No such file or directory. Stop. make: *** [pkg-rpvm] Error 2 *** Installation of rpvm failed *** Removing 'C:/R/rw2000/library/rpvm' What does this error message tell me? -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ Show us what our next emoticon should look like. Join the fun. http://www.advision.webevents.yahoo.com/emoticontest __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] (no subject)
__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] how to print error message in batch mode
Dear list, I am using R in batch mode: $ R -q --no-save prova output the input file prova has these commands: data(USArrests) x-USArrests hist(x) of course, the command hist(x) produces an error. The error message is: Error in hist.default(x) : `x' must be numeric. Is there the possibility to save this error massage in the output file? Thanks in advance, Angelo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
On Thu, 07 Apr 2005 12:45:16 +0200 Elio Mineo wrote: Dear list, I am using R in batch mode: $ R -q --no-save prova output the input file prova has these commands: data(USArrests) x-USArrests hist(x) of course, the command hist(x) produces an error. The error message is: Error in hist.default(x) : `x' must be numeric. Is there the possibility to save this error massage in the output file? You could do something like $ R -q --no-save prova prova.out 2 prova.err Best, Z Thanks in advance, Angelo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] sweave bwplot error
Hi I use sweave and have a problem with the following figure, but not with other figures: tt - data.frame(c(a, b, c), c(1.2, 3, 4.5)) names(tt) - c(x1, x2) bwplot(x2 ~x1, data = tt) ok now in sweave: \begin{figure}[H] \begin{center} echo=FALSE, fig=TRUE, height=5, width=10= lset(col.whitebg()) bwplot(x2 ~x1, data = tt) @ \caption{xxx} \end{center} \end{figure} PROBLEM: the pdf of the figure is not correctly created (neither the esp) and the error I get from sweave is: pdf inclusion: required page does not exist 0 thanks for help christoph __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
That's fine! Thanks a lot. Angelo Achim Zeileis wrote: On Thu, 07 Apr 2005 12:45:16 +0200 Elio Mineo wrote: Dear list, I am using R in batch mode: $ R -q --no-save prova output the input file prova has these commands: data(USArrests) x-USArrests hist(x) of course, the command hist(x) produces an error. The error message is: Error in hist.default(x) : `x' must be numeric. Is there the possibility to save this error massage in the output file? You could do something like $ R -q --no-save prova prova.out 2 prova.err Best, Z Thanks in advance, Angelo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Order of boxes in boxplot()
Hi Sorry for such an inane question - how do I control the order in which the boxes are plotted using boxplot() when I pass it a formula and a data.frame? It seems that the groups are plotted in alphabetical order... I want to change this Many thanks Mick __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] sweave bwplot error
On Thu, 07 Apr 2005 13:24:42 +0200 Christoph Lehmann wrote: Hi I use sweave and have a problem with the following figure, but not with other figures: tt - data.frame(c(a, b, c), c(1.2, 3, 4.5)) names(tt) - c(x1, x2) bwplot(x2 ~x1, data = tt) ok now in sweave: \begin{figure}[H] \begin{center} echo=FALSE, fig=TRUE, height=5, width=10= lset(col.whitebg()) bwplot(x2 ~x1, data = tt) @ \caption{xxx} \end{center} \end{figure} PROBLEM: the pdf of the figure is not correctly created (neither the esp) and the error I get from sweave is: pdf inclusion: required page does not exist 0 This is covered by FAQ 7.22. (you need to print() the plot) Z thanks for help christoph __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] sweave bwplot error
On Thu, 07 Apr 2005 13:24:42 +0200, Christoph Lehmann (CL) wrote: Hi I use sweave and have a problem with the following figure, but not with other figures: tt - data.frame(c(a, b, c), c(1.2, 3, 4.5)) names(tt) - c(x1, x2) bwplot(x2 ~x1, data = tt) ok now in sweave: \begin{figure}[H] \begin{center} echo=FALSE, fig=TRUE, height=5, width=10= lset(col.whitebg()) bwplot(x2 ~x1, data = tt) @ \caption{xxx} \end{center} \end{figure} PROBLEM: the pdf of the figure is not correctly created (neither the esp) and the error I get from sweave is: pdf inclusion: required page does not exist 0 From the Sweave FAQ: A.6 Why do R lattice graphics not work? The commands in package lattice have different behavior than the standard plot commands in the base package: lattice commands return an object of class trellis, the actual plotting is performed by the print method for the class. Encapsulating calls to lattice functions in print() statements should do the trick, e.g.: fig=TRUE= library(lattice) print(bwplot(1:10)) @ should work. -- --- Friedrich Leisch Institut für Statistik Tel: (+43 1) 58801 10715 Technische Universität WienFax: (+43 1) 58801 10798 Wiedner Hauptstraße 8-10/1071 A-1040 Wien, Austria http://www.ci.tuwien.ac.at/~leisch __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] sweave bwplot error
Christoph Lehmann wrote: Hi I use sweave and have a problem with the following figure, but not with other figures: tt - data.frame(c(a, b, c), c(1.2, 3, 4.5)) names(tt) - c(x1, x2) bwplot(x2 ~x1, data = tt) ok now in sweave: \begin{figure}[H] \begin{center} echo=FALSE, fig=TRUE, height=5, width=10= lset(col.whitebg()) bwplot(x2 ~x1, data = tt) @ \caption{xxx} \end{center} \end{figure} PROBLEM: the pdf of the figure is not correctly created (neither the esp) and the error I get from sweave is: pdf inclusion: required page does not exist 0 thanks for help christoph You need wrap print() round lattice functions to get them to do anything in situations like this. See the Sweave FAQ for this FAQ: http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html#x1-8000A.6 G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Gavin Simpson [T] +44 (0)20 7679 5522 ENSIS Research Fellow [F] +44 (0)20 7679 7565 ENSIS Ltd. ECRC [E] gavin.simpsonATNOSPAMucl.ac.uk UCL Department of Geography [W] http://www.ucl.ac.uk/~ucfagls/cv/ 26 Bedford Way[W] http://www.ucl.ac.uk/~ucfagls/ London. WC1H 0AP. %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] sweave bwplot error
On Thu, Apr 07, 2005 at 01:24:42PM +0200, Christoph Lehmann wrote: \begin{figure}[H] \begin{center} echo=FALSE, fig=TRUE, height=5, width=10= lset(col.whitebg()) bwplot(x2 ~x1, data = tt) @ \caption{xxx} \end{center} \end{figure} PROBLEM: the pdf of the figure is not correctly created (neither the esp) and the error I get from sweave is: pdf inclusion: required page does not exist 0 You need a print() statement around bwplot() -- see the FAQ. Dirk -- Better to have an approximate answer to the right question than a precise answer to the wrong question. -- John Tukey as quoted by John Chambers __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] build rpvm under cygwin
Read the FAQs, etc, about building R on Windows. Summary: stay away from Cygwin when it comes to R. On Apr 7, 2005 11:34 AM, Lars Schouw [EMAIL PROTECTED] wrote: I tried ot build rpvm in my own makefile. But runs into some linker errors like e.g. undefined reference to `_R_alloc' My enviornment looks like this: CYGWIN pvm 3.4 compiled under cygwin myself R installed from the rw2001.exe setup file. I guess that the R under rw2001.exe was build with some other compiler? I then tried to compile R myself under CYGWIN but runs into the following problem: building from src typing make gcc -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CON FIG_H -D__NO_MATH_INLINES -g -O2 -c dynload.c -o dynload.o In file included from dynload.c:35: ../../src/include/Defn.h:60:22: psignal.h: No such file or directory I found the psignal.h header file under gnuwin32/fixed/h/psignal.h how do incoperate this udner cygwin? I also tried to type build under src/gnuwin32 but get another error: $ make make: ./Rpwd.exe: Command not found make[1]: ./Rpwd.exe: Command not found make --no-print-directory -C front-ends Rpwd make -C ../../include -f Makefile.win version make[3]: Nothing to be done for `version'. make Rpwd.exe gcc -O2 -Wall -pedantic -I../../include -c rpwd.c -o rpwd.o rpwd.c:22:20: direct.h: No such file or directory rpwd.c: In function `main': rpwd.c:38: warning: implicit declaration of function `chdir' rpwd.c:41: warning: implicit declaration of function `getcwd' make[3]: *** [rpwd.o] Error 1 make[2]: *** [Rpwd] Error 2 make[1]: *** [front-ends/Rpwd.exe] Error 2 make: *** [all] Error 2 Ideas would be appreciated. Regards Lars Schouw --- Lars Schouw [EMAIL PROTECTED] wrote: Dear Professor Ripley The good news is that I fot PVM up and running one two Windows nodes now. I had to connect them with each other manually . for now not using rsh or ssh. Now building RPVM for Windows might not be so easy as it sounds. Did anyone try this out before successfully? Also the SNOW package but that did not look so bad. Regards Lars --- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Thu, 24 Mar 2005, A.J. Rossini wrote: Looks like you are trying to install source tarball on Windows without the relevant toolset (compiler, etc)? To save further hassle, rpvm is not going to build on Windows unless you have PVM installed and working on Windows. If that is the case, this looks like the use of the wrong make, with the wrong shell (that message is coming from a Windows shell, not sh.exe). Do see the warnings in README.packages about the MinGW make. On Thu, 24 Mar 2005 00:11:34 -0800 (PST), Lars Schouw [EMAIL PROTECTED] wrote: I am trying to install the rpvm package doing this: C:\R\rw2000\binrcmd install rpvm_0.6-2.tar.gz '.' is not recognized as an internal or external command, operable program or batch file. '.' is not recognized as an internal or external command, operable program or batch file. make: *** /rpvm: No such file or directory. Stop. make: *** [pkg-rpvm] Error 2 *** Installation of rpvm failed *** Removing 'C:/R/rw2000/library/rpvm' What does this error message tell me? -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ Show us what our next emoticon should look like. Join the fun. http://www.advision.webevents.yahoo.com/emoticontest __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- best, -tony Commit early,commit often, and commit in a repository from which we can easily roll-back your mistakes (AJR, 4Jan05). A.J. Rossini [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
It is probably easier to use BATCH as in R CMD BATCH prova output See ?BATCH. That does what Elio actually asked for. On Thu, 7 Apr 2005, Achim Zeileis wrote: On Thu, 07 Apr 2005 12:45:16 +0200 Elio Mineo wrote: Dear list, I am using R in batch mode: $ R -q --no-save prova output the input file prova has these commands: data(USArrests) x-USArrests hist(x) of course, the command hist(x) produces an error. The error message is: Error in hist.default(x) : `x' must be numeric. Is there the possibility to save this error massage in the output file? You could do something like $ R -q --no-save prova prova.out 2 prova.err Best, Z Thanks in advance, Angelo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Order of boxes in boxplot()
It is the order of the levels of the grouping factor. If `grp' is not a factor, it will be made into one (with levels in alphabetical order). On Thu, 7 Apr 2005, michael watson (IAH-C) wrote: Sorry for such an inane question - how do I control the order in which the boxes are plotted using boxplot() when I pass it a formula and a data.frame? It seems that the groups are plotted in alphabetical order... I want to change this -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
This solution is fine, too. The Achim's solution is what I have asked with this slight modification: $ R -q --no-save prova prova.out 2 prova.out Again, thanks to Achim and Prof. Ripley. Best, Elio Prof Brian Ripley wrote: It is probably easier to use BATCH as in R CMD BATCH prova output See ?BATCH. That does what Elio actually asked for. On Thu, 7 Apr 2005, Achim Zeileis wrote: On Thu, 07 Apr 2005 12:45:16 +0200 Elio Mineo wrote: Dear list, I am using R in batch mode: $ R -q --no-save prova output the input file prova has these commands: data(USArrests) x-USArrests hist(x) of course, the command hist(x) produces an error. The error message is: Error in hist.default(x) : `x' must be numeric. Is there the possibility to save this error massage in the output file? You could do something like $ R -q --no-save prova prova.out 2 prova.err Best, Z Thanks in advance, Angelo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] analyse des correspondances multiples
bonjour, Je voudrais faire une analyse des correspondances multiples avec R. avec les représentation graphiques correspondantes avec R. je ne sais pas comment procéder .. en vour remerciant par avance Faouzi __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] using command line flags with TINN-R
Philippe, I have no idea what R call-tip server means, but I will invoke it and see what happens. I will also read the FAQ more. Thanks for your help. Thanks, Roger On Apr 7, 2005 1:53 AM, Philippe Grosjean [EMAIL PROTECTED] wrote: roger bos wrote: This is a TINN-R editor question rather than an R question, but can anyone tell me how to use command line flags with TINN-R. There is a space to fill in the path to Rgui, and I have C:\Program Files\R\rw2001pat\bin\Rgui.exe. If I try to add a command line flag after that, such as --no-save or --max-mem-size then TINN-R will not open the application. No that does not work, but you can consider working in the other way: starting Tinn-R while you start R. Then you have all the flexibility to define whatever command line argument you want for R. There are many ways to do so, but I personally use the following one: 1) I define: options(IDE = c:/program files/tinn-R/bin/tinn-R.exe) (of course, the path should reflect the place you actually installed Tinn-R!) and then, I start the svGUI package (from the SciViews bundle available on CRAN). library(svGUI) Tinn-R is started (if not already running), and also, the R call-tip server (live calculation of call-tips for the syntax of R functions) is activated behind the scene. If you are happy with this, and would like to start Tinn-R and activate the R call-tip server automatically everytime you start R, just add those two lines of code in your 'Rprofile' file (the general 'Rprofile' is in /etc subdirectory of the R directory). Once it is done, do not worry about starting Tinn-R, or R from within Tinn-R, just start R with all the command line options you like, and you get Tinn-R started automatically (if it is not running yet)! Note that this tip is in FAQ 3.8, in the new version of Tinn-R FAQ to be released soon, together with the latest stable Tinn-R 1.15.1.7 next week or so ;-) Best, Philippe ..°})) ) ) ) ) ) ( ( ( ( (Prof. Philippe Grosjean ) ) ) ) ) ( ( ( ( (Numerical Ecology of Aquatic Systems ) ) ) ) ) Mons-Hainaut University, Pentagone (3D08) ( ( ( ( (Academie Universitaire Wallonie-Bruxelles ) ) ) ) ) 8, av du Champ de Mars, 7000 Mons, Belgium ( ( ( ( ( ) ) ) ) ) phone: + 32.65.37.34.97, fax: + 32.65.37.30.54 ( ( ( ( (email: [EMAIL PROTECTED] ) ) ) ) ) ( ( ( ( (web: http://www.umh.ac.be/~econum ) ) ) ) ) http://www.sciviews.org ( ( ( ( ( .. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
On Thu, Apr 07, 2005 at 01:58:55PM +0200, Elio Mineo wrote: This solution is fine, too. The Achim's solution is what I have asked with this slight modification: $ R -q --no-save prova prova.out 2 prova.out The canonical and perhaps more correct way to redirect stderr into stdout is 21, as in R -q --no-save prova prova.out 21 Best regards, Jan -- +- Jan T. Kim ---+ |*NEW*email: [EMAIL PROTECTED] | |*NEW*WWW: http://www.cmp.uea.ac.uk/people/jtk | *-= hierarchical systems are for files, not for humans =-* __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] analyse des correspondances multiples
library(MASS) ?mca On Thu, 7 Apr 2005, Faouzi LYAZRHI wrote: Je voudrais faire une analyse des correspondances multiples avec R. avec les représentation graphiques correspondantes avec R. je ne sais pas comment procéder .. en vour remerciant par avance -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] package
hello, I created a package with my functions, and i wand to hide the code of some functions. Could you help me ? Grégory -- GAZ DE FRANCE Grégory Benmenzer DIRECTION DE LA RECHERCHE Pôle Economie Statistiques et Sociologie 361 Avenue du président Wilson - BP 33 93211 La Plaine Saint Denis cedex tel : 01 49 22 55 07 fax : 01 49 22 57 10 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] how to print error message in batch mode
On Thu, 7 Apr 2005, Jan T. Kim wrote: On Thu, Apr 07, 2005 at 01:58:55PM +0200, Elio Mineo wrote: This solution is fine, too. The Achim's solution is what I have asked with this slight modification: $ R -q --no-save prova prova.out 2 prova.out The canonical and perhaps more correct way to redirect stderr into stdout is 21, as in R -q --no-save prova prova.out 21 Which as I have already pointed out, is what R CMD BATCH does. [I hope you replied to Elio as well as to the list: the headers do not show it.] -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] package
Please define what you mean by `hide'. If the functions are not to be called by users directly, just create a NAMESPACE and export only the ones you want to expose to the users. However, the users can still get to those not exported if they really want to. If you don't want the users to be able to see the code through any means, I don't know if that's possible. Andy From: Gregory BENMENZER hello, I created a package with my functions, and i wand to hide the code of some functions. Could you help me ? Grégory -- GAZ DE FRANCE Grégory Benmenzer DIRECTION DE LA RECHERCHE Pôle Economie Statistiques et Sociologie 361 Avenue du président Wilson - BP 33 93211 La Plaine Saint Denis cedex tel : 01 49 22 55 07 fax : 01 49 22 57 10 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] package
you could use a namespace, look at Writing R extensions doc, section 1.6 I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.ac.be/biostat/ http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm - Original Message - From: Gregory BENMENZER [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Sent: Thursday, April 07, 2005 2:43 PM Subject: [R] package hello, I created a package with my functions, and i wand to hide the code of some functions. Could you help me ? Grégory -- GAZ DE FRANCE Grégory Benmenzer DIRECTION DE LA RECHERCHE Pôle Economie Statistiques et Sociologie 361 Avenue du président Wilson - BP 33 93211 La Plaine Saint Denis cedex tel : 01 49 22 55 07 fax : 01 49 22 57 10 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] analyse des correspondances multiples
Also, library(ade4) Best Marwan --- Marwan Khawaja http://staff.aub.edu.lb/~mk36/ --- -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent: Thursday, April 07, 2005 3:45 PM To: Faouzi LYAZRHI Cc: R-help@stat.math.ethz.ch Subject: Re: [R] analyse des correspondances multiples library(MASS) ?mca On Thu, 7 Apr 2005, Faouzi LYAZRHI wrote: Je voudrais faire une analyse des correspondances multiples avec R. avec les représentation graphiques correspondantes avec R. je ne sais pas comment procéder .. en vour remerciant par avance -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] apply
Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) However, I would like to pass an argument to the function mean, namely na.rm=TRUE Does anyone knows how to do this? Thanks in advance, Jan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] hex format
Hello world: Has anyone used hex notation within R to represents integers? Cheers, Steve Vejcik __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Is a .R script file name available inside the script?
I think you might want 'commandArgs()' which gives you the original command line call. -roger Darren Weber wrote: Hi, if we have a file called Rscript.R that contains the following, for example: x - 1:100 outfile = Rscript.Rout sink(outfile) print(x) and then we run source(Rscript.R) we get an output file called Rscript.Rout - great! Is there an internal variable, something like .Platform, that holds the script name when it is being executed? I would like to use that variable to define the output file name. Best, Darren __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] apply
On Apr 7, 2005, at 8:27 AM, malte wrote: Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) However, I would like to pass an argument to the function mean, namely na.rm=TRUE apply(sdata,2,function(x) {mean(x,na.rm=TRUE)}) Sean __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] apply
From: malte Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) However, I would like to pass an argument to the function mean, namely na.rm=TRUE Does anyone knows how to do this? aveBehav - apply(sdata, 2, mean, na.rm=TRUE) or more efficiently: aveBehav - colMeans(sdata, na.rm=TRUE) Read ?apply and look at the ... argument. If you don't understand how it works, try the example on that page. Andy Thanks in advance, Jan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] apply
try, apply(sdata, 2, mean, na.rm=TRUE) or # assuming `sdata' is a matrix colMeans(sdata, na.rm=TRUE) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.ac.be/biostat/ http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm - Original Message - From: malte [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Thursday, April 07, 2005 2:27 PM Subject: [R] apply Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) However, I would like to pass an argument to the function mean, namely na.rm=TRUE Does anyone knows how to do this? Thanks in advance, Jan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] half-normal residual plots
Hi all, I am trying to produce a half-normal plot of residuals from a GLM. I have found the qqnorm function for producing a normal plot but can't figure out how to produce a half-normal. Can anyone help with this? Thanks Malcolm -- MJ Price, Social Medicine [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Dta structure of LOADINGS class in factanal
Hi R users, I need some help in the followings: I'm doing factor analysis and I need to extract the loading values and the Proportion Var and Cumulative Var values one by one. Here is what I am doing: fact - factanal(na.omit(gnome_freq_r2),factors=5); fact$loadings Loadings: Factor1 Factor2 Factor3 Factor4 Factor5 b1freqr2 0.246 0.486 0.145 b2freqr2 0.129 0.575 0.175 0.130 0.175 b3freqr2 0.605 0.253 0.166 0.138 0.134 b4freqr2 0.191 0.220 0.949 b5freqr2 0.286 0.265 0.113 0.891 0.190 b6freqr2 0.317 0.460 0.151 b7freqr2 0.138 0.199 0.119 0.711 b8freqr2 0.769 0.258 0.195 0.137 b9freqr2 0.148 0.449 0.103 0.327 Factor1 Factor2 Factor3 Factor4 Factor5 SS loadings 1.294 1.268 1.008 0.927 0.730 Proportion Var 0.144 0.141 0.112 0.103 0.081 Cumulative Var 0.144 0.285 0.397 0.500 0.581 I can get the loadings using: fact$loadings[1,1] [1] 0.2459635 but I couldn't find the way to do the same with the Proportion Var and Cumulative Var values. Thanks, Tamas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] 4D Plot ??
Tried to post this last night, but it doesn't seem to have appeared. Using R 2.0.1 on Windows XP + SP2. I am traveling, away from my usual references. I'm trying to make a 4-dimensional plot: a levelplot with overlaid contours, with different response variables represented by (1) colors on the levelplot and (2) the contour lines. First try was filled.contour + contour but the key printed by the first means that the scales differ. Then I tried levelplot. I couldn't figure out how to pass 3 variables to levelplot, so I duplicated all rows of the data frame and changed the z data for the second half, in order to plot one half at a time. #-- # Try at a 4D contourplot: y = x = 1:50 grid - expand.grid(x=x, y=y) grid$z = sqrt(x*y) n1 = nrow(grid) grid2 = rbind(grid,grid) grid2$z[(n1+1):(n1*2)] = log(grid2$x[1:n1] * grid2$y[1:n1] + 10) panel.4d - function(x,y,z,subscripts) { n1 = 1; n2 = length(x)/2 panel.levelplot(x[n1:n2],y[n1:n2],z[n1:n2],subscripts,region=TRUE) n1 = n2 + 1 ; n2 = length(x) panel.levelplot(x[n1:n2],y[n1:n2],z[n1:n2],subscripts,region=FALSE,contour=TRUE) } aa = levelplot(z~x*y, data=grid2, cuts = 20, panel=panel.4d) print(aa) # This gives the following error message: Error in panel.levelplot(x[n1:n2], y[n1:n2], z[n1:n2], subscripts, region = FALSE, : NAs are not allowed in subscripted assignments although that it completes when the second levelplot is set to region=TRUE, contour=FALSE (though then the second plot then hides the first). Hints or sample code will be most welcome. Once this works, the next refinement will be to replace the colored levelplot with something similar but with smooth edges produced by contouring, so advice on that is also welcome. Michael Prager, Ph.D. NOAA Center for Coastal Fisheries Habitat Research Beaufort, North Carolina, USA __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] how to analysis this kind of data set?
hi,everybody I have a *time course* data set about a CML cell line treated by two drugs and their combination.The experiment was performed on cDNA microarray platform.The green channel of all the arrays are common,the untreated cell.Here follows the experiment design: a_0hr,a_3hr,a_8hr,a_12hr,a_24hr,a_48hr,a_72hr, b_3hr,b_8hr,b_12hr,b_24hr,b_48hr,b_72hr, ab_3hr,ab_8hr,ab_12hr,ab_24hr,ab_48hr,ab_72hr. A total of 19 *cDNA microarrays*.a_0hr means* *drug *a *treament *0 hours vs. control. *And a_3hrs means drug a treatment 3 hours vs. control.So for drug *b *and their combination *ab*(drug a and drug b added together).My goal is to identify the three sets of genes,the genes differentially expressed by drug a,the genes by drug b ,and their combination. i am thinking about *ANOVA* ,but i am not sure whether it is correct. Any comments,suggestions?Any R/bioconductor packages can be used?Thanks in advance [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] apply
malte [EMAIL PROTECTED] writes: aveBehav=c(apply(sdata, 2, mean)) aveBehav= apply(sdata, 2, mean, na.rm=TRUE) and ?apply will tell you about this. + seth __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] parameterisation of Factor levels
Hi all, I am trying to fit simple 2 factor (factors A and B - 3 and 2 levels respectively) poisson regression model. Inititially a pure error model is fitted and significance tests are performed on each parameter. I wish to remove an individual parameter from the model - the interaction between the second level of factor A and factor B. However i only seem to be able to remove the AB (all interactions between A and B) term, which is no use as it also removes the interaction term between level 3 of Factor A and factor B. Can anyone help with this. Thanks Malcolm -- MJ Price, Social Medicine [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like to do the opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. On Thu, 2005-04-07 at 08:45, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] apply
On 7 Apr 2005 at 14:27, malte wrote: Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) Hallo try aveBehav=c(apply(sdata, 2, mean, na.rm=T)) Cheers Petr However, I would like to pass an argument to the function mean, namely na.rm=TRUE Does anyone knows how to do this? Thanks in advance, Jan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Petr Pikal [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] half-normal residual plots
There is a halfnorm function in the faraway package. Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: [EMAIL PROTECTED] Tel: (651) 733-6092 Fax: (651) 736-3122 MJ Price, Social Medicine [EMAIL PROTECTED] To istol.ac.uk r-help@stat.math.ethz.ch Sent by: cc [EMAIL PROTECTED] at.math.ethz.ch Subject [R] half-normal residual plots 04/07/2005 08:43 AM Hi all, I am trying to produce a half-normal plot of residuals from a GLM. I have found the qqnorm function for producing a normal plot but can't figure out how to produce a half-normal. Can anyone help with this? Thanks Malcolm -- MJ Price, Social Medicine [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Importing data into R
I have a highly formated Excel with multiple tabs. Is it currently possible to read this data into R without changing the format of the Excel file? Also, is it possible to write back to the same Excel file or at least create a new Excel file with the same formatting as before with modified data which has been processed in R. Thanks in advance for any help that you can provide. Dave __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Prof Brian Ripley [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] On Thu, 7 Apr 2005, Steve Vejcik wrote: Has anyone used hex notation within R to represents integers? Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). I see this works fine with R 2.0.0 on a Linux platform, but doesn't work at all under R 2.0.1 on Windows. as.numeric(0x1AF0) [1] NA Warning message: NAs introduced by coercion Seems to me the conversion from hex to decimal should be system independent (and makes working with colors much more convenient). Why isn't this system independent now? The prefix on hex numbers is somewhat language dependent (0x or $) perhaps but I didn't think this conversion should be system dependent. I don't remember where I got this, but this hex2dec works under both Linux and Windows (and doesn't need the 0x prefix). hex2dec - function(hexadecimal) { hexdigits - c(0:9, LETTERS[1:6]) hexadecimal - toupper(hexadecimal)# treat upper/lower case the same decimal - rep(0, length(hexadecimal)) for (i in 1:length(hexadecimal)) { digits - match(strsplit(hexadecimal[i],)[[1]], hexdigits) - 1 decimal[i] - sum(digits * 16^((length(digits)-1):0)) } return(decimal) } Example: hex2dec(c(1AF0, )) [1] 6896 65535 can be interpreted as 65535 as unsigned and -1 as signed on the same system depending on context. This isn't system dependent, but rather context dependent. I suggest as.numeric should perform the unsigned conversion on all systems. What am I missing? efg -- Earl F. Glynn Scientific Programmer Bioinformatics Department Stowers Institute for Medical Research __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] package
On Apr 7, 2005 8:43 AM, Gregory BENMENZER [EMAIL PROTECTED] wrote: hello, I created a package with my functions, and i wand to hide the code of some functions. Could you help me ? Grégory There was some discussion on the list that there is work being done on an R compiler. I don't know what the status is or whether it would indeed solve your problem but you could try googling around for it or maybe someone else on the list can provide more info. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] 4D Plot ??
On Thu, 7 Apr 2005, Mike Prager wrote: Tried to post this last night, but it doesn't seem to have appeared. Using R 2.0.1 on Windows XP + SP2. I am traveling, away from my usual references. I'm trying to make a 4-dimensional plot: a levelplot with overlaid contours, with different response variables represented by (1) colors on the levelplot and (2) the contour lines. First try was filled.contour + contour but the key printed by the first means that the scales differ. You could put the contour() call in the plot.axes argument in the filled.contour() call. This is a useful trick for getting the right scales for an overlay. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] 4D Plot ??
Hi Mike, I've done a bit of playing around with these kind of plots for visualising microarray data (to eventually go into a bioconductor package). I've attached my code for producing surfaceplots (my name for the type of plots that includes both image and contour plots) - it's all lattice based, so you'll need some familiarity with how lattice works to understand how it all works. The key function is panel.superpose.surface which you can use as follows: levelplot(surfacevar + contourvar ~ x * y, data, panel=panel.surface.smooth, asp=iso) (note that the contours are automatically smoothed using image.smooth from the fields package - you can control the amount of smoothing use contour.theta) You can supply multiple contour variables, but be advised it gets messy really quickly! You can also smooth the surface by setting panel.base = panel.surface.smooth Hope this is helpful! Hadley # Microarray surface plot # Surface plot with sensible defaults for microarrays surfaceplot.ma - function( formula, data, aspect = xy, prepanel = function() {list(dx=1, dy=1)}, allow.multiple = TRUE, as.table = TRUE, panel = panel.ma.surface, col.regions = pal(), tipRows = NULL, tipCols = 0, main = NULL, ylab = NULL, ... ) { #if plot x vs. y, draw tip grid lf - latticeParseFormula(formula, data, dim=3) if (lf$right.x.name == x lf$right.y.name %in% c(y,-y)) { try({ tipCols = seq(0,max(data$x), by=max(data$spotCol))+0.5 tipRows = seq(0,max(data$y), by=max(data$spotRow))+0.5 }) } if (missing(main)) { main - paste(Surface plot of, lf$left.name) } levelplot(formula, data, allow.multiple=allow.multiple, as.table = as.table, panel=panel, main=main, prepanel=prepanel, aspect=aspect, col.regions=col.regions, ylab=ylab, tipRows = tipRows, tipCols = tipCols, ...) } panel.ma.surface - function(..., tipRows, tipCols, grid.col = grey80) { panel.superpose.surface(...) if (!missing(tipRows)) { panel.abline(h = tipRows, col = grid.col) } if (!missing(tipCols)) { panel.abline(v = tipCols, col = grid.col) } } panel.superpose.surface - function ( x, y, z, groups, subscripts, panel.base = panel.levelplot, panel.superpose = panel.contour.smooth, contour.col = #BF2828, ..., region, contour, at ) { if (missing(groups)) { panel.base(x = x, y = y, z = z, subscripts = subscripts, contour=FALSE, at=at, ...) return() } x - x[subscripts]; y - y[subscripts]; z - z[subscripts] groups - groups[subscripts] if (is.factor(groups)) { vals - levels(groups) } else { vals - sort(unique(groups)) } nvals - length(vals) first - TRUE for (i in seq(along = vals)) { id - (groups == vals[i]) if (first) { panel.base(x = x[id], y = y[id], z = z[id], groups = groups[id], subscripts = TRUE, contour=FALSE, region=TRUE, col=col, at=at, ...) first - FALSE } else { panel.superpose(x = x[id], y = y[id], z = z[id], groups = groups[id], subscripts = TRUE, contour=TRUE, region=FALSE, col=contour.col, at = seq(min(z[id]), max(z[id]), length=10), ...) } } } panel.surface.smooth - function(x, y, z, zcol, at, subscripts, surface.theta = 3, ..., contour, region, col.regions) { x - x[subscripts]; y - y[subscripts]; z - z[subscripts]; z.smooth - smooth.grid(x, y, z, surface.theta) zcol - cut(z.smooth, at, labels=FALSE) panel.levelplot(x, y, z.smooth, zcol=zcol, subscripts=TRUE, ..., at=at, contour = FALSE, region = TRUE, col.regions=col.regions) } panel.contour.smooth - function(x, y, z, zcol, at, subscripts, contour.theta = 3, ..., contour, region, col.regions) { x - x[subscripts]; y - y[subscripts]; z - z[subscripts]; z.smooth - smooth.grid(x, y, z, contour.theta) zcol - cut(z.smooth, at, labels=FALSE) panel.levelplot(x, y, z.smooth, zcol=zcol, subscripts=TRUE, ..., at=at, contour = TRUE, region = FALSE, col.regions=col.regions) } smooth.grid - function(x, y, z, theta=2) { if (!require(fields)) {stop(Fields library required for smoothed panels)} image.smooth(matrix(z[order(y, x)], nrow=max(x)), theta=theta)[(y-1)*max(x) + x] } # Modified xy panel function that scales spots to z values panel.spot - function (x, y, z, subscripts, col, ...) { z.range - range(z, na.rm = TRUE) x - as.numeric(x)[subscripts]
Re: [R] package
You may want to think about using a package NAMESPACE. I don't know if that is what you mean, but it is something available for making some functions public and others package private, but it doesn't give a mechanism to ABSOLUTELY hide code. Sean On Apr 7, 2005, at 11:08 AM, Gabor Grothendieck wrote: On Apr 7, 2005 8:43 AM, Gregory BENMENZER [EMAIL PROTECTED] wrote: hello, I created a package with my functions, and i wand to hide the code of some functions. Could you help me ? Grégory There was some discussion on the list that there is work being done on an R compiler. I don't know what the status is or whether it would indeed solve your problem but you could try googling around for it or maybe someone else on the list can provide more info. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Assigning dates attribute
Dear List, I have a one-column data set in .csv format. I used read.csv to import the data into R, as follows: x - read.csv(data.csv, header = TRUE, sep = ,) The data points have a 'dates' attribute, which is in a separatel .csv file. I used the same command as above to import it into R. To assoicate the 'dates' attribute with the data points, I did: attributes(x)-date Which resulted in: Error in attributes-(`*tmp*`, value = date) : attributes must be in a list So then I did: attributes(x)-list(date) Again, got an error, though slightly different this time: Error in attributes-(`*tmp*`, value = list(date)) : attributes must be named Any help is appreciated. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Dta structure of LOADINGS class in factanal
At 9:44 AM -0400 4/7/05, Tamas Gal wrote: Hi R users, I need some help in the followings: I'm doing factor analysis and I need to extract the loading values and the Proportion Var and Cumulative Var values one by one. Here is what I am doing: fact - factanal(na.omit(gnome_freq_r2),factors=5); fact$loadings Loadings: Factor1 Factor2 Factor3 Factor4 Factor5 b1freqr2 0.246 0.486 0.145 ... b9freqr2 0.148 0.449 0.103 0.327 Factor1 Factor2 Factor3 Factor4 Factor5 SS loadings 1.294 1.268 1.008 0.927 0.730 Proportion Var 0.144 0.141 0.112 0.103 0.081 Cumulative Var 0.144 0.285 0.397 0.500 0.581 I can get the loadings using: fact$loadings[1,1] [1] 0.2459635 but I couldn't find the way to do the same with the Proportion Var and Cumulative Var values. Although not pretty, try colSums(fact$loading*fact$loading)/dim(fact$loading)[1] for the proportion Var and cumsum(colSums(fact$loading*fact$loading)/dim(fact$loading)[1]) to get the cumulative Var values Bill -- William Revelle http://pmc.psych.northwestern.edu/revelle.html Professor http://personality-project.org/personality.html Department of Psychology http://www.wcas.northwestern.edu/psych/ Northwestern University http://www.northwestern.edu/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Earl F. Glynn wrote: Prof Brian Ripley [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] On Thu, 7 Apr 2005, Steve Vejcik wrote: Has anyone used hex notation within R to represents integers? Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). I see this works fine with R 2.0.0 on a Linux platform, but doesn't work at all under R 2.0.1 on Windows. as.numeric(0x1AF0) [1] NA Warning message: NAs introduced by coercion Seems to me the conversion from hex to decimal should be system independent (and makes working with colors much more convenient). Why isn't this system independent now? Presumably because nobody thought it was important enough to make it so. R isn't a low level system programming language, so why should it treat hex specially? Duncan Murdoch __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Steve Vejcik wrote: Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like to do the opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. you can use chcode() to define hex.to.dec(), dec.to.hex() and sum.hex() to operate with hex numbers. Peter Wolf -- define chcode= chcode - function(b, base.in=2, base.out=10, digits=0123456789ABCDEF){ # change of number systems, pwolf 10/02 # e.g.: from 2 2 2 2 ... - 16 16 16 ... digits-substring(digits,1:nchar(digits),1:nchar(digits)) if(length(base.in)==1) base.in - rep(base.in, max(nchar(b)-1)) if(is.numeric(b)) b - as.character(as.integer(b)) b.num - lapply(strsplit(b,), function(x) match(x,digits)-1 ) result - lapply(b.num, function(x){ cumprod(rev(c(base.in,1))[ 1:length(x) ] ) %*% rev(x) } ) number-unlist(result) cat(decimal representation:,number,\n) if(length(base.out)==1){ base.out-rep(base.out,1+ceiling(log( max(number), base=base.out ) ) ) } n.base - length(base.out); result - NULL for(i in n.base:1){ result - rbind(number %% base.out[i], result) number - floor(number/base.out[i]) } result[]-digits[result+1] apply(result, 2, paste, collapse=) } @ define hex.to.dec, dec.to.hex and sum.hex= hex.to.dec-function(x) as.numeric(chcode(x, base.in=16, base.out=10)) dec.to.hex-function(x) chcode(x, base.in=10, base.out=16) sum.hex-function(x,y) dec.to.hex(hex.to.dec(x) + hex.to.dec(y)) @ quick test: define hex numbers= a-dec.to.hex(10); print(a) b-dec.to.hex(3);print(b) @ output-start decimal representation: 10 [1] 0A decimal representation: 3 [1] 03 output-end @ sum of a and b= sum.hex(a,b) @ output-start decimal representation: 10 decimal representation: 3 decimal representation: 13 Thu Apr 7 17:31:42 2005 [1] 0D output-end __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Assigning dates attribute
On Thu, 7 Apr 2005 08:26:41 -0700 (PDT) JTW wrote: Dear List, I have a one-column data set in .csv format. I used read.csv to import the data into R, as follows: x - read.csv(data.csv, header = TRUE, sep = ,) The data points have a 'dates' attribute, which is in a separatel .csv file. I used the same command as above to import it into R. To assoicate the 'dates' attribute with the data points, I did: attributes(x)-date Which resulted in: Error in attributes-(`*tmp*`, value = date) : attributes must be in a list So then I did: attributes(x)-list(date) Again, got an error, though slightly different this time: Error in attributes-(`*tmp*`, value = list(date)) : attributes must be named Any help is appreciated. The error message is pretty informative, the assignment needs a named list, e.g.: R x - 1:10 R attributes(x) - list(foo = letters[1:10]) R x [1] 1 2 3 4 5 6 7 8 9 10 attr(,foo) [1] a b c d e f g h i j Note, that this will strip off all other attributes. To add one attribute, you can do R attr(x, bar) - LETTERS[1:10] R x [1] 1 2 3 4 5 6 7 8 9 10 attr(,foo) [1] a b c d e f g h i j attr(,bar) [1] A B C D E F G H I J Furthermore, the data you describe look like a time series. So you might want to store the data as a time series. For time series with a date attribute of class Date, look at the zoo package. Z __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Steve Vejcik wrote: Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like You don't understand how R works: x - as.numeric(0x1AF0) produces an number, not its decimal representation. A number is a number is a number irrepsective of the the base of its character representation. to dothe opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. On Thu, 2005-04-07 at 08:45, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Importing data into R
I don't know if it does what you want, however you might try package RExcel. However it is not on CRAN. You can find it on http://sunsite.univie.ac.at/rcom/download/. I belive it might be obsolete and replaced by R (D)COM Server V1.35 (previously you needed this package to use RExcel) which you can find on http://cran.planetmirror.com/contrib/extra/dcom/RSrv135.html (description) or http://cran.planetmirror.com/contrib/extra/dcom/RSrv135.exe (add-on). I hope this helps! Ales Ziberna - Original Message - From: Dave Evens [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Thursday, April 07, 2005 4:47 PM Subject: [R] Importing data into R I have a highly formated Excel with multiple tabs. Is it currently possible to read this data into R without changing the format of the Excel file? Also, is it possible to write back to the same Excel file or at least create a new Excel file with the same formatting as before with modified data which has been processed in R. Thanks in advance for any help that you can provide. Dave __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Duncan Murdoch wrote: [...] If you want an integer vector to always display in hex, assign a class to it and define a print method. I don't think there's a standard library function to display in hex, but there are probably packages to do so. In R 2.1.0-to-be x - as.numeric(0x00B) # this is platform-specific x [1] 11 sprintf(0x%X, as.integer(x)) # this is not [1] 0xB -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] 2d plotting and colours
And does this work? n - 5 par(mfrow = c(2,2)) palette(default) barplot(1:25,col = 1:25) pal - rainbow(n) barplot(1:25,col = pal[(1:25-1)%%n+1]) pal - rgb((0:15)/15, g=0,b=0, names=paste(red,0:15,sep=.)) barplot(1:25,col = pal[(1:25-1)%%n+1]) Earl F. Glynn wrote: Mulholland, Tom [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Since I was only concentrating on colour issues and not on your specific problem I was just showing the possibilities. Does this code help n - 5 par(mfrow = c(2,2)) palette(default) barplot(1:25,col = 1:25) palette(rainbow(n)) barplot(1:25,col = 1:25) palette(rgb((0:15)/15, g=0,b=0, names=paste(red,0:15,sep=.))) barplot(1:25,col = 1:25) require(cluster) x - runif(100) * 8 + 2 cl - kmeans(x, n) palette(rainbow(n)) plot(x, col = cl$cluster) abline(h = cl$centers, lty = 2,col = grey ) palette(palette()[order(cl$centers)]) points(x,col = cl$cluster,pch = 20,cex = 0.4) Using Windows with R 2.0.1 this looks fine at first. But when I resize the graphic, copy the graphic to a metafile and paste it into Word, or go to an earlier graphic and come back using History, the colors ae all messed up. It's as if only the last palette is being used for all four plots in the figure. Oddly, if I copy the graphic as a bitmap, the colors are preseved in the bitmap. Is this a quirk of my machine or does this happen for others? Is it possible that the Windows palette manager is being used (which is such about obsolete) and that true color graphics are not being used (which is the easist way to avoid headaches from the Windows palette manager)? efg __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Duncan Murdoch [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Seems to me the conversion from hex to decimal should be system independent (and makes working with colors much more convenient). Why isn't this system independent now? Presumably because nobody thought it was important enough to make it so. R isn't a low level system programming language, so why should it treat hex specially? 1) While generally I'd agree with your statement, manipulating colors is one place the ability to convert to/from hex would be quite nice. rgb(1,0,0.5) [1] #FF0080 rgb returns a hex string and then R makes manipulating this string somewhat difficult. One might want to use such color values to convert to a different color space, perform some sort of manipulation in that other color space, and then convert back to rgb. 2) I would think that one of R's mathematical abilities would be to provide a way to convert from any base to base 10, and from base 10 to any base. I haven't found this general math tool yet in R. Working with base-16 (or even base 2 sometimes) could be done with such a general math tool. 3) While I may be in a minority, I would even consider exporting IEEE floating-point numbers in hex form as a way to avoid any additional conversion losses converting to/from decimal. 4) Why not make working with raw data a little easier? readbin shows hex values but they are not easy to work with inside of R. IntegerSize - 4# How do I get this value from R? i - -2:2 i [1] -2 -1 0 1 2 length(i) [1] 5 object.size(i) [1] 52 writeBin(i, big.bin, endian=big) big - readBin(big.bin, raw, length(i)*IntegerSize) big [1] ff ff ff fe ff ff ff ff 00 00 00 00 00 00 00 01 00 00 00 02 writeBin(i, little.bin, endian=little) little - readBin(little.bin, raw, length(i)*IntegerSize) little [1] fe ff ff ff ff ff ff ff 00 00 00 00 01 00 00 00 02 00 00 00 5) Does R have a hex consistency problem? The color values start with a # for hex, but the as.numeric(#FF0080) isn't allowed? Thanks for your time. efg __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 2005-04-07 at 11:06, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like You don't understand how R works: x - as.numeric(0x1AF0) produces an number, not its decimal representation. A number is a number is a number irrepsective of the the base of its character representation. as.numeric(0x1AF0) returns a decimal value for a hex string. If you prefer, substitute the word shows for returns. to dothe opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. On Thu, 2005-04-07 at 08:45, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] hex format
From: Steve Vejcik On Thu, 2005-04-07 at 11:06, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like You don't understand how R works: x - as.numeric(0x1AF0) produces an number, not its decimal representation. A number is a number is a number irrepsective of the the base of its character representation. as.numeric(0x1AF0) returns a decimal value for a hex string. If you prefer, substitute the word shows for returns. You don't seem to get the point. as.numeric() is a function that _returns_ a _value_. How you want that _value_ to be _shown_ is a different matter. Would you substitute `I gave the money to the cashier' with `I showed the money to the cashier'? Andy to dothe opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. On Thu, 2005-04-07 at 08:45, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Fitting a mixed negative binomial model
Dear list members, I want to fit a nonlinear mixed model using the nlme command. My dependent variable takes the form of event counts for different countries over a number of years, and hence I was going to fit a mixed effects negative binomial model. The problem, as far as I can glean from Pinheiro Bates 2000, is that I need a model that is not normal in the errors. All the models they discuss have linear error structures. Is there a package in the R language that fits a negative binomial mixed effects model? Thank you, Jose Aleman PhD Candidate Politics Department 130 Corwin Hall Princeton, NJ 08544 609.937.0190 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] hex format
I understand that point. Again: I would like to have numbers represented to me in hexidecimal format, not decimal format. This was my original query and I think it's clear. Let me try another variation: I would like R to recognize that I am using hexadecimal notation when I type a number at the keyboard. I would like R to have the ability to show me an integer expressed in hexadecimal format. On Thu, 2005-04-07 at 12:12, Liaw, Andy wrote: From: Steve Vejcik On Thu, 2005-04-07 at 11:06, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Thanks for your advice. Unfortunately, your answers are inconsistent: as.numeric(0x1AF0) returns a decimal value for a hex string. I'd like You don't understand how R works: x - as.numeric(0x1AF0) produces an number, not its decimal representation. A number is a number is a number irrepsective of the the base of its character representation. as.numeric(0x1AF0) returns a decimal value for a hex string. If you prefer, substitute the word shows for returns. You don't seem to get the point. as.numeric() is a function that _returns_ a _value_. How you want that _value_ to be _shown_ is a different matter. Would you substitute `I gave the money to the cashier' with `I showed the money to the cashier'? Andy to dothe opposite-use hex notation to represent a decimal. e.g. x-0x000A y-0x0001 x+y=0x00B Cheers. On Thu, 2005-04-07 at 08:45, Prof Brian Ripley wrote: On Thu, 7 Apr 2005, Steve Vejcik wrote: Hello world: Has anyone used hex notation within R to represents integers? That's a spectacularly vague question. Short answer: yes. as.numeric(0x1AF0) [1] 6896 (which BTW is system-dependent, but one person used it as you asked). PLEASE read the posting guide and try for a `smarter' question. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Notice: This e-mail message, together with any attachment...{{dropped}} __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Earl F. Glynn wrote: 1) While generally I'd agree with your statement, manipulating colors is one place the ability to convert to/from hex would be quite nice. rgb(1,0,0.5) [1] #FF0080 rgb returns a hex string and then R makes manipulating this string somewhat difficult. One might want to use such color values to convert to a different color space, perform some sort of manipulation in that other color space, and then convert back to rgb. The convertColor function in R 2.1.0 provides colorspace conversion, including hex. 5) Does R have a hex consistency problem? The color values start with a # for hex, but the as.numeric(#FF0080) isn't allowed? #ff0080 isn't a number, it's a colour (or perhaps a color). If it were converted to numeric form it would be a vector of three numbers, and which three numbers would depend on the coordinate system used for colour space. For example, R already provides both hsv() and rgb() to create colours from vectors of three numbers, but the correspondence is different in each case. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, Apr 07, 2005 at 11:58:48AM -0500, Earl F. Glynn wrote: Duncan Murdoch [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Seems to me the conversion from hex to decimal should be system independent (and makes working with colors much more convenient). Why isn't this system independent now? Presumably because nobody thought it was important enough to make it so. R isn't a low level system programming language, so why should it treat hex specially? 1) While generally I'd agree with your statement, manipulating colors is one place the ability to convert to/from hex would be quite nice. rgb(1,0,0.5) [1] #FF0080 rgb returns a hex string and then R makes manipulating this string somewhat difficult. I'd like to second this opinion. It just occasionally happens that data are available in some variant of hex format, and I've had the impression that getting such data into R is a bit less convenient than it could be. One might want to use such color values to convert to a different color space, perform some sort of manipulation in that other color space, and then convert back to rgb. 2) I would think that one of R's mathematical abilities would be to provide a way to convert from any base to base 10, and from base 10 to any base. I haven't found this general math tool yet in R. Working with base-16 (or even base 2 sometimes) could be done with such a general math tool. In fact, the ANSI C function strtol already provides conversion to any base between 2 and 36, so R's mathematical capabilities don't even need to be invoked here. An R function strtol(x, base), x being a character variable and base an integer between 2 and 36, would probably add a bit of convenience. I've never programmed that, though -- seems that I'm one of those to whom this hasn't been important enough. If it is done some day, I'd favour the strtol function over having as.numeric interpret the (rather C-ish) 0x prefix. I wasn't aware that this currently works on some platforms (and I'm glad it doesn't interpret the 0 prefix for octal, as C does, making 007 legal and 008 not. ;-) ) Best regards, Jan -- +- Jan T. Kim ---+ |*NEW*email: [EMAIL PROTECTED] | |*NEW*WWW: http://www.cmp.uea.ac.uk/people/jtk | *-= hierarchical systems are for files, not for humans =-* __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] ks.test for conditional distribution Y|x
Dear experts, Is it possible to use ks.test function to check the goodness of fit of the conditional distribution Y|X=x? For example, I would like to check that my data (Y,X) come from Norm(0.5+x,1) using KS. Thank you in advance, Victoria Landsman. [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Stratified Bootstrap question
Dear Tim, Thank you very much for taking time giving me advices on my questions. I talked with my professor about this bootstrapping question whether to resample clinic or resample clinic + resample patients within clinic. I was told that the second method might destroy the correlation structure between the patients within a clinic. So I am thinking if it is worthy that I do a simulation to compare the two kinds of bootstrapping method. I mean, is this comparision meaningful and is it worth of doing? What do you think? Thank you. Qian On 1 Apr 2005, Tim Hesterberg wrote: Qian wrote: I talked with my advisor yesterday about how to do bootstrapping for my scenario: random clinic + random subject within clinic. She suggested that only clinic are independent units, so I can only resample clinic. But I think that since subjects are also independent within clinic, shall I resample subjects within clinic, which means I have two-stage resampling? Which one do you think makes sense? This is a tough issue; I don't have a complete answer. I'd appreciate input from other r-help readers. If you randomly select clinics, then randomly select patients within the clinics: (1) by bootstrapping just clinics, you capture both sources of variation -- the between-subject variation is incorporated in the results for each clinic. (2) by bootstrapping clinics, then subjects within clinics, you end up double-counting the between-subject variation That argues for resampling just clinics. By analogy, if you have multiple subjects, and multiple measurements per subject, you should just resample subjects. However, I'm not comfortable with this if you have a small number of clinics, and relatively large numbers of patients in each clinic, and think that the between-clinic variation should be small. Then it seems better to resample both clinics and patients. I'm leery about resampling just clinics if there are a small number of clinics. Bootstrapping isn't particularly effective for small samples -- it is subject to skewness in small samples, and it underestimates variances (it's advantages over classical methods really show up with medium size samples). There are remedies for the small variance, see Hesterberg, Tim C. (2004), Unbiasing the Bootstrap-Bootknife Sampling vs. Smoothing, Proceedings of the Section on Statistics and the Environment, American Statistical Association, 2924-2930 www.insightful.com/Hesterberg/articles/JSM04-bootknife.pdf Tim Hesterberg | Tim Hesterberg Research Scientist | | [EMAIL PROTECTED] Insightful Corp.| | (206)802-23191700 Westlake Ave. N, Suite 500 | | (206)283-8691 (fax) Seattle, WA 98109-3044, U.S.A. | | www.insightful.com/Hesterberg | Download the S+Resample library from www.insightful.com/downloads/libraries *** Qian An Division of Biostatistics University of Minnesota (phone) 612-626-2263 (fax) 612-626-8892 Email: [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] ks.test for conditional distribution Y|x
Couldn't you do this by subtracting 0.5 + x from your y values and checking for normality with mean 0 and sd = 1 (using ks.test or another test of normality). If you fail, you'll have to do additional work to find out whether pairs with some particular x value (or range of x values) is causing the problem, but I think this fits the question as stated. Of course, if you have discrete x values, and enough data at each one, you could just run the check for each x. Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Vicky Landsman Sent: Thursday, April 07, 2005 3:16 PM To: R-help list Subject: [R] ks.test for conditional distribution Y|x Dear experts, Is it possible to use ks.test function to check the goodness of fit of the conditional distribution Y|X=x? For example, I would like to check that my data (Y,X) come from Norm(0.5+x,1) using KS. Thank you in advance, Victoria Landsman. [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
Thomas Lumley [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] The convertColor function in R 2.1.0 provides colorspace conversion, including hex. #ff0080 isn't a number, it's a colour (or perhaps a color). If it were converted to numeric form it would be a vector of three numbers, and which three numbers would depend on the coordinate system used for colour space. Colo(u)rs and numbers are interchangeable to me. When you look at a picture, don't you see numbers? Maybe you don't see a number here, but I do. #ff0080 is interpreted in some (non-R) contexts as a single number. In many contexts, including HTML, colors are represented as three bytes in hex with this notation and the # means hexadecimal. The RGB color componets can be discerned quite easily: hex FF is decimal 255 (red), hex 00 is decimal 0 (green), hex 80 is decimal 128 (blue). Some programs, e.g., Dreamweaver, allow specification of colors in this hex 3-byte form directly. The 16 million colors you seen on a true color display are from the 256*256*256 (or in hex FF*FF*FF) possible RGB triples. For example, R already provides both hsv() and rgb() to create colours from vectors of three numbers, but the correspondence is different in each case. Sorry if some consider this off topic: HSV as a color space is really only liked by computer scientists. Image processing and color engineers rarely if ever use HSV. There are MANY other color spaces and computations possible (see color spaces or color conversions or other color topics on this page http://www.efg2.com/Lab/Library/Color/Science.htm). Most of these color manipulations in R are not easy because the very first step, converting colors, I mean numbers g, like #ff0080 to the red, green components is hindered because one must reinvent the wheel of hex-to-decimal conversion. Perhaps R will someday introduce a pixel type that would encapsulate the three color components (for color images at least). A matrix of pixels could easily be made into an image. Some color computations such a Maxwell Triangle, or a CIE Chromaticity Chart (sorry the links are currently broken, but the image can be seen on this Chinese translation page) http://bluemoon.myrice.com/efg/color/chromaticity.htm in R is more difficult than it should be because of how R is designed now. Many image processing statistical problems could be tackled directly in R if there were an easier way to manipulate pixels and images. But the hex manipulations I'm advocating could be used for variety of other purposes. E.g, I must periodically deal with a binary data stream of flow cytometery data -- part ASCII, part binary. Reading this stream directly from R would be nice and is almost doable. Working with raw data and understanding exactly what you've got would be facilitated by better conversion capabilities within R. efg __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Adapt Function Examples
I have read the help file for Adapt, but I cannot create a functn that works. I believe this is because I do not understand how to do this, and I have not found any working examples posted in the help. I have recieved many different errors in my attempts. Please post a simple but working use of the adapt function in 2 dimensions other than the one in the help(adapt) or explain it differently. Thank you __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] /bin/exec/R: No such file or directory
I have just installed R-2.0.1 from R-2.0.1.tar.gz on SUSe 9.1 64bit. When I am trying to launch R: R_HOME_DIR/bin/R; I am getting following message: ./R: line 151: /R_HOME_DIR/bin/exec/R: No such file or directory ./R: line 151: exec: /R_HOME_DIR/bin/exec/R: cannot execute: No such file or directory I do not have exec directory in bin directory. Does anybody know what went wrong? Thank you. Jarmila. Jarmila Bohmanova University of Georgia Department of Animal and Dairy Science Athens, GA __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Earl F. Glynn wrote: For example, R already provides both hsv() and rgb() to create colours from vectors of three numbers, but the correspondence is different in each case. Sorry if some consider this off topic: HSV as a color space is really only liked by computer scientists. Image processing and color engineers rarely if ever use HSV. There are MANY other color spaces and computations possible (see color spaces or color conversions or other color topics on this page http://www.efg2.com/Lab/Library/Color/Science.htm). Most of these color manipulations in R are not easy because the very first step, converting colors, I mean numbers g, like #ff0080 to the red, green components is hindered because one must reinvent the wheel of hex-to-decimal conversion. Yes, and convertColor in R-devel does quite a few of these (XYZ tristimulus space; CIE Lab and Luv; sRGB, Apple RGB and roll-your-own RGB based on chromaticities of the primaries; and chromatic adaptation for changing the white point). The colorspace package has a more elegant implementation of a somewhat different set of color space computations, and R-devel also has hcl() for specifying colors based on hue, chroma, and luminance (polar coordinates in Luv space). Basing R graphics on these (and so making them colours rather than just data about colours) requires a further step of considering the characteristics of the output device. This might be as simple as declaring R's output to be sRGB or as complicated as worrying about ICC profiles. -thomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] vectorized approach to cumulative sampling
Hi All, I need to sample a vector (old), with replacement, up to the point where my vector of samples (new) sums to a predefined value (target), shortening the last sample if necessary so that the total sum (newsum) of the samples matches the predefined value. While I can easily do this with a while loop (see below for example code), because the length of both old and new may be 20,000, a vectorized approach will save me lots of CPU time. Any suggestions would be greatly appreciated. Thanks, Dan # loop approach old=c(1:10) p=runif(1:10) target=20 newsum=0 new=NULL while (newsumtarget) { i=sample(old, size=1, prob=p); new[length(new)+1]=i; newsum=sum(new) } new newsum target if(newsumtarget){new[length(new)]=target-sum(new[-length(new)])} new newsum=sum(new); newsum target -- Daniel E. Bunker Associate Coordinator - BioMERGE Post-Doctoral Research Scientist Columbia University Department of Ecology, Evolution and Environmental Biology 1020 Schermerhorn Extension 1200 Amsterdam Avenue New York, NY 10027-5557 212-854-9881 212-854-8188 fax deb37ATcolumbiaDOTedu __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] hex format
On Thu, 7 Apr 2005, Earl F. Glynn wrote: ... picture, don't you see numbers? Maybe you don't see a number here, but I do. #ff0080 is interpreted in some (non-R) contexts as a single number. In many contexts, including HTML, colors are represented as three bytes in hex with this notation and the # means hexadecimal. The RGB color componets can be discerned quite easily: hex FF is decimal 255 (red), hex 00 is decimal 0 (green), hex 80 is decimal 128 (blue). Some programs, e.g., Dreamweaver, allow specification of colors in this hex 3-byte form directly. The 16 million colors you seen on a true color display are from the 256*256*256 (or in hex FF*FF*FF) possible RGB triples. For example, R already provides both hsv() and rgb() to create colours from vectors of three numbers, but the correspondence is different in each case. Sorry if some consider this off topic: HSV as a color space is really only liked by computer scientists. Image processing and color engineers rarely if ever use HSV. There are MANY other color spaces and computations possible (see color spaces or color conversions or other color topics on this page http://www.efg2.com/Lab/Library/Color/Science.htm). Most of these color manipulations in R are not easy because the very first step, converting colors, I mean numbers g, like #ff0080 to the red, green components is hindered because one must reinvent the wheel of hex-to-decimal conversion. I think R has the hex to decimal OK, but might be lacking in the decimal to hex case zz-function(x){ x-as.numeric(sub(#,'0x',x)); c(x%/%256^2, x%/%256%%256, x%%256) } zz('#0f0e0d') [1] 15 14 13 zz('#ff0080') [1] 255 0 128 If you already have the 3 byte triplet in read in as a binary, the same integer arithmetic does the extraction. Perhaps R will someday introduce a pixel type that would encapsulate the three color components (for color images at least). A matrix of pixels could easily be made into an image. Some color computations such a Maxwell Triangle, or a CIE Chromaticity Chart (sorry the links are currently broken, but the image can be seen on this Chinese translation page) http://bluemoon.myrice.com/efg/color/chromaticity.htm in R is more difficult than it should be because of how R is designed now. Many image processing statistical problems could be tackled directly in R if there were an easier way to manipulate pixels and images. But the hex manipulations I'm advocating could be used for variety of other purposes. E.g, I must periodically deal with a binary data stream of flow cytometery data -- part ASCII, part binary. Reading this stream directly from R would be nice and is almost doable. Working with raw data and understanding exactly what you've got would be facilitated by better conversion capabilities within R. I'm still not sure what you mean by hex manipulations. R has string manipulations, hex-to-number manipulations, binary-file-to-number manipulations, mixed file to number manipulations, and number to number manipulations. What I think you are asking for is /displaying/ numbers. Since R's sprintf() doesn't support the %x, (or %o, or %u) formats, I'm not sure how to use R to translate the number 257 into #000101 zzinv-function(x){} # such that: zzinv(257) #or zzinv(c(0,1,1)) #000101 Is zzinv() the operation you need? Dave -- Dr. David Forrest [EMAIL PROTECTED](804)684-7900w [EMAIL PROTECTED] (804)642-0662h http://maplepark.com/~drf5n/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Zipping Rdata Files
Saving Rdata files in a zip archive form can in some cases save a considerable amount of disk space. R has the zip.file.extract function to extract files from zip archives, but appears not to have any corresponding function to save in zipped form. (At least I have not been able to find anything in the help files or through searching the mail archives.) The system function can be used to call gzip or some other utility, but perhaps there is a more direct method. Also, when I use gzip to zip a file, I get an error message when using zip.file.extract to extract the file as follows: save(trips, file=trips.Rdata) system(gzip trips.Rdata) # saves trips.Rdata in an archive named trips.Rdata.gz load(zip.file.extract(trips.Rdata, trips.Rdata.gz)) [1] trips.Rdata Warning message: error 1 in extracting from zip file Setting options(unzip=gunzip) or options(unzip=gunzip.exe) does not solve the error. load(zip.file.extract(trips.Rdata, trips.Rdata.gz)) Error in open.connection(con, rb) : unable to open connection In addition: Warning message: cannot open compressed file `trips.Rdata' Of course I could reverse the process with, system(gunzip trips.Rdata.gz) load(trips.Rdata) but perhaps there is a simpler solution. P.S. I'm running R 2.0.1 on a Windows XP computer. Brian Gregor, P.E. Transportation Planning Analysis Unit Oregon Department of Transportation [EMAIL PROTECTED] (503) 986-4120 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] off-topic question: Latex and R in industries
Wensui, I work for 'A' electronics test equipment corporation. I have been using R ( since 1.6 ) instead of MATLAB etc. as a general language for data analysis and graph generation. On they way to R I tried Python/Scipy, Scilab and others - but R wins in quality and ease of use (it just needs DSP and GPIB/HPIB libraries to be perfect ). LaTeX is also my document tool of choice .. However LaTeX generated pdfs sent out as reports are much disliked. MS Word, PowerPoint and Excel are the standards, and very importantly they offer cut and paste ability across the larger team. MS's offerings comes no where near to the quality of LaTex / R, but in world of shared authorship - it's a one sided battle. My other PC universe vs Unix/OS X problem is vector / Meta-file graphics - essential for quality reports. Postscript, PDF and MS products just don't play. The newest Office and Visio versions seem to be dropping even more of the postscript import and export filters ( which never work very well anyway ). I have never met any other colleagues who use LaTeX or R. Any one else sharing the same experiences ? Message: 37 Date: Wed, 6 Apr 2005 11:38:55 -0400 From: Wensui Liu [EMAIL PROTECTED] Subject: [R] off-topic question: Latex and R in industries To: r-help@stat.math.ethz.ch Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=ISO-8859-1 Latex and R are really cool stuff. I am just wondering how they are used in industry. But based on my own experience, very rare. Why? How about the opinion of other listers? Thanks. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] vectorized approach to cumulative sampling
Hi, sample() takes a replace argument, so you can take large samples, with replacement, like this: (In the sample() call, the 50*target/mean(old) should make it sample 50 times more than likely. This means the while loop will probably get executed only once. This could be tuned easily, and there may be better ways of guessing how much to take). old - c(1:2000) p - runif(1:2000) target - 4000 new - 0 while ( sum(new) target ) new - sample(old, 50*target/mean(old), TRUE, p) i - which(cumsum(new) = target)[1] new - new[1:i] new[i] - new[i] - (sum(new)-target) Cheers, Rich On Apr 8, 2005 9:19 AM, Daniel E. Bunker [EMAIL PROTECTED] wrote: Hi All, I need to sample a vector (old), with replacement, up to the point where my vector of samples (new) sums to a predefined value (target), shortening the last sample if necessary so that the total sum (newsum) of the samples matches the predefined value. While I can easily do this with a while loop (see below for example code), because the length of both old and new may be 20,000, a vectorized approach will save me lots of CPU time. Any suggestions would be greatly appreciated. Thanks, Dan # loop approach old=c(1:10) p=runif(1:10) target=20 newsum=0 new=NULL while (newsumtarget) { i=sample(old, size=1, prob=p); new[length(new)+1]=i; newsum=sum(new) } new newsum target if(newsumtarget){new[length(new)]=target-sum(new[-length(new)])} new newsum=sum(new); newsum target -- Rich FitzJohn rich.fitzjohn at gmail.com |http://homepages.paradise.net.nz/richa183 You are in a maze of twisty little functions, all alike __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] vectorized approach to cumulative sampling
On 07-Apr-05 Daniel E. Bunker wrote: Hi All, I need to sample a vector (old), with replacement, up to the point where my vector of samples (new) sums to a predefined value (target), shortening the last sample if necessary so that the total sum (newsum) of the samples matches the predefined value. While I can easily do this with a while loop (see below for example code), because the length of both old and new may be 20,000, a vectorized approach will save me lots of CPU time. Any suggestions would be greatly appreciated. Thanks, Dan Hi Dan, You should be able to adapt the following vectorised approach to your specific needs: old-0.001*(1:1000) new-sample(old,1,replace=TRUE,prob=p) target-200 min(which(cumsum(new)target)) ## [1] 385 This took only a fraction of a second on my medium-speed machine. If you get an Inf as result, then 'new' doesn't add up to 'target', so you have to extend it. Hoping this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED] Fax-to-email: +44 (0)870 094 0861 Date: 07-Apr-05 Time: 22:46:12 -- XFMail -- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Zipping Rdata Files
Brian.J.GREGOR at odot.state.or.us writes: Saving Rdata files in a zip archive form can in some cases save a considerable amount of disk space. R has the zip.file.extract function to I suspect you may want to read up on the compress=TRUE option to the save() function. Hth, Dirk __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Principle Component Loadings
Dear R Could you help here I'm trying to decifer what the principle component loadings are in an R output. Are they in any way related to eigen vectors or eigen values? Brett Stansfield __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] off-topic question: Latex and R in industries
Greetings, Adobe Illustrator works with PDFs, either directly or by converting them to Illustrator format. These vector graphics have infinite resolution (can be enlarged 64 fold). I find that graphics passed through MS intermediary programs lose resolution. Illustrator can also convert single-page PostScript documents (most of the time, I have seen some instrument parts diagrams with a large number of crazy loopy lines). PS documents can also be converted with Adobe Acrobat (full version). Gerard. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Donald Ingram Sent: Thursday, April 07, 2005 19:13 To: r-help@stat.math.ethz.ch Subject: Re: [R] off-topic question: Latex and R in industries Hi Bert and Jonathan, When I want a quality report - I write it with pdfLaTeX ( TexShop or TeXnicCenter) with postscript generated diagrams and R plots as pdf's - ( so I can use PC / UNIX / OS X inter-changeably with no problems ) The quality and readability of the pdf document is liked but, and it's a big but is . When someone else in the team needs to extract quality vector graphics from the report, I have to give it to them in powerpoint or word document , which means running R again on a PC to get WMF's. Not impossible just extra work. ( Is there a universal vector format I could use ? ) However, and this is probably off topic-R, when I use drawings / schematics in native postscript from a Unix box, using them is fine in LaTeX, but they can't be pasted into MS applications without first rasterizing. The other option I tried - Ghostview seems to mess up line angles and fonts in attempting conversion into WMF. ( If anyone knows a way to avoid this, I will be forever grateful ) My problems - are not R but with general UNIX - PC interoperability Thanks for the nsf links - it's good to see Latex accepted, I also think the IEEE takes LaTeX, but for the business world it's Word only. Donald On 7 Apr 2005, at 22:56, Jonathan Baron wrote: On 04/07/05 22:46, Donald Ingram wrote: However LaTeX generated pdfs sent out as reports are much disliked. Really? I don't have this problem. It may have something to do with how you make them. With TeTeX, I use either pdflatex or dvips followed by dvipdfm. The latter is required when I have figures in eps. (ps2pdf is BAD.) I believe that these meet the standards of NSF (http://www.fastlane.nsf.gov). Unfortunately, https://www.fastlane.nsf.gov/servlet/faq.Faq; jsessionid=a8301381731112910739147?areaIndex=3faqIndex=12 now recommends that you just send the dvi file. They have given up on the possibility of users getting it right, but I think this is what they do. But all my papers on http://papers.ssrn.com are done this way. Jon -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Home page: http://www.sas.upenn.edu/~baron On 7 Apr 2005, at 22:56, Berton Gunter wrote: ?? R and MS coexist quite nicely. I frequently import R graphics as .wmf's into e.g. Word and Powerpoint. So I don't understand your remarks. Of course, there's no question about R's superiority for data analysis, graphs, etc. from any MS product. Incidentally, it is possible to use R via DCOM to generate data analyses and plots within Excel -- I don't know enough to be able to do this myself, but I know it can be done. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician is to catalyze the scientific learning process. - George E. P. Box -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Donald Ingram Sent: Thursday, April 07, 2005 2:46 PM To: r-help@stat.math.ethz.ch Subject: Re: [R] off-topic question: Latex and R in industries Wensui, I work for 'A' electronics test equipment corporation. I have been using R ( since 1.6 ) instead of MATLAB etc. as a general language for data analysis and graph generation. On they way to R I tried Python/Scipy, Scilab and others - but R wins in quality and ease of use (it just needs DSP and GPIB/HPIB libraries to be perfect ). LaTeX is also my document tool of choice .. However LaTeX generated pdfs sent out as reports are much disliked. MS Word, PowerPoint and Excel are the standards, and very importantly they offer cut and paste ability across the larger team. MS's offerings comes no where near to the quality of LaTex / R, but in world of shared authorship - it's a one sided battle. My other PC universe vs Unix/OS X problem is vector / Meta-file graphics - essential for quality reports. Postscript, PDF and MS products just don't play. The newest Office and Visio versions seem to be dropping even more of the postscript import and export filters ( which never work very well anyway ). I have never met any other colleagues who use LaTeX or R. Any one else sharing the same experiences ? Message: 37 Date: Wed, 6 Apr
[R] Error in save.image
Hi, I just came across an error that I haven't seen before and hope someone can help me. When I try to save my current workspace (using save.image or on quitting) I get the error message Error in save.image() : recursive default argument reference Does anyone know what is going on and how I can quit R without losing the contents of my current workspace? I am running R2.0.1 on a windows XP Pro platform. Thanks, Angelo -- -- | Angelo J. CantyEmail: [EMAIL PROTECTED] | | Mathematics and Statistics Phone: (905) 525-9140 x 27079 | | McMaster UniversityFax : (905) 522-0935 | | 1280 Main St. W. | | Hamilton ON L8S 4K1 | __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] HTML Help Browser in R Mac OS X Aqua GUI
I'm using R 2.0.1 with the Aqua R GUI 1.0 for Mac OS X, and I would like very much to use a firefox browser window for viewing help topics. options(htmlhelp) = TRUE options(browser) = /Applications/Connections/Firefox.app/Contents/MacOS/firefox-bin Java Embedding Plugin 0.9.0 is installed (the Java Embedding Plugin (JavaEmbeddingPlugin.bundle) and the MRJ Plugin JEP (MRJPlugin.plugin), are in the /Library/Internet Plug-Ins folder, and MRJ Plugin's timestamp is more recent than the Java Embedding Plugin's timestamp) help.start() launches firefox and displays the initial html help page. however, the following error message is displayed: /Applications/Connections/Firefox.app/Contents/MacOS/firefox-bin: can't map file: /Library/Internet Plug-Ins/MRJPlugin.plugin ((os/kern) invalid argument) subsequent calls to help in the form ?help.topic do not open html help documentation for help.topic. instead, the documentation is displayed in the internal help browser for the Aqua GUI. has anyone encountered this problem and found a solution? thanks Tony __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] error in save.image (addendum)
Some further remarks: The problem is not in executing save.image as I get the same error when I try to simply print out the function. I managed to quit and save my workspace using save(list = ls(all=TRUE), file = .RData); q(no) On restarting R with the same workspace save.image works fine. Even though I have managed to solve the problem for now, I would still be interested in knowing why it happened so that I can avoid whatever it was I did! Thanks, Angelo -- -- | Angelo J. CantyEmail: [EMAIL PROTECTED] | | Mathematics and Statistics Phone: (905) 525-9140 x 27079 | | McMaster UniversityFax : (905) 522-0935 | | 1280 Main St. W. | | Hamilton ON L8S 4K1 | __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] half-normal residual plots
Dear Malcolm, I don't think that anyone fielded this question earlier today: see the halfnorm function in the faraway package. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of MJ Price, Social Medicine Sent: Thursday, April 07, 2005 8:43 AM To: r-help@stat.math.ethz.ch Subject: [R] half-normal residual plots Hi all, I am trying to produce a half-normal plot of residuals from a GLM. I have found the qqnorm function for producing a normal plot but can't figure out how to produce a half-normal. Can anyone help with this? Thanks Malcolm -- MJ Price, Social Medicine [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] axis colors in pairs plot
The following command produces red axis line in a pairs plot: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c(red, green3, blue)[unclass(iris$Species)]) Trying to fool pairs in the following way produces the same plot as above: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species,pch = +, col = c(black, red, green3, blue)[ 1+ unclass(iris$Species)]) One very kludgy work-around is to define a new level 1, say foo in the first row of iris: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[1]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) However, if any other row is redefined, the red-axis persists. For example: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[3]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) I'd appreciate suggestions for a simpler work-around. Thanks, Anne __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] NA in table with integer types
Hi -- I am having the following problem with table() when applied to vectors of type (mode) integer. When I use the table() command, I can *only obtain an entry in the table for NA values by using exclude=NULL*. Just issuing exclude=NaN will not do it. See below, where x is double at first, and then coerced to integer and notice the difference. Is this a bug or is there something that I do not understand about the integer data type? That is, is there some other value besides NA and NaN that missing integers take? Thanks -- pr -- x - c(1,2,3,3,NA) is.double(x) [1] TRUE table(x,exclude=NA) x 1 2 3 1 1 2 table(x,exclude=NaN) x 123 NA 1121 table(x,exclude=NULL) x 123 NA 1121 x - as.integer(x) x [1] 1 2 3 3 NA is.na(x) [1] FALSE FALSE FALSE FALSE TRUE is.integer(x) [1] TRUE table(x,exclude=NA) x 1 2 3 1 1 2 table(x,exclude=NaN) x 1 2 3 1 1 2 table(x,exclude=NULL) x 123 NA 1121 R.version _ platform powerpc-apple-darwin6.8 arch powerpc os darwin6.8 system powerpc, darwin6.8 status major2 minor0.1 year 2004 month11 day 15 language R -- == Paul Rathouz, Assoc. Professor ph 773-834-1970 Dept. of Health Studies, Rm. W-264 fax 773-702-1979 University of Chicago[EMAIL PROTECTED] 5841 S. Maryland Ave. MC 2007 Chicago, IL 60637 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] axis colors in pairs plot
On Thursday 07 April 2005 17:51, Anne York wrote: The following command produces red axis line in a pairs plot: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c(red, green3, blue)[unclass(iris$Species)]) Trying to fool pairs in the following way produces the same plot as above: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species,pch = +, col = c(black, red, green3, blue)[ 1+ unclass(iris$Species)]) One very kludgy work-around is to define a new level 1, say foo in the first row of iris: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[1]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) However, if any other row is redefined, the red-axis persists. For example: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[3]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) I'd appreciate suggestions for a simpler work-around. One possibility is something along the lines of pairs(iris[1:4], panel = function(...) points(..., col = c(red, green3, blue) [unclass(iris$Species)] )) Deepayan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] axis colors in pairs plot
Hi Anne, Here's one suggestion, use a simple panel function: cols - c(red, green3, blue) with(iris, pairs(iris[, -5], main = Andersons Iris Data - 3 species, panel = function(x, y, ...) points(x, y, pch = (2:4)[Species], col = cols[Species], ...) )) Bill Venables -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Anne York Sent: Friday, 8 April 2005 8:51 AM To: Help R Subject: [R] axis colors in pairs plot The following command produces red axis line in a pairs plot: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c(red, green3, blue)[unclass(iris$Species)]) Trying to fool pairs in the following way produces the same plot as above: pairs(iris[1:4], main = Anderson's Iris Data -- 3 species,pch = +, col = c(black, red, green3, blue)[ 1+ unclass(iris$Species)]) One very kludgy work-around is to define a new level 1, say foo in the first row of iris: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[1]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) However, if any other row is redefined, the red-axis persists. For example: iris2=iris iris2$Species = as.character(iris2$Species) iris2$Species[3]=foo iris2$Species = factor(iris2$Species) pairs(iris2[1:4], main = Anderson's Iris Data -- 3 species, pch = +, col = c( black,red, green3,blue)[ unclass(iris2$Species)]) I'd appreciate suggestions for a simpler work-around. Thanks, Anne __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html