similar to: Density Estimation

Dear R users: I have a dataset with two variables (>20000 observations, two samples from same subject) and I used "kernSur" from library(Genkern) to get a estimated bivariate density and corresponding plots as follows: new.data.normal<-data.normal[!is.na(data.normal[,2]),] x<-new.data.normal[,2] y<-new.data.normal[,3] op <- KernSur(x,y, xgridsize=50, ygridsize=50,

Hi all I'm trying to create a three-dimensional wireframe-type plot, much like the screenshot shown at http://www.r-project.org/screenshots/volcano-persp.jpg In my case, I'm trying to show a density estimate of two-variable system, using the GenKern package. I've used the image() and contour() to create a nice shaded plot, but I've found a 3d plot is often easier for the

Hi Spencer, I have seen your name on the web site, and perhaps you can help me with my R problem. I'm trying to use KernSur to put in evidence a substructure in a bidimensional plot. My problem is that, in order to get the density in the low density areas (in which the substructure is located) I should use different bandwidths. How I can do that? Also, I think that the best choice for my

Is there code for bivariate kernel density estimation? For bivariate kernels there is kde2d in MASS kde2d.g in GRASS KernSur in GenKern (list probably incomplete) but none of them seems to accept a weight parameter (like density does since R 2.2.0) -- Erich Neuwirth, University of Vienna Faculty of Computer Science Computer Supported Didactics Working Group Visit our SunSITE at

Hi, I have two separate questions relating to the KernSmooth package. I am using the dpik function from the KernSmooth package and receive the error Warning message: In kappam * Gcounts : longer object length is not a multiple of shorter object length I saw an earlier post , but the issue was using the bkde fxn and the person appeared to be using too small of a bandwidth.

Thank you for your reply. I know the x in dpik() means the vector. But I don't know how to import into c() with a huge metadata (>1000). Following is my some try, and the h is: [1] 0.001180569, which seems to be feasible.

Hi there and thanks in advance. Nowadays I am working on the plug-in bandwidth selection with R. Firstly, my 1010 data is the return rate from Yahoo Finance. Secondly, my code is following: > r=read.table("/Users/user/Desktop/research/a.txt",sep=",",header=TRUE) > x<-r[8:1010,] > library(KernSmooth) >

Hi I hope this is going to the right place. I am trying to write a program which uses KernSmooth library to estimate mutual information between two time series at various different lags. At the moment it’s producing negative values, which is supposed to be impossible (something is fishy). I am summing across one row of the matrix to get p(value is in bin x) and summing across the columns to get

Hi All: Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources. Best wishes, ANDREW tseries: Package for time series analysis Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6 Depends: ts, mva, quadprog Date: 2001-08-27 Author: Compiled by Adrian

Hi, I am using the KernSmooth package to estimate nonparametrically bivariate density functions. However, it seems that the bandwidths (one for each co-ordinate direction) have to be selected manually. This does not apply for the univariate case, for which dpik (included in KernSmooth) uses up-to-date plug-in rules. Does anyone know about a package, or function, which estimates bandwidths

Hi, I'm relatively new to R and R development, so please forgive me for any obvious errors. What I am trying to do is use the command dpik within the package KernSmooth to estimate bandwidth parameters for GPS telemetry data. I have been able to get this to work on a case by case basis without any problem, but would like to extend this so that I can batch process many different animals for

Epanechnikov kernel works if option canonical=TRUE, however it would be good to know why it does not for for canonical=FALSE (default). Sorry for craetaing maybe useless thread. Best regards, Ondra. [[alternative HTML version deleted]]

Are there any R function or packages that generate wind roses? These are ``star''-like plots that summarize the frequency, direction and magnigude of windspeeds. Thanks, Matt -- Matt Pocernich University of Colorado - Denver Department of Mathematics -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear R helpers I am having recovery rates as given below and I am trying to estimate the Loss Given Default (LGD) and for this I am using Kernel Density estimation method. recovery_rates = c(0.61,0.12,0.10,0.68,0.87,0.19,0.84,0.81,0.87,0.54,0.08,0.65,0.91, 0.56,0.52,0.30,0.41,0.24,0.66,0.35,0.36,0.64,0.55,0.43,0.36,0.28,0.89,0.11,0.23,0.07,

Hello, I've been using the KernSmooth package recently and think I have found a problem with it: after loading the library I can issue bkde(c(27,26,27), bandwidth=dpik(c(27,26,27)), range.x=c(4.4, 113.6), gridsize=128, truncate=T) and bkde returns an error. If I change the gridsize to 129 the function works perfectly. I have tried this on my Linux box, and on a nearby Solaris machine, both

I'm sure this question is answered in the help file, but likely I'm not reading it corrected. Running windows version 1.00.0, loading a table (35K rows by 10 columns) from Excel using the read.table command I receive the following message. Error: cons memory (350000 cells) exhausted See "help(Memory)" on how to increase the number of cons cells. >From reading the

On 27 April 2017 at 12:01, Johannes Ranke wrote: | | > so it seems to me this must affect all packages in Debian sid that were | > built before the release of R 3.4.0! | | or rather before 14 April 2017, which is when R from revision r72510 was | uploaded to sid as pre-release candidate. Another example with KernSmooth: > library(KernSmooth) KernSmooth 2.23 loaded Copyright

Hi, I want to fit a kernel density estimator by bkde of library KernSmooth. I need only the density value at the point 0. I do not understand the following behaviour: > q <- rnorm(100) > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(0,0)) Error in 0:L : NA/NaN argument > bkq <- bkde(q, bandwidth=0.11, gridsize=1, range.x=c(-1,1)) > bkq $x [1] -1 $y [1] NA > bkq

Am Dienstag, 25. April 2017, 11:21:31 schrieb Dirk Eddelbuettel: > On 25 April 2017 at 16:11, Johannes Ranke wrote: > | This looks similar to what I got this morning when I tested my > | (unreleased) > | backport of R 3.4.0 to Debian jessie. My test was > | > | library(MASS) > | example(rlm) > | > | and there was an object that was not found. I am on a train on the way

Hi, I was using a density estimation function as follows: > est <- KernSmooth::bkde(x3, bandwidth=10) When setting bandwidth less than 5, I got the error "longer object length, is not a multiple of shorter object length in: kappa * gcounts ". I wonder if there is anybody who can explain the error for me? Thanks! Hui