similar to: Oceanographic lattice plots?

Dear R listers, When using smooth.spline to interpolate data, results are generally good. However, some cases produce totally unreasonable results. The data are values of pressure, temperature, and salinity from a probe that is lowered into the ocean, and the objective is to interpolate temperature and salinity to specified pressures. While smooth.spline provides excellent values at the

R-help, I have the following file I want to import to R (some lines removed) Calibrated CTD data for station:00280001 Calibrated:23/8 2001, Salinity Unsmoothed, Fluorescence Uncalibrated Maximum observed depth: 36 m QUAL has one digit for each of pressure, temp., sal. and fluor. QUAL=1:Uncal., QUAL=2:OK, QUAL=6:Interp., QUAL=9:No data

I am exploring the use of loess for oceanographic applications and would like to plot the locations (longitude and latitude) points where the models (salinity~temperature*longitude*latitude,parametric="temperature") are fitted. Chambers and Hastie(1993) explains the locations are nodes of a k-d tree. but I cannot find anything about accessing this information. It would be useful to

In working through material on p.272 of MASS (4th ed.), I came across the following model formula: pet1.lm <- lm(Y ~ No/EP - 1, Petrol) I was at a loss to understand the use of "/" until I looked in "An Introduction [!] to R," where I found the explanation. My request is that more complete material on model formulae be lifted from "Introduction to R" (or

To follow up on my previous question, suppose a user R session wants to unload one workspace and load another within an R session. Is the following the correct sequence? 1. save.image() to save the current workspace as .Rdata in the current working directory. 2. rm(list=ls()) to remove everything from the workspace. 3. setcwd("xxx") to set the new working directory. 4.

R 2.4.1 on Windows XP. Question: In traditional graphics, is it possible to find out the height of a line of text in units that can be used in arithmetic and then in calls to text()? Context: I have written a function that draws a plot and then, depending on whether some arguments are TRUE or FALSE, draws various lines of text in the plot. The text lines may be turned on or off individually

Dear List, Is there a possibility to change the font-face for a part of the title of a plot? For example I have the following... plot(nirs, type="l", xlab="Wellenl?nge [nm]", col="darkslategray", main = "Spektrum Deschampsia caespitosa") ...and I would like to change the part of the title-string "Deschampsia caespitosa" to italics? Is

I hope to use the plotmath facility to print titles that mix math and values of R variables. The help for "plotmath" has an example, which after repeated reading, I find baffling. Likewise, I have read the help file for "substitute" (wqhich seems to be needed) without ever understanding what it does, other than being used in some magic incantations. I would like to do

I remember reading somewhere that locations on plots (in my case, arguments x and y for legend()) can be specified in several scales besides the usual data scale. I would like to set x and y as proportions of total plot size or something similar. Can anyone steer me to documentation on how to do it? -- Michael Prager, Ph.D. <Mike.Prager at noaa.gov> NOAA Center for

This is (correctly) filed under "non-reproducible." The corresponding notes say "solved by re-installing Windows (private comm)." I would like to clarify that this was solved by re-installing the patch known as "Windows 2000 Service Pack 2," not by reinstalling Windows itself. Perhaps my original message was not clear. I point this out because installing the

A minor inconvenience in updating an R installation is remembering which packages were installed previously. Has anyone written a script to inspect a previous installation, then get & install the same packages into the new installation? -- Michael Prager NOAA Center for Coastal Fisheries and Habitat Research Beaufort, North Carolina 28516 Opinions expressed are personal, not official.

I need help with replacing NaN with zero (the value '0') in my dataset. The reason is that I can't get it to graph because of the NaN in the dataset. I have tried: data[is.nan(data)] <- 0 that others have suggested in the help archives but this does nothing so I am not sure what I am doing wrong. Thanks in advance for your help. [[alternative HTML version deleted]]

library(StreamMetabolism) day <- (structure(c(10.08, 10.08, 10.05, 10.03, 10, 9.98, 9.96, 9.95, 9.95, 9.96, 9.96, 9.98, 10.01, 10.05, 10.06, 10.09, 10.11, 10.11, 10.13, 10.13, 10.15, 10.15, 10.13, 10.14, 10.11, 10.13, 10.14, 10.13, 10.12, 10.13, 10.14, 10.16, 10.18, 10.19, 10.23, 10.27, 10.33, 10.37, 10.45, 10.57, 10.6, 10.66, 10.73, 10.77, 10.84, 10.86, 10.87, 10.94, 10.98, 11.01, 11.05,

Gurus, Can I rely on the rownames() function, when applied to a matrix, always returning either NULL or an object of type character? It seems that row names can be entered as integers, but as of now (R 2.4.1 on Windows) the function returns character vectors, not numbers (which is my desired result). I am using elements of the returned vector to index the matrix. e.g., nams <-

hello, today while trying to extract data from a list for subsequent analysis, i stumbled upon this funny behavior on my system: > x<-c(0.1,0.9) > 1-x[2] [1] 0.1 > x[1] [1] 0.1 > x[1]==1-x[2] [1] FALSE > x[1]>1-x[2] [1] TRUE > x<-c(0.3,0.7) > x[1] [1] 0.3 > x[2] [1] 0.7 > 1-x[2] [1] 0.3 > x[1]==1-x[2] [1] FALSE but: >

Hi, I am thinking about IWLS vs ML estimation. When I use glm() for a 2-parameter distribution (e.g., Weibull), I can otain the MLE of scale parameter given shape parameter through IWLS. Because this scale parameter usually converges to the MLE. In this point, I am wondering: i) can you say that the direct MLE, which is obtained by maximizing a likelihood function, is equalvant to the indirect

I would like to start building R packages under Windows XP. I have programming experience and a minimal but working knowlege of many Unix (-like) programming tools. The package functions (for now) will be from R source, not C or Fortran. I've installed Rtools, Perl, the MS hhc, and so on. I am setting up a command shell with the correct PATH for building packages. (The master Windows PATH

Hello R users -- If I have a dataframe such as the following, named "frame" with the columns intended to be named col1 through col6, > frame col1 col2 cmlo3 col4 col5 col6 [1,] 3 10 2 6 5 7 [2,] 6 8 4 10 7 1 [3,] 7 5 1 3 1 8 [4,] 10 6 5 4 9 2 and I want to correct or otherwise change the

Colleagues, I am analyzing data collected during oceanographic cruises. We have conducted many cruises over the last decade. On each cruise we visit ~50 stations. At each station (termed EventNum)we lower an instrument that measures depth, temperature, salinity and oxygen every few seconds as it is lowered through the water. Data from all EventNums on all cruises are stacked, generating a data

Hi everyone, I am having a problem reading my data file in R. It only reads part of it, and stops doing it in half way through the data. I looked at the observation where reading stops, but it is really not different from the observations read it. I expanded the matrix size, but it did not help. I am using a text format to read the data file. Could someone help me? Thank you Soyoko