Re: [Rd] as.POSIXct behaviour
On 9 March 2017 at 01:29, Arunkumar Srinivasan wrote: > The time info is lost on the first index as well. And it happens *silently*. Yes, because it assumes homogeneous format on the entire vector. You may want to do two passes with different formats or run a regular expression to catch typos. > lapply(x, as.POSIXct) > A list is returned, but values are as I’d expect. Would it be possible > to retain the time info with as.POSIXct(x) directly as well? If not, This behaves differently because you are running as.POSIXct for each element separately. __ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
Re: [Rd] Error in formatDL(nm, txt, indent = max(nchar(nm, "w")) + 3)
Hello:
I tried "debug(help)" with the problem mentioned below. It
stopped with a call to "library", from which I generate the following
simple replication of this error:
> library(help = 'fda', character.only = TRUE)
Error in formatDL(nm, txt, indent = max(nchar(nm, "w")) + 3) :
incorrect values of 'indent' and 'width'
This was using R 3.3.3 inside RStudio on a Mac, as noted below.
This same "library(help='fda', character.only=TRUE)" ran fine using R
Console 3.2.1 under Windows 7.
I also tried debug(library). It died in line 298:
for (i in which(file.exists(docFiles))) pkgInfo[[i]] <-
readDocFile(docFiles[i])
with i = 1 and docFiles[i] =
"/Library/Frameworks/R.framework/Versions/3.3/Resources/library/fda/Meta/package.rds".
This looks like a bug to me. I don't know if it's new since R
3.2.1 or it works fine under Windows 7 but not a Mac.
Thanks,
Spencer Graves
*** NOTES:
1. When I tried "Q" to exit debug mode, I got, "Error retrieving
help: R code execution error." This suggests a bug in "debug", which
could be replicated by trying debug(help) with "help(package='fda',
help_type='text')". It does not want to exit from "help" using "Q".
"c" also fails.
2. For sessionInfo(), see below.
On 2017-03-08 4:06 PM, Spencer Graves wrote:
Hello:
A call to help(..., help_type='text') fails with "package='fda":
> install.packages('fda')
> help(package='fda', help_type='text')
Error in formatDL(nm, txt, indent = max(nchar(nm, "w")) + 3) :
incorrect values of 'indent' and 'width'
I have this wrapped inside "try" in PackageSum2{sos}, so it
doesn't die, but it also doesn't give me the additional information I
want from help(..., help_type='text')
Is this a bug?
Thanks,
Spencer Graves
Using RStudio 1.0.136 with
> sessionInfo()
R version 3.3.3 (2017-03-06)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierra 10.12.3
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets
[6] methods base
loaded via a namespace (and not attached):
[1] tools_3.3.3
__
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
[Rd] as.POSIXct behaviour
Dear R-devel, I have tested the code below on R v3.3.2 and v3.3.3 on
Mac and Windows.
x <- c("2017-01-01 05:00:02", "2017-01-02 03 :M:00") # note the ‘ :M’
in 2nd value
as.POSIXct(x)
# [1] "2017-01-01 GMT" "2017-01-02 GMT”
The time info is lost on the first index as well. And it happens *silently*.
On the other hand, if I do:
lapply(x, as.POSIXct)
# [[1]]
# [1] "2017-01-01 05:00:02 GMT"
#
# [[2]]
# [1] "2017-01-02 GMT”
A list is returned, but values are as I’d expect. Would it be possible
to retain the time info with as.POSIXct(x) directly as well? If not,
what’s the rationale? In addition, I think a warning message that some
malformed timestamps were found would be very useful.
Thank you,
Arun.
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
[Rd] Error in formatDL(nm, txt, indent = max(nchar(nm, "w")) + 3)
Hello:
A call to help(..., help_type='text') fails with "package='fda":
> install.packages('fda')
> help(package='fda', help_type='text')
Error in formatDL(nm, txt, indent = max(nchar(nm, "w")) + 3) :
incorrect values of 'indent' and 'width'
I have this wrapped inside "try" in PackageSum2{sos}, so it
doesn't die, but it also doesn't give me the additional information I
want from help(..., help_type='text')
Is this a bug?
Thanks,
Spencer Graves
Using RStudio 1.0.136 with
> sessionInfo()
R version 3.3.3 (2017-03-06)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierra 10.12.3
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets
[6] methods base
loaded via a namespace (and not attached):
[1] tools_3.3.3
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Re: [Rd] Bug in nlm()
{This was sent to me, MM, only, but for completeness should
have gone back to R-devel.
Further: I now *have* added Marie B to the members'of "R bugzilla"
-- M.Maechler}
I had already read the R bug reporting guide and I'm sure it is a bug.
The bug occurs when the user provides not only the analytic gradient but also
the analytic Hessian of the objective function. In that case, the algorithm
does not converge due to an erroneous implementation of the modified Cholesky
decomposition of the Hessian matrix. It is actually a bug in the C-code called
by nlm(), therefore it is hard to show that the non-convergence of the
algorithm is really due to this bug with only a MRE.
However, a short example (optimizing the Rosenbrock banana valley function with
and without analytic Hessian) is:
fg <- function(x){
gr <- function(x1, x2) c(-400*x1*(x2 - x1*x1)-2*(1-x1), 200*(x2 - x1*x1))
x1 <- x[1]; x2 <- x[2]
res<- 100*(x2 - x1*x1)^2 + (1-x1)^2
attr(res, "gradient") <- gr(x1, x2)
return(res)
}
nlm.fg <- nlm(fg, c(-1.2, 1))
fgh <- function(x){
gr <- function(x1, x2) c(-400*x1*(x2 - x1*x1) - 2*(1-x1), 200*(x2 - x1*x1))
h <- function(x1, x2){
a11 <- 2 - 400*x2 + 1200*x1*x1
a21 <- -400*x1
matrix(c(a11, a21, a21, 200), 2, 2)
}
x1 <- x[1]; x2 <- x[2]; res<- 100*(x2 - x1*x1)^2 + (1-x1)^2
attr(res, "gradient") <- gr(x1, x2)
attr(res, "hessian") <- h(x1, x2)
return(res)
}
nlm.fgh <- nlm(fgh, c(-1.2,1))
I have almost finished a more detailed bug report, which I would like to submit.
Best,
Marie Boehnstedt
> Martin Maechler
> on Fri, 3 Mar 2017 18:15:47 +0100 writes:
> Boehnstedt, Marie
> on Fri, 3 Mar 2017 10:23:12 + writes:
>> Dear all, I have found a bug in nlm() and would like to
>> submit a report on this. Since nlm() is in the
>> stats-package, which is maintained by the R Core team,
>> bug reports should be submitted to R's Bugzilla. However,
>> I'm not a member of Bugzilla. Could anyone be so kind to
>> add me to R's Bugzilla members or let me know to whom I
>> should send the bug report?
> Dear Marie,
> I can do this ... but are you really sure? There is
> https://www.r-project.org/bugs.html which you should spend
> some time reading if you haven't already.
> I think you would post a MRE (Minimal Reproducible
> Example) here {or on stackoverflow or ...} if you'd follow
> what the 'R bugs' web page (above) recommends and only
> report a bug after some feedback from "the public".
> Of course, I could be wrong.. and happy if you explain /
> tell me why.
> Best, Martin Maechler
>> Thank you in advance.
>> Kind regards, Marie B�hnstedt
>> Marie B�hnstedt, MSc Research Scientist Max Planck
>> Institute for Demographic Research Konrad-Zuse-Str. 1,
>> 18057 Rostock, Germany
>> www.demogr.mpg.de
