search for: nrd

Dear all, Sorry if I am posting to the wrong place, but I could not find the link for registration on the bug tracker, that?s why I am writing here. I think there is inconsistency between two R functions from the stats package, bw.nrd0 and bw.nrd. Consider the following vector: D <- c(0, 1, 1, 1, 1) bw.nrd(D) returns zero bandwidth for this object even without a warning. Considering the fact that in most cases, one divides something by the bandwidth, it is highly undesirable that these function return zero bandwidth without...

Good evening, I am Marta Colombo, student at Milan's Politecnico. Thank you very much for your kindness, this mailing list is really useful. I am using the function kde2d for two-dimensional kernel density estimation and I'd like to know something more about this kind of density estimator. In particular I'd like to know: what bandwidth is used ? Thank you in advance for your attention

Dear R Users, Before running ksmooth( ), a suitable bandwidth selection is needed. I use some functions for this task and receive these results for my data: width.SJ(y,nb=100,method="ste") : 40.25 bcv(y,nb=100) : 40.53 ucv(y) : 41.26 bandwidth.nrd(y) : 45.43 After implementing the function ksmooth(x,y, bandwidth= each of abovementioned bandwidths), I have some NAs in output with regard to each bandwidth. But if the bandwidth be equal to 62 or bigger, then the function ksmooth(x,y, bandwidth=62 or bigger) works without NA,i.e., the ba...

...ims = c(range(x), range(y))) { nx <- length(x) if (length(y) != nx) stop("data vectors must be the same length") gx <- seq(lims[1], lims[2], length = n) # gridpoints x gy <- seq(lims[3], lims[4], length = n) # gridpoints y if (missing(h)) h <- c(bandwidth.nrd(x), bandwidth.nrd(y)); if (missing(w)) w <- numeric(nx)+1; h <- h/4 ax <- outer(gx, x, "-")/h[1] # distance of each point to each grid point in x-direction ay <- outer(gy, y, "-")/h[2] # distance of each point to each grid point in y-direction z <-...

...is a weight vector, you must have sum(w)=length(x) nx <- length(x) if (length(y) != nx) stop("Data vectors must be the same length") gx <- seq(lims[1], lims[2], length = n) gy <- seq(lims[3], lims[4], length = n) if (missing(h)) h <- c(bandwidth.nrd(x), bandwidth.nrd(y)) if(missing(w)) w <- rep(1,nx) h <- h/4 ax <- outer(gx, x, "-")/h[1] ay <- outer(gy, y, "-")/h[2] z <- matrix(dnorm(ax), n, nx) %*% diag(w) %*% t(matrix(dnorm(ay), n, nx)) / (nx * h[1] * h[2]) return(list(x = gx,...

...nt methods for bandwidth determination. Below are the different ways in R on deciding the bandwidth. Can anyone give an idea on which ones are preferred. Also, how can I take a look at the source code for the density function? Thank you very much. Hannah x <- rnorm(1000) > bw.nrd(x) [1] 0.2688588 > bw.nrd0(x) [1] 0.2282763 > bw.ucv(x) [1] 0.2112366 > bw.bcv(x) [1] 0.2890085 Warning message: In bw.bcv(x) : minimum occurred at one end of the range > bw.SJ(x) [1] 0.2716242 > density(x, give.Rkern=T, kernel="gaussian") [1] 0.2820948 > d...

Hello, JIT compiler interferes with internal R completions: compiler::enableJIT(2) utils:::functionArgs("density", '') gives: utils:::functionArgs("density", '') Note: no visible global function definition for 'bw.nrd0' Note: no visible global function definition for 'bw.nrd' Note: no visible global function definition for 'bw.ucv' Note: no visible global function definition for 'bw.bcv' Note: no visible global function definition for 'bw.SJ' Note: no visible global funct...

We've been trying unsuccessfully to connect our Meridian Option 81C to a TE110P via PRI. We've followed the directions in asterisk-meridian-a1.pdf (link on http://www.voip-info.org/wiki/view/Asterisk+legacy+integration), but it doesn't seem to work on our 81C (even though many, many users report it works very well on Option 11's). Has anyone had any success in getting the above

Problem: Based on a n x 2 data matrix i want a kernel estimate of the bivariate density. However, i also wish to specify wich points the density should be calculated at. I can offcourse just write the full kernel density estimate as a R-code, but surely there must already exist some package for this operation? The package density(), seems to create a new matrix (depending on n), where the

Hi, This message is to let the R developers know that the pdfs for the first two volumes of the printed edition of R reference manual are available at, http://www.network-theory.co.uk/R/base/vol1.pdf http://www.network-theory.co.uk/R/base/vol2.pdf I've divided the manual into "core commands" of the base package (control flow, and related programming commands, etc) in Volume 1

Hi all, Since two days I am trying to find a solution to be able to create beanplots for my data. When I call the beanplot function the following error appears: > beanplot(y1 ~ x1, log="", what=c(1,1,1,0), ylim=c(0,1)) > Error in bw.SJ(x, method = "dpi") : sample is too sparse to find TD What is really strange: I have 32 different vectors and the problem occurs for

The following is from version 1.6.2 of R under Windows, or 1.6.1 under Mac OSX/X11 > density(1) Error in if (!(lo <- min(hi, IQR(x)/1.34))) (lo <- hi) || (lo <- abs(x[1])) || : missing value where logical needed I am not sure how this should be handled. I encountered it in connection with densityplot(). In that connection, it might be enough to modify density() so that it

...CPUB OFF DAPC NO BCOT 0 DSEL VOD PTYP PRI AUTO NO DNIS NO DCDR NO ICOG IAO SRCH LIN TRMB YES STEP ACOD 7930 TCPP NO PII NO TARG CLEN 1 BILN NO OABS INST IDC NO DCNO 0 * NDNO 0 DEXT NO ANTK SIGO STD ICIS YES TIMR ICF 512 OGF 512 EOD 13952 NRD 10112 DDL 70 ODT 4096 RGV 640 GRD 896 SFB 3 NBS 2048 NBL 4096 IENB 5 TFD 0 VSS 0 VGD 6 DRNG NO CDR NO VRAT NO MUS NO FRL 0 0 FRL 1 0 FRL 2 0 FRL 3 0 FRL 4 0 FRL 5 0 FRL 6 0 FRL 7 0 OHQ NO OHQ...

...` qq{kjs b jnmve ohq|l` rsphqrhweqjne `cemrqrbn (khv. 9 28288) Lnqjb`, Vemrp, Qbepwjnb oep., 1, nthq 443 rek.:(095) 208-0257, 208-0261, 207-9974 r/t`jq: (095) 208-0375 l."Jhr`i Cnpnd", "Whqr{e Opsd{" Hlhdf - Bq8! L{ me dna`bkel mhj`jni qnaqrbemmni m`vemjh! Dk B`q - nrd{u Dk B`xei qel|h - p`dnqr| Dk qnrpsdmhjnb B`xei jnlo`mhh - qelhm`p m` jspnpre @JRS@K\MN! Kern 2001 Rspvh Echoer Tp`mvh Ankc`ph Jp{l nr 280 s.e. nr 290 s.e. nr 500 s.e. nr 220 s.e. nr 180 s.e. Jhop Rsmhq Hr`kh J`pha{ J`bj`g nr 300 s.e. nr 435 s.e. nr 550 s.e. nr 1200 s.e. nr 150 s.e. Hqo`mh...

...hat the AIX platform seems to have similar issues. It would seem sso> that the same may be true for IRIX and Tru64 UNIX. sso> Andrew Orndorff sso> Andrew R. Orndorff sso> Technical Support Manager sso> NHTSA R&D Data Center sso> U.S. Dept. of Transportation sso> NHTSA / NRD-10 sso> 400 7th Street, SW, Rm. 2403 sso> Washington, DC 20590 sso> Phone: 202.366.8827 sso> FAX: 202.366.3986 sso> E-mail: andrew.orndorff@nhtsa.dot.gov sso> This message is not representative of the opinions, sso> practices, or policies of ISSI, the National Highway sso&g...

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML

...amba.org> List-Unsubscribe: <http://lists.samba.org/mailman/listinfo/samba>, <mailto:samba-request@lists.samba.org?subject=unsubscribe> List-Archive: http://lists.samba.org/pipermail/samba/ P`qq{kj` B`xei hmtnpl`vhh on }kejrpnmm{l `dpeq`l 30 r{qw (!) npc`mhg`vhi Lnqjb{. Jnknqq`k|m` nrd`w`. Nr $120. info_marketus@yahoo.com Return-Path: <nikolaus@erlm.siemens.de> Delivered-To: samba@lists.samba.org Received: from atlantis.erlm.siemens.de (atlantis.erlm.siemens.de [212.114.202.14]) by lists.samba.org (Postfix) with ESMTP id E7F1E5466 for <samba@lists.samba.org>; T...