Re: [R] Comparing each component of vector to each component of a Matrix.

matrix and print 1 if the component is less or equal to one and print 0 otherwise. bi <= Aij, print 1 and 0 otherwise. Any help is appreciated. outer(b,A,"<=")+0 cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS Department of Statistics University of Auckland Phone: +6

Re: [R] .Rdata files -- fortune?

you don't understand, you shouldn't be doing it. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 Home phone: +64-9-480-4619 __ R-help@r-project.org mail

Re: [R] Add the Gauss curve on histogram

* There is no such thing as "the" Gauss curve. There is an uncountably infinite number of Gauss curves. * As Jeff Newmiller told you, *do* provide a *reproducible* example when you ask a question. (Nobody but *you* has the data set "newdata".) Neither R nor the members of t

Re: [R] Add the Gauss curve on histogram

"y", main="Simulated data") x <- seq(min(simdat),max(simdat),length=100) muhat <- mean(simdat) sigmahat <- sd(simdat) y <- dnorm(x,muhat,sigmahat) lines(x,y) # ---- Imitate it. cheers, Rolf Turner On 18/03/1

Re: [R] Returning to parent function

On 18/03/15 13:22, Bert Gunter wrote: With all due respect: I don't think much respect is due. Fortune nomination! cheers, Rolf -- Rolf Turner Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 Home phone: +64-9-480

Re: [R] Fast evaluation of functions in 3D domains

is being approximated by a step function. The method seems to work well enough for the uses that we put it to in spatstat. Whether it works well enough for your purposes depends on what your purposes are. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS Department of Statistics Uni

Re: [R] A problem someone should know about

ng the error that you saw is peculiar to your system. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 Home phone: +64-9-480-4619 __ R-help@r-project.org mailing list

Re: [R] multinomial probabilities with mlogit

Reproducible example??? cheers, Rolf Turner P.S. Where does "mlogit" come from? Note fortune(182). R. T. On 31/03/15 06:46, Ingrid Charvet wrote: Hello, When fitting a logit multinomial model with "mlogit" I can retrieve the response probabilities using fit$fitted

[R] Grouped boxplots using ggplot() from ggplot2.

but thought I'd better scrounge around a bit more before posting this, and thereby managed (mirabile dictu!) to answer them myself. Can anyone help me out with questions (1) --- (4)? Please keep it simple and very explicit, for I am a bear of very little brain and long words bother me! Tha

Re: [R] Grouped boxplots using ggplot() from ggplot2.

On 28/07/18 17:03, Jeff Newmiller wrote: When you understand the strong dependence on how the data controls ggplot, using it gets much easier. I still have to google details sometimes though. Note that it can be very difficult to make a weird plot (e.g. multiple parallel axes) in ggplot becau

Re: [R] Grouped boxplots using ggplot() from ggplot2.

On 29/07/18 02:54, Jeff Newmiller wrote: 1) I don't know... it looks to me like you did not run my code. Aaaarrrgghhh. I *thought* I had, but instead left "fill=Type" inside the aes() call and neglected to add fill=NULL outside this call. Du!!! It's tough being mentally challenged, le

Re: [R] Philip Morris International - Windows10 migration assessment

. Actually, if I as a private user would ask that question, I guess I would get an answer here. So please - not alienate the business users. Cheers and good luck with Windows. On the contrary, I would say that one has a moral obligation to alienate tobacco companies. cheers, Rolf Turner

Re: [R] [FORGED] Re: Philip Morris International - Windows10 migration assessment

On 04/08/18 10:07, Dylan Distasio wrote: Nice virtue signalling, feel better now? I was going to respond to this, but then I thought that it would be better to heed the ancient wisdom: "Do not feed the trolls." cheers, Rolf Turner On Fri, Aug 3, 2018, 3:47 PM Clive Nich

[R] A slightly unorthodox matrix product.

rary(abind) xxx <- lapply(1:nrow(a),function(i,a,b){a[i,]%o%b[i,]},a=A,b=B) do.call(abind,c(xxx,list(along=3))) Is there a cleverer way? cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __

Re: [R] [FORGED] Re: A slightly unorthodox matrix product.

