Re: [R] *not* using attach() *but* in one case ....

2011-05-19 Thread Martin Maechler
over load() there is that loaded objects (and there maye be a bunch!), are put into a separate place in the search path and will not accidentally overwrite objects in the global "workspace". Of course, there are still quite a few situations {e.g. in typical BATCH use of R for simulations

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

2011-05-19 Thread Martin Maechler
s, programming (R or otherwise) is considered a creative activity to quite some extent, and creativity can be crushed by too rigid rules. I'd state that cultural history shows that human culture implicitly follows many rules, but it is (almost) only interesting because the are enough exce

Re: [R] pam() seems to ignore cluster number

2011-05-23 Thread Martin Maechler
ng performed, NAs in the computed dissimilarity matrix." But because of 'cluster.only=TRUE' *and* because of a lapsus of the 'cluster' maintainer (me), pam() returned without the error message in this case. The next release of R (or of 'cluster') will give the

Re: [R] NaN, Inf to NA

2011-05-27 Thread Martin Maechler
> Duncan Murdoch > on Fri, 27 May 2011 08:23:14 -0400 writes: > On 11-05-27 4:27 AM, Albert-Jan Roskam wrote: >> Aha! Thank you very much for that clarification! It would >> be much more user friendly if R generated a >> NotImplementedError or something similar. The 'g

Re: [R] 'matplot' for matrix with NAs: broken lines

2010-05-06 Thread Martin Maechler
orate, using your example (below)? I found Gabor's answer appropriate, I really cannot see why matplot() should behave differently here... Martin Maechler TS> >> From: ggrothendi...@gmail.com >> Date: Wed, 5 Ma

Re: [R] Cluster procedure using geographical neighborhood

2010-05-07 Thread Martin Maechler
ties, 'w' is weight factor and d~(u,v) is e.g. the geodesic distance between u and v. I'm CC'ing this to the R-help mailing list, as I think you could get more advice from there. Martin Maechler, ETH Zurich DS> However if I do that the procedure hclust () does not wo

[R] Accidental misuse of 'R-packages' list

2010-05-18 Thread Martin Maechler
good rule is to not mention a pacakge more than once a year, and much less than that for all packages with only minor changes. Martin Maechler, ETH Zurich >>>>> "MY" == Matt Young >>>>> on Wed, 12 May 2010 08:23:24 -0700 writes: MY> Lot of

Re: [R] lme4 + R 2.11.0 + mac unavailable

2010-05-18 Thread Martin Maechler
stall into the (or "a") library where Matrix is installed. For those who only use one library, this will not be a concern. Others could use the fact that MacOSX is a kind of Unix and you can simply "symlink" Matrix (from the R standard library) into the library into which you want to

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread Martin Maechler
IST. These are *PACKAGES* ! Why do you think are you talking about the function install.packages() --- To answer the question you did want to ask: Do not be afraid: Depedencies are only installed when needed, i.e., no package will be downloaded and installed if it alread

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread Martin Maechler
ill not be reinstalled. > > cheers > > milton > > On Tue, May 18, 2010 at 12:49 PM, Martin Maechler < > maech...@stat.math.ethz.ch> wrote: > >> { I've modified the subject; I can't stand it hitting square into >> my face ... } >> >

Re: [R] lme4 + R 2.11.0 + mac unavailable

2010-05-19 Thread Martin Maechler
Just a small correction to what I've said ... see below >>>>> "MM" == Martin Maechler >>>>> on Tue, 18 May 2010 10:00:20 +0200 writes: >>>>> David Winsemius >>>>> on Mon, 17 May 2010 17:44:00 -0400 writes:

[R] R-help "spam" detection; please help the moderators

2010-06-01 Thread Martin Maechler
h digits. We hereby ask those among you who use a freemail account to please no longer post via nabble. Thank you for your support of R-help, *the* "community mailing list" of the R project since even before that project existed "formally", namely since 1997-04-01, today 13 ye

Re: [R] agnes clustering and NAs

2011-01-28 Thread Martin Maechler
23371 dissimilarity matrix, stored as a double precision matrix (on a 64-bit platform) that's an object of size 4.4 GBytes, not very convenient to work with. as dissimilarity object it will only be about half of that size, but that's still ``a bit large''.. As I said above, for such

Re: [R] Positive Definite Matrix

2011-01-29 Thread Martin Maechler
> Prof Brian Ripley > on Sat, 29 Jan 2011 21:00:19 + (GMT) writes: > On Sat, 29 Jan 2011, David Winsemius wrote: >> >> On Jan 29, 2011, at 12:17 PM, Prof Brian Ripley wrote: >> >>> On Sat, 29 Jan 2011, David Winsemius wrote: >>> On Jan

Re: [R] Sparse matrix multiplication

2011-01-29 Thread Martin Maechler
rnorm(1e4), 100,100)) system.time(B. <- .bdiag(Blis)) # 1.28 sec system.time(cc <- A %*% B.) # 1.7 sec class(cc)# "dgeMatrix" .. i.e. dense ## and depending on the context you may revert to traditional unclassed matrices ## via c2 <- as(cc, "matrix") >> T

Re: [R] Positive Definite Matrix

2011-01-31 Thread Martin Maechler
>I think the bottom line can be summarized as > follows: > 1. Give up on Cholesky factors unless you have a > matrix you know must be symmetric and strictly positive > definite. (I seem to recall having had problems with chol > even with matrices that were theor

Re: [R] Meaning of pterms in survreg object?

2011-01-31 Thread Martin Maechler
reference. The reference to see what's posted or not are the mailing list archives "proper", https://stat.ethz.ch/pipermail/r-help/ , and they are available directly from the https://stat.ethz.ch/mailman/listinfo/r-help url at the bottom of every posting ... Martin Maechler, ETH

Re: [R] Preparing dataset for glmnet: factors to dummies

2011-02-01 Thread Martin Maechler
;s model.matrix() (as that produces only standard, i.e., dense matrices). The functionality of model.Matrix() has been moved out of the Matrix package into the package 'MatrixModels', and that package also provides -- somewhat experimental -- functionality for fitting GLMs with sparse mode

Re: [R] better way to iterate matrix multiplication?

2011-02-01 Thread Martin Maechler
ly when the matrix is non-singular, whereas the `%^%` in package expm should work reliably in all cases. Also for (much) smaller powers than n=10'000 the cpu time needed is more comparable. Martin Maechler, ETH Zurich EV> Am 01.02.2011 17:16, schrieb Ben Bolker: >>

Re: [R] "hubers" function in R MASS library : problem and solution

2011-02-04 Thread Martin Maechler
stimate no-longer-scale-equivariant. huberM() *has* an s argument for specifying the scale estimate, so you could use it as huberM(a, s = max(mad(a), 1e-6)) if you want. Note that your sample 'a' is constructed in a way that all scale-equivariant 50%-breakpoint robust estimate

Re: [R] Revolution Analytics reading SAS datasets

2011-02-11 Thread Martin Maechler
on R-help, please. At first, note that R is GNU R, and R-help has been about R as a Free ("Libre") Software, for all its many years and hundreds of thousands of messages. Revolutions's product may be fine for some, in some situations, but supporting non-Free parts of it, really

Re: [R] Console output

2011-02-21 Thread Martin Maechler
> Ted Harding > on Mon, 21 Feb 2011 11:08:19 - (GMT) writes: > That doesn't produce quite what Antje asked for (since each > line gets number "[1]"). The following does work: > print(cbind(NULL,(1:10))) > [,1] > [1,]1 > [2,]2 > [3,]3 >

Re: [R] Removing elements ..., NO bug in which() function

2011-02-23 Thread Martin Maechler
> "RK" == Rumen Kostadinov > on Sun, 13 Feb 2011 12:46:52 -0500 writes: RK> Thanks Sarah, RK> Yes, the function behaves Exactly as documented: RK> check this out: >> a = c(1,2,3,4,5) >> a[which(a!=6)] RK> [1] 1 2 3 4 5 >> a[!which(a==6)] RK> numeric(0)

Re: [R] The L Word

2011-02-24 Thread Martin Maechler
> "CB" == Claudia Beleites > on Thu, 24 Feb 2011 12:31:55 +0100 writes: CB> On 02/24/2011 11:20 AM, Prof Brian Ripley wrote: >> On Thu, 24 Feb 2011, Tal Galili wrote: >> >>> Thank you all for the answers. >>> >>> So if I may extend on the question - >>> W

Re: [R] extract printed value from a function

2011-02-24 Thread Martin Maechler
> "PE" == Peter Ehlers > on Thu, 24 Feb 2011 07:07:29 -0800 writes: PE> On 2011-02-24 06:32, David Winsemius wrote: >> >> On Feb 24, 2011, at 7:22 AM, Peter Ehlers wrote: >> >>> On 2011-02-24 03:26, Duarte Viana wrote: Hello all, This

Re: [R] The L Word

2011-02-24 Thread Martin Maechler
> "HW" == Hadley Wickham > on Thu, 24 Feb 2011 10:14:35 -0600 writes: >> Note however that I've never seen evidence for a *practical* >> difference in simple cases, and also of such cases as part of a >> larger computation. >> But I'm happy to see one if anyone has an

Re: [R] The L Word

2011-02-24 Thread Martin Maechler
>>>>> "MM" == Martin Maechler >>>>> on Thu, 24 Feb 2011 18:34:36 +0100 writes: >>>>> "HW" == Hadley Wickham >>>>> on Thu, 24 Feb 2011 10:14:35 -0600 writes: >>> Note however that I've ne

Re: [R] lty=NULL crashing R for x11(type="cairo")

2011-03-05 Thread Martin Maechler
> "IZ" == Ista Zahn > on Sat, 5 Mar 2011 14:07:04 + writes: IZ> I confirm this bug exists and is 100% replicable on R IZ> version 2.12.2 (2011-02-25) Platform: i686-pc-linux-gnu IZ> (32-bit) WHoa... debugging ===> it *is* a bug in R after all : > plot(1); axi

Re: [R] Teaching R: To quote, or not to quote?

2011-03-10 Thread Martin Maechler
r the fortunes package! well, yes, in particular as {R code} > maintainer("fortunes") [1] "Achim Zeileis " and BTW, Please send only *once* to R-help, and that nowadays prefers the @r-project.org address. Martin Maechler, ETH Zurich

Re: [R] Getting the source file's name where the custom function is written

2011-03-10 Thread Martin Maechler
ed when they are > installed (and is almost certainly false). and "almost certainly" is really correct and applies to "almost all" users, where my Bayesian a priori estimate of that proportion would be '> 0.'. Martin Maechler > Duncan Murd

Re: [R] Persistent storage between package invocations

2011-03-18 Thread Martin Maechler
> It would be nice to have a standard directory where R can > write things this way. A semi-standard directory is given > by Sys.getenv("R_LIBS_USER"), which defaults to ~/R/.../. > Maybe ~/R/ could serve as that convention? That way we > (various developers etc) would also n

Re: [R] RmetricsTools.R

2011-04-04 Thread Martin Maechler
It's "top-level". With the standard subversion client, you can get it, e.g., by svn export svn://svn.r-forge.r-project.org/svnroot/rmetrics/pkg/RmetricsTools.R -- Martin Maechler, ETH Zurich __ R-help@r-project.org mailing list htt

Re: [R] [R-pkgs] Boost library in R package.

2011-04-23 Thread Martin Maechler
Wrong mailing list!! This got posted accidentally to R-packages -- which is exclusively for (semi-)important announcements by package *authors*. Please do only reply to R-help --- and do use R-help for *ALL* questions about R packages (unless you ask the package maintainer directly ..). Martin

Re: [R] 3D Bar Graphs in ggplot2?

2011-08-04 Thread Martin Maechler
a >> 2D bar graph with a 3D shaped bard. See attached excel file >> for an example. >> >> Before anyone asks I know that 3D looking bars don't add >> anything except "prettiness". > If you want graphs like that, you sho

Re: [R] Possible bug of QR decomposition in package Matrix

2011-08-04 Thread Martin Maechler
> C6H5NO2 > on Thu, 4 Aug 2011 16:03:54 +0800 writes: > Thank you very much, Josh! > As you suggested, I will contact the developers of "Matrix". > PS, C6 are just initial characters of my email account :-) > Best wishes, > C6 well, as the posting guide (http:/

Re: [R] Can glmnet handle models with numeric and categorical data?

2011-08-05 Thread Martin Maechler
27;Matrix' as part of your R installation). For exactly this reason, 'glmnet' has supported the use of sparse matrices for a long time, and we have provided the convenience function sparse.model.matrix() {package 'Matrix'} for easy construction of such matrices.

[R] [R-pkgs] New Matrix and lme4: Must reinstall lme4 if got new Matrix

2011-08-10 Thread Martin Maechler
or possibly R-devel@...), so we can prepare a shiny sparkling Matrix_1.0-0 in time. Thank you in advance for such a contribution to the Free Software universe. Martin Maechler and Doug Bates (and Ben Bolker for lme4). -- -- -- - Why must lme4 be re-installed as well ? [Technical !] This is

