similar to: variable object naming

1) When adding additional data sets to a plot using "plot" followed by "lines", is there a way to automate the scaling of the axes to allow for all data sets to fit within the plot area? 2) I attempted to solve this by setting xlim=c(min(c(data1,data2,data3)),max(c(data1,data2,data3))) however, there are some NAs and Infs in these data sets, and min(data1) and max(data1) both

When reading a data set into R using readLines, na.strings="-99.99" is ignored. Is there an equivalent command for readLines? Alternatively, either immediately after reading into R or once the data set has been converted to a data frame, what is the appropriate command (or appropriate use of is.na) to convert my -99.99s to NAs? Thanks, Ben Osborne -- Botany Department University of

How can I extract K (number of parameters) from an AIC calculation, both to report K itself and to calculate AICc? I'm aware of the conversion from AIC -> AICc, where AICc = AIC + 2K(K+1)/(n-K-1), but not sure of how K is calculated or how to extract that value from either an AIC or logLik calculation. This is probably more of a basic statistics question than an R question, but I thank

What is the difference between the two R-squareds returned for a linear regression by summary(lm...)? When might one report multiple vs. adjusted R-squared? Thank you, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405 benjamin.osborne at uvm.edu phone: 802-656-0297 fax: 802-656-0440

I'm attempting to do model selection with AIC, using a glm and a lognormal distribution, but: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian(link="log")) ## gives the same result as either of the following: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian) fit1<-lm(BA~Year,data=pdat.sp1.65.04) fit1 #Coefficients: #(Intercept) Year2004 # -1.6341

I have a set of time-series climate data with missing entries. I need to add rows for these missing entries to this data set. The only way I know to do this is unsing a foreloop, but this won't work on a list. I've tried to convert the list to a data frame, but that won't happen, either. I want to fill rows in this table: > newtest[10:15,] yrmos yearmo snow.sum snow.mean

I want to organize outputs from several regressions into a handy table. When I try the following, each of my "fit_s" is replaces instead of read. Is there a way to read from the regression summaries that does not require writing separate lines of code for each? -Ben Osborne > fit1<-lm(dBA.spp16$sp2.dBA.ha~dBA.spp16$sp1.dBA.ha) >

Does anyone know of a list of functions that R already "knows?" I can't seem to find this anywhere in the help documentation. For example, I want to count the number of occurences of a certain value in a column of a data frame: What do I have to do to tell R to "Count?" Thanks, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405

Hello! Thanks for reading this request for assistance. I have a question regarding creating a histogram-like figure from data that are not currently in the correct format for the "hist" command. Specifically, my data have been processed and are in a matrix with columns containing the variables of interest and separate columns containing the number of times this variable was

I have one data frame with a column of dates and I want to fill another data frame with one column of dates, one of years, one of months, one of a unique combination of year and month, and one of days, but R seems to have some problems with this. My initial data frame looks like this (ignore the NAs in the other fields): > mans[1:10,] date loc snow.new prcp tmin snow.dep tmax 1

Hello, I have a question regarding how to subset/select parts of a data frame (matrix) in order to plot data associated only with this subset. Specifically I have a large data frame in which one column contains ID values (dates), and other columns contain data I would like to plot (temperature, light, etc.). I would like to break up this large matrix so as to plot data associated

Message: 63 Date: Wed, 26 Jan 2005 04:28:51 +0000 (UTC) From: Gabor Grothendieck <ggrothendieck at myway.com> Subject: Re: [R] chron: parsing dates into a data frame using a forloop To: r-help at stat.math.ethz.ch Message-ID: <loom.20050126T052153-333 at post.gmane.org> Content-Type: text/plain; charset=us-ascii Benjamin M. Osborne <Benjamin.Osborne <at> uvm.edu>

I've downloaded and installed R 1.9.0 under OS X 10.3.3. I then removed my previous R 1.8.0 installation by following the instructions in the readme file, e.g., 1. drag the "StartR" icon into the Trash; 2) go inside your Home directory and search for the Library folder. Inside the Library folder in your Home directory you'll find a RAqua directory: drag the RAqua

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I downloaded the package source onto my Desktop, unzipped it, and then typed: R CMD INSTALL /Users/brianbeckage/Desktop/Hmisc * Installing *source* package 'Hmisc' ... ** libs g77 -fno-common -g -O2 -c cidxcn.f -o cidxcn.o g77 -fno-common -g -O2 -c cidxcp.f -o cidxcp.o g77

Dear List, I'd like to add text to a plot where a text symbol has both a superscript and subscript. For example, the variable S with a superscript 2 and a subscript t. I have been able to accomplish this using either expression( paste(S,atop(scriptstyle(2),scriptstyle(t))) ) or expression( {S[t]}^2 ) but the spacing isn't quite right (or rather what I'd like) using either of

Hi all, I'm having difficulty converting a 'dates' object to a POSIXct object: testDATES<-c(35947,35971,36004,36008,36053,36066) testDATES<-chron(dates=testDATES, format = c(dates = "m/d/y"), origin=c(month = 12, day = 30, year = 1899)) >[1] 06/01/98 06/25/98 07/28/98 08/01/98 09/15/98 09/28/98 > as.POSIXct(testDATES) [1] NA NA NA NA NA NA

I'm having trouble loading Ryacas. I've downloaded and extracted Ryacas 0.2-9 (also tried Ryacas ) and yacas 1.0.63, have the latest version of R and have tried the following (this works for installing other packages): install.packages("Ryacas") library(Ryacas) install.packages("yacas") library(yacas) and get this: Warning message: package ‘yacas’ is not available

Hi all, Earlier today I posted this question on s-news, so apologies to some for the duplication. > Please put aside your snobbery about Microsoft products for a moment. > > I am fitting population models to annual survey data for trout. For those of you familiar with ecological models, I am working in the Lefkovitch matrix framework; for those unfamiliar with that shorthand, the