Displaying 20 results from an estimated 9000 matches similar to: "Anyone know a forum for stats advice?"
2011 Mar 10
2
Not sure how to handle hazard in my survival model
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
2011 Jan 18
1
Circular variables within a GLM, GLM-GEE or GAM
Hi,
I have a variable (current speed direction) which is circular (0=360 degrees), and I'd like my GLM to include the variable as a circular variable. Can I do this? And what is the code?
I'm actually doing a GLM-GEE using the 'geepack' package, so want to use it in that, but also interested in whether it can also be used in GLMs and GAMs (I use the 'mgcv' package for
2009 Oct 07
2
Second y-axis --- alternative to par(new=
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
2006 Sep 06
4
problem with putting objects in list
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,
2009 Oct 08
2
intersection of 2 density curves
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
2009 Oct 30
2
Ubuntu packages compiled with 2.9
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.
2009 Dec 11
1
Creating bibtex file of all installed packages?
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),
2010 Jan 21
1
superimpose histogram and fitted gamma pdf
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)
2011 Feb 24
2
debugging
Dear R user
How to make the program stop at the spot where the error occurs? I mean
inside the iterations,.
Many thanks
yan
**********************************************************************
This email and any files transmitted with it are confide...{{dropped:10}}
2006 Sep 20
4
Calculating mean together with split
Hi
I have a table called npl containing results of simulations.
It contains about 19000 entries and the structure looks like this:
NoPlants sim run year DensPlants
1 6 lng_cs99_renosterbos 1 4 0.00192
.
.
.
it has 43 different entries for sim and year goes from 1 to 100, and run
from 1 to 5.
I would like to calculate the mean of DensPlants for each
2011 Mar 04
1
Time series analysis for a daily series
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
2011 Jan 10
2
write.table equivalent for lists?
-----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
2006 Oct 25
3
simplification of code using stamp?
Hi
I have the following code which I would like to simplify. Id does linear
regressions and returns the r-squares, and the coefficients.
It runs slow, as it is doing the regressions for each - is it possible
to get the values in a dataframe which looks as follow:
expert | xx | seeds | r.squared | slope | intercept
Thanks in advance,
Rainer
library(reshape)
rsqs <- as.data.frame(
2012 Aug 22
3
Question concerning anova()
Hi
I am comparing four different linear mixed effect models, derived from updating the original one. To
compare these, I want to use anova(). I therefore do the following (not reproducible - just to
illustration purpose!):
dat <- loadSPECIES(SPECIES)
subs <- expression(dead==FALSE & recTreat==FALSE)
feff <- noBefore~pHarv*year # fixed effect in the model
reff <-
2011 Feb 01
2
How to change size of xlab inplot?
-----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
2008 Aug 01
1
bug in readRAST6 function in package spgrass6
Hi
if I try to import a raster layer which consists only of NULL values
from grass by using the readRAST6, I get an error message:
> readRAST6("HSericea_seedsDisperse_2007")
ERROR: Invalid value for null (integers only)
Error in readBinGrid(rtmpfl11, colname = vname[i], proj4string = p4,
integer = to_int) :
no such file:
2010 Jul 29
2
Replace last element in a vector - elegant solution?
-----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
2009 Sep 30
2
Compress (gzip) a pdf device
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
2011 Mar 11
2
How to get all combinations between two character vectors?
-----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
2006 Dec 01
2
specify point shape for ggplot (equivalent to pch)?
Hi
is it possible to specify the shape of the point to be used in ggplot
(as with pch in plot)? I couldn't find anything in the help.
Thanks
Rainer
--
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)
Department of Conservation Ecology and Entomology
University of Stellenbosch
Matieland 7602
South Africa
Tel: +27 - (0)72 808 2975 (w)
Fax: +27 - (0)86 516 2782
Fax: