similar to: debugging

Hi R experts :) I'm trying to carry out a survival model on my data, but I am unsure of whether it's appropriate or if I should do something specific in regards to hazard. My data is time to death by predator where I have 8 prey and one predator in the setting. This means that two prey can't possibly die at the same time and I can't quite get my head around how to include this in

Hi is there an easy and fast way, to generate a BibTeX file of all installed / loaded packages and R? I know about toBibtex(citation()) to extract the BibTeX for a single package, but how can I generate a file containg citations for all installed / loaded packages? Cheers, Rainer -- NEW GERMAN FAX NUMBER!!! Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT),

Hi is there an alternative to par(new), for ading data to a plot for a different y-axis? My problem with par(new=TRUE) is, that it re-defines all axis and labels (as in example 1) and one has to use xlim=... to fix the x-axis. I am looking for something, which simply resets the y-axis, so that a new plot() (or points()/lines()) keeps the x-axis, but re-defines the y-axis. Is there something

Hi, I would like to find out the coordinates of the intersection points of 2 density curves. I did a search but i didn't get any significant results. I really hope some of you have some ideas. here it is an example: set.seed(123) x1 <- rnorm(100, 1, 1) x2 <- rnorm(100, 0, 1) d1 <- density(x1) d2 <- density(x2) plot(d1, col = 2) lines(d2, col = 4) Now i would like to know

Hi r-users,   I try to draw histogram (in terms of probabilty) and superimpose with the gamma pdf. Using this code below I can get the plots BUT the y scale for the density is out of scale.  How do I change the y-axis scale to max 1?  Another thing, how I do draw smooth line rather that points?   Nota that my observed data is hume_pos and the fitted data is rgam1.   hist(hume_pos,prob=TRUE)

Hi I upgraded to R 2.10.0, and I encounter a problem with some packages on R, installed via apt-get.. They are still compiled by 2.9 and give me an warning when I load the packages. Will these packages be recompiled with R 2.10, or should I rather switch to self-compiling? Cheers, Rainer -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys.

Hi, I know this forum is for R-related issues, but the question I have is a statistical question & I was wondering if anyone could recommend a good statistics forum where I can ask the question? My question is relating to bootstrapping of binary data (ecology data) - I can give more detail, but wasn't sure I could address the question here as it is more statistical based than R based

Hi everyone I am trying to do some time series analysis with daily temperature data (40 years). I have created a zoo object and ts object but can't apply stl function. It says the series is not periodic or has less than two periods. I've searched through google and found a lot of messages about this problem but not a solution/example to look for trend and seasonal component of a

Here are what i am going to accomplish: I have 400 files named as xxx.txt. the content of the file looks like the following: name count 1. aaa 100 2. bbb 2000 3. ccc 300 4. ddd 3000 ........ more that 1000 rows in each files. these are the areas i need help: 1. how can i only read in the files with the string patterns ggg or fff as part of the file names? for

Can I save R graphs as a R objects ? IF yes then if I click this R object can I edit my chart as in Excel. Please suggest your views -- View this message in context: http://n4.nabble.com/Can-I-save-R-graphs-as-a-R-objects-tp1014621p1014621.html Sent from the R help mailing list archive at Nabble.com.

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I am writing simulations in R, and quite regularly, I have to save lists and objects to HDD and load it later again. So I am wondering: why is there no function to write lists (and S3, S4 objects) onto HDD WITHOUT keeping the name? What I mean is: For data.frames I can use x <- data.frame(x = runif(10)) write.table(x, "x.txt) rm(x) y

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I want to increase the size of the text for xlab in the plot command: plot(1,1, xlab="Label") I tried cex, cex.caption and others, but none worked. Is there a way of using cex=2 for the xlab? Thanks, Rainer - -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centre of

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I have to replace the last element of a vector with 0, but the length of x is not known at coding time. So I do the following: x <- 1:10 x x[length(x)] <- 0 x it is working nicely, but I am wondering: is there a more elegant solution( like tail(x, 1) <- 0 , which does not work)? Thanks, Rainer - -- Rainer M. Krug, PhD (Conservation

I have not found an easy way to compress a file on filesystem. Especially I'd like to compress a pdf from pdf() function/device. Is it possible to compress It on the flight? I'd like to do something like: pdf(gzipconnection()) dev.off() If It is not possible, how can I create a gzip with the pdf? Thanks Daniele A. ________________________________ ORS Srl Via Agostino Morando 1/3

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I know there is a function - I have used it before - but I always forget what it is called... I need the combination of two character vectors, i.e: x <- c("a", "b") y <- c("x", "y") z <- THEFUNCTION(x, y) z == c("ax", "ay", "bx", "by") I promise I will

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Assuming the following: > x <- data.frame(a=1:10, b=runif(10)) > str(x) 'data.frame': 10 obs. of 2 variables: $ a: int 1 2 3 4 5 6 7 8 9 10 $ b: num 0.692 0.325 0.634 0.16 0.873 ... > write.csv(x, "x.csv") > x2 <- read.csv("x.csv") > str(x2) 'data.frame': 10 obs. of 3 variables: $ X:

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi rgl crashes my R session, when resizing the rgl graphic window. I am using Ubuntu Maversick, with dual monitor setup. If I disconnect one monitor, I can resize it "a little bit", but it still craches if I enlarge it to much. I assume that the problem has to do with allocated graphic memory in the kernel, but why is R crashing

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi rgl crashes my R session, when resizing the rgl graphic window. I am using Ubuntu Maversick, with dual monitor setup. If I disconnect one monitor, I can resize it "a little bit", but it still craches if I enlarge it to much. I assume that the problem has to do with allocated graphic memory in the kernel, but why is R crashing

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I know this has been asked before, but I can not find it - so my apologies. I want to compile R myself, to optimize it for speed. There is a small section (a few lines actually) in the R-admin manual (Compilation-flags) which states that the LDFLAGS "on recent systems ?'-Bdirect,--hash-style=both,-Wl,-O1'? is recommended". But

Hi I use the following code and it stores the results of density() in the list dr: dens <- function(run) { density( positions$X[positions$run==run], bw=3, cut=-2 ) } dr <- lapply(1:5, dens) but the results are stored in dr[[i]] and not dr[i], i.e. plot(dr[[1]]) works, but plot([1]) doesn't. Is there any way that I can store them in dr[i]? Thanks a lot, Rainer -- Rainer M. Krug,