similar to: Cluster Analysis: Density-Based Method

Hi All , For an academic project I am trying to do the following Step 1 ) Draw and cluster a N ( lets say 3 ) column dataset by dbscan algorithm using R-project’s fpc package ( let say they are "training clusters" ) , Using dbscan as number of clusters are not know before hand Step 2 ) Once that is done i want to spread some new data points in the

R 1.7.1 on Windows XP The "normal Q-Q plot" produced by plot.lm() mislabels points when the model is fitted using na.action=na.exclude. Example: x <- 1:50 y <- x + rnorm(50) y[c(5,10,15)] <- NA # insert some NA's y[40] <- 50 # add an outlier plot(lm(y ~ x, na.action=na.omit)) # outlier correctly labeled in all # four plots

What would be nice is to have something like a "robust" task view... Andy From: Berton Gunter > > Many thanks for this Martin. There now are several packages > with what appear to be overlapping functions (or at least > algorithms). Besides those you mentioned, "robust" and > "roblm" are at least two others. Any recommendations about > how or

Dear R-help, Can some one tell me what could be the problem? I downloaded the R-1.7.1 source and try to compile it on WinXP. "make", "make bitmapdll", "make tcl", "make recommended" all ran w/o problem. However, when I ran "make docs", it failed at creating refman.pdf, with no apparent error that I could see. I have attached the log file

Hello everyone, When looking for information about clustering of spatial data in R I was directed towards DBSCAN. I've read some docs about it and theb new questions have arisen. DBSCAN requires some parameters, one of them is "distance". As my data are three dimensional, longitude, latitude and temperature, which "distance" should I use? which dimension is related to

Hi, I am running Samba 3.6.19 and making my way through all the Active Directory documentation. I have joined various Linux boxes to a Windows Server 2008 R2 Active Directory domain. Windows 7 users who log on with their domain accounts can mount the correct shares exported by the Linux boxes without having to supply any additional credentials. All working perfectly! However, I have a

> From: Thomas Lumley [mailto:tlumley at u.washington.edu] > > Please don't send these messages to r-announce. r-announce is for > announcements. Isn't r-announce supposed to be moderated? I posted to r-announce about a minor update in randomForest and was rejected. I was told that it's meant for new or significant upgrades in packages, etc. Andy

Dear R-help, A colleague of mine was running some code on two of our boxes, and noticed a rather large difference in running time. We've so far isolated the problem to the difference between R 1.7.1 and 1.8.0, but not more than that. The exact same code took 933.5 seconds in 1.7.1, and 3594.4 seconds in 1.8.1, on the same box. Basically, the code calls boot() to bootstrap fitting mixture

R-listers, In using the adapt function, I am getting the following warning: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(ndim = 2, lower = lower.limit, upper = upper.limit, functn = pr.set, Would someone explain what the 'lenwrk' value indicates in order to help diagnose this issue. Also, what are the possible codes for Ifail, so I can set

Hello Bug people, I have an unexpected behavior and am unsure whether the problem is in my thinking, my implementation or the program R. Basically I get two different answers depending on how I call a function which takes other functions as arguments as indicated below. To me it should make no difference if f is a function that returns the function g then z(f(x)) whould give the same as y<-

Hi, I'm trying to implement the DBSCAN algorithm to get O(N*LogN) complexity and I'd need a spatial tree of some sort (kd,r,bd..), or a function that computes radius-based search on spatial data, i.e. given a radius eps finds ALL the points which fall in the corresponding hypersphere centered on the current examined point. Is there a package with this features? So far I found RANN

Dear People, I'm trying to use the function adapt, from the adapt library package, which does multidimensional numerical integration. I think I must be using the wrong syntax or something, because even a simple example does not work. Consider foo <- function(x){x[1]*x[2]} and adapt(2, lo = c(-1,-1), up = c(1,1), functn = foo) This simply hangs. A more complicated example crashes R,

I am trying calculate a probability using numerical integration. The first program I ran spit out an answer in a very short time. The program is below: ## START PROGRAM trial <- function(input) { pmvnorm(lower = c(0,0), upper = c(2, 2), mean = input, sigma = matrix(c(.1, 0, 0, .1), nrow = 2, ncol = 2, byrow = FALSE)) } require(mvtnorm) require(adapt) bottomB <- -5*sqrt(.1) topB <-

I tried to integrate numerically a function wich is similar to the following: > dummy <- function(x) { exp(-1*x) * dnorm(x) } > dummy(-100) [1] 0 > dummy(-1000) [1] NaN > dummy(-10000) [1] NaN If I choose the lower boundary to be too small integrate causes a segmentation fault: > library(integrate) > integrate(dummy, -100, 0)$value [1] 1.387143 > integrate(dummy, -1000,

Hi there. Thanks for your time in advance. I am using R 2.2.0 and OS: Windows XP. My final goal is to calculate 1/2*integral of (f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes: $\frac{1}{2}\int^{{\infty}}_{\infty} (\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx $.) where f1(x) and f2(x) are two marginal densities. My problem: I have the following R codes using "adapt" package. Although "adapt"

[Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets created. Also, is it possible to make that an option that can be set at the R prompt, or perhaps even better, as the default no matter how/where R is started? For those of

[Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets created. Also, is it possible to make that an option that can be set at the R prompt, or perhaps even better, as the default no matter how/where R is started? For those of

Hi all I was wondering if someone could help me. I have to estimate some parameters, so I am using the function nlm. Inside this function I have to integrate, hence I am using the function adapt. I don't understand why it is giving the following warnings: At the beginning: Warning: a final empty element has been omitted the part of the args list of 'c' being evaluated was:

Dear all, I get the following error message. And I cannot quite work out what is wrong. I think the optim gets infinite values. Certainly my data do not have any infinite values. How can I solve this? fitdistr(A1, "weibull") Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) : non-finite value supplied by optim I am using R version 1.9.1 on RedHat Linux, Kernel 2.6.8.

Hi I have a FORTRAN version of the L-BFGS-B algorithm and I was comparing it to the code in the lbfgsb.c file available at R-4.0.0.tar.gz Everithing looks the same, except for those two lines that must be printed by the prn3lb function in case of an error (lines 3559 and 3561 in lbfgsb.c): case -5: Rprintf("l(%d) > u(%d). No feasible solution", k, k); break; case -7: