similar to: define subset argument for function lm as variable?

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 <-

Hi I know how to referenc R in a scientific paper - but is there a standardised way to reference packages? Thanks Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Plant Conservation Unit Department of Botany University of Cape Town Rondebosch 7701

Hi I am sure this is simple - but how can I convert one charecter vector into one string? example: x <- c("This ", "is ", "one ", "sentence.") should become "This is one entence" Thanks Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys.

Hi Dirk, Thanks for the explamnayion - Debian is running in a VM (and nothing really installed on it) , and I think the easiest is for me just to install Debian 10 and to use that one. Thanks a lot, Rainer > On 22 Sep 2020, at 15:14, Dirk Eddelbuettel <edd at debian.org> wrote: > > > On 22 September 2020 at 14:49, Rainer M Krug wrote: > | I know this is likely

Hi > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 2 minor 5.1 year 2007 month 06 day 27 svn rev 42083 language R version.string R version 2.5.1 (2007-06-27) I want to create a pdf withe three graphs on a page and with two pages:

Hi this might be a little bit off topic, but here it goes: lets assume I have the following: set.seed(13) dat1 <- rnorm(2000, mean=10, sd=10) dat2 <- rnorm(100, mean=10, sd=20) d.all <- density(dat, n=1024) d.co <- density(x[[v]], , from=min(d.all$x), to=max(d.all$x), n=1024) d.diff <- list( x = d.all$x,

Hi I am trying to install the package "pheno", but it needs the package "nprq" by Roger Koenker et al. which I can I find this package? It does not seem to be on CRAN and googling also doesn't give me an URL - is it still somewhere available? Thanks, Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF

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 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 am testing a package, and after I make changes, I have to close R and open R again to load the new version (same version number) of the package I am working on. So my question: is there a function which removes a package, i.e library(myPackage) > Package is loaded unlibrary(myPackage) > package is not loaded any more Thanks, Rainer -- Rainer M. Krug, PhD (Conservation Ecology,

Hi I have a plist t which contains size measurements of individual plants, identified by the field "plate". It contains, among other, a field "year" indicating the year in which the individual was measured and the "height". The number of measurements range from 1 to 4 measurements in different years. My problem is that I would need the LAST measurement. I only

Hi I know this is likely documented somewhere, but I can?t find it. How can I install R > 4.0.0 on Debian 9 Stretch? Thanks a lot, Rainer -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Orcid ID: 0000-0002-7490-0066 Department of Evolutionary Biology and Environmental Studies University of Z?rich Office Y34-J-74

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),

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Just to be sure: I couldn't find a way of creating an interactive 3d graph which can be embedded in a pdf - is it correct that this can not be done with R? Cheers, Rainer - -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centre of Excellence for Invasion Biology Stellenbosch

-----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

-----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

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

-----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 have 10.000 simulations for a sensitivity analysis. I have done a few sensitivity analysis for different response variables already, but now, as most of the simulations (if not all) show some cyclic behaviour, see how the independent input parameter influence the frequency of the cyclic changes and "how cyclic" they actually are. So