Re: [R] round() a data frame containing 'character' variables?

2011-08-11 Thread Martin Maechler
> Liviu Andronic > on Wed, 10 Aug 2011 11:46:55 +0200 writes: > Hello On Wed, Aug 10, 2011 at 11:41 AM, Dimitris > Rizopoulos wrote: >> One approach is the following: >> >> numVars <- sapply(iris, is.numeric) iris[numVars] <- >> lapply(iris[numVars], round, d

Re: [R] question on silhouette colours

2011-08-29 Thread Martin Maechler
a matter of fact, I'm pretty sure you found a bug. Note that it would be better in such cases (a function in an R package) to first contact the package maintainer, in this case > maintainer("cluster") [1] "Martin Maechler " but I did see your message on R-help

Re: [R] The elegant way to test if a number is a whole number

2011-09-09 Thread Martin Maechler
; >>> On Sep 8, 2011, at 2:42 PM, David Winsemius wrote: >>> >>>> >>>> On Sep 8, 2011, at 3:35 PM, Alexander Engelhardt wrote: >>>> >>>>> Am 08.09.2011 20:48, schrieb Marc Schwartz: >>>>> There was a post from Mar

Re: [R] dgamma in jags within r

2011-09-10 Thread Martin Maechler
> Marc Girondot > on Sat, 10 Sep 2011 11:43:20 +0200 writes: > Marc Girondot > on Sat, 10 Sep 2011 11:43:20 +0200 writes: > I define priors in jags within r using a gamma distribution. I would > like to control the shape but I have problems. Any help will be usef

Re: [R] Very slow using S4 classes

2011-09-10 Thread Martin Maechler
> Martin Morgan > on Sat, 10 Sep 2011 10:18:11 -0700 writes: > On 09/10/2011 08:08 AM, André Rossi wrote: >> Hi everybody! >> >> I'm creating an object of a S4 class that has two slots: ListExamples, >> which is a list, and idx, which is an integer (as the code be

Re: [R] Boxplot with Log10 and base-exponent axis

2012-06-22 Thread Martin Maechler
ackage "sfsmisc" does (something like) this (and a bit more) even more nicely : install.packages("sfsmisc") require("sfsmisc") x <- list(One=10^(sin(1:10)+5), Two=10^(cos(1:30)*2)) boxplot(x, log="y", yaxt="n") eaxis(2) -- Martin Maechler,

Re: [R] scatter.smooth() line colour

2012-06-22 Thread Martin Maechler
nd to use do.call(lines, *) As a little exercise for you to see, I'll add that to R-devel .. even though I wonder if its worth the extra work {{but then, as otherwise someone gets the idea to file it in the bug tracker ..}} soon to be visible at https://svn.r-project.org/R/trunk/src/libra

Re: [R] R code for to check outliers

2012-07-18 Thread Martin Maechler
> Bert Gunter > on Wed, 18 Jul 2012 07:14:31 -0700 writes: > checkforoutliers <- function(series) NULL > Cheers, Bert > *Explanation: There is no such thing as a statistical > outlier -- or, rather,"outlier" is a fraudulent > statistical concept, defined arbitrari

Re: [R] computer name

2011-06-14 Thread Martin Maechler
deed! BTW: Learn something more by considering the slight variation Sys.info()[['nodename']] which is what you typically really want (and which will be very very slightly more efficient in typical use cases). Martin Maechler, ETH Zurich > On 6/13/2011 7:38 AM, David Scott wr

Re: [R] Tricky (?) conversion from data.frame to matrix where not all pairs exist

2011-06-22 Thread Martin Maechler
> Marius Hofert > on Wed, 22 Jun 2011 01:00:21 +0200 writes: > Thanks a lot! Works like charm :-))) Cheers, well, I'd still *strongly* vote for Bill Dunlap's version. Martin > Marius > On 2011-06-22, at 24:51 , Dennis Murphy wrote: >> Ahhh...you want a matrix. xt

Re: [R] working with sparse matrix

2011-06-22 Thread Martin Maechler
rspMatrix() functions which the Matrix package provides. >> M2 <- as.matrix(M) Why would you have to make your sparse matrix dense? (It won't work for much larger matrices anyway: They can't exist in your RAM as dense ..)... Yes, I see you work with further code that do

Re: [R] ?plot: Add an example on how to plot functions to the help of `plot`.

2011-07-28 Thread Martin Maechler
he help page would start with > Generic X-Y Plotting > > Description: > > Generic function for plotting of R objects. For more details > about the graphical parameter arguments, see ‘par’. > > For simple scatter plots, ‘plot.default’ will be used. However, >

Re: [R] bug in dev.copy2pdf output?

2011-07-29 Thread Martin Maechler
e bug that has been on R's ?pdf help page for years [] --> search for "q" (including the quotes). Why oh why oh why ... are people always thinkg of bugs in R when they notice something and why oh why oh why do they then "google around" instead of just once care

Re: [R] Communative Matrix Multiplcation

2012-08-16 Thread Martin Maechler
should) update packages also if for some reason you cannot update your version of R! Martin Maechler, ETH Zurich > -Original Message- From: > r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of > Doran, Harold Sent: Tuesda

Re: [R] reshape cast to a sparse matrix?

2012-08-16 Thread Martin Maechler
mostly advertize sparseMatrix() which has a large 'Examples' section. If you want, you can contact me offline for more help, but you'd definitely have to provide more details, and a selfcontained example. Martin Maechler, ETH Zurich, one of the two 'Matrix' author

Re: [R] Communative Matrix Multiplcation

2012-08-16 Thread Martin Maechler
out this thread! > {in addition to the fact that he should have used all.equal() > and not identical() as someone had pointed out already ..} > And at last: Yes, some of the methods for Diagonal() matrices in the > 'Matrix' package did have bugs,

Re: [R] using apply with sparse matrix from package Matrix

2012-09-04 Thread Martin Maechler
>> my.sum<-apply(sm, 1, sum) > Error in asMethod(object) : > Cholmod error 'problem too large' at file ../Core/cholmod_dense.c, line 106 So, actually it would have worked (though not efficiently) if your sm matrix would have been much smaller. However,

Re: [R] arrow egdes and lty

2011-11-14 Thread Martin Maechler
> Ted Harding > on Mon, 14 Nov 2011 10:39:35 + writes: > On 14-Nov-11 09:57:52, Matthias Gondan wrote: >> Dear R developers, >> >> I want to draw an arrow in a figure with lty=2. The lty >> argument also affects the edge of the arrow, which is >> unfortuna

Re: [R] How to compute eigenvectors and eigenvalues?

2011-11-14 Thread Martin Maechler
,] -3.743456e+07 9.358640e+06 9.358641e+06 it isn't... > My question is: how can I fix it? No need to fix anything, as nothing is broken ;-) Martin Maechler, ETH Zurich __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] modify the name of axis of an R function

