All: Kindly take this offline please.
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 Mon, Aug 3, 2020 at 12:39 PM Abby Spurdle wrote:
> > Sorry, Abby, I do
> Sorry, Abby, I do disagree here ((strongly enough as to warrant
> this reply) :
Which part are you disagreeing with?
That unambiquous names/references should be used, or that there are
many R functions for GLMs.
The wording of your post, suggests (kind of), that there is only one R
function for
> Abby Spurdle
> on Sun, 2 Aug 2020 15:13:51 +1200 writes:
> That's a bit harsh. Isn't the best advice here, to post a
> reproducible example... Which I believe has been
> mentioned.
> Also, I'd strongly encourage people to use
> package+function name, for this
That's a bit harsh.
Isn't the best advice here, to post a reproducible example...
Which I believe has been mentioned.
Also, I'd strongly encourage people to use package+function name, for
this sort of thing.
stats::glm
As there are many R functions for GLMs...
On Sun, Aug 2, 2020 at 12:47
On 2/08/20 5:39 am, Paul Bernal wrote:
Dear friends,
Hope you are doing great. I want to fit a logistic regression in R, where
the dependent variable is the covid status (I used 1 for covid positives,
and 0 for covid negatives), but when I ran the glm, R complains that I
should make the
On 22/07/20 2:26 pm, Paul Bernal wrote:
Dear friends,
I have a sample dataset, which is basically the number of transits through
a particular waterway, and is on a daily basis.
MyDat <- dataset$DailyTransits
What I´d like to do is to test whether MyDat follows a poisson distribution
or
On 15/07/20 7:45 pm, Sam H wrote:
Hi,
I am trying to download some data using read.csv and it works perfectly in
RStudio and fails in the R console in the terminal in Ubuntu 18.04 after
upgrading from R 3.6.3 to 4.0.2. Before upgrading this worked in the R
console in the terminal also
On 2/07/20 11:24 am, Alexey Shipunov wrote:
I need to clarify. This part of example(dotchart) does not show group
labels in my case:
===
dotchart(VADeaths, main = "Death Rates in Virginia - 1940")
Ah! I see. I didn't know what I should be looking for in this
instance. Sorry for the
I need to clarify. This part of example(dotchart) does not show group
labels in my case:
===
dotchart(VADeaths, main = "Death Rates in Virginia - 1940")
===
My session is
===
> sessionInfo()
R version 4.0.2 (2020-06-22)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.6 LTS
On 2/07/20 2:58 am, Alexey Shipunov wrote:
Dear colleagues,
There is a new problem with dotchart(), and it is very simple to reproduce.
Just run example(dotchart).
On R versions < 4, group labels ("Urban Female" and so on) were
visible. Now they are not visible.
I just tried
On 3/06/20 2:21 pm, p...@philipsmith.ca wrote:
Thanks Bert. That did it.
Philip
On 2020-06-02 22:02, Bert Gunter wrote:
In a function you must explicitly print/plot the ggplot() object, I
assume. i.e. plot(ggplot(...)) etc.
I do not use ggplot, so if I'm wrong, sorry. But try it.
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
Hi
Something like this might do what you want ...
## From example(manhattan)
manhattan(gwasResults, annotatePval = 0.0001)
library(gridGraphics)
grid.echo()
## grid.ls()
## The annotations have "text" in their name
labels <- grid.grep("text", grep=TRUE)
grid.edit(labels, gp=gpar(cex=1))
On 3/05/20 11:30 am, Yousri Fanous wrote:
library (car)
aa <- data.frame(x=c(2, 5, 6, 7, 8),
+ y=c(5, 10, 9, 12, 11),
+ ch=c("N", "Q", "R", "S", "T"),
+ stringsAsFactors=FALSE)
scatterplot(aa$x,aa$y,smooth = FALSE, grid = FALSE, frame = FALSE,regLine=F)
Both x and y boxplots are
It may be of some comfort to some readers, that I received emails from
two moderators, telling me that I was
I acknowledge that my initial email was a mistake.
However, I need to note recent public comments from the author of the
spatstat package:
Your internet skills are pathetic.
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
> Rolf Turner
> on Thu, 30 Apr 2020 10:42:46 +1200 writes:
> 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
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
Thank you
Yousri
On Wed, Apr 29, 2020 at 5:12 PM Paul Murrell
wrote:
> Hi
>
> This behaviour is as expected.
>
> The layout.show() function is just there to help visualise what the
> layout will look like.
>
> So for testing purposes, you would do something like ...
>
> layout(...)
>
Hi Abby,
Once again I must say your form of communication puzzles me. Do you
believe it benefits the open source community to communicate like this?
What about contacting the package maintainer in this way:
"Hi, I was helping someone else using your package to estimate a pair
correlation
Dear all,
I see two issues here:
1. A new user has a hard time finding and using a specific function in
spatstat. As package authors we are always interested in such reports
and we then try to improve documentation, which is indeed a very
important part of any software project. The package is
Hi
This behaviour is as expected.
The layout.show() function is just there to help visualise what the
layout will look like.
So for testing purposes, you would do something like ...
layout(...)
layout.show(n)
Then to actually use the layout, you would do something like ...
layout(...)
> suggests that the package isn't fit for CRAN, which I consider a
> direct insult to Adrian and all the hard work he has done
This is my last post on this subject.
I just ran R check on the source package.
After 40 minutes, R check wasn't complete.
And I note the CRAN check results show the
NOTE CITATIONS USE BRIEF EXCERPTS
> It's a complaint that no-one responded to your query within 24 hours.
Correction, it wasn't my query.
I was replying to someone else's query.
> Finally, you cast doubt on whether the function pcf3est actually does
> calculate...
You've taken what I said out
I am the main author of spatstat, and the author of the code and documentation
for pcf3est.
This is the first time in 25 years that I can remember anyone complaining about
the documentation for the spatstat package.
The available documentation for spatstat includes:
- welcome
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 developers.
Best regards
Eric LEROY
On 2020/4/29 11:55, Richard M. Heiberger wrote:
I pasted it into Emacs ESS on Macintosh, Emacs ESS on Windows, and Rgui
on Windows.
All 4.0.0 . It runs fine.
I test the code using ESS on Windows. It also give error, but with
different position:
> gz <- within(gz,
+ {
+
I should have noted that my comments weren't directed towards the main
authors, but to all people listed in the description file, which is
many, including some R core members.
Also, overall, I'm impressed by the effort here. It's just I strongly
feel that good documentation is crucial (especially
I pasted it into Emacs ESS on Macintosh, Emacs ESS on Windows, and Rgui on
Windows.
All 4.0.0 . It runs fine.
On Tue, Apr 28, 2020 at 11:36 PM Rolf Turner
wrote:
>
> 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
Hackles down, Rolf... most documentation can benefit from the perspective of a
new user. It would be helpful to link the mention of pp3 to the the pp3
function via hyperlink to help clarify what this argument is supposed to be.
Abby, FWIW I tend to recommend reading the vignettes before trying
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 Zhao wrote:
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
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
A different solution:
grr2 <- function( rn ) {
f <- function( i ) {
if ( 0 == i %% 2 ) seq.int( i )
else seq( i, 1 )
}
L <- lapply( seq.int( rn - 1 ), f )
do.call( c, L )
}
On April 17, 2020 5:11:40 PM PDT, Jeff Newmiller
wrote:
>A useful help page:
>
>?Syntax
>
>On April 17,
A useful help page:
?Syntax
On April 17, 2020 4:26:19 PM PDT, Rolf Turner wrote:
>
>
>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
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.
cheers,
Rolf Turner
On 18/04/20
On 10/04/20 12:09 pm, Bernard Comcast wrote:
I want to create a Monte Carlo simulation with 4 input parameters
that are correlated with each other. The parameters have normal
distributions and the variance/covariance matrix is known. Are there
any R functions available to generate such
On 10/04/20 10:59 am, petr smirnov wrote:
I am having trouble parsing the documentation for sapply and vapply,
and I cannot understand if it explains the different behaviour of
USE.NAMES between the two.
I noticed the following different behaviour between the two functions:
sapply(c("1"=1,
On 30/03/20 7:55 am, David wrote:
I’m trying to write a function that will add items to two vectors, and
then to create a third vector that is of the form 1, 2, 3, …, length of
one of the newly modified vectors. My problem is that what I’ve written
doesn’t seem to return any of those
Or even split -> lapply -> unsplit, in cases where you want the results put
back in the original order.
(Doesn't matter here, but it would, had it been, say, ordered 1,2,3,1,2,2,3).
-pd
> On 22 Mar 2020, at 08:44 , Deepayan Sarkar wrote:
>
> Another possible approach is to use split ->
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)
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)
dlist <- lapply(dlist, within,
{
Serial_test
On 22/03/20 4:01 pm, Thomas Subia via R-help wrote:
Colleagues,
Here is my dataset.
Serial Measurement Meas_test Serial_test
1 17 failfail
1 16 passfail
2 12 passpass
2 8
On 29/02/20 8:11 am, Abby Spurdle wrote:
There was clearly a non-linear growth pattern such that an additive
mixed effects model was proposed to model the behavior of biomass
as a function of time and treatments.
The presence of non-linearity (in x) does not necessarily mean that
you can't
On 23/02/20 11:42 am, greg holly wrote:
Dear R-Help members;
I have the following error messages when I would like to create
training and testing data for Random Forest.
Your help is highly appreciated.
Regards,
Greg
inTrain <- createDataPartition(a, p = 0.7, list = FALSE)
Error: cannot
Dear Yawo,
I would suggest that you learn to use R, rather than thrashing around
blindly and expecting or hoping to get others to do your work for you.
To get you started, the characteristics that you call "labels" are
stored as *attributes* of the columns of your tibble/data frame. E.g.
> should have been read_excel whose parameters need to be adjusted to not look
> for a header line.
And Jeff's probably correct (from the package documentation):
> read_excel(path, sheet = NULL, range = NULL, col_names = TRUE,
Set col_names to FALSE...???
But still be careful with the
My assessment was and is likely correct.
My error was in writing the wrong function name... should have been read_excel
whose parameters need to be adjusted to not look for a header line.
On February 5, 2020 3:48:36 PM PST, Rolf Turner wrote:
>
>On 6/02/20 11:01 am, Thomas Subia wrote:
>
>>
On 6/02/20 11:01 am, Thomas Subia wrote:
Jeff,
You wrote: " Pay attention to whether the read_csv call is configured to expect
first line as header."
Here is the code I'm using to extract one cell from a series of Excel files
having the same physical format.
library(plyr)
library(readxl)
Thanks all for your solutions.
I like the setup best where I startup an xserver connection and then
halting the script with readlines.
I want to write and run small scripts and that seems to me the nearest
to a leightweight solution to the problem in R.
Kind regards,
Felix
On 1/30/20 7:15 AM,
> .prompt = function (pkg, fun)
> { fbody = character (1)
And I've already noticed a bug.
This isn't going too well...
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
> I am using Linux and device X11 worked for me in the past. I will check
> your script later and tell you if it works.
The previous script isn't very good.
I've improved it.
The new version has been tested on Linux and Windows, and with lattice
and ggplot2.
I note that the script is more
On Tue, 28 Jan 2020 at 13:03, Felix Blind wrote:
>
> I am using Linux and device X11 worked for me in the past. I will check
> your script later and tell you if it works.
[...]
I am using R version 3.6.2 on ubuntu 18.04.3 LTS and the X11-device
when I make plot in the terminal works without
Thanks Abby,
I am using Linux and device X11 worked for me in the past. I will check
your script later and tell you if it works.
On 1/28/20 1:27 AM, Abby Spurdle wrote:
> Here's a solution for Windows, for use with graphics::plot.
> Using, Linux and grid or ggplot2, may require some
Here's a solution for Windows, for use with graphics::plot.
Using, Linux and grid or ggplot2, may require some modification.
(However, I assuming the modifications would be easy to make).
You can create a file, say quasi_interactive.r
quasi_interactive.r
#for use with windows only
It is possible to produce interactive-style plots, while running
scripts from the terminal.
(Or at least, it was, because I've done before).
Can I verify, ***which OS*** are you using...???
I'm offsite now.
If you're using Windows or Linux, I can check later.
(If no one else posts the answer,
Thanks a lot you guys,
you all were right; I used most of your suggestions as part of the solution.
It is a considerable effort to get persistent plots when running an R
script.
The problem is even when I open an interface with X11() it closes
immediately afterwards, when the script is done
Second.
Bert
On Sun, Jan 26, 2020 at 2:20 PM Rolf Turner wrote:
>
> 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!!! :-)
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
On 27/01/20 9:47 am, Sorkin, John wrote:
Felix, I suggest you consider using an IDE such as RStudio as you
develop and run R code. An integrated development environment will
allow you to concentrate on learning R rather on the mechanics of
running R in a non-standard environment.
I
Dear All
I could solve the problem: I removed my .Rprofile and I could install
the package via the usual
install.packages("rmgarch")
My .Rprofile contains few settings, therefore, it was not too difficult
to understand what was conflicting with the installation. The issue in
.Rprofile was
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
On 31/12/19 12:29 am, Pietro Coretto wrote:
On 30/12/2019 12.11, Duncan Murdoch wrote:
On 30/12/2019 5:59 a.m., Pietro Coretto wrote:
[...]
You didn't show us the command you used to install it.
Duncan Murdoch
Sorry for this,
I used the following:
install.packages("rmgarch")
from
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
--
Honorary Research
On 8/12/19 9:35 pm, gabriele pallotti wrote:
Dear John,
thank you for your prompt reply. An inexperienced user like me tends to see
the .Rdata folder like the document folder for other programs, and, as one
doesn't have to delete the document folder when updating Libreoffice, tends
to think
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="*.xls", full.names = FALSE)
Agree, especially there is an "Urgent" on the title. He must be too
"urgent" to think about your answer. I will wonder if your effort will be
in vain.
Best,
Jiefei
On Tue, Nov 5, 2019 at 4:52 PM Rolf Turner wrote:
>
> Richard: I know that you mean well, but *please* don't do people's
>
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.
Let's approach it
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
Thank you very much everyone. All fine now!!
On Fri, Oct 18, 2019 at 11:42 AM Rolf Turner
wrote:
>
> On 18/10/19 2:43 PM, ani jaya wrote:
>
> > Dear R-Help,
> >
> > I have a list of data frame that I import from excel file using read.xlsx
> > command.
> >
> > sheets <-
On 18/10/19 2:43 PM, ani jaya wrote:
Dear R-Help,
I have a list of data frame that I import from excel file using read.xlsx
command.
sheets <- openxlsx::getSheetNames("rainfall.xlsx")
test <- lapply(sheets,function(i) read.xlsx("rainfall.xlsx", sheet=i,
startRow=8, cols=1:2))
names(test)
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), "
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
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
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.
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
On 3/08/19 3:51 PM, Bert Gunter wrote:
Spencer:
Sorry, but I'll be blunt. IMO, you are misusing this list (see the
posting guide). You clearly don't know what you're doing statistically
and need to consult with your advisors. This list cannot and is not
meant to serve that purpose -- it is
Hey Anamika,
I only caught the tail end of what became an off-topic thread, but backed up a
bit to your original q. If I'm duplicating anything previous, apologies.
If you are going to ship your "product" to end users directly (vs provide via
an API or web application) I'm not sure how you get
I most firmly do not agree with you, Rolf.
Over the years that I've been a member of this community I've noted the actions
of a small group who have seemed to feel that it is ok to be an asshole to
other people, and a group of hangers-on who have applauded and egged on this
shameful behaviour.
Hi Rolf:
As they say, do read the posting guide:
> Good manners: Remember that customs differ. Some people are very direct.
> Others surround everything they say with hedges and apologies. Be tolerant.
> Rudeness is never warranted, but sometimes `read the manual’ is the
> appropriate
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 example, the search string
"proprietary use GPL" produces one hit that's clearly relevant to you:
This method is more neutrally worded. It doesn't insult the original
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,
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
Who let the editor in here? :-)
Anyway, I stand corrected.
On July 6, 2019 3:37:52 PM PDT, Rolf Turner wrote:
>
>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
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 1/07/19 9:41 PM, Marvin Kiene wrote:
Hello dear helpers,
I am currently running a small R-crash course for beginners at my
university, since I believe that there a far too few lectures about how to
use R.
Thereby, I showed the '??' or 'help.search()' function to the students as
well as the
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
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 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
I will take a look but the hdrcde package appears to give me the plots I was
looking for.
Thanks
Bernard
Sent from my iPhone so please excuse the spelling!"
> On Mar 28, 2019, at 10:16 PM, Jeff Newmiller wrote:
>
> I just found the "ks" package which looks promising... [1]
>
>
I just found the "ks" package which looks promising... [1]
https://cran.r-project.org/web/packages/ks/vignettes/kde.pdf
On March 28, 2019 12:16:12 PM PDT, David Winsemius
wrote:
>
>On 3/28/19 11:33 AM, Bernard McGarvey wrote:
>> pracma meets my needs - thanks - this is an amazing package with
On 3/28/19 11:33 AM, Bernard McGarvey wrote:
pracma meets my needs - thanks - this is an amazing package with lots of very
very useful functions. Thanks for bringing it to my attention.
Excellent. I'll mention also the hdrcde package which addresses some of
your questions as well. It also
pracma meets my needs - thanks - this is an amazing package with lots of very
very useful functions. Thanks for bringing it to my attention.
Bernard McGarvey
Director, Fort Myers Beach Lions Foundation, Inc.
Retired (Lilly Engineering Fellow).
> On March 28, 2019 at 1:40 PM David Winsemius
On 3/27/19 3:43 PM, Bernard Comcast wrote:
To follow on Jeff, is there a function to do 2-D (double) numerical integration
in R?
Packages pracma and cubature offer a variety of solutions to that task.
--
David.
Bernard
Sent from my iPhone so please excuse the spelling!"
On Mar 27,
ing with Flow Cytometry data.
>
> Especially flowViz package is designed to visualise such data.
>
> Cheers
> Petr
>
>
>> -Original Message-
>> From: R-help On Behalf Of Bernard McGarvey
>> Sent: Wednesday, March 27, 2019 9:55 PM
>> To: Paul Mur
R-help On Behalf Of Bernard McGarvey
> Sent: Wednesday, March 27, 2019 9:55 PM
> To: Paul Murrell ; John Kane
>
> Cc: R. Help Mailing List
> Subject: Re: [R] [FORGED] Re: Quantile Density Contours
>
> If I understand correctly the ContourLines function gives you the contour
> li
https://cran.r-project.org/web/views/Multivariate.html
https://cran.r-project.org/web/views/
On March 27, 2019 4:05:31 PM PDT, Bernard Comcast
wrote:
>No - how do I access that?
>
>Bernard
>Sent from my iPhone so please excuse the spelling!"
>
>> On Mar 27, 2019, at 6:57 PM, Jeff Newmiller
>
You might be wishing for a contour plot of the density, labeled by the
probability mass outside of each contour, but there is no general simple
connection between density contours and the mass inside of them. You can work
it out (I think) for elliptically contoured distributions, but I suspect
No - how do I access that?
Bernard
Sent from my iPhone so please excuse the spelling!"
> On Mar 27, 2019, at 6:57 PM, Jeff Newmiller wrote:
>
> I don't know. Have you looked at the Multivariate Task View?
>
>> On March 27, 2019 3:43:52 PM PDT, Bernard Comcast
>> wrote:
>> To follow on Jeff,
I don't know. Have you looked at the Multivariate Task View?
On March 27, 2019 3:43:52 PM PDT, Bernard Comcast
wrote:
>To follow on Jeff, is there a function to do 2-D (double) numerical
>integration in R?
>
>Bernard
>Sent from my iPhone so please excuse the spelling!"
>
>> On Mar 27, 2019, at
1 - 100 of 637 matches
Mail list logo