[R] Incorporate a shapefile with an package animation

Hello I'm working with NetCDF files in an animation and am trying to superimpose a shapefile on the image as it is generates the html pages. If I take out the lines for the shapefile, it works correctly, however I'm having difficultly including the shapefile. If anyone has any ideas I would

Re: [R] automating a script to read a file

Petr, R- Users The scheme you provided yesterday worked very well. I am now trying to add more information to the graph(s) to make them more informative for my purposes. Essentially, I am trying to use a double Y axis graphic that will include the skew normal curves that I wrote about

[R] automating a script to read a file

Hi, The following script (which I did not develop) is used to calculate and plot a skewed normal curve. The script currently requires the user to input six parameters, rather than reading these directly from a file. I've been spinning wheels here, trying to figure out how to modify the script

Re: [R] automating a script to read a file

Petr, Thank you very much this works. A little more tweaking and I'll have what I need. Thanks Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224

Re: [R] introducing R to high school students

Christopher, I originally thought about writing off list to avoid a plethora of babble. However here goes. I don't see any reason why good students can 't learn the fundamentals of R. It has lots of advance methods that perhaps are too complex to handle for younger - less experienced people.

Re: [R] scatter3d: problem with spheres-color

I do something similar my script is structured as follows: perhaps it will help library(car) library(rgl) library(mgcv) scatter3d(MADep2009 ~ MADwet2009 + MedAD2009 | as.factor(VegClass), data = SFCN, sphere.size=1.5, surface=FALSE, parallel=FALSE, elliposiod=TRUE, surface.col=c(green, red,

[R] Calculating the overlapping area of ellipsoides

I'm generating several scatter3d visualizations of vegetation community distribution against a number parameters using the scatter3d function car package, The visualizations are fantastic and they provide great support for the analysis. As far as I can determine, scatter3d does not save the

[R] Help with rJava

I am trying to invoke an in house developed java program from R and am running into some confusing issues. I've written this function: DoNormalSkew - function() { .jaddClassPath(U:/development/software/untilities/ELVeSkew/dist/ELVeSkew.jar)

[R] Is there a function for scatter3d with Categorical responses?

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

Re: [R] Is there a function for scatter3d with Categorical responses?

John, I would like ellipsoids to represent the distributions relative to each of the vegetation types. I tried your approach thinking that it would work as well, but I'm getting an error scatter3d(MADep2004 ~ Dhydro2004 + MedWet2004 | VegtypeID, ellipsoid=TRUE. data= CSSS) Error in sprintf

Re: [R] Contingency tables example

This is a elementary problem concerned with tests of association. You should read up on categorical statistics, chi square tests, and general linear models. Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor)

Re: [R] Upgrade R?

Why don't you just download the latest release from CRAN - R. It is the recommended approach to installing R. Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov

[R] Plotting skewed normal distribution with a bar plot

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.

Re: [R] Random Forest Classification

Explore the ModelMap package. It might offer some useful tools for your application. Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax

Re: [R] Logistic Regression - Variable Selection Methods With Prediction

Try the glm package Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

This is the first time I've seen an R Style Guide. I will admit that I haven't looked for one previously, but nevertheless I still haven't seen one. My code style simply evolved (perhaps, chugged along) by reading posts from other users who post to the r-help community. I regularly program with

Re: [R] join tables in R

Look at the merge command ?merge Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] Syntax coloring in R console

RStudio is a new interface. I just started using it and find it very good. You might want to check it out. It does have Ubuntu capabilities, but you have be sure that your graphic drivers are the same as is specified in the RStudio system. Hope this helps Steve Friedman Ph. D. Ecologist /

[R] filled contour plot with contour lines

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

Re: [R] filled contour plot with contour lines

Hi Peter, Thanks for taking the time to consider the problem. I realize the procedure is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in size to pull data into this process. Nobody in their right mind should work with such things. It's a spatial temporal database with 10 years of

Re: [R] filled contour plot with contour lines

Hi Peter, Thanks for taking the time to consider the problem. I realize the procedure is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in size to pull data into this process. Nobody in their right mind should work with such things. It's a spatial temporal database with 10 years of

[R] MultiDimensional arrays

hello, I'm working with a 4 dimensional NetCDF file, trying to calculate some summary statistics on a parameter. The dimensions are number of rows and columns, time and days. The attribute in question in this case happens to be a probability of woodstork foraging success. But it could be

[R] Setting maximum value of the legend on an image.plot and animation

Hello, I'm doing the following: library(ncdf) library(fields) library(animation) saline - open.ncdf(salinity_1990.nc) salt = get.var.ncdf(nc=saline, varid=Salinity) # create an animation of the complete temporal domain in the ncdf file. saveHTML({ for (i in 1:364) {

Re: [R] Tinn-R 2.3.7.0 released

I am also finding the link between TINN - R (2.3.7.0) and R (2.12.0 2010 - 10 - 15) to be problematic. source(.trPaths[5], echo=TRUE, max.deparse.length=150) Error in source(.trPaths[5], echo = TRUE, max.deparse.length = 150) : object '.trPaths' not found Steve Friedman Ph. D. Ecologist

[R] Setting default path to library

Hello I recently upgraded from 2.11.1 to 2.12.0 on a windows machine. When I launch R via TINN - R,(2.3.7.0) most things appear correct. The exception is the path to the library. I store all of the packages in C:\Program_Files \R\R-2.12.0\library Last week when I upgraded I rec'd an error:

[R] repeated measure test

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

[R] Clustering with ordinal data

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

Re: [R] Clustering with ordinal data

Thanks Phil, I'll do so now. Much appreciated. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] lower triangle of the correlation matrix with xtable

try lower.tri and see ??lower.tri Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] Data format question for triangle.plot package ade4

hello, I am trying to develop a triangle plot but am having difficultly assigning the row.names to the 3 columns in the data.frame Here is what I've done, attach(SoilVegHydro) dim(SoilVegHydro) 129239 # now take 3 variables from main data.frame for plotting dat - cbind.data.frame(TP,

[R] Ctree Question

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,

[R] variograms and kriging

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,

Re: [R] variograms and kriging

Please disregard. I've posted to the wrong site. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] ctree

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)

Re: [R] ctree

Joris, Thank you, your solution did what I needed. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] conditional inference trees extracting terminal node probabilities

The ctree function (package party) provides a method for running conditional inference trees. Plotting results of ctree returns a binary map of the tree and for each terminal node a barplot of the probabilities of the response categories. For example: iris.ct - ctree(Species ~ . , data =

[R] script development for Unconditional Density and Probability estimation

Hello, I'd like to automate this script a bit more and cycle several parameters(both the species and the metric). For example where AnnualDepth occurs, I need to process about 12 metrics so instead of writing this entire script 12 times once for each metric I'd like to be able to automatically

[R] Plotting Question

Hello, I would like to produce a series of graphs comparing the probability distributions for 8 factors against a continuous metric. The kind of graph I'm hoping to produce would look like the density comparison graphs (library sm) using the function sm.density.compare. However, instead of

[R] restructuring by output for use with write.table

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

[R] General Question regarding Lists

Hello, Here is a general question that I think many of you will have insight to. I have 40 independent dataframes. Each has a dimension of 1359 x 15. For example dim(meanAnnualWaterDepth) 1539 15 I need to calculate the average value for each dataframe but only for columns 7: 15. I can do

[R] Help with RandomForest

I'm working with the randomForest package and have successfully build a model. I'd like to go one step further however, and use the output from the model to construct a map using the output. My input data are spatial, and I have an independent set of rasterize maps for each of the predictor

[R] Installing randomForest on Ubuntu Errors

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'

Re: [R] mlogit is not an S4 object error

Achim - Everyone Yes, I did exactly as you've suggested. I did the analysis over on both my windows XP machine and on my Linux machine. It worked correctly. I don't recall which other packages I installed the first time, but there was a conflict. Thanks you for the attention and the great

Re: [R] cluster analysis

Without know what your data set really looks like, I'd look to decision trees - specifically package rpart and use method = classify. Your problem may not be appropriate in that environment, but it is hard to say with limited explanation of issues. good luck Steve Friedman Ph. D. Spatial

[R] mlogit: Error reported using sample dataset

I've been working on a multinomial logit model, trying to predict vegetation types as a function of total phosphorus. Previous responses to my postings have pointed me to the mlogit package. I'm now trying to work examples and my data using this package. data(Fishing, package = mlogit) Fish -

Re: [R] number of Excel worksheets

Kevin This function is very valuable to me as well. Can you explain what the function is looking for with regards to the arguments sheet and condition? Thanks Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead,

[R] glm models with more than one response

Hi everyone, I am trying to construct a glm and am running into a couple of questions. The data set I am using consists of 6 categories for the response and 6 independent predictors representing nutrient concentrations at sample point locations. Ultimately I'd like to use the probabilities for

[R] Install R 2.10.1 on Windows XP Errors

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,

Re: [R] OT: Software for specific visualisation of data...ideas?

Whats wrong with Power Point or anyone of its equivalents? Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] Formula for normal distribution with know mean and standard error and n terms

Hello, I am searching for a method to calculate a normal distribution. For example this equation is used to calculate the normal curve when the mean and standard deviation are know. p(x) = (1/σ*sqrt(2π)) x exp (- (x-μ)2/2σ2) or (Embedded image moved to file: pic27350.jpg)Normal Probability

Re: [R] Drought Severity Index (DSI)

A few years ago, I work with Stuart Gage who had developed a Heat / Precipitation Index as a measure of drought severity. It works just as well if not better than the Palmer Drought Index. You can find the formula in this only pdf report: Climate Variability in the North Central Region:

[R] Upgrading To 2.10 from 2.6.2

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

Re: [R] Upgrading To 2.10 from 2.6.2

Ok, that seems to work. Thanks for the help. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] how to use different ylim scales on a lattice bwplot?

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

Re: [R] How to plot an image in R

You are going to have to subset the necdf dataset by dates or some other logical aggregation scheme, and then plot. Since the data is spatial, I would think that plotting could be achieved with the image( ) function. Good Luck Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and

[R] Help with time series

Hello I'm working with a bunch of time series data. The data are downloaded from a server and stored as ascii files prior to reading them into R. After reading the data sets read into R with no problem and I can us the ts function to coerce them to time series, sometimes this works and

[R] Time Series Format Error

I have rec'd the following error: P34annual - read.table(A:\\Data\\Output\\Sparrow\\Hydro_Data\\P34_Annual.txt, header=TRUE, sep=,, stringsAsFactors= FALSE, skip=1) P34annual$GS - rep(1.86, dim(P34annual)[1]) P34annual$Depth - as.numeric(P34annual$P34_stage) - as.numeric(P34annual$GS)

Re: [R] Tinn-R setup

Kevin, I did the same installation a week or two ago. The same downloading process occurred, but not so many iterations. I don't know if it was normal or not , but the system seems to work very well. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National

[R] Problem with R.2.9.2 Plot functions

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

[R] Preparing for multi-core CPUs and parallel processing applications

Hello I am fortunate (or in really big trouble) in that the research group I work with will soon be receiving several high end dual quad core machines. We will use the Ubuntu OS on these. We intend to use this cluster for some extensive modeling applications. Our programming guru has

[R] help with superscripts in simple plots

Hello All, When I use the following lines of code to create a plot and add labels with R-square values the labels have a superscripted R2. library(lattice) xyplot(PropHatchedNests$Phatched + PropHatchedNests$PropNests + PropHatchedNests$meanHSI + PropHatchedNests$RelMeanEggsNest ~

Re: [R] off topic but need your pointers about statistics

Chapter 10 A Twentieth-Century Tour of Categorical Data Analysis in Agresti's smaller book with this title An Introduction to Categorical Data Analysis includes this history. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor)

[R] dropping empty panels in xyplot

I note that someone asked for help a few weeks ago regarding the ability to drop empty panels on an xyplot. I did not see a reply to that inquiry so I am asking for assistance with the same problem To prepare the data for the plotting routine, I did the following to restrict the content of the

Re: [R] loglinear analysis

Look to the glm function then pass the output to the step function Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] Superscripts and subscripts in trellis graphics

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,

Re: [R] error in importing text files

Eric, The error suggests that the file has an inconsistent number of attributes. You should inspect the file to make sure that each line has the same number of expected attributes (columns). If one or more is missing, it suggests that the data collection process from the Palm Pilot is not what

Re: [R] Superscripts and subscripts in trellis graphics

Gabor's solution has worked very nicely. Thank you very much. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] I'm offering $300 for someone who know R-programming to do the assignments for me.

Shame on you! Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] I'm offering $300 for someone who know R-programming to do the assignments for me.

Wacek raises a very valid perspective. I have previously not responded to these requests, and I've seen a bunch. I have thought they were a waste of my time to respond to. In the future I think I'll simpy disregard these requests, because I simply can not crystal ball the true needs of the

[R] General Question about Graphics

Greetings I have encountered a situation with regards to plotting barcharts with associated error bars. My search for clues on how to accomplish this turned up some interesting information. Basically, I found that including error bars with barplots is not desirable and hence there appears that

[R] setting trellis auto.key color values

I'm working with Lattice graphics and I would like very much to color code the auto.key fill color with the same corresponding colors that I use in the panels. I've looked on the web for clues, and on the CRAN-R help sites searching on trellis auto.key color and variations, unfortunately the

Re: [R] setting trellis auto.key color values

Sundar, Thanks for the assistance. Much appreciated. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] setting key boxes in xyplot

David and the list Yes the graph is produced with the code I submitted. However, I was not clear in my initial post. and it is totally my mistake for taking up everyone's valuable time. My main question focuses on the key which prints in black as opposed to using the colors (black, red and

Re: [R] odbcConnectAccess function

Feilipe, You should look at the RODBC package It works great ! Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] setting key boxes in xyplot

I thought I mimicked the coded correctly, but find an error/omission somewhere in the xyplot function shown below. eggs # alligator egg fate dput(round(eggs, 2), file = ) structure(list(Year = c(1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000),

Re: [R] help working with date values

Gabor, These options will work just fine. Thank you Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] tiff graphics

My Linux installation of R is x86_64-redhat-linux-gnu-library When I tried to write a tiff file the system responded that tiff was not supported with this version. Which version of R on a linux environment does support tiff graphic files ? Thank you Steve Steve Friedman Ph. D. Spatial

[R] help working with date values

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

[R] Setting lattice par parameters

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,

Re: [R] Setting lattice par parameters

Sundar, Thank you very much. I see my mistake. Much appreciated. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

Re: [R] Reading in a large number of dbf files

Duncan, Thank you very much. I did not have a chance to examine this approach today. I will because I need to have this functionality. I appreciate you time and insights Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd

[R] Reading in a large number of dbf files

good morning This question is not a stats question per say but a data management and lattice plotting problem. I apologize now if I'm asking an inappropriate question to this gracious group. I'm need to bring in approximately 100 *.dbf files into R but I'm having difficultly understanding

[R] numbers not behaving as expected

If someone can explain this odd behavior I'd appreciate it. I have a data.frame Cell.ave (attached and created via dput(Cell.ave, Cell.ave) which contains three columns of parameters year, month and AveRain. I need to subset the data.frame by months such that DrySeaonMonths are 1,2,3,4, 11,

Re: [R] numbers not behaving as expected

My apologies for sending a binary file. I was following advice from someone (from this list) who insisted I send data via dput. I guess that is frown upon. Anyway. here is a Cellave.txt file. Stavros, here is the output following your suggestion. WetMonths - Cell.ave[Cell.ave$month = 5

Re: [R] help with ggplot2 -- ggpoint function missing?

I tried the example in ggplot as well and rec'd the same error message. ggpoint is not found in the ggplot2 library. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office

Re: [R] two monitors

Hello, I work with Windows XP and do what you are hoping to do. Instead of using the built in script dialogue editor I use TINN-R. I can drag either this editor or the R Console to different monitors. I often have to make sure that the TINN-R editor panel is not maximized as that will hide any

[R] labeling panels in lattice plots

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,

Re: [R] labeling panels in lattice plots

Sundar, Your suggestion to convert year to a factor worked like a charm. Thank you very much. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282

[R] Requesting help with lattice again

Hello, this is a request for assistance that I submitted earlier, this time with the dataset. My mistake for taking up bandwidth. I've also rephrased the question to address an additional concern. I'm working on a windows XP machine with R 2.8.1 1). I'd like a barchart (or other lattice type

[R] Help with lattice

Hello. I'm working on a windows XP machine with R 2.8.1 I'm working with lattice and I don't understand how to set the size of the text for the labels. Can anyone suggest a solution? I'm establishing a melted dataframe first and then setting a key and calling the plot as follows;

[R] Requesting assistance installing ESS for R on Redhat

I'm moving my R applications to a Redhat OS and want to install ESS. My sys admin has downloaded the rpm (emacs-common-ess-5.3.8-1.fc8.src.rpm), but when he tried to do the install he rec'd a number of warnings. We are not sure if the warnings are telling us that the installation did not occur

Re: [R] Using very large matrix

I'm very interested in the bigmemory package for windows 32-bit environments. Who do I need to contact to request the Beta version? Thanks Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034

[R] Trouble - Installing R on RedHat el5

I am lucky and am now working with a new Redhat Linux 64-bit OS But I am getting the following error can anyone provide some graciously needed assistance: [r...@bluebird system-files]# rpm -i R-2.8.1-1.rh5.x86_64.rpm warning: R-2.8.1-1.rh5.x86_64.rpm: Header V3 DSA signature: NOKEY, key ID

[R] Setting initial starting conditions in scripts

Hello, I'm writing a variety of R scripts and want to code the loadhistory and workspace from within the script. I found the loadhistory function but do not see a comparable function for load workspace. Is there one ? Working with R 2.8.1 (2008-12-22) on a windows platform. Thanks for any

[R] Changing Ylab and scale in hclust plots

Hello, Running R 2.8.1 (2008-12-22) on Windows. I running a series (25) of clustering procedures using the hclust function and would like each of the plots to have the same yaxis label and scale in all of the plots. Is there a procedure to change the scale on these plots? Or is there an

[R] Reading Binary Files

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

[R] More help with Binary Files

Does anyone else have any insights to this issue: Henrick, thank you for your very quick response. I've examined the readBin help file with respect to endian and I'm still not sure I'm getting this correct. Here is what I'm coding: con - file(file.choose(), open=rb) Year66 - readBin(con,

[R] Upgrading to TINN - R 2.1.1.6

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:

[R] help with time series

Hello everyone I'm working with R 2.8.1 on a windows machine I have a question regarding time series analysis The first question is how does R expect the input file to be structured? I'm working with a *.txt file similar to the abbreviated one here: Date,stage 4/2/1953,7.56 4/3/1953,7.56

Re: [R] Useful books for learning the R software and the S programming language

Look here to start. http://www.r-project.org/doc/bib/R-books.html Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

  1   2   >