.
I'll post R code later, being busy for the next few hours.
Martin
>>>>> "MM" == Martin Maechler <[EMAIL PROTECTED]>
>>>>> on Tue, 15 Apr 2008 18:13:43 +0200 writes:
MM> I am looking for an algorithm (written in R (preferably) or C,
d as signCycLength(cycLengths2(p))
e.g.,
p <- sample(50)
signCycLength(cycLengths2(p))
Thank you, again for all the feedback, and suggestions,
Martin Maechler, ETH Zurich
>>>>>> Peter Dalgaard wrote yesterday
>> Martin Maechler wrote:
>> > Ok,
>>
; (I would use R in batch mode but I find the startup time prohibitive,
unless
ML> there is a way to speed it up)
there definitely is, and I wonder why you haven't found it.
E.g.,
Rscript --vanilla --slave --default-packages=stats -e 'pnorm(1.96)'
[1] 0.9750021
i
Hmm,
> "KeBe" == Beck, Kenneth (STP) <[EMAIL PROTECTED]>
> on Thu, 24 Apr 2008 10:12:19 -0500 writes:
KeBe> OK I've spent a lot of time with the core
KeBe> documentation, and I never found anything as simple as
KeBe> their table 2.1, which elucidated the difference
KeB
t understand why you want to throw away the dimnames of
the resulting numeric matrix. Rather, in general you'd want
to keep them.
Martin Maechler, ETH Zurich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEAS
R-help is not at all the place to ask such a question.
Please use ess-help (or an Xemacs mailing list).
> "WS" == Werner Stanzl <[EMAIL PROTECTED]>
> on Fri, 2 May 2008 20:41:21 -0400 writes:
WS> Hi,
WS> When I start Xemacs version 21.5.28 (in Unix), using the init.el:
, df=1.66)? Where should I put m and s when I do simulation?
m + s * rt(n, df= df)
[I still hope this isn't a student homework problem...]
Martin Maechler, ETH Zurich
k> m s df
k> 0.2340746 4.0447124 1.6614823
k> (0.3430796) (0.4158891) (0.263870
e at the moment to look at the example data and
the model, and show you how to use it for nlrob();
if you find a way to you it for nls() , then the same should
work for nlrob().
I'm CCing this to the specialists for "Robust Stats with R"
mailing list, R-SIG-robust.
Best re
> "RT" == Rolf Turner <[EMAIL PROTECTED]>
> on Wed, 14 May 2008 13:12:02 +1200 writes:
RT> I have written a function which reads data from files using read.delim
RT> ().
RT> The names of these files are complicated and are built using arguments
RT> to the wrapper funct
Oops, there where some typos / thinkos in my post,
in the "more sophisticated" tryCatch() example.
I append the full message with corrected example.
--
Martin
---
> "RT" == Rolf Turner <[EMAIL PROTECTED]>
> on Wed, 14 Ma
le with R 2.7.0, David?
I'll not really rewrite it (to make it "nice" in my eyes),
but I'm fixing the problems with the fact that the 'Matrix'
package you were using back then has been much changed in the
mean time.
Expect a corrected benchmark R script within a day or s
>>>>> "MM" == Martin Maechler <[EMAIL PROTECTED]>
>>>>> on Wed, 14 May 2008 12:05:00 +0200 writes:
>>>>> "PhGr" == Philippe Grosjean <[EMAIL PROTECTED]>
>>>>> on Tue, 13 May 2008 16:10:15 +0200 w
t;R"
> gc()
used (Mb) gc trigger(Mb) max used(Mb)
Ncells1570998.4 3518.7 35 18.7
Vcells 537027103 4097.2 1409694707 10755.2 1342332915 10241.2
## after GC, we are back to 4 GB :
> sfsmisc::Sys.ps()
pid pcpu t
:
Do use data.matrix(dd) instead of
as.matrix(dd)
if dd is a data frame ... it will also produce a numeric matrix
when dd contains factor (and similar) columns.
Martin Maechler, ETH Zurich
EW> - Edward
EW> On Wed, May 21, 2008 at 11:31 AM, Peter Alspach
E
ntioned below,
I found it useful to add
dDA()
diagDA()
to package 'sfsmisc'
---> install.packages("sfsmisc")
library(sfsmisc)
?diagDA
and that help page mentions that these are just improvements of
earlier code of Sandrine Dudoit and Jane Fri
;, what's your name?
or equivalently
> cat("My name is \"Luis\", what's your name? \n")
My name is "Luis", what's your name?
Pay close attention to both single and double quotes
used above and to the use of '\&#x
added the following
to the 'ChangeLog' of my package "sfsmisc" :
2007-10-13 Martin Maechler <[EMAIL PROTECTED]>
* R/prettylab.R (eaxis): new function for nice (log) axis
labeling. (pretty10exp): drop.1: -10^k instead of -1*10^k
Here is the code -- which a
> "JH" == Johannes Huesing <[EMAIL PROTECTED]>
> on Thu, 22 Nov 2007 22:14:57 +0100 writes:
JH> Antony Unwin <[EMAIL PROTECTED]> [Thu, Nov 22,
JH> 2007 at 12:43:07PM CET]:
>> There have been several constructive responses to John
>> Sorkin's comment, but none of them ar
> "GaGr" == Gabor Grothendieck <[EMAIL PROTECTED]>
> on Mon, 26 Nov 2007 23:25:45 -0500 writes:
GaGr> See ?package.skeleton
GaGr> ?prompt
Yes, indeed. And I think that was the only pertinent answer to
Edna's question.
GaGr> On Nov 26, 2007 6:59 PM, Edna Bell
GaG
hc <- hclust(...)
X <- as.dendrogram(hc)
and then "work with that" (add 'nodePar' properties to certain nodes).
example(plot.dendrogram)
gives you a glimpse of the possibilities.
Someone else will have to give you details,
or provide you with utility functions
.r-project.org/incoming
That's a place NOBODY can get it - but the CRAN maintainers,
and they are not very pleased to see a package there which does
not pass R CMD check.
Please use a different place!
If you cannot, please
copy-paste one of the *.Rd files into an e-mail (reply to this).
Martin M
other color
space functionality in R all of which deserve more usage
in my oppinion and also in my own code ! ;-)
Package 'vcd' (and others) use package 'colorspace',
and I have wondered in the past if these color space computations
should not be merged into to standard R (package &
D() in package 'Matrix' (based on a donation from
Jens Oehlschlaegel) provides a more sophisticated algorithm for
this problem.
If you really only need the eigenvalues of the "corrected"
matrix, you might want to abbreviate the nearPD() function by
just returning the final &#x
s all of
approxfun(), splinefun() or ecdf()
}
You can also use
ecdf(x)(x)
and indeed check that it is identical to the convoluted
percentrank() function above :
> ecdf(x)(0.48742905)
[1] 0.52
> ecdf(x)(x)
[1] 0.20 0.44 0.12 1.00 0.48 0.16 0.56 0.72 0.60 0.28 0.96 0.52 0.24 0.04 0.92
ve just posted about to R-help as well].
Please say
example(colorRamp)
in R
and slowly watch the output, and I expect you will never ever
want to use the horrible "Matlab-like" color palette again..
Regards,
Martin Maechler, ETH Zurich
Earl> efg
Earl&g
> "MS" == Marc Schwartz <[EMAIL PROTECTED]>
> on Wed, 05 Dec 2007 12:43:50 -0600 writes:
[]
MS> Martin,
MS> Thanks for the corrections. In hindsight, now seeing the intended use of
MS> ecdf() in the fashion you describe above, it is now clear that my
MS> a
is one imprecise citation).
And that's not counting the hacker's attempt to break into mail
servers, notably if they are at "prestigious" places, DNS
attacks, etc.
Currently, we *must* use blacklist services, or our servers
would crumble.
The exact set of blacklist services o
Dir)
for(i in grep(paste("^", pfnStub,".*\\.RData$", sep = ""), cvListFiles)) {
load(paste(fnDir, cvListFiles[i], sep = "/"))
myFunction(rliObject)
rm(rliObject)
}
Something many readers on this list will much prefer to your
original, and hen
2) ## using plot() method for functions
or simpler to understand:
xx <- seq(0,1, length=100)
lines(xx, myF(xx), col=2)
t> Thank you very much
t> TDB
you're welcome.
Martin Maechler, ETH Zurich
__
R-help@r-project.org mail
> "KK" == Knut Krueger <[EMAIL PROTECTED]>
> on Sat, 15 Dec 2007 08:39:02 +0100 writes:
KK> Bob Green schrieb:
>>
>> > dates <- read.csv("c:\\dates.csv",header=T)
>> > dates
>> v1 v2
>> 1 12/12/1978 12/12/2005
>> 2 23/01/1965 23/09/2001
>> 3 24/
Its definition is currently
inv.seq <- function(i) {
## Purpose: 'Inverse seq': Return a short expression for the 'index' `i'
##
## Arguments: i: vector of (usually increasing) integers.
## --
>>>>> "TP" == Tony Plate <[EMAIL PROTECTED]>
>>>>> on Fri, 21 Dec 2007 18:17:18 -0700 writes:
TP> Martin Maechler wrote:
>>>>>>> "MS" == Marc Schwartz <[EMAIL PROTECTED]>
>>>>>>
ith sparse matrices:
> library(Matrix)
> m <- Matrix(0,6,6)
> diag(m[-1,]) <- 1:5
> m
6 x 6 sparse Matrix of class "dgCMatrix"
[1,] . . . . . .
[2,] 1 . . . . .
[3,] . 2 . . . .
[4,] . . 3 . . .
[5,] . . . 4 . .
[6,] . . . . 5 .
- - - - - - - - - - -
M
> "ML" == Mark Leeds <[EMAIL PROTECTED]>
> on Fri, 21 Dec 2007 19:38:19 -0600 (CST) writes:
ML> I was playing around with a simple example using solve.qp ( function is
in the quadprog package ) and the code is below. ( I'm not even sure there if
there is a reasonable solution bec
> "SG" == Sven Garbade <[EMAIL PROTECTED]>
> on Fri, 11 Jan 2008 09:47:39 +0100 writes:
SG> Hi list,
SG> from time to time I got an "Error: bad value" and must restart R. The
SG> mail archives suggests memory corruption, but I do not run "special" C
SG> code, only base
> "GS" == Gavin Simpson <[EMAIL PROTECTED]>
> on Thu, 10 Jan 2008 14:16:36 + writes:
GS> On Thu, 2008-01-10 at 10:48 +, Marc Moragues wrote:
>> Hi,
>>
>> I am trying to calculate a distance matrix on a binary
>> data frame using dist.binary() {ade4}. This i
aster, the Fedora-maintainer of its
ESS package) that it was only for Fedora-devel; maybe my
misunderstanding.
For that reason I had only updated the web page
http://ESS.r-project.org/
without making a "big" announcement
[ which would have been to the ESS-help mailing list,
--&g
95)/2
RM> [1] 0.15
>> sd(c(1.25, 0.95))
RM> [1] 0.2121320 # why it is different from 0.15?
because 1 is different from 2 !
If 2 was 1, than sqrt(2) == 1 as well, but actually I don't
think the universe and we all would exist in that case
Martin Maechler, ETH
ed spam rate on r-help.
I have to tell you (or rather the R-help audience) that I've
spent too much time (again!) during this week trying to fight
spam.
Maybe the time of "free" (unauthorized) posting to mailing lists
is approaching its end...
Martin Maechler, ETH Zurich
_
ed spam rate on r-help.
I have to tell you (or rather the R-help audience) that I've
spent too much time (again!) during this week trying to fight
spam.
Maybe the time of "free" (unauthorized) posting to mailing lists
is approaching its end...
Martin Maechler, ETH Zurich
_
g/
-
which is slightly preferred.
BTW: I've been looking for a (or a few) volunteer to rewriting
that web page using things like CSS and/or PHP but not using a
CMS with data base back end).
Martin Maechler, ETH Zurich
ChuckB>
er seeing a few particularly nasty orbits.
As we know ``from Chaos theory'', there can be delicate
inhereent and numerical problems in ODE solving..
Regards,
Martin Maechler, ETH Zurich
__
R-help@r-project.org mailing list
https://stat.eth
etc etc (because the so-called *raw* periodgram is
*in*consistent as an estimate fo the underlying true spectrum).
The Time-Series chapter/section of the MASS book is very helpful
here, IIRC.
Martin Maechler, ETH Zurich
>>>>> "OB" == Oliver Bandel <[EMAIL PROTECTED]>
&
n/listinfo/r-help
{please, Ram, also do read the above web page carefully ..}
Regards,
Martin Maechler,
ETH Zurich (provider of all the @r-project.org mailing lists)
ss> On Sat, Sep 6, 2008 at 4:48 AM, Ram Kumar Basnet <[EMAIL PROTECTED]>
wrote:
>> hi
>> I am expecti
nge the code
to the following but it didnt work, any suggestions?
>>
>> labels<-as.character(levels(list(sex,Age)))
>> ybar<-as.vector(tapply(ConditionIndex,sex,Age,mean))
>> error.bars(ybar,se,labels)
FEH> See http://biostat.mc.vanderbilt.edu/Dy
al to the 2.7.2 one:
https://svn.r-project.org/R/trunk/src/library/stats/R/t.test.R
Note that most people would want to *search* in the R sources,
and therefore still get the complete R sources (either by
subversion or from a *.tar.gz file, from CRAN).
Yes, Uwe Ligges wrote an R News (2006
functions as ASCII files or
YX> even write a package are also good solutions :-)
^
"also good..." baahh... *THE* good!
I'd strongly second Adai's view here:
*DO* use *.R source files --- together with an R-syntax aware
editor, there are many
hank you again
TS> for the help!
So you now have understood that R's behavior of handling +/- Inf
in this respect is rather excellent than bogous ?
Martin Maechler, ETH Zurich (and R-core team)
__
R-help@r-project.org mailing list
https
> "JIV" == Jorge Ivan Velez <[EMAIL PROTECTED]>
> on Fri, 12 Sep 2008 18:39:04 -0400 writes:
JIV> Dear Amanda,
JIV> Try this:
JIV> # Data set
JIV> set.seed(123)
JIV> DF=data.frame(A=rnorm(10),B=rpois(10,10),C=sample(1:3,10,replace=T))
JIV> attach(DF)
JIV>
r robust
"location+scatter" [aka "(mu , Sigma)"] from packages
'rrcov', 'robustbase' and/or 'robust', e.g.,
robustbase::covMcd
If you need to know more, please divert to the
SIG (special interest group) mailling list R-SIG-robust
to which I CC this.
M
DR> a
Yes, indeed, it needs the t(.) trick.
Note that 'Matrix' package has a function forceSymmetric(.) to
do this for you (faster, using C code):
A <- forceSymmetric(Matrix(rnorm(36), 6))
is all you'd need {if can afford to trash half of the random
numbers
notably since in R we had introduced the explicit long (= long
integer) constants, using the 'L' suffix,
i.e., 7L is "integer"
7 is "double"
Note however that for both, is.numeric(.) is fulfilled and
class(.) and mode(.) return "numeric".
Only t
tatistical Computing,
Martin Maechler
ETH Zurich
Note: As an exception, we are announcing this through R-announce.
___
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-announce
_
0 0
[3,]001100110 0 1
[4,]01010101 0 1 0
But really look at the examples in help(digitsBase)
for more.
Martin Maechler, ETH Zurich
>> Since I have to teach number base conversion within 2 wee
> "AR" == Arthur Roberts <[EMAIL PROTECTED]>
> on Fri, 26 Sep 2008 01:54:59 -0700 writes:
AR> Hi, all,
AR> Is there an R-command that will convert shell("cd") to a string, so
AR> that I can use it in the paste() command?
I'm pretty sure that for the specific case, you sh
I'm replying to this on the R-devel mailing list and the
Subject "On modes, types and R documentation"
Martin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/po
gs as normal.
The mails should be kept in a secondary/substitute mail server
and be forwarded to ours as soon as that "gets up" again.
But you will notice that your message is not posted to the lists
for several hours.
We are sorry for the inconvenience to the R community.
Martin Maech
> <[EMAIL PROTECTED]>
> on Thu, 2 Oct 2008 10:02:04 -0700 writes:
> Sorry, I must be looking at a different section but when I look at 3.4 in
r-intro.pdf I see:
> 3.4 The class of an object
> All objects in R have a class, reported by the function class. For simple
vect
Yes masking var() is a bad idea;
a much better (and much less error-prone) idea would be to install
R 2.8.0 alpha even now.
It will become 'beta' early next week.
We are asking the R community to please install and use
pre-release versions of R (if you can / are allowed to)
at least fro
> "AMT" == Arnau Mir Torres <[EMAIL PROTECTED]>
> on Tue, 14 Oct 2008 17:13:01 +0200 writes:
> "AMT" == Arnau Mir Torres <[EMAIL PROTECTED]>
> on Tue, 14 Oct 2008 17:13:01 +0200 writes:
AMT> Hello.
AMT> I need to know how can R compute AIC when I study a regression
c(-1,2,-3,-4,5,6,-7,8,-9,10)
f <- cbind(c,d)
dat <- data.frame(f)
rm(c,d,f) ## << check to make sure 'c' or 'd' are not picked up from GlobalEnv
[MM]
##
daN.lm <-lm(c ~ d, data = dat, subset = d < 0)
summary(daN.lm)
Does this really not work for
"")) : cannot open 'pdf' file argument
SK> 'Pictures/full_map.pdf'
SK> Just close the Acrobat Reader or whatever.
or even better: Close your OS (MS Windows) and start using a
better one.
E.g., on all Linux versions I know I can write to files open in
ac
at this
would still not be equivalent to "kmeans with manhattan
distance".
In particular, pam() should be more robust than any kmeans
incantation, since pam() does not use any (non-robust) mean.
Martin Maechler,
ETH Zurich
TG> On Wed, Dec 17, 2008 at 1:25 PM, Christian Hennig
a co-maintainer of the package, I'm biased of course, but
also extremely interested in feedback / questions from users
like you, about missing features, unclear documentation,
possible bugs, ...
Regards,
Martin Maechler, ETH Zurich
>>>
>>>
> "BertG" == Bert Gunter
> on Tue, 6 Jan 2009 10:13:29 -0800 writes:
BertG> Folks: I am getting duplicate messages on
BertG> posts. Please correct my details if I'm wrong, but I
BertG> believe it's because folks are posting to **both**
BertG> the addresses, r-help@r-p
illing to fulfill a New Year's wish for me,
I'd be very grateful for a (digital scan) copy of that page.
Martin Maechler,
ETH Zurich
[..]
>> Thought you might be interested in reading this article,
>> which appears in the 1/6/9 online edition of T
>>>>> "MM" == Martin Maechler
>>>>> on Wed, 7 Jan 2009 14:18:24 +0100 writes:
[..]
MM> At the bottom of the page it says
MM> __ A version of this article appeared in print on
MM> January 7, 2009, __ on page B6 of
or those in the S / R tradition you already go a
traditional one...
but to finally help you:
use
ag <- agnes(...)
dg <- as.dendrogram(as.hclust(ag))
plot(dg)
and if you read and look at the examples of
help(plot.dendrogram)
you'll see that you have many more option
do rocket science (:-)
I've solved this exercise for you:
require("Matrix")
repRow <- function(m, i, times)
{
## Purpose: return a sparse matrix containing row m[i,] 'times' times
## --
## Arguments: m: sparseMatrix; i: row index; time
r text above could be read differently) this
bug ... quite long standing! ... seems to be much wider spread
than just (Open)SuSE. I have seen and still see it on Redhat
(RHEL 5) and many Ubuntu versions too.
And indeed, the workaround recommended (somewhere) in the above
URL is working for me too:
A
s) maintainable by SVN"
but *not* restricting it to "no-server-side modules required".
Martin Maechler, ETH Zurich
Ao> - I volunteer in both cases :)
Ao> Winner of Design Contest should get
Ao> some bragging rights in a small hyperlink (with nofollow t
precision x now use a long
double accumulator where available and so more closely match
sum() and prod() in potentially being more accurate.
and indeed, in R-devel,
sum(x) - cumsum(x)[length(x)]
gives 0 for your example.
Martin Maechler, ETH Zurich and R-core team
GaGr&g
probably yes, platform-dependently;
However I vaguely remember that I didn't see one such case in the few
experiments I did.
Martin
SM> On Wed, Feb 18, 2009 at 8:03 AM, Martin Maechler
SM> wrote:
SM> ...
>> o cumsum(x) and cumprod(x) for double precision x no
onality of
such a function, bandMatrix(.) say ?
It probably should have an option for specifying symmetric
vs. "general" band matrices, something that the "corresponding"
Matlab function [ spdiags() ] does not.
Martin Maechler, ETH Zurich
TL> -thom
# Arguments: (n,m) : Matrix dimension
##k : integer vector of "diagonal numbers", with identical
##meaning as in band(*, k)
## listDiag: list of (sub)diagonals
## symmetric: if TRUE, specify only upper or lower triangle;
7;d seen that they
really are supposed to be binary only.
A very small change in the sources does accomplish this, passes
the standard checks (and I cannot imagine reasonable code that
would have relied on the more lenient behavior), so
this will have changed in one of the next versions of R-deve
>>>>> "WK" == Wacek Kusnierczyk
>>>>> on Mon, 23 Feb 2009 19:29:31 +0100 writes:
WK> Martin Maechler wrote:
>>>>>>> "WK" == Wacek Kusnierczyk
>>>>>>> on Mon, 23 Feb 2009 12:06:32 +0
me to provide something
"uniformly better" (i.e., for those who want, it *must* work via
e-mail exclusively, not via compulsory web-browser clicking).
Martin Maechler,
ETH Zurich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
'else' to avoid a syntax error in entering a 'if ... else'
> construct at the keyboard or via 'source'. For that reason, one
> (somewhat extreme) attitude of defensive programming is to always
> use braces, e.g., for 'if' clauses.
Regar
>>>>> "WK" == Wacek Kusnierczyk
>>>>> on Wed, 25 Feb 2009 09:46:19 +0100 writes:
WK> Martin Maechler wrote:
>>
>> I think this is FAQ (or should become one):
>>
>> ?if [the help page you really shou
http://cran.uk.r-project.org/bin/windows/contrib/2.8/Matrix_0.999375-21.zip'
>> cannot open: HTTP status was '404 Not Found'
"Further" (well, maybe the reason for the whole thread):
Matrix has recently been upgraded, so the above URL indeed
should no
n negative values for x and y axis
>>
>> i can plot point using points(0,0) but if i want to draw lines
>> along it, can't seem to get it right
PaCo> ?abline
PaCo> maybe abline(h = 0)
PaCo> abline(v = 0
for large x with only few NAs.
But this now has *REALLY* changed into a topic belonging to
R-devel, not R-help
~~~ --> hence I've diverted the thread to there.
Martin Maechler, ETH Zurich
CAPE> Charles Annis, P.E.
CAPE> charles.an...@statisticalengineering.com
[.]
the footer here tells us that you still don't follow the posting
guide (namely, not sending HTML) even though you may have
started doing so by trying to produce reproducible code.
Regards,
Martin Maechler, ETH Zurich
aw>
ommenges is pretty active in this area.
Further, Marloes Maathuis is active in this field, too.
Her (CRAN) R package is 'MLEcens'
which computes the (nonparametric!) MLE for (bivariate or univariate)
interval censored data, which -- I think -- is what the
KM estimator is for the one-
> "GS" == Gavin Simpson <[EMAIL PROTECTED]>
> on Tue, 29 Jan 2008 00:09:05 + writes:
GS> hits=-2.6 tests=BAYES_00
GS> X-USF-Spam-Flag: NO
GS> On Mon, 2008-01-28 at 12:17 -0700, Michelle DePrenger-Levin wrote:
>> I was asked for the following information and hope it
cribers, these latter decision coming
later.
Consequence:
Please subscribe to R-help __ at least in digest mode __
now, and ``spread the word''
so that you will be able to continue asking questions
and getting advice here.
Martin Maechler, ETH Zurich
... a sad list administrato
>>>>> "DM" == Dieter Menne <[EMAIL PROTECTED]>
>>>>> on Thu, 31 Jan 2008 08:58:54 + (UTC) writes:
DM> Martin Maechler stat.math.ethz.ch> writes:
>>
>> The consequence of this will be that .. before the weekend
> "NS" == Ng Stanley <[EMAIL PROTECTED]>
> on Thu, 31 Jan 2008 16:56:26 +0800 writes:
NS> Hi,
NS> Is this a bug from hist() ? The total density is greater than 1.
no, the bug is in your thinking...
Hint: 2 * 1/2 = 1
>> test
NS> [1] 0.05077802 -0.50585520 -0.980536
>>>>> "DM" == Duncan Murdoch <[EMAIL PROTECTED]>
>>>>> on Thu, 31 Jan 2008 08:07:06 -0500 writes:
DM> On 1/31/2008 6:48 AM, Martin Maechler wrote:
>>>>>>> "DM" == Dieter Menne <[EMAIL PROTECTED]&g
> "BB" == Ben Bolker <[EMAIL PROTECTED]>
> on Fri, 1 Feb 2008 14:03:15 + (UTC) writes:
BB> Paul Smith gmail.com> writes:
>>
>> Dear All,
>>
>> Take this code:
>>
>> > f <- function(x) exp(-x)*x-0.05
>> > g <- function(x) 0
>> > curve(f,0,5)
ole R-help list)
if you are volunteering for the above,
not later than Monday, Feb.4 (any timezone).
Martin Maechler,
for the R-core team
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
#-> large (!) 3000 x 3000 / 2
sk <- silhouette(km$cl, dissE)
plot(sk)
## but kmeans is rather equivalent to work with {D_ij}^2,
## hence this better corresponds:
dE2 <- dissE^2
sk2 <- silhouette(km$cl, dE2)
plot(sk2)
LS> Many thanks!
LS> Linda
you're welco
s / functions in R,
I'm very strongly convinced that using GNUplot for plotting is
``the wrong'' approach by almost all definitions of "wrong".
Martin Maechler, ETH Zurich
__
R-help@r-project.org mailing list
https://stat.ethz
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]>
> on Mon, 3 Mar 2008 14:16:15 + (GMT) writes:
[]
BDR> In a later message Louise mentioned the desire to use TeX fonts for
BDR> annotation, to match a LaTeX document. Paul Murrell has pointed out
his
BDR>
overlooked,
I'll add it to 'See also' as well.
Adding an example to 'Examples' is a bit more awkward,
since it needs to load the methods package;
but I'll add a "don't run" example.
Martin Maechler, ETH Zurich
__
s = )
and setAs(.)
much more frequently. That may well depend on the kind of your
classes of course.
For the 'Matrix' package at least, our experience and use
testing has resulted in almost no use of setIs()
but rather hierarchical class definitions, i.e., the above
setClass(
> "GaGr" == Gabor Grothendieck <[EMAIL PROTECTED]>
> on Tue, 18 Mar 2008 23:50:28 -0400 writes:
GaGr> Try:
GaGr> eapply(.GlobalEnv, class)
GaGr> or perhaps
GaGr> unlist(eapply(.GlobalEnv, class))
GaGr> or
GaGr> str(eapply(.GlobalEnv, class))
nice! {and we'
x27;) not an "error" !)
> Attaching package: ''
>
>
> The following object(s) are masked from package:stats :
>
> var
Are you?
In any case, package 'circular' is one CRAN package example that
uses a NAMESPACE and redefines var(),
for this listserv indicated
JG> that PDFs were fine):
Yes, indeed, PDFs *are* fine
(and I don't know exactly what Dieter's problem was with it;
he was right that you should include self-contained code).
Martin Maechler, ETH Zurich
__
R-h
701 - 800 of 837 matches
Mail list logo
