similar to: logit regression, test among groups

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 =

Windows 2000 (v.5.00.2195), R 1.5.1 I have an lme object, fm0, which I examine with anova() and intervals(). The anova output indicates one of the interaction terms is significant, but the intervals output shows that the single parameter for that term includes 0.0 in its 95% CI. I believe that the anova is a conditional (sequential) test; is intervals marginal or approximate? Which should I

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

Windows 2000 . 5.00.2195 with Service Pack 1. R 1.5.1 Output from my split-split plot aov "alerted" me that I have done something wrong. I designed an experiment with all combinations of all levels of each treatment, but lost a little data (3 out of 192 plots). With the following data, I run the following model: > collim[c(1:6,187:192),c(1,3:6,9)] plot Litter Fert

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

How do I control the text in strips? Specifically, I want to put in the ranges generated in shingle(x) where x is continuous. Thanks, Hank Martin Henry H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.muohio.edu/~botcwis/bot/henry.html

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

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

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

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

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

Hi All! As a new useR, I am currently working in R 2.5.1 to produce Ripley's L (the linear version of Ripley's K) for a spatial point dataset (originally a ESRI shapefile). I have coerced the data to a ppp object (because I was originally using 'kest' in the spstat package to produced Ripley's K) but I need to convert my spatial data to a spp object to be able to use

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

After checking the usual archives (Baron search engine) I came up nothing. It appears that ylab=expression() doesn't work in lattice? How do I put superscripts in ylab? Many Thanks, Hank Martin Henry H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243

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