2011-12-16 Thread Martin Maechler
disappointed. Not at all much more convoluted, just different. I think you should make a 2nd attempt, starting with things like showMethods(plot) with versions showMethods(plot, class = "...") showMethods(plot, include = TRUE) showMethods(plot, class = "...

Re: [R] nlrob problem

2011-12-23 Thread Martin Maechler
‘weights’. Can also be a list or an >environment, but not a matrix. > I think it would be good to mention in the documentation of 'nls' that > weights can be passed as 'w' as well. You have overlooked that nlrob uses a different *formula* for nls() --- nowadays

Re: [R] stop R from rounding

2011-10-20 Thread Martin Maechler
> David Winsemius > on Thu, 20 Oct 2011 01:51:28 -0400 writes: > On Oct 19, 2011, at 11:29 PM, Alyse wrote: >> Hello, >> >> I have a column in a data frame that need to be 10 digits long. As >> such: >> >> Decimal.Year >> 1 1994.25997 >> 2 1

Re: [R] Building package/DESCRIPTION file not existing?

2011-10-26 Thread Martin Maechler
ntly taught a course on R Package building and on Windows, the user had problems because of an 'ä' (a-Umlaut) in one of the directories in her 'path'. So if you work from another place than 'C:/Users/Propri�taire

Re: [R] Call the Standard Error and t-test probability in linear regression

2012-03-02 Thread Martin Maechler
> On 02-Mar-2012 IOANNA wrote: > > Hello, > > I run a linear regression I get the summary, e.g.: > > > >> > >Call: > >lm(formula = signal ~ conc) > >Residuals: > > 12 3 4 5 > > 0.4 -1.0 1.6 -1.80.8 > >Coefficients: > >

Re: [R] max.print

2012-03-12 Thread Martin Maechler
e user more clearly! I propose to implement that it will produce an *error*, if the result is not a positive integer. Another question: Would anyone (on R-core ?) know why we print two closing "]]" but only one initial "[" in the following ? [ reached getOption("max.

Re: [R] what does "rlm" do if it fails to converge within iteration limits?

2012-03-13 Thread Martin Maechler
da after save(., file="mydata.rda") for download 2) What's the result of using lmrob() {package 'robustbase'} instead of rlm() {package 'MASS'} ? Best regards, Martin Maechler, ETH Zurich > On Fri, Mar 9, 2012 at 1:21 PM, Berend Hasse

Re: [R] Urgente: retiro de la lista

2012-03-16 Thread Martin Maechler
> "JB" == Jaime Bernal-Hadad > on Thu, 15 Mar 2012 19:11:19 -0700 writes: JB> Urgente JB> Dedido al gran n�mero de correos recibidos que afectan seriamente mi trabajo, es urgente que me reitren de esta lista. JB> Alguien me puede decir que debo hacer JB> �

Re: [R] Problems with CRAN package Ryacas

2012-04-09 Thread Martin Maechler
Apropos: I don't have the problems, the OP had, but on my ubuntu notebook, Ryacas does not return expressions (just the strings), and hence as.expression( ) always gives NULL and e.g. the demo(Ryacas-Function) also fails: > yacas(expression(deriv(BurrCDF(x,c,k k*c*x^(c-1)*(x^c+1)^

Re: [R] Schedule R script using cron

2012-04-12 Thread Martin Maechler
> winie > on Wed, 11 Apr 2012 18:51:28 -0700 writes: > I am trying to schedule my R script using cron, but it is not > working. It seems R can not find packages in cron. Anyone can help me? > Thanks. > The following is my bash script > # source my porfile >

[R] #!/usr/bin/env Rscript --vanilla ??

2012-04-13 Thread Martin Maechler
#! /usr/bin/env Rscript --vanilla the error message simply is /usr/bin/env: Rscript --vanilla: No such file or directory I have tried a few variations on the theme, using quotes in different places, but have not succeeded till now. Any suggestions? Martin Maechler, ETH Zurich __

Re: [R] Problems with plotCI

2012-01-14 Thread Martin Maechler
pts text/plain If your e-mail software has changed, and no longer uses (or allows you to use) text/plain for plain text such as R source code, then you should blame the provider of your e-mail software ... or alas provide text inline, as you did. Martin Maechler, ETH Zurich (and R-help &q

Re: [R] fUtilities removed --> use fBasics

2012-01-14 Thread Martin Maechler
one of rmetrics-c...@r-project.org, as having been an active co-maintainer of some parts of the Rmetrics collection, but I have not been part of the "management" nor the foundation. Martin Maechler, ETH Zurich > You need to send your questions to the maintainers. They > ap

Re: [R] R equivalent of Python str()?

2012-02-08 Thread Martin Maechler
> "R Michael Weylandt " > > on Tue, 7 Feb 2012 20:23:57 -0500 writes: > Possibly as.character() is what the OP was seeking > Michael or format() which is closer for numeric data > On Feb 7, 2012, at 7:15 PM, jim holtman wrote: >> ?dump >> ?dput >>

Re: [R] Huge difference btw system time and elapsed time

2012-02-21 Thread Martin Maechler
ystem.time()[[1]] etc That's typically good enough [as long as you don't use internally multithreaded BLAS, or parallelize / distribute your job, or other things you won't be doing (yet)]. Martin Maechler, ETH Zurich > I also tried to improve the code by using 'cmp

Re: [R] How to "flatten" a multidimensional array into a dataframe?

2012-04-20 Thread Martin Maechler
> baptiste auguie > on Fri, 20 Apr 2012 09:01:35 +1200 writes: > library(plyr) > adply(my.array,1:3) Hmm, as Marc Schwartz has already said: There's the base R as.data.frame() method for "table" objects which also works for arrays, so as.data.frame(as.table(my.a

Re: [R] PCA sensitive to outliers?

2012-04-23 Thread Martin Maechler
lable in an R package is cov.rob() from MASS ('Recommended' package available with every R). *And* you can use standard R's princomp(x, ... , covmat = (x)) to get robust PCA. I'll add a note with that to the 'Robust' CRAN task view. Martin Maechler, ETH Zurich

Re: [R] Selecting columns whose names contain "mutated" except when they also contain "non" or "un"

2012-04-27 Thread Martin Maechler
> David Winsemius > on Mon, 23 Apr 2012 12:16:39 -0400 writes: > On Apr 23, 2012, at 12:10 PM, Paul Miller wrote: >> Hello All, >> >> Started out awhile ago trying to select columns in a >> dataframe whose names contain some variation of the word >> "mutant"

Re: [R] masked by GlobalEnv ???

2012-05-12 Thread Martin Maechler
nce the latter overwrites objects in .GlobalEnv, I do recommend using the former in this case. Martin Maechler ETH Zurich and R core Team > On Sat, May 12, 2012 at 5:54 AM, Jhope > wrote: >> Hi R Listers, >> >> I am trying to upload a data file and I r

Re: [R] Job posting - Statistical Consultant ...

2012-05-21 Thread Martin Maechler
ael Sorry, but that's ABSOLUTELY *not* appropriate! {for one: R-help is international, so many job offers do not make sense anyway}. We have an R-SIG-jobs (slight misnomer, yes) mailing list --> http://stat.ethz.ch/mailman/listinfo/r-sig-jobs for job postings, if you really must use a

Re: [R] fitting t copula with fixed dof

2012-05-21 Thread Martin Maechler
ll appreciate any help/comments. I searched wrongly. The copula package has been able to do that for a very long time. --> fitCopula() Martin Maechler, ETH Zurich __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

Re: [R] Shared Library Error

2008-05-24 Thread Martin Maechler
my source code is written in f90. But in that case, I'd almost bet that .Fortran("merge_xtabs_patterns_file", ) will work. The idea being that you use the same name as in the fortran code, namely with*OUT* trailing "_". Regards, Martin Maechler, ETH Zurich

Re: [R] need help for building R in Ubuntu 8.04

2008-05-28 Thread Martin Maechler
r that it did not contain *really* all the tools I'd wanted, but almost all. e.g., you may also want the two packages tcl8.4-dev tk8.4-dev which probably are not in the dependencies of 'r-base-dev'. --- Martin Maechler, ETH Zurich EH> On Sun, May 25, 2008 at 1:20 PM

Re: [R] need help for building R in Ubuntu 8.04

2008-05-28 Thread Martin Maechler
> "DE" == Dirk Eddelbuettel <[EMAIL PROTECTED]> > on Wed, 28 May 2008 12:46:50 + writes: DE> On Wed, May 28, 2008 at 02:29:10PM +0200, Martin DE> Maechler wrote: >> > "EH" == Erin Hodgess <[EMAIL PROTECTED]> >> > on Sun, 25 May 2008 13:27:04 -0500 writes:

Re: [R] probleme with R CMD check

2008-06-03 Thread Martin Maechler
> "c" == cgenolin <[EMAIL PROTECTED]> > on Tue, 03 Jun 2008 00:32:13 +0200 writes: c> Hi the list. c> I got a probleme with R CMD check. I run a package.skeleton on a very c> simple code (four toy function, no probleme with the R CMD build). But c> on the R CMD check

Re: [R] probleme with R CMD check

2008-06-03 Thread Martin Maechler
> "c" == cgenolin <[EMAIL PROTECTED]> > on Tue, 03 Jun 2008 09:55:15 +0200 writes: >> provide commented, minimal, self-contained, reproducible code. >> Please do -- both -- c> Sorry... Usualy, I do... c> --- 8< c> userA <- function(x)x+1 c>

Re: [R] probleme with R CMD check

2008-06-03 Thread Martin Maechler
> "c" == cgenolin <[EMAIL PROTECTED]> > on Wed, 04 Jun 2008 00:30:08 +0200 writes: c> I find the problem: the name of my base directory (the c> one that contain R, man and test) was not the same than c> the name of the package... Aha! c> I thaught this was allowed, I

Re: [R] Bug/Error in formatC? (Was: Why doesn 't formatC( x, digits=2, format= " g")...)

2008-06-10 Thread Martin Maechler
", flag = "#") ends up calling sprintf("%#.2g", x) where the is carefully determined from x, in your example sprintf("%#3.2g", 0.005999) For me, this correctly gives > sprintf("%#3.2g", 0.005999) [1] "0.0060" and I assume that for

Re: [R] Bug/Error in formatC? ....

2008-06-10 Thread Martin Maechler
>>>>> "MM" == Martin Maechler <[EMAIL PROTECTED]> >>>>> on Tue, 10 Jun 2008 15:05:05 +0200 writes: >>>>> "PD" == Peter Dunn <[EMAIL PROTECTED]> >>>>> on Tue, 10 Jun 2008 09:25:07 +1000 writes:

Re: [R] Output of silhouette (cluster package)

2008-06-14 Thread Martin Maechler
ilhouette CV> where the fourth group is composed by units from 46 to 48 instead of CV> units from 58 to 60 (belonging to the third cluster), see You are right. Indeed, I see that this has been a bug in sortSilhouette() ever since that had been introduced in 2002. It will be fixed in c

Re: [R] R web site-Useability & finding varous bits of documentation

2008-06-19 Thread Martin Maechler
;> on CRAN <<< Note that the R-project page (www.r-project.org) and CRAN are two "things", albeit closely related. CRAN is for "DOWNLOAD"ing, including free contributed docs. So that is the main reason, "contrib.docs" are not there in the www.r-project.o

Re: [R] highest eigenvalues of a matrix

2008-06-19 Thread Martin Maechler
> "KateM" == Katharine Mullen <[EMAIL PROTECTED]> > on Thu, 19 Jun 2008 11:50:22 +0200 (CEST) writes: KateM> On Thu, 19 Jun 2008, Simon Wood wrote: >> > >> > I happily use eigen() to compute the eigenvalues and >> eigenvectors of > a fairly large matrix (200x200, say),

Re: [R] error related to approxfun in R 2.7.0

2008-06-21 Thread Martin Maechler
tats". If f <- approxfun(...) in R 2.0.1 and you save f , load it in R 2.7.0, you can do body(f)[[2]][["PACKAGE"]] <- "stats" and then use f(..) in R 2.7.0. Martin Maechler, ETH Zurich PS: Do I guess correctly, that you are Swiss? ___

Re: [R] transform a table into a matrix

2008-06-24 Thread Martin Maechler
JohnF> from your question where the "table" is located -- in JohnF> a file perhaps? If so, you can use read.table() to JohnF> read the data into a dataframe. Yes, and then possibly use data.matrix(.) to produce a numeric matrix where the "categorical variables

Re: [R] as.Date() clarification

2008-07-03 Thread Martin Maechler
p pages on the PD> web, but I can't remember how. http://stat.ethz.ch/R-manual/ (I think one of the oldest still working links in the R-project) where, for as.Date , as in our question, you may end up here: http://stat.ethz.ch/R-manual/R-patched/library/base/html/format.Date.h

Re: [R] Can't install package Matrix on solaris.

2008-07-03 Thread Martin Maechler
e nice feature that it *will* check if you are using GNU make and complain in English (instead of giving a pretty unhelpful error message later) if you don't use GNU make. Hoping, you'll enjoy "the Matrix", Martin Maechler, ETH Zurich BDR> On Wed, 2 Jul 2008, Scott Br

Re: [R] S4 : dumpMethod

2008-03-19 Thread Martin Maechler
bug that dumpMethod() and dumpMethods() are currently not working at all. Do you want to submit a (legitimate (!)) bug report to R-bugs ? Thank you for your reproducible example and subsequent questions! Martin Maechler, ETH Zurich __ R-

Re: [R] Inverse t-distribution(TINV function in Excel)

2008-03-25 Thread Martin Maechler
> "FC" == Felipe Carrillo <[EMAIL PROTECTED]> > on Mon, 24 Mar 2008 18:42:36 -0700 (PDT) writes: FC> Thanks everyone for your help the qt function worked FC> like a dream. I ended up using it like this; T <- FC> qt((1-0.05/2),7) to estimate 95% intervals divided by FC>

[R] [R-pkgs] "The Matrix" is approaching version 1.0-0

2008-03-28 Thread Martin Maechler
e many (mostly small) problems to announce that release 1.0-0 is imminent. In the DESCRIPTION of the package we say -- Title: A Matrix package for R Author: Douglas Bates <....> and Martin Maechler <> Maintainer: Doug and Martin

Re: [R] S4 : the list of all the object

2008-03-31 Thread Martin Maechler
> "CG" == Christophe Genolini <[EMAIL PROTECTED]> > on Sat, 29 Mar 2008 16:20:41 +0100 writes: CG> Hi the list, Is it possible to get the list of all the CG> S4 user define classes? do you mean those in .GlobalEnv? If you have written a package using S4 classes, these are als

Re: [R] hsv2rgb in R

2008-03-31 Thread Martin Maechler
> "DM" == Dieter Menne <[EMAIL PROTECTED]> > on Mon, 31 Mar 2008 08:06:39 + (UTC) writes: DM> I noted that there is a hsv2rgb in the C-API, but no DM> corresponding function in R, while rgb2hsv is available DM> in R. DM> Does the functionality of hsv2rgb hide under

Re: [R] Thinking about using two y-scales on your plot?

2008-04-02 Thread Martin Maechler
> "ASR" == Allen S Rout <[EMAIL PROTECTED]> > on Wed, 02 Apr 2008 10:39:31 -0400 writes: ASR> "hadley wickham" <[EMAIL PROTECTED]> writes: >> Please read this first: >> http://www.perceptualedge.com/articles/visual_business_intelligence/dual-scaled_axes.pdf >> >>

Re: [R] How to pack my stuff into a package (library, collection)?

2008-04-08 Thread Martin Maechler
think "Matlab" One R source file typically contains several function (and maybe class and method) definitions. Often on Rd file contains documentation for more than one object. Look at the R sources , at any package sources (i.e. not the *.zip, but the *.tar.gz files). Martin Maechler, ETH Z

[R] sign() in R ?

2008-04-15 Thread Martin Maechler
y, a random permutation is simply found by 'sample(n)'} Thank you in advance for further pointers, or even working R code. Best regards, Martin Maechler, ETH Zurich __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

<    2   3   4   5   6   7   8   9   >