similar to: Making a group membership matrix

I have a vector of land cover class data from a GIS: > landcov[1:10,2] #the first ten examples of a large (100k+) object 1 2 3 4 5 6 7 8 9 10 12 12 01 12 01 15 15 15 15 15 etc. I also have a lookup table for the class data that gives the cover type as a string: > cnd.names #the look up table, i.e., landcov[3,2] == 1 == "Montane Meadow" CndVal

Dear List, I am looking to calculate two things from my data frame and was after some advice. For the example below i want to know. 1. How many unique Orders/Families and Genera there are per eco-name 2. How many incidences are there for each Order/Family and Genus there are per eco-region I have 650 econame. I.e for Biak-Numfoor rain forests there are 2 orders, 2 families and two genera.

Dear R_help list members!   I’m studding forest carnivores and I have data on resting site selection and use!   I’m trying to model the probability that a forest carnivore might have in be located in a tree hollow (1) (the main resting site selected) rather than elsewhere (0) (other resting site types, dens, nests, etc.). The model should be related with several variables such as for instance:  

Hello, I have a data set that needs to be combined so that rows are summed by a group based on a certain variable. I'm pretty sure rowsum() or rowsums() can do this but it's difficult for me to figure out how it will work for my data based on the examples I've read. My data are structured like this: Plot SizeClass Stems 12 Class3 1 12 Class4

I'm trying to use tapply to output means and SD or SE for my data but seem to be limited by how many times I can subset it. Here's a snippet of my data > stems353[1:10,] Time DataSource Plot Elevation Aspect Slope Type Species SizeClass Stems 1 Modern Cameron 70F221 1730 ESE 20 Conifer ABCO Class1 3 2 Modern Cameron 70F221 1730

Hello all, This may be a simple question to answer, but I'm a little bit stumped with respect to the calculation of the F statistics in nested anovas with unbalanced design in R. In my case, I have 11 vegetation transects (with 1000 10cmx10cm squares), where we estimated shrub cover. We have two different treatments: wildfire (4 transects) and prescribed burning (7 transects) and we want to

I am a beginner with R and have managed to use vegan to run a CCA on my data. I am evaluating the habitat relationships of select species in riparian areas around creeks before and after logging across three sampling intervals (pre-logging, two-years post-logging,. and 10 years post-logging). I sampled from 18 sites of 3 treatments (n = 6 sites for each treatment). When I plot my

I am interested in estimating the variance components from different groups of variables in a linear regression. For example, I want to model dependent variable ?Y? based on six covariates as follows: (these are land use types measured at three scales for each Y, "Local", "Riparian", and "Catchment" - these are not nested as there is a unique area for each Y -

Hello useRs, I have some problems when I try to plot a lattice. The lattice has 200 columns (X value) and 5 rows (Y value), with a total of 1000 squares. I want to plot this lattice with the squares having different colours according to a variable value. The variable I want to plot is cover value of a shrub, and it has values between 0 and 1. I want to plot the squares with value 1,

CRAN (and crantastic) updates this week New packages ------------ Updated packages ---------------- New reviews ----------- This email provided as a service for the R community by http://crantastic.org. Like it? Hate it? Please let us know: cranatic at gmail.com.

I'm trying to update a package and would like to crush a WARNING message for a clean build. I've been struggling with this question and haven't gotten any traction on the web either. I've got a document file (Rd) that contains the following \usage statement: \name{sample.data} \alias{sample.data} \title{CONIFERS forest growth model sample data} \description{ A list object of

Hello list, I have installed 4 Etch, with dovecot 1.0.rc15-2, and I wanted that when the users make IMAP/POP3, the dovecot verifies the user/passwd in a Active Directory of Windows 2003. Greetings -- Usemos el Software Libre "Con todos y para el bien de todos" Michel Vega Fuenzalida Usuario Linux: 353763 Coordinador del?Grupo Linux Pinero T?cnico de Laboratorio Centro

---NOTICE--- This is a Notice that an iDetect 2002 user is currently doing research on you. This is not a real Detective. You can do a little investigating for yourself by getting iDetect 2002. You will also get A list of people who have done Investigations on you and your family. ---Read Below--- Now you can learn EVERYTHING about your friends, neighbors,

Hello useRs! I have % cover data for different plant species in 300 plots, and I use the ARCSINE transformation (to deal with % cover data). When I use a GLM I do not have any problem. But when I am trying to use a GAM model using mgcv package, to account for non-linearity I get an ?error message?. I use the following model: sp1.gam<-gam(asin(sqrt(0.01*SP1COVER))~

R developers, I've been happily building packages, under windows, for some time now and upgraded to R 2.1.0 and now when I attempt to build a package, I get the following errors... C:\conifers>rcmd build Rconifers * checking for file 'Rconifers/DESCRIPTION' ... OK * preparing 'Rconifers': * checking DESCRIPTION meta-information ... OK * cleaning src * checking whether

I've upgraded laptops and have put a new version of R on the new machine as well as the tools for building packages. I know R has required a "non-broken" version of tar for some time now and was hoping to not have to install any cygwin tools. Since building packages requires the use of MinGW and MSYS, why can't we simply use the tar that's part of msys? After looking at

I have long axis labels (8units) on the x axis and 2 unit number labels on the y axis. I need to move the x axis title without moving the y axis title. I have used mgp=c(title, label and line), but mgp moves both the x and y titles. The following puts the y axis title way way out to the left. Anyone know how to move just the x axis. X11() conif<-c("Control", "Switched

Hello. I am sorry for posting a (seemingly) simple question, but I have just spent 2 hours trying to find the answer, without success. I want to make a histogram with conditioning on a factor, using Trellis graphics. However, I do not want any colours (only black and white) either in the histograms or in the strip. There must be some simple argument but I can’t find it. Here is my code so

Is there an R function that does canonical correspondence analysis. Can it be done using the VR function corresp()? If not, how hard it be to write R code to do it? I am a population biologist with long but patchy programming experience in C, Smalltalk, Java and other languages. Thanks, Patrick Foley patfoley at csus.edu

I'm adding a couple of vignettes to an existing package. When I make a change to the sweave file, and run the check command, c:\conifers\trunk>R CMD check rconifers I get the following message(s) in the 00check.log file: * checking for unstated dependencies in vignettes ... OK * checking package vignettes in 'inst/doc' ... WARNING Package vignette(s) without corresponding PDF: