similar to: simplifying a GLM-removing categorical variables

Hello, I'm working with a series (30+) of hydrologic metrics and 10 vegetation communities and I need to determine which of the metrics provide the best separability for each of the vegetation communities. The hydrologic metrics are highly correlated, therefore the need to reduce the number of them considered is critically important. I've been looking at the scatter3d function in the

If I submit a form with lots of checkboxes using IE, rails, lighttpd (server on OSX Server, client on XP), then I get a "This page cannot be displayed" message from IE after some delay and nothing shows up in the rails or lighttpd logs. Cutting the number of checkboxes fixes the problem. Safari works fine. I''ve googled for anwers and apparently there is an upload problem

Hello list, is it necessary to properly answer() an incoming call ? I don't want to answer a call because the caller has to pay even if the attached SIP-phones do not answer the phone call. Because I answer() the incoming call, the caller has to pay for 60 seconds of 'ringtone'. On the other hand, sometimes an incoming call is send to a macro where the caller is given the

Dear list, I have data like this: dat1 <- data.frame(subject=rep(1:10,2), cond1=rep(c("A","B"),each=5), cond2=rep(c("C","D"),each=10), choice=sample(0:1,10,replace=TRUE)) I would like to compare subjects' "choice" for (cond1=="A" & cond2=="C") vs

Virt SIG folks, please join me in welcoming our latest SIG member, Marianne Lombard (username: jehane). She's been a Fedora contributor for quite sometime already, and will now be helping to ensure I'm not being a bottleneck RE: docker on CentOS virt :) Welcome aboard, Marianne. -- Lokesh Freenode: lsm5 GPG: 0xC7C3A0DD -------------- next part -------------- A non-text attachment was

I have two lists of the same shape, like this: x <- list() x[[1]] <- c("one","two") x[[2]] <- c("three","four","five") y <- list() y[[1]] <- c(TRUE,FALSE) y[[2]] <- c(FALSE,TRUE,TRUE) I would like to index x "by" y, that is, the result in this case should be: z [[1]] [1] "one" [[2]] [1] "four"

I'm not sure if this is a bug, or if I'm doing something wrong. =20 =46rom the worms dataframe, which is at in a file called worms.txt at =20 http://www.imperial.ac.uk/bio/research/crawley/therbook <http://www.imperial.ac.uk/bio/research/mjcraw/therbook/index.htm>=20 =20 the idea is to extract a subset of the rows, sorted in declining order of worm density, with only the maximum

Hi list, I have data in a dataframe t1, with a column for different amounts spent, a column what it was spent on, and a column with dates, from which I create a new column with months. Example: amount <- rep(c(10,20,30),3) what <- rep(c("food","books","cycling"),3) when <- c(rep("2008-09-05",5),rep("2008-10-07",4)) t1 <-

Dear list, I've written a function that plots subjects. Something like: myplot <- function(subject) { plot(subject) } Subjects are vectors, e.g. ... s1 <- c(200,200,190,180) ... and plotting them works fine, e.g. ... myplot(s1) Now I want to have "s1" etc appear in the plot title, but I don't know how to refer to this generically (the object "name"? I tried

Dear list, I have a file that is comma delimited but contains some erroneous non-delimiter commas. I would like to replace these commas with semicolons and then read the correct file into R as a data frame. I want to do this from within R, without changing the original data file. My current idea of how to do this would be to use system("sed ...") and feed the result to read.csv(), but

Dear list, I have two vectors: x <- c("one","two") y <- paste(rep(x,2),"blah") I want to replace all occurrences of each element of x in y with something else, so that y looks like this: y [1] "something else blah" "something else blah" "something else blah" [4] "something else blah" I can do this using a loop: for (

Dear Rers, I have a function to barplot() a matrix, eg myfun <- function(x, ...) { barplot(x , ... )} (The real function is more complicated, it does things to the matrix first.) So I can do: m1 <- matrix(1:20,4) myfun(m1) myfun(m1, main="My title") I'd like to be able to add the number of rows of the matrix passed to the function to the "..." argument, eg

First I would like to say thank you for taking the time to read it.Here is my problem. I am running a lmer analysis for binary longitudinal (repeated measures) data. Basically, I manipulated fruits and vegetation to two levels each(present and absent) and I am trying to access how these factors affect mice foraging behavior. The design consist of 12 plots, divided in 3 blocks. So each block

Dear R users, I create a graphic with two plots side by side using layout(), like this: layout(matrix(c(1,2),1)) plot(1:10,main="left plot") lines(c(3:7,7:3),col="red") plot(10:1,main="right plot") The lines() obivously get added to the "left plot" plot. Now, I'm trying to write a function that builds up a plot bit by bit to then include it in a

Hi R listers, I am trying to group distances of nests to the vegetation into classes that are define by (0-5m, 6-10m, 11-15m, 16-20m, 21-25m, 26-30m, 31-35m, 36-40m, 41-45m, 46-50m, 51-55m, 56-60m). Each row is a nest and all the distances to the vegetation is in a column. In plyr, I have tried - below script but I think I am going about this the wrong way and am not successful. #Veg index

Would anybody be able to advise on which package would offer the best approach for producing a model able to predict the probability of species occupation based upon a range of variables, some of them catagorical (eg. ten soil types where the numbers assigned are not related to any qualitative/quantitative continuum or vegetation type) and others continuous such as field size or vegetation height.

Dear list, I would like to start some R workshops at King's College London, and to do so, I would like to use the "Use R!" logo at http://www.agrocampus-ouest.fr/math/useR-2009//useR%21%202008_fichiers/useR-middle.png Since it seems to be difficult to get a shell account at KCL, I also went ahead and registered use-r.org.uk and am starting to put together a website at

Hi, I have a qplot like the one in the minimal example below, except I also have faceting like this: qplot(jitter(Goodall),jitter(Better.adapt),colour=Second.adapt,facets=~Pol,data=d1) and with the real data I get quite a lot of overplotting, so I would like to add an alpha channel. In addition, I would like to be able to control which colors are used for each value of Second.adapt (which

Hello all, I am going slightly mad trying to create a table for running co-correspondence analysis. What I have is seed bank and vegetation data, and my aim is to see if the vegetation found in a site (containing several seed bank samples) can predict the composition of a seed bank sample within that site. So for this I need two tables with matching rows. I have created an empty matrix,

Hi, This is a question regarding technique rather than an R specific issue. I have been asked to evaluate a 30+ year long term continuous survey of bird presence/absence data that has an associated ocular estimate of the vegetation community percent coverage. The data are organized by subpopulations (5), and by year ( 1991 - present). We are interested in gaining understanding on whether bird