similar to: Plotting from a list

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

Hello, This is a re-submittal of question I submitted last week, but haven't rec'd any responses. I need to extract the probabilities used to construct the barplots displayed as part of the graph produced by plot("ctree"). For example, library(party) iris.ct <- ctree(Species ~ . , data = iris) plot(iris.ct) Instead of a simple example with only 4 terminal nodes, my

Hello I have a Linux machine (Ubuntu 8.04 hardy, Gcc version 4.2.4 (i486-linux-gnu) currently running R 2.6.2. I'd like to upgrade to 2.10. First Question): What is the appropriate way to remove the old version of R? Part 2. After downloading r-base_2.10.0.orig.tar.gz and opening the archive. I ran the ./configure routine. It failed claiming that it could not find the F77 compiler. My

Hello I've been asked to help evaluate a vegetation data set, specifically to examine it for community similarity. The initial problem I see is that the data is ordinal. At best this only captures a relative ranking of abundance and ordinal ranks are assigned after data collection. I've been trying to find a procedure in R that can handle ordinal based classification and so far have

Hello, I've been using ctree and have developed a 55 node - 28 terminal solution. As can be imagined, the plot is difficult to travel down each of the major branches. I've read the help files for ctree I saw where terminal nodes can be color coded. plot(airct, type = "simple") > plot(airct, terminal_panel = node_boxplot(airct, col = "blue", + fill = hsv(2/3,

I am using windows XP with R 2.8.1 I am generating a lattice plot of annual rain patterns using the following function: > xyplot(rain.stats$min+ rain.stats$max + rain.stats$ave ~ rain.stats$month |rain.stats$year, lty = 1, data = rain.stats, type = c("l","l", "l"), col = c("red", "blue", "green"), distribute.type =

I know this has been addressed before, but I'm still confused by the solution. When running TINN - R (ver 2.1.1.6) with R (2.8.0) I have lost the functionality of TINN-R It returns an error source(.trPaths[5], echo=TRUE, max.deparse.length=150) I modified the RProfile.site file as follows: options(pager="interal") .trPaths=c( 'C:/Documents and

Hello, I've been getting the following error when using the mlogit function from the mlogit package This is one of the examples provided in the Package "mlogit" January 27, 2010 description data("Fishing", package="mlogit") Fish <- mlogit.data(Fishing, varying = c(4:11), shape="wide", choice="mode") summary(mlogit(mode ~ pr + ca - 1,

Hello, I've tried to install randomForest on a Ubuntu 8.04 Hardy Heron system. I've repeatedly rec'd the error: > install.packages("randomForest", dependencies = TRUE) ERROR: compiliation failed for package 'randomForest' ** Removing '/home/admuser/R/i486-pc-linux-gnu-library/2.6/randomForest' The downloaded packages are in

I'd like to annotate an xyplot with an R^2 value, but can not find the syntax to define the superscript format for the text. I'd appreciate suggestions, xyplot(SharkSloughEggs.df$Sharkeggs.rel + SharkSloughEggs.df$SharkHatched.rel + SharkSloughEggs.df$SharkFlooded.rel + HSI.shark$MEAN ~ Year, pch=c(1,9,5,4), lty=c(1,6,3,4), lwd = c(2), col = c("black",

Hi, I need to create a plot (type = "h") and then overlay a skewed-normal curve on this distribution, but I'm not finding a procedure to accomplish this. I want to use the plot function here in order to control the bin distributions. I have explored the sn library and found the dsn function. dsn uses known location, scaling and shape parameters associated with a given input

Hello Trying to develop variograms and kriged surfaces from a point file. Here is what I've done so far. library(gstat) # also loads library(sp) library(lattice) soilpts$x <- soilpts$UTM_X soilpts$y <- soilpts$UTM_Y soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN, TC, Total_Mg)) dim(soil.dat) [1] 1292 7 coordinates(soil.dat) <- ~ x+y

I'm stumped, can anyone find my error in this sequence. for(j in 1:(varsize[4]-1)) temp <- get.var.ncdf(nc=input, varid="p_foraging",c(1,1,j),c(varsize[1],varsize[2],1)) filled.contour(x, y, temp, color = terrain.colors, plot.title = title(main = paste("Everglades Wood Stork Foraging Potential \nYear", (2000+j)), xlab =

Hello I'm plotting a large suite of barcharts and need to modify the size of the text for both the yaxis and xaxis labels. I've tried using the following: > trellis.par.set(list(par.ylab.text = list(cex = 0.65)), trellis.par.set(list = par.xlab.text = list(cex = 0.65)))) On inspection, however after I invoke this line, > trellis.par.get("par.ylab.text") >

Hello, Last week I installed ver 2.9.1 and it worked fine. This morning I have been working on some simple schemes, plot(PropHatchedNests$Phatched, PropHatchedNests$MEAN) abline(lm(PropHatchedNests$Phatched, PropHatchedNests$MEAN)) This oddly produces a box and whiskers plot. I uninstalled 2.9.1 and installed 2.9.2 and ran the same commands. Same output. Is anyone else experiencing this

A colleague is interested in modifying the ylim definition for individual panels of a common bwplot plotting statement. Is there an approach to modifying the bwplot function to allow for a dynamic ylim range given different panel factors ? He is using R 2.6.2 on a Linux distribution running from a CD. Thanks for the assistance. Steve Steve Friedman Ph. D. Spatial Statistical Analyst

I have just upgraded from 2.9.2 to 2.10.1 on my XP machine. I rec'd the following error message: Error in strsplit(x[ok], "[.-]") : 5 arguments passed to .Internal(strsplit) which requires 6 I also tried update.packages(checkBuilt=TRUE, ask = FALSE) > update.packages(checkBuilt=TRUE, ask=FALSE) Error: could not find function "update.packages" Has anyone else

Hello, vegMeans <- by(SoilVegHydro[3:37] , SoilVegHydro$Physiogomy, mean) vegSD <- by(SoilVegHydro[3:37] , SoilVegHydro$Physiogomy, sd) write.table(vegMeans, file="A:\\Work_Area\\Steve\\Hydrology_Data\\data\\vegMeans.txt") Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class '"by"' into a

My data contains a variable "observation_date" and it contains values as: 1985-09-02 1985-09-15 1985-07-31 1985-09-02 I need to process data annually rather than daily, therefore I'm trying to 1) either extract the first 4 digits from this field and use them as a new variable "year" or 2) keep the variable as it is and process the analysis using the first 4 digits of

Hello I'm encountering some difficulty correctly reading binary files. The binary files store data as "short" rather than "double" , "int", or any of the other modes of the vector being read. The data represents a regular grid of size 419 rows by 264 columns, to make it more interesting, the data are daily records, for a total of 37 years. The file size is