Thanks to Eric Berger, Jeff Newmiller and Chuck Berry for their help with this. Thanks especially to Eric for catching the bugs in my proposed solution. I can *never* keep the indexing straight in multidimensional arrays! It's tough being a -wit! I wouldn't have figured out the "ans0b"

Re: [R] A slightly unorthodox matrix product.

On 05/08/18 16:41, Thomas Jagger wrote: > Date: Sat, 4 Aug 2018 17:52:37 +1200 > > From: Rolf Turner <mailto:r.tur...@auckland.ac.nz>> > To: r-help mailto:r-help@r-project.org>> > Subject: [R] A slightly unorthodox matrix product. > Message-ID: <mail

[R] Mysterious seg fault.

and the Fortran code. The package is a bit complicated, so giving more detail would be cumbersome. Also I have no idea what aspects of detail would be relevant. If anyone would like more info, feel free to ask. I would really appreciate it if someone could give me some suggestions before

Re: [R] Mysterious seg fault.

On 12/08/18 09:52, Göran Broström wrote: Rolf, have you tried to run R with a debugger, as in $ R -d gdb the gnu debugger? No, I haven't. Did not know about the gnu debugger. I shall investigate. Thanks cheers, Rolf -- Technical Editor ANZJS Department of Statistics University of A

Re: [R] Mysterious seg fault.

On 12/08/18 17:42, Eric Berger wrote: Hi Rolf, When faced with such a situation I take the following approach which often helps. Use the same setup that caused the seg fault (you need a reproducible problem.) Start your R session using valgrind. e.g. in linux I do: $ valgrind R Assuming t

Re: [R] Mysterious seg fault.

On 13/08/18 12:03, Duncan Murdoch wrote: So I did: /usr/local/bin/valgrind R I believe on your system R is a script, so you can't run valgrind this way.  It's just debugging bash, not R.  You need to use R -d valgrind (though with your weird path problems, you might need a fully qualifi

Re: [R] Mysterious seg fault --- SOLVED

OK everybody! You can relax. :-) I managed to spot the loony. After mucking around with valgrind, and before trying gdb, I had one more look at my code and *finally* saw the stupid thing that I had been doing. In the call to .Fortran() I had a line nphi=as.integer(nphi), but "nphi"

Re: [R] Mysterious seg fault --- SOLVED

On 13/08/18 20:45, Henrik Bengtsson wrote: On Mon, Aug 13, 2018 at 3:51 AM Rolf Turner wrote: OK everybody! You can relax. :-) I managed to spot the loony. After mucking around with valgrind, and before trying gdb, I had one more look at my code and *finally* saw the stupid thing that

Re: [R] Mysterious seg fault --- SOLVED

On 13/08/18 23:39, peter dalgaard wrote: It's odd, possibly a bug, that you don't get Error: object 'nphi' not found but I can't offhand see where the evaluation of args to .C/.Fortran is supposed to take place. If it is indeed a bug then it would be nice an it were fixed. If that is possi

Re: [R] [FORGED] Spline function

for which you want the corresponding x value, where spf() is monotone on [a,b], and where there *exists* an x value between a and b such that spf(x) = y0. Good luck. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-75

Re: [R] Mysterious seg fault --- SOLVED

