similar to: Point of intersection

Hi, I like to generate two normal densities such that the overlapping area between them is 0.35. Is there any code/package available in R to do that?? Regards, Lavan -- View this message in context: http://www.nabble.com/densities-with-overlapping-area-of-0.35-tp19384741p19384741.html Sent from the R help mailing list archive at Nabble.com.

Hi, I have a vector or random variables and I'm estimating the density using "bkde" function in the KernSmooth package. The out put contains two vectors (x and y), and the R documentation calls y as the density estimates, but my y-values are not exact density etstimates (since these are numbers larger than 1)! what is y here? Is it possible to get the true estimated density at each

Hi, I posted this message earlier in "Rmetrics" and I don't know whether I posted in the wrong place, so I'm posting it again in Rhelp. I have a function in x and y and let's call it f(x,y). I need to get the Hessian matrix. i.e I need (d^2f/dx^2), (d^2f/dxdy), (d^2f/dydx), (d^2f/dy^2).I can get these using the D function. now I need to evaluste the hessian matrix for

Dear R users, I used sm.density function in the sm package and kde2d() in the MASS package to estimate the bivariate density. Then I calculated the Kullback leibler divergence meassure between a distribution and the each of the estimated densities, but the asnwers are different. Is there any difference between the kde2d and sm.density estimates? if there is a difference, then which is the best

Hi there, I'm trying to find the KL divergence measure between a prior and it's posterior distributions, and I'm using the KLdiv method in the flexmix package. plese see the example below: require(flexmix) x=seq(-4,4,length=100) d1=dnorm(x,0,1) d2=dunif(x,-3,3) y=cbind(d1,d2) kl=KLdiv(y) but let say, x1=seq(-5,5,length=100) d3=dunif(x1,-3,3) y1=cbind(d1,d3) kl1=KLdiv(y1) Notice

Is there a function, or code snippet that anyone is aware of that will give the number of elements in a MECE intersection of x vectors where the total number of vectors is n? I have several vectors, A, B, C, ....,n of varying length. I am trying to figure out a way to systematically calculate the number of MECE elements in all possible combinations of 2 vectors only, then do the same

Hi, I have two matrices as follow: matrix A = a=matrix(c(c("abc","abc","bcd","bcd","bce","bce"),c("a1","d2","d1","d2","a1","a2")),6,2) and matrix B which contains pair of values : b=matrix(c(c("a1","a1"),c("a2","d2")),2,2) In short,

v1<-c("a","b","c","d") v2<-c("a","b","e") v3<-c("a","f","g") I want to get the intersection of v1,v2,v3,ie "a" How can I do then? What I know is only for 2 vectors via "intersect" function,but don't know how to deal with multiple vectors. Many thanks

Hi there, I am working on a project which needs the value of the interaction area under two distributions( eatimated by kernel density estimators). For example: x<-rnorm(100,0,1) y<-rnorm(100,0.2,1) density(x) # This produces the summary of dependent variable and independent variable. How can I get the individual values of variables and reform a curve to calculate the area under

Hi all, As far as I can tell, the only canned way to do an intersect between two vectors of ints is the intersect(vec1, vec2) function -- is there another function I'm missing for intersecting more than two vectors?? TIA, Ken

Hello sir: group1:1,2,3,4 group2:1,3,4,5 group3:2,4,8,9 .. ... group1000:9,3,8,2 I wanna get the "intersection set" and "union set". I've tried command "setdiff"and "union",but only two groups is permited.How can I deal with multi groups to find the "intersection set" and "union set"? Thanks a lot! My best

Hi, Previously I posted a thread asking help on how to best-fit (in the least squares sense) a straight line through a set of data points. Thanks a lot to all replying to it. I managed it in Matlab using a function 'fit_3D_data' (link: http://webscripts.softpedia.com/script/Scientific-Engineering-Ruby/Statistics-and-Probability/Orthogonal-Linear-Regression-in-3D-space-35532.html). But to

I was wondering if there were any R packages supporting fast set intersection. I am aware of base::intersection, bit::`&`, and set::set_intersection, each of which has its advantages. base::intersection is space-efficient and accepts arbitrary unsorted lists of arbitrary base types. set::set_intersection is very fast but requires space and time proportional to the size of the universe, a

How can I create intersections of vectors? a <- c(1,2,3) b <- c(1,5,6) the intersected list c should contain c(1)... Ralf

Hi, I'm a beginner with R. I have two different  tables with the variable name dmr1 and tp2 a given following. v1 is the common  column  field of the both tables. In the first table column v3 is always v2+1 while in the second table v2 and v3 hold the range. I want to know which rows of tb1 intersect with the range between v2 and v3 of tb3. I need some one's kind help to solve this

Dear All, I'm looking an appropriate way in R to compute/estimate? points of intersection between a line and a curve and will really appreciate for any suggestion or ideas? Thank you, Fir

Hi all,     I have matrix A of 67420 by 2 and another matrix B of 59199 by 2. I would like to find the number of rows of matrix B that I can find in matrix A (rows that are common to both matrices with or without sorting). I have tried the "intersection" and "is.element" functions in R but it only working for the vectors and not matrix i.e,    intersection(A,B) and

Hi, I want to merge multiple chromosomal regions based on their common intersecting regions. I tried couple of things using while and if loops but did not work out. I would appreciate if anyone could provide me a small piece of code in R to get the intersection of following example: chr1: 100-150 chr1: 79-250 chr1: 100-175 chr1: 300-350 I want the intersection of all four regions as follow:

Hello, This code builds a simple example of 2 wireframes : require(lattice) x <- c(1:10) y <- c(1:10) g <- expand.grid(x = 1:10, y = 1:10, gr = 1:2) g$z <- c(as.vector(outer(x,y,"*")), rep(50,100)) wireframe(z ~ x * y, data = g, groups = gr, scales = list(arrows = FALSE)) However, the intersection between the wireframes is not properly drawn. Is there a way to fix this

I've got two tables.... first one(table1): ID chrom start end Ex1 2 152 180 Ex2 10 2000 2220 Ex3 15 3000 4000 second one ( table2): chrom location name 2 160 Alv 2 190 GNN 2 100