On 30/04/19 12:01 PM, Marc Schwartz wrote:
Hi Guys,
I suspect that this entry from news() for 3.6.0 is relevant:
"When loading namespaces, S3 method registrations which overwrite
previous registrations are now noted by default (using
packageStartupMessage())."
Always a good idea to read th
On 30/04/19 1:31 PM, Duncan Murdoch wrote:
On 29/04/2019 9:25 p.m., Duncan Murdoch wrote:
On 29/04/2019 8:43 p.m., Rolf Turner wrote:
On 30/04/19 12:19 PM, Jeff Newmiller wrote:
Since you now have this indirect dependency, you should make sure you
have updated this gaggle of packages. I
On 14/06/19 5:30 AM, Tina Chatterjee wrote:
Hello everyone!
I have the following dataframe(df).
a<-c("a1","a2","a2","a1","a1","a1")
b<-c("b1","b1","b1","b1","b1","b2")
c<-c("c1","c1","c1","c1","c1","c2")
time <- c(runif(6,0,1))
I know that it's basically harmless, but it drives me nuts
On 22/06/19 5:21 AM, Steven Yen wrote:
Unhappy but thanks.
Steven
Well, if you're *really* unhappy (???) you could always do
suppressWarnings(RNGversion("3.5.3"))
Does that make you happier?
cheers,
Rolf
On 6/22/2019 1:13 AM, Uwe Ligges wrote:
On 21.06.2019 19:09, Steven Yen wr
On 30/06/19 8:08 PM, Eric Berger wrote:
Nothing came through. Note that this is a plain text list and content that
is not plain text may be stripped off (or mangled).
On Sat, Jun 29, 2019 at 10:14 PM Don or Charlotte Smith
wrote:
Sent from my iPhone
Eric: This is most likely some sort
. (But it sounds to me more like a bug than a feature.)
I would just like to point out that when I do
help.search()
(with no argument) I get:
Error in help.search() : do not know what to search
This seems reasonable, apart from being ungrammatical.
If I do, e.g.,
help.search("elli
On 6/07/19 11:15 PM, Jeff Newmiller wrote:
I have tried both Mathematica and Wolfram Alpha at times and found neither
of them so compelling that I felt the slightest bit jealous of people who
have licensed access to Wolfram's tools and services.
You mean "... *envious* of people ... "!!!
On 7/07/19 11:59 AM, Jeff Newmiller wrote:
Who let the editor in here? :-)
Sorry! Can't help it! :-)
Anyway, I stand corrected.
Sad to say, you are far from being alone in your breach of this
particular rule of usage.
cheers,
Rolf
--
Honorary Research Fellow
Department of Statistics
On 13/07/19 10:54 AM, Spencer Graves wrote:
Hello:
What do you suggest I do about modeling random truncation?
Good question! Probably the best answer is "Give up and go to the pub!" :-)
But seriously, there is a package DTDA on CRAN which purports to analyse
randomly truncated da
On 13/07/19 3:31 PM, Abby Spurdle wrote:
> The distribution of the randomly truncated variable has thus four
> parameters: a, b, mu and sigma. I was able to write down the likelihood
> and attempted to maximise it
I read the Wikipedia article more carefully.
The formula is relatively simp
Jeff: Your comments are often (almost always?) a bit rough about the
edges, and on the recipient, but are always cogent. Although John M.
has a valid point, I tend to agree with you. I would say that if you
want to make money, trying to sell your own software is a bad way to go
about it,
ipley into people! If the recipient can't take the heat, he or
she should get out of the kitchen!
See also fortunes::fortune(87).
cheers,
Rolf Turner
P.S. Jeff makes a huge and extremely useful contribution to R-help. He
gives generously of time and effort to solve beginners'
waste time discussing such matters on the list. Ad hominem
comments are absolutely out of place.
-Roy
On Jul 24, 2019, at 2:49 PM, Rolf Turner wrote:
On 25/07/19 4:36 AM, Weiwen Ng, MPH wrote:
Here's one way to phrase your reply:
"I'd recommend you search Google. For exam
questions he asks are generally ill-conceived and ill-posed.
He appears to be attempting advanced statistical analysis without having
mastered the elementary basics.
Spencer should focus on remedying his misconceptions and filling in the
lacunae in his knowledge.
cheers,
Rolf Turner
--
Hon
I'm obviously not understanding something here, but it seems to me that
the conjecture
It appears to me that the cause is buried in effects:::Analyze.model
in or close to the the lines
if (is.numeric(xlevels) & length(xlevels) == 1L) {
levs <- xlevels
for (name in focal.predictors)
On 6/09/19 5:30 PM, pusuluri madhu wrote:
Please unsubscribe me
Go unsubscribe yourself! :-)
See the footer at the bottom of every r-help posting:
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
HTH
cheers,
Rolf Turner
t gets.
IMHO this is a precious PC quibble, taking offence where no offence is
intended.
If you are really concerned about literal slavery --- as everyone should
be --- then join/contribute to an appropriate activist organisation
(e.g. Amnesty International, which is my personal choice
On 22/09/19 11:19 PM, Richard O'Keefe wrote:
Whenever you want a vector that counts something,
cumsum of a logical vector is a good thing to try.
Fortune nomination.
cheers,
Rolf
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 882
maxtemp ~ sampdate, data=watertemp, type="h", col="red",
)
points(maxtemp ~ sampdate, data=watertemp, type="h",col="blue")
(Not tested since your example is not reproducible.)
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of
On 27/09/19 11:08 AM, David Winsemius wrote:
Instead of trying to mix lattice and base functions, you might try using the
formula:
maxtemp+mintemp ~ sampdate
And then: col= c(“red”, “blue”)
Sent from my iPhone, so make sure those quotes are ordinary double quotes.
Ah-ha! I've learned so
On 28/09/19 9:16 AM, Jim Lemon wrote:
On Sat, Sep 28, 2019 at 2:04 AM Martin Maechler
wrote:
For back compatibility reasons, the old command line option will
continue to work so the many shell and other scripts that use
it, will not stop working.
That's a relief. I was getting worried tha
? :-)
This appears to be an Rstudio question and not an R question. It should
be directed at the proper resource.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
__
R
60
.
.
.
Should I extract the data frame from the list one by one?
Try (something like):
revisedTest <- lapply(test,function(x){x$RR[x$RR==] <- NA; x})
Or, using Bert Gunter's slick is.na() trick:
revisedTest <- lapply(test,function(x){is.na(x$RR) <- x$RR==; x})
On 18/10/19 2:58 PM, Thevaraja, Mayooran wrote:
Hello
You can use the following function,
##
replace_missings <- function(x, replacement) {
is_missings <- is.na(x)
x[is_missings] <- replacement
message(sum(is_missings), " miss
I have the vague impression that "chi-squared" is more common in British
usage and "chi-square" more common in American usage. I'm pretty sure
that either is acceptable, although "chi-squared" sounds much better to
my ear.
Of course within a given document (or collection of related documen
- as.data.frame(lX)
This seems to work, but is rather kludgy. Is there a better way?
Thanks for any pointers.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
__
R-help@r-p
On 20/10/19 3:00 PM, Richard M. Heiberger wrote:
Look at
methods(as.data.frame)
Define your specialized columns to have a newly defined class, say "myclass".
Then write as.data.frame.myclass
It will be similar to the function you already have in the lapply statement.
Now your statement
X <- X[ok,
On 21/10/19 1:15 AM, Rui Barradas wrote:
Hello,
Richard's idea is good but shouldn't it be `[.myclass` instead?
Yes, I kind of thought that, and cobbled together something on that
basis that seemed to work. However my code was rather a hodge-podge. I
kept having to work around errors th
On 21/10/19 11:07 AM, Rui Barradas wrote:
Hello,
Sorry, you're right, in the method it's x, X is the test dataframe.
Repost:
`[.myclass` <- function(x, i, j, drop = if (missing(i)) TRUE else
length(cols) == 1){
SaveAt <- lapply(x, attributes)
x <- NextMethod()
lX <- lapply(names(x
Sigh. See FAQ 7.31. (As someone else has remarked on this list, 7.31
is by far the most frequently asked of all frequently asked questions.)
An aside: you should seriously consider upgrading your R installation;
the current version is 3.6.1.
cheers,
Rolf Turner
--
Honorary Research
Richard: I know that you mean well, but *please* don't do people's
homework for them!!! (They are *cheating* by asking R-help to do their
homework.)
cheers,
Rolf Turner
On 6/11/19 4:27 AM, Richard O'Keefe wrote:
This looks vaguely like something from exercism.
Le
s something like the plot in the attached pdf file
test.pdf (which I managed to produce using lattice).
Can anyone point me in the right direction?
Thanks.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 8
like to get the same thing as if I said vars(clyde) which gives
>
[[1]]
expr: ^clyde
env: global
There *must* be some magic arcane incantation that I can apply to
rowName (and colName) to get what I want. Mustn't there?
I tried things like vars(as.name(rowName)) --- nope, no h
ncol = ncols)+
ylab("Success fraction") +
theme_bw()
print(ggplot(X) +
geom_point(aes(y = Sfrac , x = x)) +
facet_grid_paginate(facets = Species ~ LifeStage:degC,
page = page,
nrow = nrows, ncol = ncols)+
ylab("Success fract
]]))
Thanks again.
cheers,
Rolf
On Saturday, November 16, 2019, Rolf Turner <mailto:r.tur...@auckland.ac.nz>> wrote:
I need to call ggplot() from another function with the names of the
faceting variables supplied as arguments to the calling function.
These names (which are
base graphics to
illustrate roughly what I am after.
I would be grateful if someone could point me in the right direction.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
#
# Reproducible example.
#
library
On 2/12/19 3:03 am, Rui Barradas wrote:
Hello,
See if this is it. The standard trick is to reshape the data from wide
to long, see the SO post [1]. Then add a scale_shape_* layer to the plot.
yyy <- cbind(xxx, y3 = y3)
long <- reshape2::melt(yyy, id.vars = c("x", "y1", "grp"))
ggplot(lon
On 2/12/19 10:45 am, Rui Barradas wrote:
Hello,
Here are two ways.
The first is an adaptation from your code. It uses facet_wrap_paginate,
not *_grid_*.
plotObj2 <- vector("list",2)
for(pg in 1:2) {
plotObj2[[pg]] <- ggplot(egDat) +
geom_point(aes(y = obsd, x = x),
On 2/12/19 5:08 pm, Rui Barradas wrote:
Hello,
Here are two ways of drawing the lines black and at the same time
removing the lines in the legend. The second way is more idiomatic.
1. Override the colour setting in the ggplot call when drawing the lines:
geom_line(aes(y = y1), colour = "bla
The best advice that anyone could give:
See fortunes::fortune("Friends") .
cheers,
Rolf Turner
On 6/12/19 4:39 am, Thomas Subia wrote:
Colleagues,
I'm trying to extract a cell from all Excel files in a directory.
library(readxl)
files <- list.files(pattern="*.x
ontains" (or may contain) one or
may R objects which similarly are "like" individual files. However
there is a distinction and it can be important to keep this distinction
in mind in order to understand subtle phenomena that may appear.
cheers,
Rolf Turner
--
Honorary
ore the = sign, I seem to be successful only at grabbing
> the stuff after the equal sign.
>
> Pointers to the obvious fix? Thanks...
Perhaps names(test) ???
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Sta
I have just changed my email address with the R-help mailing lists, and
I want to check that things are working. Sorry for the noise.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373
Eglhmm"). The source package and a Windoze binary
are available. I have not yet submitted my package to CRAN, although I
hope to do so in the fairly (???) near future.
My package may well have a fair few bugs lurking in its innards. (But
then, don't they all?) If you do experiment with
On Mon, 3 Jul 2023 13:40:41 -0700
Bert Gunter wrote:
> I am not going to try to sort out your confusion, as others have
> already tried and failed.
Fortune nomination!!!
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats.
On Fri, 18 Aug 2023 12:17:51 +0200
Martin Maechler wrote:
> I think it would be nice to provide the average R user with a
> (possibly super small) R package that allows to turn on (and off)
> such CNR reproducibility.
Would it be possible to effect this on/off via options()?
cheers,
Rolf
advice: graphics can be very revealing, and are easy to
effect in R. Relevant method: plot.function(); relevant utility:
abline(). Look at the help for these.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
Df Sum Sq Mean Sq F value Pr(>F)
> Treatment 3 200046668 9.553 0.000246 ***
> Expression1 24902490 3.568 0.071050 .
> Treatment:Expression 3 2386 795 1.140 0.353142
> Residuals24 16751
to your enquiry so
far may be relied upon to give sound advice.)
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373-7599 ext. 89622
Home phone: +64-9-480-4619
___
e maintainer of UMR has been notified, and a revision of
UMR should become available sometime in the near future.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of AucklandZhiwei Zhang
Stats. Dep't. (secretaries) pho
ch leaves me about 150 points short!) to be
able to cope with its intricate and unintuitive syntax.
I believe that what you want to do is triv in base graphics, but I must
confess that I have not gone through the details.
cheers,
Rolf Turner
--
Honorary Research Fellow
Departmen
On Sat, 13 Jan 2024 17:59:16 -0500
Duncan Murdoch wrote:
> My guess is that one of the bootstrap samples had a different
> selection of countries, so factor(Country) had different levels, and
> that would really mess things up.
>
> You'll need to decide how to handle that: If you are trying t
On Sat, 13 Jan 2024 16:54:01 -0800
Bert Gunter wrote:
> Well, this would seem to work:
>
> e <- data.frame(Score = Score
> , Country = factor(Country)
> , Time = Time)
>
> ncountry <- nlevels(e$Country)
> func= function(dat,idx) {
>if(length(unique(dat[idx,'Count
Please see fortunes::fortune(285).
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373-7599 ext. 89622
Home phone: +64-9-480-4619
__
R-h
d love to share my script so you guide me where I did wrong .
(1) This post is far too vague to be appropriate for this list.
(2) You should learn some statistics; probably linear modelling.
(3) You should talk to your thesis advisor.
(4) Please see fortunes::fortune(285).
cheers,
Rolf Tur
<- lm(response~time*treatment,data=xmpldata)
yhat <- fitted(mod)
xtb <- with(xmpldata,xtabs(yhat ~ time+treatment))
print(xtb)
> treatment
> time Control Treatment
> Post 94.10501 201.99112
> Pre 101.63792 210.04248
Is that (something li
See fortunes::fortune(36).
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373-7599 ext. 89622
Home phone: +64-9-480-4619
__
R-help@r-projec
On Tue, 7 May 2024 06:34:50 -0400
Duncan Murdoch wrote:
> On 07/05/2024 6:31 a.m., Iago Giné Vázquez wrote:
> > Thanks Duncan.
> >
> > I am currently on Windows. Is there any solution for it?
>
> Switch to Linux or MacOS?
Fortune nomination!
cheers,
Rolf Tu
ing wrong with my (tested) solution.
>
>
> Hope this helps,
>
> Rui Barradas
I never cease to be impressed and indeed amazed by Rui's apparently
inexhaustible patience.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats.
Apologies for the noise.
R.T.
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Stats. Dep't. (secretaries) phone:
+64-9-373-7599 ext. 89622
Home phone: +64-9-480-4619
__
R-help@r-project.org mailing list -- To UNS
that are tested.
The last point is focussed on in:
https://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf
I don't agree with everything that Prof. Venables says, but "Exegeses"
is *well* worth a read.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
Univer
On Mon, 16 Sep 2024 09:28:14 +0200
Francesca wrote:
> Dear Contributors,
> I hope someone has found a similar issue.
I hope *not*! 😊️
> I have this data set,
You may have, but we haven't. The data you provided have an
incomprehensible (to me at least) structure. Please use dput()
to include
On 20/12/19 1:30 am, Bert Gunter wrote:
"But the important point is:
If you know the structure of the data you want to
parse, then it is best to tell R (or any other language)
this structure explicitly. "
Fortune nomination!
Second the nomination!
cheers,
Rolf Turner
--
Honorar
the "P" in "Point type" to be
directly above the "o" in "obsd" and "f" in "fitted".
Is there any way to effect this? Thanks.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone:
On 24/12/19 2:29 pm, Jim Lemon wrote:
Hi Rolf,
Following the docs back to draw.key, It looks like the ellipsis
argument is ignored. I was hoping for a brilliant solution along the
lines of:
adj=0
that could be passed down the functions like a hot potato, but was disappointed.
Thanks for giv
ource "rmgarch_1.3-7.tar.gz" must be place in the same
directory as that from which you opened the terminal window.
Something like this should also work under Windoze, but there are
probably extra "gotchas" under that system, and I can give you no
guidance there.
cheers,
211,1))
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.ethz.ch/mailman/listinfo/r-help
PLEA
On 2/01/20 9:51 pm, Pietro Coretto wrote:
No problem Rolf. Thanks for you interest. But the problem is still
unsolved!
I experimented and found that I too could not install rmgarch. However
the string of error messages that I got was quite different from yours.
I did some scrounging ar
mportant to get your
terminology right and to understand the concepts that you are dealing
with. Slap-dash hammer and hope is a recipe for disaster, especially
in R.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9
accomplish simple actions and
impossible to accomplish complex actions.
--- Doug Gwyn
(22/Jun/91 in `comp.unix.wizards')
But as I say it's a matter of taste and inclination. It is likely that
Felix is among the overwhelming majority of users who prefer GUIs.
On 27/01/20 11:06 am, Jim Lemon wrote:
Hi Puja,
Three things:
3) You should probably change the subject line of your message to
"Would anyone care to do my work for me?"
Fortune nomination!!! :-)
cheers,
Rolf
--
Honorary Research Fellow
Department of Statistics
University of Auckland
f the
problematic *.xls files (note that you *cannot* attach such a file to
your email; it will get stripped; you'll need to provide a URL from
which the file can be obtained). Given such a file someone on the list
may be able to help you.
cheers,
Rolf Turner
--
Honorary Research Fellow
ould start by doing something like
CatsDogs$Animal <- factor(CatsDogs$Animal,levels=c(0,1),
labels=c("Cat","Dog"))
and similarly for the other columns. When you have learnt a bit about
R, doing your frequency tabulations and barplots will then be
can't find anything
helpful. I've fiddled about with trellis.par.set() and cannot seem to
get any effect.
Could someone please give my poor feeble brain some guidance? Ta.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ex
On 13/02/20 6:16 pm, Deepayan Sarkar wrote:
On Thu, Feb 13, 2020 at 10:39 AM Richard M. Heiberger wrote:
It works as anticipated for me
xyplot(1 ~ 1,
+ main="The quick brown fox jumped\n over the lazy dog.")
xyplot(1 ~ 1,
+ main="The quick brown fox jumped over
e to say more.
Members of this list do not, as far as I know, have psychic powers.
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 --
numerals come out as being less than 0.
Can anyone explain the rationale to me? Not that it matters a damn.
Just idle curiosity.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
___
Thanks to Marc and Duncan for setting me straight. I guess the piece of
the puzzle that I was overlooking is the fact that lexicographic
ordering for string comparison depends on locale.
It would also have helped me a bit if I'd done the RTFM thing and looked
at ?"<" !!!
Thanks again.
c
se a point that a
lot of people don't seem to get. A linear model is
*linear in the parameters*
not in the predictors. For example
y = beta_0 + beta_1 * x + beta_2 * x^2 + E
is a linear model. It is linear in (beta_0,beta_1,beta_2). It is not
of course "linear in x
pass
5 passpass
6 failfail
7 passfail
If you want input and output combined, as in the way that you presented
your data use cbind(X,foo(X)).
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64
On 22/03/20 8:44 pm, Deepayan Sarkar wrote:
Another possible approach is to use split -> lapply -> rbind, which I
often find to be conceptually simpler:
d <- data.frame(Serial = c(1, 1, 2, 2, 2, 3, 3),
Measurement = c(17, 16, 12, 8, 10, 19, 13))
dlist <- split(d, d$Serial)
d
I have written a function that calls xyplot() from the lattice package.
In this function (which is actually a method for plot(), I assign a value
plotObj <- xyplot(fmla,groups=...)
and then either either do print(plotObj) --- to display a plot as well
as returning plotObj, to be kept f
triangle).
The data from data2.txt have an extremes column with all three levels
and the corresponding plot is as desired.
I'm pretty sure that the solution to my problem is quite simple, but I
am too stupid to figure it out. Can anyone guide me in the right direction?
Thanks.
ompelled to examine something that I had not previously
considered --- and the penny dropped!!! I saw what I was doing wrong,
and it was easy to fix.
Thanks everybody!!! :-)
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88
,
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.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
n
Introduction to R", readily available from the R web page (under "Manuals").
HTH
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
__
R-help@r-project
an error in the documentation here. Clearly USE.NAMES has
a different impact on vapply() than it has on sapply() and the
documentation does not indicate this, in fact quite the opposite.
It might be appropriate to submit a bug report (see
https://www.r-project.org/bugs.html) but it's
correlated normal random
variables?
?MASS::mvrnorm
?mvtnorm::rmvnorm
There may be others!
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
On 12/04/20 6:09 am, Martin Maechler wrote:
Rolf Turner
on Fri, 10 Apr 2020 12:23:49 +1200 writes:
> IMHO there is an error in the documentation here. Clearly
> USE.NAMES has a different impact on vapply() than it has
> on sapply() and the documentation
The answer is very simple: parentheses. (Also think about "operator
precedence".) If you assign rn <- 3, then 1:rn-1 is:
[1] 0 1 2
The "-" operator is applied *after* the ":" operator.
You want 1:(rn-1) which gives
[1] 1 2
and the desired result.
ch
)}
fit <- optim(par0,fn=foo,method="L-BFGS-B",lower=lwr,topn=5,x=x)
then optim() returns a result without complaint.
Am I somehow messing up the syntax for fitdistr()?
cheers,
Rolf Turner
P. S. I've tried supplying the "method" argument, method="L-BFGS-B&quo
...)
I can't see any reason why it should call dens with parm[1] < lower[1].
On Sun, Apr 26, 2020 at 5:50 PM Abby Spurdle wrote:
I haven't run your example.
I may try tomorrow-ish if no one else answers.
But one question: Are you sure the "x" and "i" are corr
On 27/04/20 3:01 am, J C Nash wrote:
:
Peter is correct. I was about to reply when I saw his post.
It should be possible to suppress the Hessian call. I try to do this
generally in my optimx package as computing the Hessian by finite differences
uses a lot more compute-time than solving the op
On 29/04/20 10:07 am, Abby Spurdle wrote:
I haven't attempted this.
(Mainly because I'm not familiar with the theory surrounding it).
However, I looked at the documentation for the spatstat package.
There are are several functions prefixed with pcf, including one named pcf3est.
According to i
The (excellent) MWE that you provided runs just fine for me (in R 4.0.0
under Ubuntu 18.04) either by sourcing "aaa.R" or by using copy-and-paste.
Must be Windoze thing. Switch to Linux!
cheers,
Rolf Turner
n 29/04/20 2:25 pm, Jinsong Zhao wrote:
On 2020/4/29 8:05, Jinsong
On 30/04/20 12:28 am, Eric Leroy wrote:
Dear all, I am sorry to see all the reactions I provoked from a
newbie user. Anyway, thank you for the answer I think that the
pcf3est function responds to my question.
Indeed the spatstat is a very impressive library and I am very grateful to the
the d
On 2/05/20 4:46 am, Ogbos Okike wrote:
Dear Contributors,
I am trying to do a plot with multiple y-axis on a common x-axis.
I would strongly advise you *not* to. This, although often done, is a
bad practice, and can sometimes (often) give misleading impressions. See
https://blog.d
dents. (See posting guide.)
(2) Please do not post in html. (See posting guide.)
(3) I agree that this looks like a bug.
(4) Email about issues like this should be sent to the maintainer of the
package in question (see maintainer("car")) and not to r-help.
cheers,
Rolf Turner
Hi Rui. Doesn't happen to me under Ubuntu 18.04:
install.packages('cowplot',lib=.Rlib)
trying URL 'https://cloud.r-project.org/src/contrib
/cowplot_1.0.0.tar.gz'
Content type 'application/x-gzip' length 1275585 bytes (1.2 MB)
==
downloaded 1.2
There is a function excel_sheets() in the readxl package which will tell
you the names of the sheets.
Using that you should probably be able to take the appropriate evasive
action.
cheers,
Rolf Turner
On 27/05/20 2:59 pm, Ravi Jeyaraman wrote:
I’ve already tried that and doesn’t work
901 - 1000 of 2111 matches
Mail list logo