On 14/08/18 23:01, peter dalgaard wrote: Hmm, .Fortran(stats:::C_setsmu, as.integer(0)) [[1]] [1] 0 .Fortran(stats:::C_setsmu, as.integer(fumble)) Error: object 'fumble' not found .Fortran(stats:::C_setsmu, fumble=as.integer(fumble)) Error: object 'fumble' not found .Fortran(stats:::C_se

Re: [R] Mysterious seg fault --- SOLVED

On 15/08/18 13:00, Richard M. Heiberger wrote: There is no explanation other than gremlins and the malevolence that the computer gods hold towards me. fortune nomination. I demur. I already have a fortune with gremlins in it attributed to me (fortune(213)). cheers, Rolf -- Technical Editor

Re: [R] bar plot add space to group data

OP did not supply the data, I cannot experiment. Perhaps someone else will have some insight. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org mailing

Re: [R] [FORGED] Re: bar plot add space to group data

Jim: (a) There's no legend. (b) I am still curious as to why the OP's code didn't work, in that the "space=c(0,2)" argument seemed to have no effect. cheers, Rolf On 18/08/18 20:45, Jim Lemon wrote: Hi citc, Try this: geac<-matrix(c(9,9,8,8,8,23,23,23,23,22,27,27,27,25,24, 19,19,19,20,2

Re: [R] [FORGED] Re: bar plot add space to group data

On 19/08/18 16:58, Peter Langfelder wrote: My guess is that space has no effect because (1) the first element is zero and (2) the code in OP's message has barplot(gceac[,3], ... i.e. barplot does not see a matrix, only a vector. To the OP, try formatting the data to be plotted as a matrix, n

Re: [R] [FORGED] Re: bar plot add space to group data

On 19/08/18 16:12, Jim Lemon wrote: Hi Rolf, That's what comes of being in a hurry. legend(4.1,30,c("A+","A","B","C","D","E","U"), fill=c("white","lightblue","blue","orange","green","red","pink")) and I thank you for alerting me to the fact that the legend arguments in barp don't position the

Re: [R] differing behavior of mean(), median() and sd() with na.rm

E > [1] NaN Seems quite consistent/coherent to me. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and m

Re: [R] importing .v8x file in R

On 08/23/2018 08:35 PM, Rui Barradas wrote: Hello, Sorry but I don't believe this is a question for r-help. r-help is meant for questions about R code, you should find out what type of file do you have. Maybe open it and see its contents. There is really nothing we can do. Indeed. But the

[R] R shared library (/usr/lib64/R/lib/libR.so) not found.

4/R/lib/libR.so? It bothers me that /usr/lib/R/lib/libR.so was not "refreshed" from my most recent install of R. I plead for enlightenment. cheers, Rolf Turner P.S. I'm running Ubuntu 18.04. And the previous install of R was done under Ubuntu 18.04. R. T. -- Technical Editor

Re: [R] R shared library (/usr/lib64/R/lib/libR.so) not found.

On 08/23/2018 11:09 PM, Jan T Kim via R-help wrote: Hi Rolf & All, I haven't built R in a while, but my general expectation of an autotools based build & install would be that the default prefix is /usr/local, rather than /usr. So I'd expect the shared libs in /usr/local/lib, /usr/local/lib64

Re: [R] R shared library (/usr/lib64/R/lib/libR.so) not found.

See in-line below. On 08/25/2018 01:10 AM, Ista Zahn wrote: On Thu, Aug 23, 2018 at 6:57 AM Rolf Turner wrote: I *think* that this is an R question (and *not* an RStudio question!) I think this is actually and Ubuntu question, and probably belongs on R-sig-debian. Well, it's

Re: [R] R shared library (/usr/lib64/R/lib/libR.so) not found.

On 08/26/2018 02:59 AM, Berwin A Turlach wrote: G'day Rolf, On Sat, 25 Aug 2018 11:20:09 +1200 Rolf Turner wrote: I was pretty sure that the foregoing was a complete red herring. And I was right. I am not sure whether I agree. :) Huh? Black is white and up is down??? I d

Re: [R] [FORGED] installing R in ubuntu

key? Any general help on installing R in ubuntu AWS instance? very many thanks for your time and effort... You may find the following link helpful. I did. https://www.digitalocean.com/community/tutorials/how-to-install-r-on-ubuntu-18-04-quickstart cheers, Rolf Turner -- Technical Editor

Re: [R] [FORGED] Main label on Cullen and Frey

from source. This way you will get your modified version of descdist() every time you load the package. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org

Re: [R] [FORGED] Main label on Cullen and Frey

Fortune nomination. cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listin

Re: [R] sink() output to another directory

n sink("gorp/mung.txt") will put the sink() output into the file "mung.txt" in the directory "gorp". cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 ___

Re: [R] [FORGED] Question on Binom.Confint

"c") but it does not have an argument "method". Please do not expect those whom you are addressing to be telepathic. Point 3: Having "method"=lrt in the call is decidedly weird. Perhaps you meant method="lrt"; this is entirely different. cheers,

Re: [R] sink() output to another directory

ary(estacada_se_wx)) sink() system("cat stat-summaries/estacada-se-precip.txt") I get: Min. 1st Qu. MedianMean 3rd Qu.Max. -1.7292 -0.7621 0.5428 0.2808 1.2603 1.7702 OMMM!!! cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of

