search for: turlach

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 files, but the options seem numerous. I haven't...

...st a warning that there are new updates available, just no new binaries. Best, Philipp ? ? ? Gesendet:?Mittwoch, 14. Februar 2024 um 18:44 Uhr Von:?"Duncan Murdoch" <murdoch.duncan at gmail.com> An:?"Martin Maechler" <maechler at stat.math.ethz.ch>, "Berwin A Turlach" <berwin.turlach at gmail.com> Cc:?"gernophil--- via R-help" <r-help at r-project.org> Betreff:?Re: [R] Packages sometimes don't update, but no error or warning is thrown On 14/02/2024 5:50 a.m., Martin Maechler wrote: >>>>>> Berwin A Turlach >&gt...

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 programming problem : min(c'*x + 1/2 * x' * H * x) subject to: b <= A * x <= b + r l <= x <= u But my pro...

Hi I want to write a little function that takes a matrix X of size m-by-n, and a list L of length "m", whose elements are matrices all of which have the same number of columns but possibly a different number of rows. I then want to get a sort of dumbed-down kronecker product in which X[i,j] is replaced by X[i,j]*L[[j]] where L[[j]] is the j-th of the "m" matrices. For

Hi look at the following session, in which I have a dataframe, and I want to extract the second row, without the first column. Everything works as expected until the last line, where I set the names of x to NULL, and get a non-desired object (I want c(4,3).). Three questions: (1) why is as.vector(a[2,-1]) not a vector? (2) How come setting names to NULL gives me bad weirdness? (3) Can I

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 > > Sweave vignette source files in @file{vignettes}. If > &gt...

I have found a strange scoping problem in Sweave. The following Rnw file doesn't produce the same output in Sweave as it does if I produce an R file using Stangle and execute that: \documentclass[12pt]{article} \begin{document} <<R>>= election <- data.frame(A=1:3, B=9:7, C=rep(0,3)) partytotal <- rep(0, ncol(election)) for (i in 1:ncol(election)) { partytotal[i] <-

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 development (unstable) (2020-01-25 r77715)" on 27 January. Since then the script always fails as a regression test seems to fail.

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 Sweave vignette source files in @file{vignettes}. If @file{Makefile} is found in the vignette source directory, then @code{R CMD build} will try to run @command{make} after the @code{Sweave} runs, otherwise @code{texi2pdf} is run on

Hi, Forgive a very basic question... I need to take two lists-of-lists, and apply a function to each pair of elements in the lists to return a single list... For example l1 <- list(1:5,6:10,2:15) l2 <- list(1:8,4:12,1:19,4:20) I could easily do an lapply across each of them, but is there a function that does a sort-of pairwise-apply across both together? Does anybody know of a good

.... I hope that's OK.) Regards, John -------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -------------------------------- > -----Original Message----- > From: Berwin A Turlach > [mailto:berwin at bossiaea.maths.uwa.edu.au] On Behalf Of Berwin > A Turlach > Sent: Tuesday, July 18, 2006 9:28 PM > To: Andrew Robinson > Cc: Ulrich Keller; John Fox > Subject: Re: [R] Test for equality of coefficients in > multivariate multipleregression > > G...

G'day Federico, On Wed, 17 May 2023 10:42:17 +0000 "Calboli Federico (LUKE)" <federico.calboli at luke.fi> wrote: > sexsnp = rep(NA, 1750) > for(i in 1:1750){tryCatch(sexsnp[i] = fisher.test(table(data[,3], > data[,i + 38]))$p, error = function(e) print(NA))} Error: unexpected > '=' in "for(i in 1:1750){tryCatch(sexsnp[i] =" Try: R> for(i in

...ata[,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 8:22?AM Berwin A Turlach <berwin.turlach at gmail.com> wrote: > > G'day Federico, > > On Wed, 17 May 2023 10:42:17 +0000 > "Calboli Federico (LUKE)" <federico.calboli at luke.fi> wrote: > > > sexsnp = rep(NA, 1750) > > for(i in 1:1750){tryCatch(sexsnp[i] = fisher.test...

G'day Thomas, On Sat, 12 Aug 2023 04:17:42 +0000 (UTC) Thomas Subia via R-help <r-help at r-project.org> wrote: > Here is my reproducible code for a graph using geom_smooth The call "library(tidyverse)" was missing. :) > I'd like to add a black boundary around the shaded area. I suspect > this can be done with geom_ribbon but I cannot figure this out. Some >

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 (8) 9380 3383 (self) The University of Western Australia FAX : +61 (8) 9380 1028 35 Stirling Highway Crawley WA 6009 e-mail: berwin at maths.uwa.edu.au Au...

I need to make a bunch of PDF files of histograms. I tried gatelist = unique(mdf$ArrivalGate) for( gate in gatelist) { outfile = paste("../", airport, "/", airport, "taxiHistogram", gate, ".pdf", sep="") pdf(file = outfile, width = 10, height=8, par(lwd=1)) title=paste("Taxi time for Arrival Gate", gate, "by

Hi, R users, Here is an example. k <- c(1,2,3,4,5) i <- c(0,1,3,2,1) if k=1, then j=0 from i if k=2, then j=0, 1 from i if k=3, then j=0, 1, 2, 3 from i if k=4, then j=0, 1, 2 from i if k=5, then j=0, 1 from i so i'd like to create a list like below. > list k j 1 1 0 2 2 0 3 2 1 4 3 0 5 3 1 6 3 2 7 3 3 8 4 0 9 4 1 10 4 2 11 5 0 12 5 1 I tried expand.grid, but I

...uality 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 version of quadprog, Version: 1.4-11, Date: 2007-07-12 and my problem is solved correctly. I have tried to contact Berwin A. Turlach <Berwin.Turlach at gmail.com> (maintainer for quadprog package) a week ago, with no success. ############################################################## load(file='quadprog.Rdata') # solve QP using quadprog require(quadprog) sol = solve.QP(Dmat, dvec, Amat, bvec, meq) x...

...rectories as soon as they are no longer needed and/or that 'R --arch=name CMD INSTALL --libs-only' installs only libraries for the specified architecture (as an unwary user might expect). Cheers, Berwin ========================== Full address ============================ A/Prof Berwin A Turlach Tel.: +61 (8) 6488 3338 (secr) School of Maths and Stats (M019) +61 (8) 6488 3383 (self) The University of Western Australia FAX : +61 (8) 6488 1028 35 Stirling Highway Crawley WA 6009 e-mail: Berwin.Turlach at gmail.com Australia...

hello I tried to use lars() but neither with my own data nor with the sample data it works. I get in both cases the following error prompt: > data(diabetes) > par(mfrow=c(2,2)) > attach(diabetes) > x<-lars(x,y) Error in one %*% x : requires numeric matrix/vector arguments > x<-lars(x,y, type="lasso") Error in one %*% x : requires numeric matrix/vector arguments