search for: lgcgroup

> On 10 Jan 2013, at 15:56 , S Ellison <S.Ellison at lgcgroup.com> wrote: > > > >> I am working with large numbers and identified that R looses >> precision for such high numbers. > Yes. R uses standard 32-bit double precision. Well, for large values of 32... such as 64. -- Peter Dalgaard, Professor, Center for Statistics, C...

...] "p" I think that's what I fall over mostly: that named, empty arguments behave entirely different from omitting them (", ,") And I definitely agree we need a guru to explain it all to us ( Cheers, Emil Bode ?On 30/11/2018, 15:35, "S Ellison" <S.Ellison at LGCGroup.com> wrote: > Yes, I think all of that is correct. But y _is_ missing in this sense: > > plot(1:10, y=) > > ... > Browse[2]> missing(y) Although I said what I meant by 'missing' vs 'not present', it wasn't exactly what missing()...

On Thu, Nov 29, 2018 at 1:10 PM S Ellison <S.Ellison at lgcgroup.com> wrote: > > > > > plot(x=1:10, y=) > > > plot(x=1:10, y=, 10:1) > > > > > > In both cases, 'y=' is ignored. In the first, the plot is for y=NULL (so not > > 'missing' y) > > > In the second case, 10:1 is positionally ma...

...*********************************** This email and any attachments are confidential. Any use, copying or disclosure other than by the intended recipient is unauthorised. If you have received this message in error, please notify the sender immediately via +44(0)20 8943 7000 or notify postmaster at lgcgroup.com and delete this message and any copies from your computer and network. LGC Limited. Registered in England 2991879. Registered office: Queens Road, Teddington, Middlesex, TW11 0LY, UK

I'm looking for a histogram variant in which data points are plotted as labelled rectangles 'piled up' to form a histogram. I've posted an example at https://www.dropbox.com/s/ozi8bhdn5kqaufm/labelled_histogram.png?dl=0 It seems to have a long pedigree, as I see it (as in this example) in documents going back beyond the '80s. But I've not seen it in recent textbooks. So it

...List of 3 $ : symbol ? $ : symbol a $ : language b <- c > str(as.list(parse(text="a <- b ? c")[[1]])) List of 3 $ : symbol ? $ : language a <- b $ : symbol c Bill Dunlap TIBCO Software wdunlap tibco.com On Fri, Aug 30, 2019 at 4:41 AM Stephen Ellison <S.Ellison at lgcgroup.com> wrote: > > From: R-devel [mailto:r-devel-bounces at r-project.org] On Behalf Of Ant F > > Sent: 29 August 2019 12:06 > > To: r-devel at r-project.org > > Subject: [Rd] ?Syntax wrong about `?`'s precedence ? > > ... > > See the following example : &gt...

...uot;, "var2", "mortality_probability") X[, "mortality_probability"] <- qpois(X[, "mortality_probability"], 0.9) hist(X[, "mortality_probability"]) Thanks for your time Marine ________________________________ De : S Ellison <S.Ellison at LGCGroup.com> Envoy? : lundi 7 ao?t 2017 14:36 ? : Marine Regis; r-help at r-project.org Objet : RE: Latin hypercube sampling from a non-uniform distribution > How can I draw a Hypercube sample for the variable mortality_probability so > that this variable exhibits the same pattern as the observ...

...39;s what I fall over mostly: that named, empty arguments behave entirely different from omitting them (", ,") > > And I definitely agree we need a guru to explain it all to us ( > > Cheers, Emil Bode > > > ?On 30/11/2018, 15:35, "S Ellison" <S.Ellison at LGCGroup.com> wrote: > > > Yes, I think all of that is correct. But y _is_ missing in this sense: > > > plot(1:10, y=) > > > ... > > Browse[2]> missing(y) > > Although I said what I meant by 'missing' vs 'not present', it wasn&...

...the best option. That doesn?t rule out exceptions. Best, John ------------------------------------------------- John Fox, Professor Emeritus McMaster University Hamilton, Ontario, Canada Web: http::/socserv.mcmaster.ca/jfox > On Dec 17, 2018, at 7:49 AM, S Ellison <S.Ellison at LGCGroup.com> wrote: > > > >> From: Thomas Yee [mailto:t.yee at auckland.ac.nz] >> >> Thanks for the discussion. I do feel quite strongly that >> the variables should always be a part of a data frame. > > This seems pretty much a decision for R core, and I thi...

...*********************************** This email and any attachments are confidential. Any use, copying or disclosure other than by the intended recipient is unauthorised. If you have received this message in error, please notify the sender immediately via +44(0)20 8943 7000 or notify postmaster at lgcgroup.com and delete this message and any copies from your computer and network. LGC Limited. Registered in England 2991879. Registered office: Queens Road, Teddington, Middlesex, TW11 0LY, UK

...ead. Though wondering how 53 bits were supposed to fit into 32 might just warrant revivification. -- Sent from my phone. Please excuse my brevity. On July 20, 2017 5:33:34 AM PDT, peter dalgaard <pdalgd at gmail.com> wrote: > >> On 10 Jan 2013, at 15:56 , S Ellison <S.Ellison at lgcgroup.com> wrote: >> >> >> >>> I am working with large numbers and identified that R looses >>> precision for such high numbers. >> Yes. R uses standard 32-bit double precision. > > >Well, for large values of 32... such as 64.

> However, my variable is simulated from the cumulative distribution function > of the Poisson distribution. Then I am afraid I don't know what you're trying to achieve. Or why. However, the principle holds; write a function that maps [0,1] to the 'pattern' you want, do that and apply it to the result from randomLHS. It happens that for generating data that follow a given

Hello, I am performing a sensitivity analysis using a Latin Hypercube sampling. However, I have difficulty to draw a Hypercube sample for one variable. I?ve generated this variable from a Poisson distribution as follows: set.seed(5) mortality_probability <- round(ppois(seq(0, 7, by = 1), lambda = 0.9), 2) barplot(mortality_probability, names.arg = seq(0, 7, by = 1), xlab = "Age

...ice. Beginners ought to >>> learn to keep data in data frames and not to overuse attach(). Experts can >>> do otherwise at their own risk, but they have less need of explicit >>> examples. >>> >>> On Fri, 14 Dec 2018 at 14:51, S Ellison <S.Ellison at lgcgroup.com> wrote: >>> >>>> FWIW, before all the examples are changed to data frame variants, I think >>>> there's fairly good reason to have at least _one_ example that does _not_ >>>> place variables in a data frame. >>>> >>>>...

On Thu, Nov 29, 2018 at 10:51 AM S Ellison <S.Ellison at lgcgroup.com> wrote: > > > When trying out some variations with `[.data.frame` I noticed some (to me) > > odd behaviour, > > Not just in 'myfun' ... > > plot(x=1:10, y=) > plot(x=1:10, y=, 10:1) > > In both cases, 'y=' is ignored. In the first, the plot...

...close and built in. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Jun 21, 2017 at 9:01 PM, S Ellison <S.Ellison at lgcgroup.com> wrote: > > I'm looking for a histogram variant in which data points are plotted as labelled rectangles 'piled up' to form a histogram. I've posted an example at https://www.dropbox.com/s/ozi8bhdn5kqaufm/labelled_histogram.png?dl=0 > > It seems to have a long pedi...

...bability"], 0.9)?, I don?t obtain a LHS value for each age class and the distribution is not a cumulative distribution as in the attached figure. So, I am afraid I don?t know how to do that. Thanks so much for your time Marine ________________________________ De : S Ellison <S.Ellison at LGCGroup.com> Envoy? : mardi 8 ao?t 2017 14:48 ? : Marine Regis; r-help at r-project.org Objet : RE: Latin hypercube sampling from a non-uniform distribution > However, my variable is simulated from the cumulative distribution function > of the Poisson distribution. Then I am afraid I don't kn...

Hi, I have one basic doubt. Suppose X ~ N(50,10). I need to calculate Probability X = 50. dnorm(50, 50, 10) gives me [1] 0.03989423 My understanding is (which is bit statistical or may be mathematical) on a continuous scale, Probability of the type P(X = .....) are nothing but 1/Infinity i.e. = 0. So as per my understanding P(X = 50) should be 0, but even excel also gives 0.03989422. Obviously

...;a <- b ? c")[[1]])) >> List of 3 >> $ : symbol ? >> $ : language a <- b >> $ : symbol c >> >> Bill Dunlap >> TIBCO Software >> wdunlap tibco.com >> >> >> On Fri, Aug 30, 2019 at 4:41 AM Stephen Ellison <S.Ellison at lgcgroup.com> >> wrote: >> >>>> From: R-devel [mailto:r-devel-bounces at r-project.org] On Behalf Of Ant F >>>> Sent: 29 August 2019 12:06 >>>> To: r-devel at r-project.org >>>> Subject: [Rd] ?Syntax wrong about `?`'s precedence ? >>&...

FWIW, before all the examples are changed to data frame variants, I think there's fairly good reason to have at least _one_ example that does _not_ place variables in a data frame. The data argument in lm() is optional. And there is more than one way to manage data in a project. I personally don't much like lots of stray variables lurking about, but if those are the only variables out