Re: [R] [FORGED] modify the supposed return value of a function during evaluation?

uot; syntax that Micro$oft has inflicted upon the world. But never mind.) I am however very sure that your proposed approach is Not A Good Idea. If you want an interactive structure for your function, make it interactive. In a well thought out and well organised manner. cheers, Rolf Turner --

Re: [R] [FORGED] Unable to update R 3.4 to R 3.5 in Ubuntu 18.04 LTS

how-to-install-r-on-ubuntu-18-04-quickstart cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see ht

Re: [R] [FORGED] Re: tibble question with a mean

Please see inline below. n 09/21/2018 01:07 PM, Peter Langfelder wrote: I don't know tibble, so I'll do the same with a plain data frame: a = data.frame(x=LETTERS[1:4],y=1:4,z=rnorm(4),a=c("dog","cat","tree","ferret")) a x y z a 1 A 1 -0.08264865dog 2 B 2 0.32344426

Re: [R] [FORGED] Re: Wavelet Package for Time Series

On 10/14/2018 01:26 PM, Jeff Newmiller wrote: Any question about "best" implies an optimization function that will almost certainly vary depending on what you are interested in ... please avoid asking such questions. Fortune nomination. cheers, Rolf -- Technical Editor ANZJS Department of

Re: [R] [FORGED] Re: Using apply

On 10/31/18 3:47 PM, jim holtman wrote: s2 <- apply(x*x, 2, sum) s2 [1] 55 330 It seems to me to be more "natural" (and perhaps more amenable to generalisation) to do: s2 <- apply(x,2,function(v){sum(v^2)}) But it's probably just a matter of taste. cheers, Rolf -- Technical Ed

[R] Problem with the matrix() function.

of ncol is not an integer scalar. I have also discerned that if ncol is not an integer, it is replaced by its floor value. Is this a Good Thing? What do others think? cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-3

[R] The Suggests field in a DESCRIPTION file.

checking R code for possible problems ... NOTE foo: no visible binding for global variable ‘fortunes’ Undefined global functions or variables: fortunes What am I doing wrong? Thanks for any insight. cheers, Rolf Turner P. S. I tried putting an "imports(fortunes)" in the NAMESPACE f

Re: [R] The Suggests field in a DESCRIPTION file --- never mind!!!

e. Trap for young players. Psigh. cheers, Rolf Turner On 11/18/18 11:16 AM, Rolf Turner wrote: I am building a package which contains a function from which I wish to call the fortune() function from the fortunes package --- if that package is available. I have place the line    Sugg

Re: [R] The Suggests field in a DESCRIPTION file.

On 11/18/18 11:22 AM, Fox, John wrote: Dear Rolf, "fortunes" needs to be quoted in requireNamespace("fortunes", quietly=TRUE). I hope this helps, John Thanks. I actually figured this out myself, just before getting your message! cheers, Rolf -- Technical Editor ANZJS Department of Sta

Re: [R] [FORGED] Re: R plot split screen in uneven panels

sible* to fight your way through the help. Experimenting is useful! And when you succeed in fighting your way through, it works *well*. (Not "good"!!!) cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-

Re: [R] [FORGED] SE for all fixed factor effect in GLMM

