search for: campden

...andbook/pri/section5/pri54.htm>. I can't believe that R doesn't have facilities to design and analyse such experiments, but I haven't been able to find them (I have looked quite hard!). Can anyone point me in the right direction? Thanks in advance, Keith Jewell mailto:k.jewell at campden.co.uk telephone (direct) +44 (0)1386 842055 Released by Mr. K. Jewell _____________________________________________________________________ The information in this document and attachments is given after the exercise of all reasonable care and skill in its compilation, preparation and issue, b...

...path(.libPaths(), "*/etc/menus.txt")) fails to find files which are there b) update.packages(ask='graphics') fails when copying the updates into the destination folders In Renviron.site I define the site library with forward slashes, not backslashes thus... R_LIBS_SITE=//campden/shares/workgroup/stats/R/library/%v ... but the startup process seems to replace them with forward slashes. I guess because .libPaths with a 'new' argument calls normalizePath which changes leading slashes to backslashes, even with winslash="/" > normalizePath("//campden...

...log, rel.tol = rel.tol, ...)$value }, x, meanlog, sdlog, ... ) } --------------------------------------------------------------- Best regards, Keith Jewell mailto:k.jewell at campden.co.uk telephone (direct) +44 (0)1386 842055 > -----Original Message----- > From: Jewell, Keith > Sent: 15 February 2008 16:48 > To: 'r-help at r-project.org' > Subject: Behaviour of integrate (was 'Poisson-lognormal probability > calculations') > > Hi again...

...;t fix it; if the factorial (=gamma) is inside the integrand I get silly answers, if it's outside I get non-finite integrand. I'm tearing my hair. Can anyone suggest where I may be going wrong? Any suggestions at all will be appreciated. Thanks in advance, Keith Jewell mailto:k.jewell at campden.co.uk telephone (direct) +44 (0)1386 842055 _____________________________________________________________________ The information in this document and attachments is given after the exercise of all reasonable care and skill in its compilation, preparation and issue, but is provided without liab...

...it, and I still get the marked difference between the presence and absence of those lines (and I'm increasingly unsure whether either answer is correct!). I'm still trying. but I really think I'm going nowhere. Has anyone any ideas? Thanks in advance, Keith Jewell mailto:k.jewell at campden.co.uk telephone (direct) +44 (0)1386 842055 > -----Original Message----- > From: Jewell, Keith > Sent: 15 February 2008 11:16 > To: 'r-help at r-project.org' > Subject: Poisson-lognormal probability calculations > > Hi, > > just for the record, although I don...

Hi All, This example code ---------------- dDF <- structure(list(y = c(4.75587, 4.8451, 5.04139, 4.85733, 5.20412, 5.92428, 5.69897, 4.78958, 4, 4), t = c(0, 48, 144, 192, 240, 312, 360, 0, 48, 144), Batch = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ), T = c(2, 2, 2, 2, 2, 2, 2, 2, 2, 2), pH = c(4.6, 4.6, 4.6, 4.6, 4.6, 4.6, 4.6, 4.6, 4.6, 4.6), S = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0), N = c(0, 0, 0, 0,

...n a known or unknown LOD? -- and/or is there positive mass on zero? In either case, using logs to calculate a geometric mean may not make sense. Paraphrasing Greg Snow, what is the scientific question? What is the model? Cheers, Bert On Mon, Jan 17, 2011 at 9:13 AM, Keith Jewell <k.jewell at campden.co.uk> wrote: > Just in case some of x are negative (the desired median still exists, as > long as the two middle values are non -ve), how about: > > x <- runif(20, -1, 100) > exp(median(log(pmax(0,x)))) > > It'll give -Inf if the two middle values are negative, when...

Dear Keith, >>>>> <Keith.Jewell at campdenbri.co.uk> >>>>> on Thu, 16 Jul 2015 08:58:11 +0000 writes: > Dear R Core Team, > Last week I made a post to the R-help mailing list > ?predict.poly for multivariate data? > <https://stat.ethz.ch/pipermail/r-help/2015-July/430311.html> &gt...

...)) to anyDuplicated(); thanks Benjamin Tyler * getRs: new function to interact with https://github.com/harrelfe/rscripts * knitrSet: new function to setup knitr with nice defaults for books etc. * rcorr: fixed sensing of NAs and diagonal elements of n matrix; thanks: Keith Jewell, Campden BRI Group; similar for hoeffd -- ------------------------------------------------------------------------ Frank E Harrell Jr Professor and Chairman School of Medicine Department of *Biostatistics* *Vanderbilt University* [[alternative HTML version deleted]] ______________________________...

