similar to: R keeps crashing when executing 'rlogspline'

Dear R-list, I have made a simple kernel density estimation by x <- c(2,1,3,2,3,0,4,5,10,11,12,11,10) kde <- density(x,n=100) Now I would like to know the estimated probability that a new observation falls into the interval 0<x<3. How can I integrate over the corresponding interval? In several R-packages for kernel density estimation I did not found a corresponding function. I

Hi all, I've successfully compiled R-1.2 on a Linux box (Mandrake 7.1). However, when I installed packages from sources, I run into problems with the packages logspline and tseries. The error messages are as follows. Can anyone help? The compiler is gcc 2.95.3, if that helps. Andy ================================================ Installing source package `logspline' ... libs gcc

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

So after about 4 hours struggling with Try and TryCatch I am throwing in the towel. I have a more complicated function that used logspline through iterative distributions and at some point the logspline doesnt function correctly for some subsets but is fine with others so I need to be able to identify when the error occurs and stop curtailing the distribution and I think this Try or TryCatch

Hi interested readers, I have a function that creates several functions within a loop and I would like them to be returned for further use as follows: Main.Function(df,...){ # df is a multivariate data funcList<-list(NULL) for (i in 1:ncol(df)){ temp<-logspline(df[,i],...) # logspline density estimate funcList[[i]]<-function(x){expression(temp,x)} } return(funcList) } I have tried

I using oldlogspline (from logspline package) to model data distributions, and having a problem. My data are search area sizes. They are based on circular search radii from random points to the nearest edge of the nearest grass tussock. Search area sizes are distributed from 0 (the random point intercepts a tussock) and upwards (as points are further from any tussocks). The density of all my

Given the following, one of the things I am trying to see is what % of draws are below a certain number: lambda <- 3 rate <- 5 n <- 5 set.seed(123) v <- replicate(n, rexp(rpois(1,lambda), rate)) vv <- unlist(v) cat("% of draws below 0.1:", round(length(subset(vv, vv < 0.1))/length(vv)*100,0), "%\n") In actuality, my lambda, rate, and n are 26, 10, 1000000,

I have a model as follows: x <- replicate(100, sum(rlnorm(rpois(1,5), 0,1))) y <- quantile(x, 0.99)) How would one go about estimating the boundaries of a 95% confidence interval for y? Any pointers would be greatly appreciated. > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 5.1 year 2007 month 06 day 27 svn rev 42083 language R

Full_Name: Emmanuel Charpentier Version: 1.2.0 OS: Linux 2.2.18 (Debian 2.2) Submission from: (NULL) (193.251.31.31) When trying "mass" installation of available packages, the package logspline does not compile, at least when installed through install.packages("logspline"). It seems that it is a small bug, such a missing #include <math.h> or -libm linking switch ?

Hi All: Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources. Best wishes, ANDREW tseries: Package for time series analysis Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6 Depends: ts, mva, quadprog Date: 2001-08-27 Author: Compiled by Adrian

Hi everyone, Using R 2.10.1 on Windows Vista. DOWNLOADED PACKAGES DO NOT INSTALL. I expect to see them in C:\Program Files\R\R-2.10.1\library These files download (see below). And they are all in zip format. What am I doing wrong? Please help. All suggestions appreciated. trying URL 'http://cran.cnr.Berkeley.edu/bin/windows/contrib/2.10/RColorBrewer_1.0-2.zip' Content type

Dear List, I have read that a lognormal mixture model having a pdf of the form f(x)=w1*f1(x)+(1-w1)*f2(x) fits most data sets quite well, where f1 and f2 are lognormal distributions. Any pointers on how to create a function that would produce the 5 parameters of f(x) would be greatly appreciated. > version _ platform i386-pc-mingw32 arch i386 os

Hi folks, I know that density function will give a estimated density for a give dataset. Now from that I want to have a percentage estimation for a certain range. For examle: > y = density(c(-20,rep(0,98),20)) > plot(y, xlim=c(-4,4)) Now if I want to know the percentage of data lying in (-20,2). Basically it should be the area of the curve from -20 to 2. Anybody knows a simple

Friends, I changed one line of a package at the source level and then rebuilt it. When I run R CMD check, I get an error: installing R.css in C:/polsplineRS.Rcheck ---------- Making package polsplineRS ------------ adding build stamp to DESCRIPTION making DLL ... making hareall.d from hareall.c making heftall.d from heftall.c making lsdall.d from lsdall.c making lspecall.d from lspecall.c

Can anyone recommend a good nonparametric density approach for data bounded (say between 0 and 1)? For example, using the basic Gaussian density approach doesn't generate a very realistic shape (nor should it): > set.seed(1) > dat <- rbeta(100, 1, 2) > plot(density(dat)) (note the area outside of 0/1) The data I have may be bimodal or have other odd properties (e.g. point mass

The following are now on CRAN: KernSmooth: version 2.2 of the code for Wand & Jones book on kernel smoothing. logspline: spline fits to log denisites, with automatic choice of smoothing. ppr: projection pursuit regression. rpart: recursive partitioning (CART-like) VR: Venables & Ripley libraries 5.3pl021 for 0.62.1 and in the devel section tree: a clone

The following are now on CRAN: KernSmooth: version 2.2 of the code for Wand & Jones book on kernel smoothing. logspline: spline fits to log denisites, with automatic choice of smoothing. ppr: projection pursuit regression. rpart: recursive partitioning (CART-like) VR: Venables & Ripley libraries 5.3pl021 for 0.62.1 and in the devel section tree: a clone

Dear all, This is a simple probability problem. I want to know, How to generate a normal distribution with mean=1, min=0.2 and max=0.8? I know how the generate a normal distribution of mean = 1 and sd = 1 and with 500 data point. rnorm(n=500, m=1, sd=1) But, I am confusing with how to generate a normal distribution with expected min and max. I expect to hear your directions. Thanks in

I am aware of the nice R package "logspline", which does smooth density estimation from interval-censored data (that is, values that are known to lie in a specified interval rather than known exactly). Function logspline.fit uses a maximum penalized likelihood method, with the penalty related to the number of knots used in a cubic regression-spline fit. I need to be able to do some

Dear List I just downloaded and installed R 2.12.0 and then installed R Commander . First it got RCmdr and Car, and then suggested for other packages for utilizing the full functionality- I clicked yes! I got 140 packages installed!!! Cran Mirror was UCLA... Here is the list. Is this intentional- I can see some packages like snow and multicore which are desirable but quite optional.(see list