search for: densitiy

Hi, I am about to write functions for multivariate kernel densitiy estimation with mixed categorical and continuous date (accoring to Jeff Racine and Qi Li), and the leave-one-out window esitmation needs a lot of computation. I am now optimizing the code performance and therefore fhe following questions: As R uses call-by-value for functions, is it computational...

Hi My simulation with the followin R code works perfectly: sim <- replicate(999, sum(exp(rgamma(rpois(1,2000), scale = 0.5, shape = 12)))) But now I do not want to have values in object "sim" exceeding 5'000'000, that means that I am just using the beginning of densitiy function gamma x < 15.4. Is there a possibility to modify my code in an easy way? Thanks for any help! Regards, Brigitte [[alternative HTML version deleted]]

...,kernel="epanechnikov",n=n) #evaluation print(f$y[5]) Here I can only evaluate the estimation at given points. These points are determined by the parameter n. By default they are equidistant distributed on the interesting interval. But I need to evaluate the estimation (the estimated densitiy function) at manually specified points. For example I want to compute f(z[i]). This means I am interested in the estimated density at a the observation z[i]. Does anyone know how I can compute this? I think this is an ordinary task so I would be surprised if R can not manage this. But even after a...

Hello, I'm writing some genetic simulations in R where I would like to place genes along a chromosome proportional to the density of markers in a given region. For example, a chromosome can be presented as a list of marker locations: Chr1<-c(0, 6.5, 17.5, 26.2, 30.5, 36.4, 44.8, 45.7, 47.8, 48.7, 49.2, 50.9, 52.9, 54.5, 56.5, 58.9, 61.2, 64.1) Where the numbers refer to the locations of

Dear list: I am attempting to use what I thought would be a pretty straightforward practical application of Cox regression. I figure users of the survival package must have come across this problem before, so I would like to ask you how you dealt with it. I have set up an illustrative example and included it at the end of this post. I took a sample of 100 data points from each of two populations

...mulation with the followin R code works perfectly: > sim <- replicate(999, sum(exp(rgamma(rpois(1,2000), scale = 0.5, shape = > 12)))) > > But now I do not want to have values in object "sim" exceeding 5'000'000, > that means that I am just using the beginning of densitiy function gamma x < > 15.4. Is there a possibility to modify my code in an easy way? > > Thanks for any help! > > Regards, Brigitte > > > > > > > > [[alternative HTML version deleted]] > > > > ------------------------------ > > Messa...