On 12/30/18 5:31 PM, Marc Girondot via R-help wrote: Dear members, Let do a example of simple GLMM with x and G as fixed factors and R as random factor: (note that question is the same with GLM or even LM): x <- rnorm(100) y <- rnorm(100) G <- as.factor(sample(c("A", "B", "C", "D"), 100,

Re: [R] SE for all fixed factor effect in GLMM

.factor(sample(c("A", "B", "C", "D"), 100, replace = TRUE)); G <- relevel(G, "A") m <- glm(y ~ x + G) summary(m)$coefficients No SE for A level in G category is calculated. * Here is a synthesis of the answers: 1/ The first solution

Re: [R] SE for all fixed factor effect in GLMM

Please keep communications on-list. On 1/2/19 10:57 AM, Marc Girondot wrote: Le 01/01/2019 à 22:35, Rolf Turner a écrit : On 1/2/19 9:35 AM, Marc Girondot wrote: Hello members of the list, I asked 3 days ago a question about "how to get the SE of all effects after a glm or glmm&quo

Re: [R] [FORGED] Tukey Test

On 1/25/19 4:51 PM, reichm...@sbcglobal.net wrote: R-Help There is an R library that will perform a Tukey test ... Surely you mean *package*. cheers, Rolf Turner P.S. You are probably thinking of the agricolae *package*. There is also the TukeyC package, which might be relevant

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

here are a number of good online tutorials. cheers, Rolf Turner -- Honorary Research Fellow Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat

Re: [R] [FORGED] (no subject)

(a) Read the posting guide. (b) This looks like homework (perhaps a homework question that you have completely misunderstood) and this list has a no-homework policy. (c) As Berend Hasselman has pointed out, your question makes no real sense anyway. cheers, Rolf Turner -- Honorary Resea

[R] data.frame() versus as.data.frame() applied to a matrix.

h more natural for converting a matrix to a data frame, yet it doesn't get it quite right, sometimes, in respect of the names. Is there some reason that as.data.frame() does not apply make.names()? Or was this just an oversight? cheers, Rolf Turner -- Honorary Research Fellow Department o

Re: [R] data.frame() versus as.data.frame() applied to a matrix.

pending on whether one uses data.frame(X) or as.data.frame(X). I'll spare you the details. :-) cheers, Rolf On February 5, 2019 2:22:24 PM PST, Rolf Turner wrote: Consider the following: set.seed(42) X <- matrix(runif(40),10,4) colnames(X) <- c("a","b"

Re: [R] [FORGED] Recreate for loop without using for loop

p function"? qprob <- function(pp) { qq <- 1 - pp -1 stotal <- 0.0 for (i in 1:length(pp)) stotal <- stotal + pp[i] * prod(qq[-i]) return(stotal) } cheers, Rolf Turner -- Honorary Research Fellow Department of Statistics University of Auck

Re: [R] [FORGED] I'm Working On A Data Security Article - Quick Question

joke! Apparently some idiotic Windoze antiviral software once identified the R executable for Windoze as being a virus. This is not to taken seriously. cheers, Rolf Turner -- Honorary Research Fellow Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 __

Re: [R] [FORGED] How to fit the following data and get the equation describing them?

y of these things? Thank you very much. You seem to be way out of your depth. You would be well advised to seek local help from a qualified statistician. cheers, Rolf Turner -- Honorary Research Fellow Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 _

Re: [R] particle count probability

On 2/21/19 12:16 AM, PIKAL Petr wrote: Dear all Sorry, this is probably the most off-topic mail I have ever sent to this help list. However maybe somebody could point me to right direction or give some advice. In microscopy particle counting you have finite viewing field and some particles coul

Re: [R] [FORGED] Q re: logical indexing with is.na

So (1:4)[1:3] just (sensibly) gives [1] 1 2 3 and *not* [1] 1 2 3 1 Perhaps a bit subtle, but it gives what you'd actually *want* rather than being pedantic about rules with a result that you wouldn't want. cheers, Rolf Turner P.S. If you do y[1:3][!is.na(y[1:3])] i.e. if y

Re: [R] [FORGED] Q re: logical indexing with is.na

On 3/10/19 6:07 PM, Jeff Newmiller wrote: Regarding the mention of logical indexing, under ?Extract I see: For [-indexing only: i, j, ... can be logical vectors, indicating elements/slices to select. Such vectors are recycled if necessary to match the corresponding extent. i, j, ... can also

Re: [R] [FORGED] Re: Issue with t.test

On 13/03/19 9:06 AM, Greg Snow wrote: The only time I have seen t.test give a p-value of 1 is when the data mean exactly equals the null hypothesis mean and the alternative is the default of two.sided. Doesn't have to be *exact* equality. Just close! E.g.: set.seed(42) x <- runif(10) me

Re: [R] [FORGED] Re: Potential Issue with lm.influence

On 4/04/19 5:34 AM, Richard M. Heiberger wrote: fortune nomination. The lesson to me here is that if you fit a sufficiently unreasonable model to data, the computations may break down. I second the nomination! cheers, Rolf -- Honorary Research Fellow Department of Statistics University o

Re: [R] analyzing qualitative data sets

n. You are encouraged to think more about the end result before you seek magic to rescue inappropriate data. -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

Re: [R] a knitr question

7;s post was the way that Duncan interpreted it. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-

Re: [R] how to process multiple data files using R loop

is to use incorrect terminology --- and the former has the advantage of not misleading all concerned. Dave Winsemius has already told you how to solve your immediate problem, using get(). cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS ___

Re: [R] C.D.F

? The list has a no homework policy. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] A basic statistics question

]-1) then you would get the same result as you get from cor(Data) (to within about 1e-15). cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] A basic statistics question

On 16/08/14 01:29, Joshua Wiley wrote: On Wed, Aug 13, 2014 at 7:41 AM, Rolf Turner mailto:r.tur...@auckland.ac.nz>> wrote: On 13/08/14 07:57, Ron Michael wrote: Hi, I would need to get a clarification on a quite fundamental statistics property, hope expeRt

Re: [R] Generating a polygon around points

out the vertices of the polygon (in the appropriate order) in some "by hand" manner. It's possible that clickpoly() might be of help to you. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list

Re: [R] Generating a polygon around points

e. I think that the Universe is picking on Dave (instead of on me, for once!). :-) cheers, Rolf -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] Generating a polygon around points

On 16/08/14 11:01, David Winsemius wrote: On Aug 15, 2014, at 3:11 PM, Rolf Turner wrote: I tried RSiteSearch("convex hull") just now, as an experiment, and it seemed to work just fine. I think that the Universe is picking on Dave (instead of on me, for once!). :-) R

Re: [R] Generating a polygon around points

On 17/08/14 01:57, Greg Snow wrote: On Fri, Aug 15, 2014 at 4:06 PM, Rolf Turner wrote: OTOH R is still lacking a mind_read() function so it probably would NOT give you *exactly* what you want. We can try the (very pre-alpha) esp package: source('g:/R/esp/esp.R') esp() [1] &

Re: [R] "no visible binding for global variable" and with() vs. within()

line: utils::globalVariables("x") cheers, Rolf -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/post

Re: [R] loading saved files with objects in same names

ou mean plot(envA$object, envB$object) ??? Or am I misunderstanding something? cheers, Rolf -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

Re: [R] spatstat rmh problem

ct "im.pop" is present in the global environment. cheers, Rolf Turner On 16/09/14 02:30, Sebastian Schutte wrote: Dear R and spatstat developers, Thanks so much for the time and effort that you invest into this awesome software. I have a problem simulating from a Point Process Model in spat

Re: [R] spatstat rmh problem

