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 <jdnew...@dcn.davis.ca.us> wrote: > > I don't know. Have you looked at the Multivariate Task View? > >> On March 27, 2019 3:43:52 PM PDT, Bernard Comcast >> <mcgarvey.bern...@comcast.net> 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 6:38 PM, Jeff Newmiller >> <jdnew...@dcn.davis.ca.us> wrote: >>> >>> Regardless of how many dimensions you have for independent variables, >> the density is one-dimensional, and if you assume the density function >> has been determined (e.g. by kernel estimation or by a Gaussian copula) >> then if you integrate the density function along that dimension there >> will be unique slices of the multivariate input domain determined by >> those slices. They might in general be disjoint regions of the >> independent variable space, but that is what the contour function does. >>> >>> I am not seeing your point, Bert, unless you are unwilling to assume >> a density function model? >>> >>>> On March 27, 2019 2:18:18 PM PDT, Bert Gunter >> <bgunter.4...@gmail.com> wrote: >>>> You are missing a crucial point. The reals are well ordered; higher >>>> dimensions are not. Therefore 2d quantile contours are not unique. >>>> >>>> Of course assuming I understand your query correctly. >>>> >>>> >>>> Bert >>>> >>>> On Wed, Mar 27, 2019, 13:55 Bernard McGarvey >>>> <mcgarvey.bern...@comcast.net> >>>> wrote: >>>> >>>>> If I understand correctly the ContourLines function gives you the >>>> contour >>>>> lines when you put in the data. But before this I need to data to >> put >>>> into >>>>> that function. I think this is something like a 2D CDF of the data >>>> that >>>>> then leads to the 2D quantiles but I am not 100% sure. What I am >>>> basically >>>>> looking for is the 2D curve that encloses say 95% of the data in a >>>> similar >>>>> fashion to a 1D quantile where the quantile represents the value >> that >>>> x% of >>>>> the data is below. I think what I am looking for is the 2D >> bivariate >>>>> version of the 1D quantile plot (where the quantile value is >> plotted >>>> vs the >>>>> % value). >>>>> >>>>> I hope this makes some sense. >>>>> >>>>> Bernard McGarvey >>>>> >>>>> >>>>> Director, Fort Myers Beach Lions Foundation, Inc. >>>>> >>>>> >>>>> Retired (Lilly Engineering Fellow). >>>>> >>>>> >>>>>> On March 27, 2019 at 3:57 PM Paul Murrell >>>> <p...@stat.auckland.ac.nz> >>>>> wrote: >>>>>> >>>>>> >>>>>> >>>>>> Are you looking for the contourLines() function ? >>>>>> >>>>>> Paul >>>>>> >>>>>>> On 28/03/19 8:37 AM, Bernard McGarvey wrote: >>>>>>> John, I have attached a pdf of the plot. Hopefully you can read >>>> this. >>>>>>> >>>>>>> If I understand correctly, this plot is basically the 2-D version >>>> of >>>>> the 1-D quantile plot. >>>>>>> >>>>>>> Thanks >>>>>>> >>>>>>> Bernard McGarvey >>>>>>> >>>>>>> >>>>>>> Director, Fort Myers Beach Lions Foundation, Inc. >>>>>>> >>>>>>> >>>>>>> Retired (Lilly Engineering Fellow). >>>>>>> >>>>>>> >>>>>>>> On March 27, 2019 at 7:44 AM John Kane <jrkrid...@gmail.com> >>>> wrote: >>>>>>>> >>>>>>>> >>>>>>>> The figure did not get through. Perhaps try a pdf? >>>>>>>> >>>>>>>> On Tue, 26 Mar 2019 at 13:41, Bernard McGarvey >>>>>>>> <mcgarvey.bern...@comcast.net> wrote: >>>>>>>>> >>>>>>>>> I want to see if I can reproduce the plot below in R. If I >>>>> understand it correctly, i takes my bivariate data and creates >>>> quantile >>>>> density contours. My interpretation of these contours is that they >>>> enclose >>>>> a certain % of the total data. I am using the bkde2D function in >>>> library >>>>> KernSmooth which gives density values that can be plotted on a >>>> contour plot >>>>> but I would like the curves that enclose a given % of the data, if >>>> that is >>>>> possible >>>>>>>>> >>>>>>>>> >>>>>>>>> Thanks >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> Bernard McGarvey >>>>>>>>> >>>>>>>>> Director, Fort Myers Beach Lions Foundation, Inc. >>>>>>>>> >>>>>>>>> Retired (Lilly Engineering Fellow). >>>>>>>>> >>>>>>>>> >>>>>>>>> ______________________________________________ >>>>>>>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, >>>> see >>>>>>>>> https://stat.ethz.ch/mailman/listinfo/r-help >>>>>>>>> PLEASE do read the posting guide >>>>> http://www.R-project.org/posting-guide.html >>>>>>>>> and provide commented, minimal, self-contained, reproducible >>>> code. >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> -- >>>>>>>> John Kane >>>>>>>> Kingston ON Canada >>>>>>>> >>>>>>>> ______________________________________________ >>>>>>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, >>>> see >>>>>>>> https://stat.ethz.ch/mailman/listinfo/r-help >>>>>>>> PLEASE do read the posting guide >>>>> http://www.R-project.org/posting-guide.html >>>>>>>> and provide commented, minimal, self-contained, reproducible >>>> code. >>>>>> >>>>>> -- >>>>>> Dr Paul Murrell >>>>>> Department of Statistics >>>>>> The University of Auckland >>>>>> Private Bag 92019 >>>>>> Auckland >>>>>> New Zealand >>>>>> 64 9 3737599 x85392 >>>>>> p...@stat.auckland.ac.nz >>>>>> http://www.stat.auckland.ac.nz/~paul/ >>>>> >>>>> ______________________________________________ >>>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>>>> https://stat.ethz.ch/mailman/listinfo/r-help >>>>> PLEASE do read the posting guide >>>>> http://www.R-project.org/posting-guide.html >>>>> and provide commented, minimal, self-contained, reproducible code. >>>>> >>>> >>>> [[alternative HTML version deleted]] >>>> >>>> ______________________________________________ >>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>>> https://stat.ethz.ch/mailman/listinfo/r-help >>>> PLEASE do read the posting guide >>>> http://www.R-project.org/posting-guide.html >>>> and provide commented, minimal, self-contained, reproducible code. >>> >>> -- >>> Sent from my phone. Please excuse my brevity. > > -- > Sent from my phone. Please excuse my brevity. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.