similar to: evaluating density objects

Hi there! Both on windows and Solaris platforms I run into troubles when trying to draw two plots on top of each other. Here a short code example: > x1 <- c(2,3,4,5,6,7,8) > y1 <- sin(x1) > x2 <- c(1,2,3,4,5,6,7,8,9,10) > y2 <- cos(x2) > d1 <- cbind(x1,y1) > d2 <- cbind(x2,y2) > plot(d1) > par(new=T, xaxs="d") > plot(d2) Error in plot.new()

Hi all. I would like to highlight overlapping regions of two densities and I could not find a way to do it. Here is the sample code: myd <- c(2,4,5, 4,3,2,2,3,3,3,2,3,3,4,2,4,3,3,3,2,2.5, 2, 3,3, 2.3, 3, 3, 2, 3) myd1 <- myd-2 plot(range(density(myd)$x, density(myd1)$x), range(density(myd)$y, density(myd1)$y), type = "n") lines(density(myd), col=1, lwd=4)

Are there any functions in R for fitting distributions? In particular, I would like to fit the weibull and extreme-value distributions when: i) the data are given and the maximum likelihood estimates are required; ii) the sample moments are available but the complete data record is not (so moment estimates are required). Any suggestions will be gratefully received.

Hi all, this is OT, but maybe someone can give me a clue. I've got data from eye tracker experiments (750 data points). These figures show how the data lock like |** * ***** |*** ****** | ** **** | * * | | * * | | * |-------------- |--------------- Y axes display velocity, x axes display time, * are the data points.

Hi all, I'm very new to R. I have installed R 1.9.0 on Linux (Fedora). Now I got an self-made package comprising R functions as well as C-Code which are used in several R functions. I installed the package without any error (see install log below). Then, I checked in /usr/lib/R/library if the package izbi exists and it exists. But whenever I try to load the library on the command line I

sorry for cross-posting Dear all, I have tow (several) bivariate distributions with a known mean and variance-covariance structure (hence a known density function) that I would like to compare in order to get an intersect that tells me something about "how different" these distributions are (as t-statistics for univariate distributions). In order to visualize what I mean hear a little

Hi, I apologize for this often/old question. I found some hints but couldn't solve the problem so far. I have C functions (incl. the header files) as well as the R wrapper functions which I want to use for faster calculations. These functions are included in a R package. The installation process seems to be ok (no errors). I also can load the package without errors. But when I call the

Dear All, I wonder whether anyone has an experience with FAME package written by Jeff Hallman. All my attempts to send him the following problem report did not succeed (the mail system says that my e-mail could not be delivered), so I turn for help to this list. I tried to use your FAME package written for R, but somehow I cannot get it working. I am using Windows XP and the newest R

Hello all, I'm a Wine newbie who tried to get Quicken 2009 (german) installed on a openSUSE 11.1 box. This is what I did according to the guy who succeeded in getting it up and running: 1) Start with a fresh Wine Profile 2) Follow the Gecko trick (not sure if needed was done for another app) http://wiki.winehq.org/Gecko 3) From winetricks install corefonts and vcrun6 4) From winetricks

Dear friends, Do you know how to calculate the CONCAVE hull of a set of points (2- dimensional or n-dimensional)? is that possible in R? (With a "smoothing" parameter of course). Best, -- Corrado Topi Global Climate Change & Biodiversity Indicators Area 18,Department of Biology University of York, York, YO10 5YW, UK Phone: + 44 (0) 1904 328645, E-mail: ct529 at york.ac.uk

Dear all, I have an array of strings a <- c("2008q3","2005q1","2004q3") I would like to add to each a[i], with i = 1,2,3 the following string "IMT", such that in the end I could get b <- c("IMT2008q3","IMT2005q1","IMT2004q3") Is it possible to accomplish this without a loop command? Thank you. Best,

Hi, When R generates density plots and these are exported to postscript( a=matrix(1:100,10,10);image(a,col=rainbow(100);dev.copy2eps(file="image.eps") ) The result often looks bad when rendered on screen. The help page states that this is because programs use anti-aliasing. That seems to be true - turning off anti-aliasing for gs (-dGraphicsAlphaBits=1) of in OSX's preview makes

Dear developers, the current implementation of hist.default() calculates 'density' (and 'intensities') as dens <- counts/(n*h) where h has been calculated before as h <- diff(fuzzybreaks) which results in 'fuzzy' values for the density, see e.g. > tmp <- hist(1:10,breaks=c(-2.5,2.5,7.5,12.5),plot=FALSE) > print(tmp$density,digits=15) [1]

For the hist estimate >par(mex=1.3) >dens<-density(q) >options(scipen=4) > ylim<-range(dens$y) > h<-hist(q,breaks="scott",freq=FALSE,probability=TRUE, +? right=FALSE,xlim=c(9000,16000),ylim=ylim,main="Histogram of q(scott)") > lines(dens) >box() ? For the kernel estimate>options(scipen=4) > d <- density(q, bw =

Hello, I am try to make a density plot where plots are stacked like the one found here: http://dsarkar.fhcrc.org/lattice/book/images/Figure_14_03_stdBW.png I am facing problems, however. Using the code example below, I'd like to generate a separate panel for each val of id2. Within each panel, I'd like to have individual histograms each on separate lines based on the value of id1. ?Note

Hi all, I'm looking for an efficient solution (speed and memory) for the following problem: Given - a data.frame x containing numbers of type double with nrow(x)>ncol(x) and unique row lables and - a character vector y containing a sorted order labels Now, I'd like to sort the rows of the data.frame x w.r.t. the order of labels in y. example: x <- data.frame(c(1:4),c(5:8))

Hi, I have a couple of observations, a question or two, and perhaps a suggestion related to the plotting of density on the y-axis within the hist() function when freq=FALSE. I was using the function and trying to develop an intuitive understanding of what the density is telling me. After reading through this fairly helpful post:

Is there a way to do text plots in R? I'd like to do some simple XY plots in R with the output in text (ascii). For example, with gnuplot I can do the following: echo 'set terminal dumb ; plot sin(x)' | gnuplot To generate a simple sin wave. Since I connect to a remote Linux machine using SSH, being able to generate a rough idea of what a plot will look like in text would be

Hi, I am wondering if there is a way to call Gnuplot from R and/or if anyone can recommend a package on CRAN capable of doing this? Thanks, Chris PS - Please cc me on the response. [[alternative HTML version deleted]]

Hello, I'm new to this and am trying to teach myself some R by plotting biological data. The growth curve in question is supposed to be fitted to the Verhulst equation, which may be transcribed as follows: f(x)=a/(1+((a-0.008)/0.008)*exp(-(b*x))) - for a known population density (0.008) at t(0). I am trying to rework the example from "An Introduction to R" (p. 72) for my case and