ity(rmh(mod,new.coef=c(1,200 (2) However, even when the correct call is given you still wind up with identical densities!!! Hm. I think this may be a bug; I'll will check with the other authors of spatstat and report back. cheers, Rolf Turner On 16/09/14 16:30, Sebastian Schut

Re: [R] spatstat rmh problem

There was indeed a bug in rmh() w.r.t. the "new.coeff" argument. The bug has been fixed and will not be present in the next release of spatstat. cheers, Rolf Turner On 16/09/14 16:30, Sebastian Schutte wrote: Thanks so much for your comments. Sorry for not having sent a runni

Re: [R] Error in predict.lm() for models without intercept?

seeing to make you think that there is an error, but I am not seeing anything untoward. cheers, Rolf Turner On 17/09/14 10:51, isabe...@ghement.ca wrote: Hi everyone, It appears my R code didn't come through the first time (thanks for letting me know, Ista). Here is my message again

Re: [R] Error in predict.lm() for models without intercept?

ata=data.frame(v = seq(0,to=max(v),by=100), x = mean(x),z=mean(z)), interval="confidence") conf.band.v <- data.frame(lwr=conf.band.v[,"lwr"], fit=conf.band.v[,"fit"], upr=conf.band.v[,"upr"]) matplot(seq(from=0,to=max(v),by=

Re: [R] plot

,y2,y3,y4 against x) including title and key y1<-c(0.84,1.03,0.96) y2<-c(1.30,1.46,1.48) y3<-c(1.32,1.47,1.5) y4<-c(0.07,0.07,0.07) x<-c(500,1000,2000) -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https

Re: [R] Fwd: Dynamic regex/sub changes to function

On 23/09/14 12:18, Bert Gunter wrote: Daniel: Do you understand what a parse tree is? Paraphrasing a quote attributed to Feynman: "If you think you understand what a parse tree is, then you don't understand what a parse tree is." :-) cheers, Rolf -- Rolf Turner T

Re: [R] division of col by the sum of the col

On 24/09/14 03:10, carol white wrote: Hi, If I want to divide the column of a matrix by the sum of the column, should I loop over the columns or can I use apply family? m1 <- apply(m,2,function(x){x/sum(x)}) should do what you want IIUYC. cheers, Rolf Turner -- Rolf Turner Technical Edi

Re: [R] Median of streaming data

ptotic properties as the median. See: * U. Holst, Recursive estimators of location. Commun. Statist. Theory Meth., vol. 16, 1987, pp. 2201--2226. and * Murray A. Cameron and T. Rolf Turner, Recursive location and scale estimators, Commun. Statist. Theory Meth., vol. 22, 1993, pp. 2503--

Re: [R] Median of streaming data

On 24/09/14 20:16, Martin Maechler wrote: Rolf Turner on Wed, 24 Sep 2014 18:43:34 +1200 writes: > On 24/09/14 17:31, Mohan Radhakrishnan wrote: >> Hi, >> >> I have streaming data(1 TB) that can't fit in memory. Is >> there a

Re: [R] Median of streaming data

f the old (Splus?) code in my chaotic archives. Should be able to get it done "real soon now". cheers, Rolf -- Rolf Turner Technical Editor ANZJS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

Re: [R] Median of streaming data

xt" so that they will get through the mailing list processor without being stripped. With a bit of luck. If they *don't* get through, anyone who is interested should contact me and I will send them to you "privately". cheers, Rolf -- Rolf Turner Technical Editor ANZJS r

Re: [R] Histogram from a single column of a data frame

n ***really***; it *IS* obvious) is hist(oded$Breed) --- which tells hist() where to find "Breed". (4) You could also use with(oded,hist(Breed)) which is a useful syntax in some circumstances. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS

Re: [R] Median of streaming data

On 26/09/14 21:48, Martin Maechler wrote: Rolf Turner > I have coded up the algorithm from the Cameron and Turner > paper. Dunno if it gives exactly the same results as my > (Splus?) code from lo these many years ago (the code that > is lost in the mists o

Re: [R] Regarding R package spatstat

tried a test just now and managed to do an install of version 1.26-0 with no problems and no hangup. cheers, Rolf Turner P. S. I am puzzled as to how you are getting 'Target "all" is up to date' given that you are installing from the tarball. This should only happen w

Re: [R] Best Distribution

attached script should do what you want. cheers, Rolf Turner -- Rolf Turner Technical Editor ANZJS # # Script botto.txt # require(MASS) require(FAdist) x <- c(42.2, 45.2, 46, 48, 54, 54.1, 59.4, 61, 62.2, 63.5, 65.024, 71.9, 73.4, 76.6, 76.708, 77.5, 77.724, 78, 81.3, 84.7, 84.

Re: [R] Anderson-Darling gof test (was "no subject")

not appropriate for this list. Moreover it would appear that you are well out of your statistical depth. You should seek local statistical consultation rather than trying to scramble your way through topics of which you have no proper understanding by peppering the r-help list with ill-posed

Re: [R] Bug in myintegrate in Elliptic package

t;x" as the name of one of the "..." arguments you wind up getting a call of the form "f(x,x=0.1,psi_measure=psi)" --- no wonder the poor damned thing gets confused and gets it wrong. Try: melvin <- function(x,...){x^2} melvin(2) melvin(2,x=5) and you'

Re: [R] grep won't work finding one column

vector. If you want a (one-column) data frame you need to use "drop=FALSE" in your subscripting call. You need to study up on R and learn how it works (read the Introduction to R) and stop going off half-cocked. cheers, Rolf Turner P.S. It is a ***bad*** idea to use "df

Re: [R] Ternary Plots Do Not Display Ellipses in PDF

In all probability this is FAQ 7.22. I.e. use: print(ellipses(mean=mn, var=vr, r=r, steps=72, thinRatio=NULL, aspanel=FALSE,col='red', lwd=2)) cheers, Rolf Turner P. S. See also fortune(123). R. T. On 15/10/14 11:00, Rich Shepard wrote: A rather strange situ

<    1   2   3   4   5   6   7   8   9   10   >