similar to: Find x value of density plots

Dear All, I was wondering whether someone has created a kernel density evaluator that estimates the density at given specified points. The regular density() function evaluates the kernel at equidistant points, but I am interested in doing such evaluation along a list of values existing in a pre-specified vector. (Similar to the option at() in the kdensity command in Stata). This question has

Hello, I have a data x with normal (or very close to normal) distribution, I can plot a density distribution with density(x,...). My question is is there any way to calculate an area under this distribution (=probability) for particular range of x values, let's say for x from 0 to 2? I was not able to find any kind of simple procedure to do this. Thanks in advance for your help, Evgeny.

Hi all, I am looking for a command for doing 2D filtering (rectangular or Gaussian) in R... I have looked at ksmooth, filter and convolve but they seem to be 1D... Any thoughts? Thanks a lot! [[alternative HTML version deleted]]

Folks, I have a small dataset of counts of recoveries on defaulted loans: recoveries<-structure(c(0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 0, 0, 0, 0, 0, 4, 0, 1, 2, 2, 12), .Dim = c(11L, 2L), .Dimnames = list( NULL, c("pcts", "counts"))) Here is the data in columnar form: pcts counts [1,] 0.0 0 [2,] 0.1 0 [3,] 0.2 0 [4,] 0.3

Dear All, would you please provide your thoughts on the following: let us say I have: a <-c(1,5,8,15,32,69) b <-c(8.5,33) and I would like to extract from "a" the two values that are closest to the values in "b", where the length of this vectors may change but b will allways be shorter than "a". So at the end based on this example I should have the result

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130513/fe848ce7/attachment.pl>

Hello, I am brand new in this site as well as in R. Please accept my early apology if it's not the right place to ask such a question. I have the following equation: b=[2(1-a)(1+2a)^0.5]/(1+3a) How can I solve the above equation for 'a' using R? Thanks in advance. KD [[alternative HTML version deleted]]

Following the beginning tutorial, I typed help.start() and was given a choice of browsers. I picked chrome but got back that chrome is not found. I cannot seem to change or get rid of it now. I looked at /etc/R/Renviron but don't see anywhere to set browser (R_BROWSER = ${R_BROWSER ...etc) - nothing about chrome. I tried starting as "R_BROWSER= R" - same thing. How do I change the

Hi all, I am trying to estimate the mode of a 4-dimensional nonparametric density estimator (any) using a sample of size n=10,000. I have tried using the package 'ks' and 'np' but they are extremely slow; this is related to the estimation of the bandwidth matrix. I also checked the package 'modeest' but it contains only methods for univariate distributions. I am only

Dear all, I am a R beginner, and I am looking for a way to do the same thing for all levels of a column in a table. Basically, I have a bunch of protein sequences composed of different amino acid residues, and each residue is represented by an uppercase letter. I want to calculate the ratio of different amino acid residues at each position of the proteins. Here is an example table: Proteins

Hi, Instead to put colour in my histogram, I want to put symbols like lines, dots etc. Do you know the function that does it? Thank you in advance -- View this message in context: http://r.789695.n4.nabble.com/Not-colour-but-symbols-tp4490030p4490030.html Sent from the R help mailing list archive at Nabble.com.

Hi all, I am using "locator" to select the points from a scatter plot... This is all fine. But the problem is that the locator only returns the axis values of the selected points. Instead, I would like to get the index of these select points... The axis values are real-values so it's a bit hard for me to directly reverse-engineer the index nubmers.. How to do that? Thank you!

Dear all I have a question about quantiles standard error, partly practical partly theoretical. I know that x<-rlnorm(100000, log(200), log(2)) quantile(x, c(.10,.5,.99)) computes quantiles but I would like to know if there is any function to find standard error (or any dispersion measure) of these estimated values. And here is a theoretical one. I feel that when I compute median from given

I have a data in the following form : CIN TRN_TYP 9079954 1 9079954 2 9079954 3 9079954 4 9079954 5 9079954 4 9079954 5 9079954 6 9079954 7 9079954 8 9079954 9 9079954 9 . . . . . . there are 100 types of CIN (9079954,12441087,15246633,...) and respective TRN_TYP first of all, I want this data to be

suppose I have two sets of vectors: x1,x2,...,xN and y1,y2,...,yN. I want N correlations: cor(x1,y1), cor(x2,y2), ..., cor(xN,yN). my sets of vectors are arranged as data frames x & y (vector=column): x <- data.frame(a=rnorm(10),b=rnorm(10),c=rnorm(10)) y <- data.frame(d=rnorm(10),e=rnorm(10),f=rnorm(10)) cor(x,y) returns a _matrix_ of all pairwise correlations: cor(x,y)

Hello all, I cannot seem to figure out this seemingly simple procedure. I want to sort a data frame by a specified character vector. So for : df.. <- data.frame(Season=rep(c("Summer","Fall","Winter","Spring"),4),Obs= runif(length(rep(c("Summer","Fall","Winter","Spring"),4)))) I want to sort the data frame

I have a stupid problem that is currently driving me crazy... Let us suppose that I want to draw a big red square in the middle of my current device (say X11) I tried the following code : pushViewport(viewport(xscale=c(0,1), yscale=c(0,1)), just=c("center", "center")) vp1 <- viewport(x=unit(0.5, "native"), y=unit(0.5, "native"), width=unit(0.4,

Hello, Does the (core) designation in the CRAN Task View indicate packages that provide basic/recommended functionality in that domain? I don't see it documented anywhere on the site. Thanks for your help.

Hi, is there a way to replicate a data.frame like you can replicate the entries of a vector (with the repeat-function)? I want to do this: x <- data.frame(x, x) (where x is a data.frame). but n times. And it should be as cpu / memory efficient as possible, since n is pretty big in my case. thanks for any suggestions!

Hello :) I'm just starting out with R and would appreciate your help with a couple of problems I am running into. I have used Sys.glob to get a list of all filenames having a particular file extension (in my case, *.txt) I would now like to use this list in the following manner: I would like to use each filename from the list and open that file into a tab separated matrix and proceed. How can