similar to: Convert from ppp to spp objects

Good day to the R community, I am interested to run the plot.count() function in the "untb" package. My script is as follows:- > library(untb) > Community1 <- > structure(c(371,167,119,78,74,53,50,31,28,25,20,19,19,17,13,12,12,10, > 9,9,8,8,7,7,7,7,6,6,6,6,5,5,5,5,4,4,4,3,3,3,2,2,2,2,2,2,2,1,1,1,1,1, > 1,1,1,1,1,1,1,1), .Dim = 60, .Dimnames =

##I am trying to test for fixed factor main effects in an unbalanced mixed effects model but when I fit the reduced model for "mic" factor effects, the extra degrees of freedom are being allocated to a nested term rather than the residuals. The model has inc, mic and spp are independent variables and vial nested within spp. inc and spp are already coded as factors since they were

Hello, I'm trying to convert a character column in several dataframes to lower case. ### # # Sample data and 'spp' column summaries: # dput(ban.ovs.1993[sample(row.names(ban.ovs.1993), 20), 1:4]) ban.ovs.93 <- structure(list(oplt = c(43L, 43L, 38L, 26L, 35L, 8L, 39L, 1L, 34L, 50L, 10L, 29L, 31L, 24L, 18L, 12L, 27L, 49L, 28L, 51L), rplt = c(NA_integer_, NA_integer_,

Hello, I am unable to figure out how to code a new column in a data frame based on an existing column that matches a column in a reference data frame, in a relational-db fashion. I would like this to maintain a minimum set of reference tables that may be reused over several similar datasets. Specifically, I have two data frames as listed below, 'Bos' and 'tree.' For each case in

Hi, I would need some help with the splans package in R. I am using a Shapefile (downloadable at) http://rapidshare.com/files/215206891/Redlands_Crime.zip and the following execution code setwd("C:\\Documents and Settings\\Dejan\\Desktop\\GIS\\assignment6\\DataSet_Redlands_Crime\\Redlands_Crime") library(foreign) library(splancs) auto_xy<-read.dbf("Auto_theft_98.dbf")

Dear all, I have a big data frame: str(data1) 'data.frame': 18272 obs. of 11 variables: $ tag : int 100001 100002 100003 100005 100007 100008 100009 100011 100012 100014 ... $ sp : Factor w/ 18 levels "acassp","acocar",..: 13 5 7 14 14 18 3 11 13 10 ... $ gx : num 20 10 35 68 88 63 123 115 137 136 ... $ gy : num 30 25 24 1 10 40 45 25 23 45 ...

Dear R-help folks, Can someone guide me to a source where I can learn more about the / operator in model formulae? I found a reference to it in Venables and Ripley's MASS, p. 142, where it says, in reference to ANCOVA: "Terms of the form a/x, where a is a factor, are best thought of as 'separate regression models of type 1 + x within the levels of a.'..." This seems

Dear List, Im new to R - making a transition from SAS. I have a space delimited file with the following structure. Each line in the datafile is identified by the first letter. A = Inventory (Inventory) X = Stratum (Stratum_no Total Ye=year established) P = Plot (Plot_no age slope= species) T = Tree (tree_no frequency) L = Leader (leader diameter height) F = Feature (start_height finish_height

Hello: I am using R mainly on windows XP, version 2.5. I?m a biologist, with a medium level statistics background. I have a problem stating a two-way factorial design where both factors are random. I?m using the lmer() function implemented in the Matrix package version 0.99. My design is as follows: Two species were randomly selected from a total of 4 species. This species are present

[+CC llvm-dev this time] Hi, As discussed in the review thread in http://reviews.llvm.org/D16439, the future plan around statepoints, deopt bundles, PlaceSafepoints etc. is to "constant fold" -spp-no-statepoints and -rs4gc-use-deopt-bundles to true. We (Azul) have moved to a representation of safepoint polls, deopt state etc. that enables us to do the above; and at this point I'm

Hi all! I've written a handy script that uses a for loop to allow me to generate a large number of figures and statistical outputs for a large dataset. I am using indexing to retrieve a species name for the title of my graphs- which worked fine. However, I need to italicize these species names. I originally used the paste function, and had no problems with indexing: *main=paste("Yield

I am trying to rasterize a SpatialPolygon object by a RasterBrick object. The documentation of the raster::rasterize function explicitly says this is allowed. Here's what I am doing # load the raster package library("raster") # create a raster brick object using the example from the brick function documentation b <- brick(system.file("external/rlogo.grd",

Dear all of you, I would like to get some help from you. Here I attach you a model, that I would like to be stochastic so I would need each time the value of x changed, the values of parameters (ss, emrg, gf, spp, sr) did too according to a normal distribution, with mean its value and standard deviation given by the parameter name preceded by SDV (In the case of parameter ss, mean is 0.3 and

Dear all: My logistic regression model includes one qualitative and one quantitative predictor variable, aes <- glm(p.a ~ spp * log(light), family=binomial(link=logit)), where spp is abundance of 3 species and light is subcanopy light availability varying from 0 1. I want to test differences among levels of the quantitative variable at a value of x other than the current log(light)=0.

I'm stuck on a seemingly simple problem. I'm trying to subset the data by several numbers and it cuts out half of the rows. Here is the sample code: test <- as.matrix(c(1,1,1,1,3,3,7,7,7,7)) Count <- tapply(test[,1], test[,1], length) # count for each value spp <- unique(test[,1]) Count1 <- as.data.frame(cbind(Count,spp)) Max <- max(Count) Count1$sppMax <-

I am trying to create a heap of boxplots, by looping though a series of factors and variables in a large data.frame suing paste to constrcut the facto and response names from the colnames I thought I could do this using get() however it is not working what am I doing wrong? thanks Nevil Amos > sp.codes=levels(data.all$CODE_LETTERS) > > for(spp in sp.codes) { > > >

Dear All I would really appreciate some help with a script which a colleague wrote for me (attached), but I am having problems running (and have not been able to contact my colleague). The script is designed to compare the area of suitable habitat in binary projections of a large number of species current and future distributions, and create an excel file detailing the total area suitable in the

Hi Listers, I am trying to reflect a PCA biplot in the x-axis (i.e. PC1) but am not having much success. In theory I believe all I need to do is multiply the site and species scores for the PC1 by -1, which would effectively flip the biplot. I am creating a blank plot using the plot command and accessing the results from a call to rda. I then use the calls to scores to obtain separate site and

Hello I have downloaded the spp_1.10.tar.gz and it contains compiled code. I installed the Rtools &lt;http://www.murdoch-sutherland.com/Rtools&gt; and I tried to install the package (spp_1.10.tar.gz) but it not worked bellow it is the error message. Could somebody please help me with this. Thank you * installing *source* package 'spp' ...

Hello All, I am trying to open all files within a folder and create multiple histograms from each file, *and* have it save with the original file name plus some new information. The way I have it set up right now, I keep saving over each new graph. I can turn the history on and see them all, but I want them all to save as unique files as well. Idealy they would be "filename CHN 1mm Length