[Rd] print.citation, small bug?

e goat were married and lived together until the end of their days. Either the goat went mad or the bear became sane." Nicholas Lewin-Koh Genentech [[alternative HTML version deleted]] __ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel

Re: [Rd] raster support in graphics devices

Hi Charlotte, I think I would have to respectfully disagree with you on many points. While we all sigh when we see heatmaps used as data to explain everything, they do have their uses. But beyond heat maps, the nature and demands of modern statistical graphics has changed. Plotting huge data sets r

Re: [Rd] buglet in is.na?

, I have a member of an object that can be an integer, a function returning an integer or NA, perhaps it would be better to specify it as an integer, a function or NULL. Nicholas On Wed, 30 Sep 2009 14:41 -0400, "Duncan Murdoch" wrote: > On 9/30/2009 2:29 PM, Nicholas Lewin-Koh

[Rd] buglet in is.na?

Hi the following example I think demonstrates the inconsistency > f<-function(x) x > length(f) [1] 1 > is.na(f) [1] FALSE Warning message: In is.na(f) : is.na() applied to non-(list or vector) of type 'closure' The documentation states: Arguments x an R object to be tested. value a suitabl

Re: [Rd] Non-GPL packages for R

With header this time On Thu, 10 Sep 2009 16:08 -0700, "Nicholas Lewin-Koh" wrote: > Hi, > I would also be in favor of a stronger stance on licenses. In > industry, where we can really get in big trouble for violating a > license, > we often maintain internal repositori

[Rd] (no subject)

Hi, I would also be in favor of a stronger stance on licenses. In industry, where we can really get in big trouble for violating a license, we often maintain internal repositories, or need to be careful about filtering what is used from CRAN. I think that is should really be a requirement the packa

Re: [Rd] proposed simulate.glm method

ors from a zipf distribution, and I did have a predict method (which is a fair assumption) would that fail because the rzipf function would not be accessed? Nicholas On Sat, 14 Feb 2009 20:10:26 +0100, "Martin Maechler" said: > >>>>> "NicLK" == Nicholas Lew

Re: [Rd] proposed simulate.glm method

Hi, For extended glms such as gams, gnm or other distributions such as negative binomial, would there need to be a separate simulate method? Or, could the current framework, rather than stopping with an error look for the appropriate model matrix, coefficients, distribution function and family obje

Re: [Rd] NY Times article

Hi Frank, Thank you for the correction, I don't work in clinical, so good you keep me honest. Nicholas On Sun, 11 Jan 2009 09:21:23 -0600, "Frank E Harrell Jr" said: > Nicholas please note that FDA does not require type III sums of squares, > LOCF, and the use of any particular software. It is j

Re: [Rd] NY Times article

t; lists, where a common refrain in response to bug reports is that "you > aren't using the latest and greatest version", without an appreciation > of the fact that there can be good reasons for not changing something > that you know works). > > Best, > K

Re: [Rd] NY Times article

Hi, Kudos, nice exposure, but to make this more appropriate to R-devel I would just like to make a small comment about the point made by the SAS executive about getting on an airplane yada yada ... 1) It would seem to me that R has certification documents 2) anyone designing airplanes, analyzing c

Re: [Rd] Unix-like touch to update modification timestamp of file?

Wow, this has to win the prize for some of the most obscure documentation ever. Kudos on your M$-archaeology. Nicholas Message: 9 Date: Wed, 27 Feb 2008 15:23:51 -0500 From: "Gabor Grothendieck" <[EMAIL PROTECTED]> Subject: Re: [Rd] Unix-like touch to update modification timestamp of file?

Re: [Rd] [R] color palette from red to blue passing white (shifted from R-help)

Hi, The move to sRGB is nice, is there any interest in adding an interface to lcms, http://www.littlecms.com, to allow gamut matching? I can think of a lot of instances where I would like to render a figure as it would appear on my printer. This is probably best done as a separate package though, a

Re: [Rd] nmle: gnls freezes on difficult case

sentence to "optim method 'L-BFGS-B' will not work since there is way to specify the upper and lower constraints" Thanks Nicholas Nicholas On Thu, 18 Oct 2007 15:46:01 +0100 (BST), "Prof Brian Ripley" <[EMAIL PROTECTED]> said: > On Thu, 18 Oct 2007, Nicholas

[Rd] nmle: gnls freezes on difficult case

Hi, I was just able to try this on a windows xp machine using R-2.5.1. gnls exits gracefully with an error so maybe this is a problem with my R installation. If so sorry to bother everyone. Does anyone have an idea which installation parameters to tweek that might cure this? Thanks Nicholas ___

[Rd] nmle: gnls freezes on difficult case

Hi, Following up on my own post, if in gnlsControl I specify opt='optim' gnls exits quite nicely, with an error, which is what I would expect. Is this a bug in nlminb? Nicholas tt<-gnls(response~SSllogis(conc,A,B,xmid,scal),tdat,start=start,weights=varPower(),verbose=TRUE,control=gnlsControl(opt

[Rd] nmle: gnls freezes on difficult case

Hi, I am not sure this is a bug but I can repeat it, The functions and data are below. I know this is nasty data, and it is very questionable whether a 4pl model is appropriate, but it is data fed to an automated tool and I would have hoped for an error. Does this repeat for anyone else? My details

[Rd] Optimization in R

Hi, In my earlier post I eluded to a situation where that would be useful. In nlme, there is a choice of optimizers, minpack.lm has Levenberg-Marquardt, while nlminb has the port routines. For the same starting values, different optimizers will present different solutions, having a common interfac

[Rd] Optimization in R

Hi Andrew, I have been working quite a bit with optim and friends on automated nonlinear fitting, mainly for calibration. All of the optimizers seem to have trouble in the tricky situation of fitting a log-logistic model when the upper asymptote is not well defined, and trying to estimate a vari

[Rd] Coloring leaves, twigs and labels in plot.dendrogram()

Hi Walton, Sorry I don't follow R-help very much, but here is a nice solution to what you are looking for in the A2R package. I don't think it is on cran yet. http://addictedtor.free.fr/graphiques/graphcode.php?graph=79 Nicholas > > Message: 3 > Date: Sat, 22 Oct 2005 20:12:52 -0400 (EDT) > Fr

Re: [Rd] evaluating variance functions in nlme

Those results plus "str(fm1)" looked like they might help answer your > question. However, I don't understand enough of what you are asking to > say. > > If an answer might still be worth pursuing to you, I suggest you read > the posting guide and submit

[Rd] evaluating variance functions in nlme

Hi, I guess this is a final plea, and maybe this should go to R-help but here goes. I am writing a set of functions for calibration and prediction, and to calculate standard errors and intervals I need the variance function to be evaluated at new prediction points. So for instance fit<-gnls(Y~SSl

Re: [Rd] extracting and manipulating components from nlme and gnls objects

at x0 using predict(fit, newdata=list(covar=x0)) to get the fitted value but how to I plug that in to the varStruct stored in fit$modelStruct? Thanks Nicholas On Tue, 19 Jul 2005 11:48:18 -0500, "Douglas Bates" <[EMAIL PROTECTED]> said: > On 7/19/05, Nicholas Lewin-Koh <[EMAIL P

[Rd] extracting and manipulating components from nlme and gnls objects

Hello, I am writing a set of functions to do prediction and calibration intervals for at least the subset of selfstarting models if not some more general ones. I need to be able to extract the varFunction from a fit object and evaluate it at a predicted point. Are there any examples around? Also