...)) to anyDuplicated(); thanks Benjamin Tyler * getRs: new function to interact with https://github.com/harrelfe/rscripts * knitrSet: new function to setup knitr with nice defaults for books etc. * rcorr: fixed sensing of NAs and diagonal elements of n matrix; thanks: Keith Jewell, Campden BRI Group; similar for hoeffd -- ------------------------------------------------------------------------ Frank E Harrell Jr Professor and Chairman School of Medicine Department of *Biostatistics* *Vanderbilt University* [[alternative HTML version deleted]] ______________________________...

...h is not a successive integers sequence") Lval <- list(xy$y) names(Lval) <- mCall["Coeffs"] Lval } , parameters = c("Coeffs") ) Sorry for wasting anyones time. Keith Jewell ------------------------------------------------- "Keith Jewell" <k.jewell at campden.co.uk> wrote in message news:... > Hi everyone, > > I thought that for a selfStart function, these two should be exactly > equivalent >> nls(Aform, DF) >> nls(Aform, DF, start=getInitial(Aform, DF)) > but in this example that is not the case in R (although it is in S-...

Hi, I'd like to know which functions in a package call one specific function. I think I've seen a tool for identifying such dependencies, but now I can't find it :-( Searches of help and R site search for keywords like function, call, tree, depend haven't helped :-( Can anyone point me in the right direction? Thanks in advance, Keith Jewell

Hi, I seek to identify those points in/outside a multidimensional convex hull (geometry::convhulln). Any suggestions? Background just in case I'm going down a really wrong road: Given an observed data set with one dependent/observed variable (Y) and multiple (3 to 10) independent/design variables (X1, X2, ...) I want to increase the number of points by interpolating. I'm using

Hi, > update.packages(ask='graphics') gives me multiple warning (one per updated package?) similar to ... Warning: unable to move temporary installation '\\Server02\stats\R\library\2.10\file3de56e0d\locfit' to '\\Server02\stats\R\library\2.10\locfit' The final, updated, folders do not end up where they should be. I can move them 'by hand', but it is an

I'm working with 3 different data sets and applying this non-linear regression formula to each of them. nls(Y ~ (upper)/(1+10^(X-LOGEC50)), data=std_no_outliers, start=list(upper=max(std_no_outliers$Y),LOGEC50=-8.5)) Previously, all of the regressions were calculated in Prism, but I'd like to be able to automate the calculation process in a script, which is why I'm trying to move to

Hello This request asks something beyond the technicalities of the R language, I would like to ask you wonderful people if you have ever suffered as programmers ( or de facto programmers like myself though I am a 'research assistant') from tendonitis and how you coped with it, i have golfer's elbow on both sides. Any resources? Pancho Mulongeni Namibia

Hi All, Another nls related problem (for background, I'm migrating a complicated modelling package from S-plus to R). Below I've reduced this to the minimum necessary to demonstrate my problem (I think); the real situation is more complicated. Two similar selfStart functions, ssA and ssB. The 'initial' function for ssB modifies its arguments a little and then calls getInital

Hi Everyone, I posted this a couple of weeks ago with no responses. My interface (via gmane) seemed a bit flakey at the time, so I'm venturing to repost with some additional information. I'm trying to write selfStart non-linear models for use with nls. In these models some combinations of parameter values are illegal; the function value is undefined. That's OK when calling the

...help > > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. ------------------------------ Message: 47 Date: Mon, 4 Aug 2008 09:58:14 +0100 From: "Keith Jewell" <k.jewell at campden.co.uk> Subject: Re: [R] Unexpected nls behaviour: Solved To: r-help at stat.math.ethz.ch Message-ID: <g76gbc$h01$1 at ger.gmane.org> Hi Everyone, I'd omitted the non-optional 'parameters' argument to selfStart. Making this change to SSbatch gives the same (successful) result...

...4 1865 272595 > > > > ------------------------------ > > Message: 2 > Date: Sat, 10 Jan 2009 08:42:46 -0500 > From: David Winsemius <dwinsemius@comcast.net> > Subject: Re: [R] rpart with interval censored data crashes R > To: "Keith Jewell" <k.jewell@campden.co.uk> > Cc: r-help@stat.math.ethz.ch > Message-ID: <5467EAED-927D-4C5C-A656-FB1D510C63BC@comcast.net> > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > On a Leopard Mac with the Urbanek compiled 64 bit R, one sees this: > > > library(rpart...