. The bug is still there with some
scaling of the matrix (from about divide by 2 to multiply by 8). The
bug disappears if the matrix is dumped and sourced back in again.
The only clue that I can offer is that it is vectors 43:50 (I think) that
are NaN's and the matrix is logically of r
Doing
help.search('precedence')
comes up empty. A fix would be to have the title:
Operator Syntax and Precedence
instead of
Operator Syntax
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the U
the proper
length of the answer at the outset since that is known.
It could then be used as an example (somewhere) of
how not to fragment memory.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Prof
On page 10 (section 1.1.4) of Writing R Extensions
version 2.1.1 the following two phrases appear:
define these function in a file
after the packages is
The 's' from 'functions' fell down the page and attached
itself to 'package'.
Patrick Burns
[EMAIL PROTE
I'm not so sure that it is a typo or ungrammatical as much as
an awkward way of stating it. Maybe something along the
lines of:
The resolution of the times will be system-specific; it is
common for the resolution to be on the order of 1/100
second, ...
Patrick Burns
[EMAIL PROTECTED]
+44
You can look at the definition of 'corner' in the public
domain area of the Burns Statistics website. It uses
'do.call' on '[' to achieve (sort of) what you want.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and &qu
12
day 07
svn rev 36656
language R
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
_
____
R-devel@r-project.org <mailto:R-devel@r-project.org> mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
--
Patrick Burn
ow you to
specify the version of a dependency so that this doesn't happen. You
can do this with R, but because it's hard to have multiple versions of
the same package installed at the same time, it's not as useful.
Hadley
--
Patrick Burns
pbu...@pburns.seanet.com
http://www.burns-
he "surprise" section of the examples
in the 'sample' help file
sample(x[x > 9])
and
sample(x[x > 9], max=FALSE)
have different behaviours.
By the way, I'm certainly not convinced that
'max' is the best name for the argument.
--
Patrick
Is there a reason that 'which.min' and
'which.max' don't have an 'arr.ind'
argument?
The context in which I wanted that was
a grid search optimization, which seems
like it would be reasonably common to me.
--
Patrick Burns
pbu...@pburns.seanet.com
http://www.
On 05/07/2010 10:56, Martin Maechler wrote:
"PatB" == Patrick Burns
on Sun, 04 Jul 2010 09:43:44 +0100 writes:
PatB> Is there a reason that 'which.min' and
PatB> 'which.max' don't have an 'arr.ind'
PatB> argument?
w
E' as argument to
arrayInd() but have which() call arrayInd(..., useNames=TRUE).
This is a back-compatibility change in arrayInd() -- which has
existed only since 2.11.0 anyway, so would seem ok, to me.
Opinions ?
I find it hard to believe that would cause
too much trauma.
Pat
I think it makes sense to have
'symbol' in the See Also of 'symbols'
and vice versa.
--
Patrick Burns
pbu...@pburns.seanet.com
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of 'Some hints for the R beg
The 'extended' argument to 'strsplit'
has been removed, but it is still mentioned
in the argument items in the help file
for 'fixed' and 'perl'.
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.com/blog
http:/
ist DIGESTED
https://stat.ethz.ch/mailman/listinfo/r-devel
End of R-devel Digest, Vol 95, Issue 14
***
[[alternative HTML version deleted]]
______
R-devel@r-project.org mailing list
https://stat.ethz.c
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
La version française suit le texte anglais.
--------
This email may contain privi
a Beleites wrote:
On 01/18/2011 10:53 AM, Patrick Burns wrote:
I'm not at all a fan of thinking
of the examples as being tests.
Examples should clarify the thinking
of potential users. Tests should
clarify the space in which the code
is correct. These two goals are
generally at odds.
Patric
packages:
[1] xts_0.8-0 zoo_1.6-5
loaded via a namespace (and not attached):
[1] grid_2.13.0 lattice_0.19-23 tools_2.13.0
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of 'Some hints for the R beginner'
uot;))
}
Pat
On 29/04/2011 15:25, Martin Maechler wrote:
Paul Johnson
on Thu, 28 Apr 2011 00:20:27 -0500 writes:
> On Wed, Apr 27, 2011 at 12:44 PM, Patrick Burns
>wrote:
>> Here are some data frames:
>>
>> df3.2<- data.frame(1
break
}
if (inherits(res, "try-error")) {
if (ThrowError) {
stop("file connection failed")
}
else {
warning("file connection failed")
}
}
}
People have done network traces and such, but so far no
in mind uses 'is.unsorted' in
'stopifnot'.
Pat
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @burnsstat @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of:
'Impatient R'
'The R Inferno'
'Tao Te Programming
x27;m missing some subtlety, but it makes
things work in my case.
Pat
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @burnsstat @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of:
'Impatient R'
'The R Inferno'
'Tao Te Programming
kage.
This is being done in RStudio on Windows
with R version 3.0.2.
The package used to work. There was a very
minor change to one function when the package
started not working.
The zip files that don't work are the same size
as the one that does work.
Thanks for any suggestions.
Pat
project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of 'Some hints for the R beginner'
and 'The R Inferno')
__
.15.0
>>>
>>> __
>>> R-devel@r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-devel
>>
>
> -BEGIN PGP SIGNATURE-
> Version: GnuPG v1.4.10 (GNU/Linux)
> Comment: Using GnuPG with Mozilla - http
but
to be fair it is a pretty gruelling process
to find the Task Views.
May I suggest that there be a "Task Views" item
on the left sidebar of the R website in the
Documentation section?
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.
ertain amount
of brevity.
Pat
On 12/02/2012 10:43, Barry Rowlingson wrote:
On Sat, Feb 11, 2012 at 5:55 PM, Patrick Burns wrote:
Now it could be that people are not trying
very hard to solve their own problems, but
to be fair it is a pretty gruelling process
to find the Task Views.
May I sugge
on my own?
Thanks for your help.
Kind Regards,
Thorn
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(ho
on purpose?
Pat
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of 'Some hints for the R beginner'
and 'The R Inferno')
__
R-devel@r-project.org mai
Both 'exists' and 'get' silently ignore all but the
first element.
My view is that 'get' should do what it currently does
except it should warn about ignoring subsequent elements
if there are any.
I don't see a reason why 'exists' shouldn't be
LC_CTYPE=English_United Kingdom.1252
[3] LC_MONETARY=English_United Kingdom.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United Kingdom.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @
I think the rule is that you can do anything as long as you don't
complain. If you want to complain, you must follow the instructions.
-- Jari Oksanen in
Re: [Rd] Keeping up to date with R-devel
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @burnsstat @portfolioprobe
Brian Ripley:
If things are not readily available in R it is always good to pause and
reflect if there might be a good reason.
In the R-help thread: How to get the t-stat for arima()?
Pat
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @burnsstat @portfolioprobe
http
st(a, b, c) else list(x, y, z)
myvar
[[1]]
[1] 4
[[2]]
[1] 5
[[3]]
[1] 6
--
Patrick Burns
pbu...@pburns.seanet.com
twitter: @burnsstat @portfolioprobe
http://www.portfolioprobe.com/blog
http://www.burns-stat.com
(home of:
'Impatient R'
'The R Inferno'
'Tao Te Prog
'as.name' doesn't recognize a name with
its namespace extension as a name:
> as.name("lm")
lm
> as.name("stats::lm")
`stats::lm`
> as.name("stats:::lm")
`stats:::lm`
Is there a reason why it shouldn't?
Pat
--
Patrick Burns
pbu...
Okay, that's a good reason why it shouldn't.
Why it should is that I want to substitute
the first element of a call to be a function
including the namespace.
Pat
On 23/04/2013 18:32, peter dalgaard wrote:
On Apr 23, 2013, at 19:23 , Patrick Burns wrote:
'as.name' does
7;s third method, but it seems like
it assumes the solution. Moving functions around
is unappetizing for my use -- this is for testing
and keeping things as faithful to real use is a
good thing.
Pat
On 23/04/2013 21:18, Duncan Murdoch wrote:
On 13-04-23 3:51 PM, Patrick Burns wrote:
Okay, that
with different levels of detail. A
sophisticated version of this scheme
might even assume different levels of
knowledge of the user so that the least
detailed level might be longer (but
easier) than a more detailed level.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-sta
Under the system of development we
now have, I agreee with Seth's
assertion. But if there were
people dedicated to documentation,
then I think something like what I
described could be workable.
Pat
Seth Falcon wrote:
Writing good documentation is hard. I can appreciate the desire to
find tec
is missing, with no default
Then that would be satisfactory, I think.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno" and "A Guide for the Unwilling S User")
__
R-devel@r-
ars of typing the "1".
Duncan Murdoch
In the spirit of sour grapes, the proposed default
might discourage some users from vectorizing their
thinking.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno&quo
I'm confused about how to change a repository
so that it doesn't hit the bug in 'available.packages'
in 2.10.0 that was fixed in 2.10.1.
I presume it involves adding fields to the
PACKAGES file.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat
ist(lapply(factors, as.character))
factor(char, levels = levels)
}
c(factor("a"), factor("b"), factor(c("c", "b","a")), factor("d"))
# [1] a b c b a d
# Levels: a b c d
Hadley
--
Patrick Burns
pbu...@pburns.seanet.
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
--
Patrick Burns
pbu...@pburns.seanet.com
http://www.burns-stat.com
(home of 'Some hints for the R beginner'
and 'The R Inferno')
___
e error quickly.
Also, using '<<-' is generally a bad idea.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Tong Wang wrote:
>HI, all
>I am debugging an R code with dynamical
If the sentence in question were amended to:
Values of trim outside that range ...
then I think it would rule out the misinterpretation of
the sentence.
Pat
Prof Brian Ripley wrote:
>There is only one _range_ mentioned, (0, 0.5). I don't see how you can
>construe 'that range' to be a refere
ter removing NAs. The
first is clearly impossible. The second has some actual
value, but we don't have enough information to have an
estimate of the value.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User&
les
should have a Side Effects section (to make it clear to
new users what side effects are and that they are not a
good thing for most functions to have). I can understand
there could be another point of view on that. However, I
definitely think that there should be a Side Effects section
in the he
Patches to the help files sound like a good idea. However,
it isn't something I'm likely to get to immediately. I'm
hoping that some other nice person will volunteer.
Pat
Duncan Murdoch wrote:
> On 12/17/2007 6:00 AM, Patrick Burns wrote:
>
>> I recently had a d
der.max)[xaic == 0]
else order.max
}
I presume something similar can be done in 'ar.burg',
but I haven't looked. 'ar.ols' and 'ar.mle' look problematic
for fixes.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20
Martin Maechler wrote:
> [ ... ]
>
>But actually, part of the changed behavior may be considered
>undesirable:
>
> rnorm(2, mean = NA)
>
>which gives two NaN's would now produce a warning,
>where I could argue that
> 'arithmetic with NAs should give NAs without a warning'
>since
> 1:2 + NA
>
Peter Dalgaard wrote:
> Duncan Murdoch wrote:
>
>> I haven't done it, but I suspect we could introduce special behaviour
>> for ??foo very easily. We could even have a whole hierarchy:
>>
>> ?foo, ??foo, ???foo, foo, ...
>>
>>
>>
> Heh, that's rather nice, actually. In words, that
xx
$A
changed b
1 2
This is observed in 2.4.0 on Linux as well as 2.7.0 and
2.8.0 on Windows XP.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
___
seldom enough that adding a warning isn't going
to slow things down appreciably.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Basically the only thing in the thread that was clear
to me was Brian's phrasing. So I'd suggest basing
any changes on that.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
S Ellison
If the positive definiteness of the covariance
is the only issue, then you could base a penalty on:
eps - smallest.eigen.value
if the smallest eigen value is smaller than eps.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide fo
re well defined.
This is the reason why I wonder if setting by myself a $-\infty$ in
the composite likelihood function is appropriate here.
However, you might be right in setting a tolerance value 'eps' instead
of the theoretical bound eigen values > 0.
Thanks for your tips,
Bes
upon a time 'by' gave 3 as the answer in both
situations.
'by.default' used to be a one-liner, but now decides
what to do based on 'length(dim(data))'.
This specific problem goes away if the line:
if (length(dim(data)))
is replaced by:
if(length(dim(data)) > 1)
There is a minor typo in the help file for assignOps:
There is no space in "operator<-" in the second
sentence of the second paragraph of the Details
section.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for th
27;d prefer
an error like the matrix case.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
markle...@verizon.net wrote:
could someone explain what is happening below ? I was trying to solve
a re
ay the 'names' attribute really is
'dimnames[[1]]'.
Thanks for this explanation. One could then argue that [1:10,] is
somewhat better than [,1:10], but that is just polish.
Perhaps it could be:
[1:10(,)]
That is weird enough that it should not lead people to
believe
'The R Inferno' page 93 and page 99.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno" and "A Guide for the Unwilling S User")
twoutop...@gmail.com wrote:
The length() of a POSIXlt object is given as 9 r
between
the items.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno" and "A Guide for the Unwilling S User")
ivo welch wrote:
hi gabor: this would be difficult to do. I don't think you want to
read my progra
ly to mind
is 'read.table'.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Original Message
Subject:Re: [R] a strange problem with integrate()
Date: Wed, 01 Mar 2006 11
ed
> set.seed(2)
> r2 <- .Random.seed
> all.equal(r1, r2)
[1] "Mean relative difference: NA"
Warning message:
NAs produced by integer overflow in: target - current
A small change to 'all.equal.numeric' would fix that if it
is felt to be worthwhile.
Patrick Burns
[EM
sed not to encounter that in R and
hope it continues.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Prof Brian Ripley wrote:
>Looking at the code it occurs to me that there is another case yo
s.matrix(x)
if (diag)
ans <- row(x) <= col(x)
else ans <- row(x) < col(x)
if(value) x[ans] else ans
}
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
_
Gabor came close to the situation I had yesterday
that prompted me to write a local version of 'lower.tri'.
It was approximately:
x[sub, sub][lower.tri(x[sub,sub])]
Pat
Gabor Grothendieck wrote:
> On 8/6/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
>
>> Is there a case to be made for this?
Kingdom.1252;LC_CTYPE=English_United
Kingdom.1252;LC_MONETARY=English_United
Kingdom.1252;LC_NUMERIC=C;LC_TIME=English_United Kingdom.1252
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "utils" "datasets"
I rather like it.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Brian Ripley wrote:
>Unix versions of R CMD BATCH have reported proc.time() unless the script
>ends in q(). E.g. if the
71 matches
Mail list logo
