similar to: Summary --- testing for existence of file.

Hey everyone, Thanks for all the input. It's happening again. This time for the packages "DBI", "parallelly", "segmented", "survival", "V8". So, RStudio shows updates for those and updating them via RStudio leads to this output: ``` > install.packages(c("DBI", "parallelly", "segmented", "survival",

Dear Berwin, Simply stacking the problems and treating the resulting observations as independent will give you the correct coefficients, but incorrect coefficient variances and artificially zero covariances. The approach that I suggested originally -- testing a linear hypothesis using the coefficient estimates and covariances from the multivariate linear model -- seems simple enough. For

Dear all, I noticed the following bug in points.formula > library(DAAG) > data(roller) > fm <- lm(depression ~ weight, data=roller) > plot( depression ~ weight, data=roller, type="n") > abline(fm) > attach(roller) > points( depression~weight, subset=1:7) > points( depression~weight, subset=8:10, col="blue") Error in if (length(x) == l) x[s] else x :

--fupGvOGOQM Content-Type: text/plain; charset=us-ascii Content-Description: message body and .signature Content-Transfer-Encoding: 7bit There's a simple bug in the handling of the start and offset arguments in glm and glm.fit. The bug exists in the latest development version of R (version information below), but it appears that glm.R has not been touched much lately, so the bug affects at

Forgot to send one copy to R help. Sorry Megh Dal <megh700004@yahoo.com> wrote: Date: Wed, 7 May 2008 02:45:09 -0700 (PDT) From: Megh Dal <megh700004@yahoo.com> Subject: Re: [R] Solution of function To: Berwin A Turlach <berwin@maths.uwa.edu.au> Hi Berwin, Thanks for having look on my problem. However on ipop() function I see following: ipop solves the quadratic

... or just put the R expression inside curly brackets, e.g. tryCatch({ sexsnp[i] = fisher.test(table(data[,3], data[,i+38]))$p }, error=function(e) print(NA)) Exercise: Compare > list(a = 2) $a [1] 2 with > list({ a = 2 }) [[1]] [1] 2 and > list(b = { a = 2 }) $b [1] 2 BTW, note how the latter two assigned a <- 2 to the global environment. /Henrik On Thu, May 18, 2023 at

Hi all, people who don't follow Linux Today regularly may want to check out: http://linuxtoday.com/developer/2003082000626OSSVDV My apologies if this is considered spam. Cheers, Berwin ========================== Full address ============================ Berwin A Turlach Tel.: +61 (8) 9380 3338 (secr) School of Mathematics and Statistics +61

G'day Duncan, On Tue, 12 Jan 2016 07:32:05 -0500 Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > On 11/01/2016 11:59 PM, Berwin A Turlach wrote: > > G'day all, > > > > In Chapter 1.4 (Writing package vignettes) the Writing R Extensions > > manual states: > > > > By default @code{R CMD build} will run @code{Sweave} on all > >

G'day all, I guess it is still early enough in the year to wish everybody a happy and successful new year. I thought I should report that the installation of the CRAN package rstan regularly fails on my machine (a 64 bit linux box running Xubuntu 15.10). The reason being that I have the 32-bit and the 64-bit architecture of R installed, and my /tmp file is on a partition with about 1Gb

Dear all, I would like to suggest changing the pairs.formula command such that a command like pairs(GNP ~ . - Year - GNP.deflator, longley) would behave in a similar fashion as lm(GNP ~ . - Year - GNP.deflator, longley) i.e., make a pairwise scatterplot of GNP and all other variables in the (longley) dataframe except for Year and GNP.deflator. The above command, with the

Dear R developers, is there a preferred format or strategy for making a patch to contribute to a package that is maintained by R-core? Berwin Turlach and I have written a very minor extension to lmeControl to allow it to pass an argument to nlminb for the maximum number of evaluations of the objective function. I've edited the nlme/R/lme.R and nlme/man/lmeControl.Rd files. I can diff the

On 17.04.2016 11:01, Prof Brian Ripley wrote: > On 17/04/2016 07:25, Berwin A Turlach wrote: >> G'day all, >> >> probably you have noticed this by now, but I thought I ought to report >> it. :) > > Already fixed for Unix by the time this reached me. Since that version > of Survival has been put into 3.2 patched, that also needed its > Makefile.in

--MP_/kvy20nVajVG/n.8m=_ZjLAX Content-Type: text/plain; charset=US-ASCII Content-Transfer-Encoding: 7bit Content-Disposition: inline On Tue, 7 Oct 2008 19:31:03 +0800 Berwin A Turlach <berwin at maths.uwa.edu.au> wrote: > The attached patch (against the current SVN version of R) implements > the latter strategy. With this patch applied, "make check > FORCE=FORCE" passes

Thanks a lot, Berwin. Unfortunately, pK1 may well be negative and as I understand the literature it may be poorly defined as such, and also seems to be at a boundary, since when lower is set to say rep(-4,3) pK1 is returned as -4 while pK2 and pK3 are undisturbed. Perhaps the point is that pK1 is not carrying any information at the pH around 5. Fair enough, I guess. Only, I believe I need

The only observation I can make is that the change to round() was made in r77727 whereas your R-devel appears to be r77715 (so would not exhibit the fixed behaviour). My guess is that there was a perpetual installation failure after r77715 but that the test folder was still retrieved and used. On Sat, 8 Feb 2020 at 19:27, Berwin A Turlach <berwin.turlach at gmail.com> wrote: > >

>>>>> Hugh Parsonage >>>>> on Sat, 8 Feb 2020 21:12:43 +1100 writes: > The only observation I can make is that the change to > round() was made in r77727 whereas your R-devel appears to > be r77715 (so would not exhibit the fixed behaviour). My > guess is that there was a perpetual installation failure > after r77715 but that

I forgot to attach the problem data, 'quadprog.Rdata' file, in my prior email. I want to report a following error with quadprog. The solve.QP function finds a solution to the problem below that violates the last equality constraint. I tried to solve the same problem using ipop from kernlab package and get the solution in which all equality constraints are enforced. I also tried an old

G'day Erin, On Sat, 16 Dec 2017 08:00:38 -0600 Erin Hodgess <erinm.hodgess at gmail.com> wrote: > I'm in the process of writing a package, and I'm using the lovely "R > Package" book as a guideline. > > However, in the midst of my work, I discovered that I had omitted a > function and am now putting in it the package. Not a problem. But > the

On Sat, Feb 8, 2020 at 9:27 AM Berwin A Turlach <berwin.turlach at gmail.com> wrote: > > G'day all, > > I have daily scripts running to install the patched version of the > current R version and the development version of R on my linux box > (Ubuntu 18.04.4 LTS). > > The last development version that was successfully compiled and > installed was "R Under

I?disagree?with?what?is?assessed?as?"correct"?in?vignette?of?package?'round'. With x?<-?9.18665 ,?what?is?actually?stored?in?'x'?is?a?number?that?is?slightly?larger?than?9.18665.?So,?as?said?in?the?vignette,?it?is?closer?to?9.1867?than?9.1866. sprintf("%.4f",?x) giving "9.1867" is?correct,?as?it?is?a?string?representing?the?_exact_?number?9.1867.