[R] Learn R in a Day - new ebook

2014-02-22 Thread Steve Murray
ader is encouraged to replicate on their computer. Each chapter also includes exercises (with solutions) to practice key skills and empower the reader to build on the essentials gained during this introductory course. Steve Murray [[alternative

[R] Summing data frame columns on identical data

2011-01-17 Thread Steve Murray
Dear all, I have 9 data frames, and I'm simply trying to sum the values of column 3 (on a row-by-row basis). However, there are a slightly different number of rows in each data frame, so I'm receiving the following error: "Error in Ops.data.frame(mrunoff_207101[3], mrunoff_207102[3]) :   + on

[R] Grouped bars in barplot

2011-01-12 Thread Steve Murray
Dear all, I am trying to make a barplot with clustered pairs of bars, using class=numeric data and the following command: barplot(c(bline_precip[10,9], bline_runoff[10,9], cccma_precip[10,9], cccma_runoff[10,9], csiro_precip[10,9], csiro_runoff[10,9], ipsl_precip[10,9], ipsl_runoff[10,9], m

[R] Alignment of lines within barplot bars

2010-09-09 Thread Steve Murray
Dear all, I have a barplot upon which I hope to superimpose horizontal lines extending across the width of each bar. I am able to partly achieve this through the following set of commands: positions <- barplot(bar_values, col="grey") par(new=TRUE) plot(positions, horiz_values, col="red", pch="

Re: [R] Unwanted boxes in legend

2010-04-21 Thread Steve Murray
e: [R] Unwanted boxes in legend > > On 2010-04-21 4:35, Peter Ehlers wrote: >> The 'border' argument was added in 2.1.10. > > Egad! Did I really type that? > I meant 'in R 2.10.0'. > > -Peter Ehlers > >> >> >> On 2010-04-21 1:53, St

Re: [R] Unwanted boxes in legend

2010-04-21 Thread Steve Murray
t, >>> >>> Steve >>> >>> >>> >>>> Date: Thu, 15 Apr 2010 12:13:40 -0600 >>>> From: ehl...@ucalgary.ca >>>> To: smurray...@hotmail.com >>>> CC: r-help@r-project.org

Re: [R] Unwanted boxes in legend

2010-04-19 Thread Steve Murray
r-project.org > Subject: Re: [R] Unwanted boxes in legend > > On 2010-04-15 11:10, Steve Murray wrote: >> >> Dear all, >> >> I am using the following code to generate a legend in my plot (consisting of >> both bars and points), but end up with boxes around

[R] Unwanted boxes in legend

2010-04-15 Thread Steve Murray
Dear all, I am using the following code to generate a legend in my plot (consisting of both bars and points), but end up with boxes around my points: legend(10, par("usr")[4], c("A", "B", "C", "D"), fill=c(NA,NA, "grey28", NA), pch=c(16,4,NA,18), col=c("red","blue","grey28","yellow"), lty=FALS

[R] Alignment of x-axis labels

2010-04-15 Thread Steve Murray
Dear all, I'm having trouble getting the correct spacing between x-axis labels on a barplot. This is the command I'm using to generate the plot: temp <- barplot(precip, beside=TRUE, xaxt="n", las=1, xpd=FALSE, col="grey28", ylim=c(0, max(precip))) Here is the structure of temp: > str(temp)  n

Re: [R] Summing data based on certain conditions

2010-04-02 Thread Steve Murray
ct.org > Subject: Re: [R] Summing data based on certain conditions > > ?by may also be helpful. > > Stephan > > > Steve Murray schrieb: >> Dear all, >> >> I have a dataset of 1073 rows, the first 15 which look as follows: >> >>> data[1:15,] >>

[R] Summing data based on certain conditions

2010-03-31 Thread Steve Murray
Dear all, I have a dataset of 1073 rows, the first 15 which look as follows: > data[1:15,]     date year month day rammday thmmday 1   3/8/1988 1988 3   8    1.43    0.94 2  3/15/1988 1988 3  15    2.86    0.66 3  3/22/1988 1988 3  22    5.06    3.43 4  3/29/1988 1988 3  29  

[R] Legends for line and point data

2010-03-29 Thread Steve Murray
Dear all, I have a plot which contains 4 data series displayed using solid, dashed and dotted lines, and also points. How do I use lty and pch together to signify that the first legend item is a solid line, the second is point data (pch=16), the third is dashed and the fourth is dotted? Many

[R] Creation of a new data frame

2010-03-25 Thread Steve Murray
Dear all, I have a data frame of 18556 rows and 19 columns and wish to create a new grid from these data of dimensions 360 rows and 720 columns. The existing data frame is set up so that every 38 rows makes up one row of the new data frame, with 2 NA values at the end of each 'block' that shou

Re: [R] Filling a grid based on existing data

2010-03-25 Thread Steve Murray
-- > From: dwinsem...@comcast.net > To: smurray...@hotmail.com > Subject: Re: [R] Filling a grid based on existing data > Date: Wed, 24 Mar 2010 20:50:40 -0400 > > > On Mar 24, 2010, at 2:34 PM, Steve Murray wrote: > >> >> Dear all, >> >&

[R] Filling a grid based on existing data

2010-03-24 Thread Steve Murray
Dear all, I currently have a data frame of dimensions 18556 rows by 19 columns. I want to convert this into a grid of dimensions 720 rows by 360 columns. The problem in this case is that not all rows in the initial data frame are complete (there are gaps). Therefore I am perhaps looking for a

[R] Calculating mean for a number of columns

2010-02-10 Thread Steve Murray
Dear all, I am attempting to perform what should be a relatively simple calculation on a number of data frame columns. I am hoping to find the average on a per-row basis for each of the 50 columns. If on a particular row a 'NA' value is encountered, then this should be ignored and the mean for

Re: [R] Resampling a grid to coarsen its resolution

2010-02-10 Thread Steve Murray
That sounds like a sensible way of dealing with the - values... ...but doesn't solve the more important question of how to perform the resampling. Are there are functions in R which have been designed to achieve this? Or is there a standard way of going about this? Many thanks for any advi

[R] Resampling a grid to coarsen its resolution

2010-02-09 Thread Steve Murray
Dear all, I have a grid (data frame) dataset at 0.5 x 0.5 degrees spatial resolution (720 columns x 360 rows; regular spacing) and wish to coarsen this to a resolution of 2.5 x 2.5 degrees. A simple calculation which takes the mean of a block of points to form the regridded values would do the

Re: [R] Counting by rows based on multiple criteria

2010-02-08 Thread Steve Murray
Spot on, thanks very much indeed. Steve _ Got a cool Hotmail story? Tell us now __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

[R] Counting by rows based on multiple criteria

2010-02-08 Thread Steve Murray
Dear all, I have a data frame of 6 columns and ~6 rows which I hope to perform the following calculation on. For each row, I wish to determine whether there are a greater number of positive or negative numbers. Then, if there are more positive numbers in the row, count how many occur - bu

[R] ifelse on a series of rows for multiple criteria

2010-02-05 Thread Steve Murray
Dear all, I am attempting to perform a calculation which counts the number of positive (or negative) values based on the sample mean (on a per-row basis). If the mean is>0 then only positive values should be counted, and if the mean is <0 then only negative values should be counted. In cases w

[R] Finding ranges on a per-row basis from several objects

2010-02-05 Thread Steve Murray
Dear all, I'm attemping to find the overall range of values from column 5 of a series of data frames, on a per-row basis, and assign the results to a new object. At present, I'm only able to receive the overall range of all values, whereas I'm intending to get the results of the range for each

Re: [R] Isolating a column within a loop

2010-02-01 Thread Steve Murray
Thanks! I've learnt something there! Steve _ We want to hear all your funny, exciting and crazy Hotmail stories. Tell us now __ R-help@r-project.

[R] Isolating a column within a loop

2010-02-01 Thread Steve Murray
Dear all, I am trying to calculate the mean of one column for many data frames. The code I am using is as follows: > for (i in seq(nrow(index))) {         assign(paste(model, "_mean_",index$year[i], index$month[i], sep=''),   mean(get(paste(model, index$year[i], index$month[i], "[,3

Re: [R] R Output and ArcGIS

2010-01-26 Thread Steve Murray
Dear all, Just to let you know that thanks to your help, I've managed to solve it. For future reference, if anyone's interested (!), if you're having problems reading R-generated data from a Mac, into ArcMap on a PC, then ensure that you're using eol="\r\n" in the write.table command and that

Re: [R] R Output and ArcGIS

2010-01-26 Thread Steve Murray
Dear all, Thanks for the replies so far. Just to emphasise, I'm not using Excel in any way. I have many many files to output, so it'd take considerable time to export from R, reprocess in Excel, then load into Arc! On a PC I'm able to go directly from R to ArcMap (9.3) without having to go vi

[R] R Output and ArcGIS

2010-01-25 Thread Steve Murray
Dear all, I've been using R on a Mac to process some data for export to ArcMap GIS (which only runs on Windows). ArcMap seems to require tab-delimited data (my data are in 3 columns), so I've been using the sep="\t" argument. However, this resulted in strange end-of-line characters when displa

[R] Looping multiple dimensions

2010-01-22 Thread Steve Murray
Dear all, I have 30 arrays, each with dimensions 720,360,12. The naming format for each of these 30 objects is: mrunoff_5221, mrunoff_5222... mrunoff_5250. For example: > str(mrunoff_5221)   num [1:720, 1:360, 1:12] NA NA NA NA NA NA NA NA NA NA ...  (the initial NA's are nothing to worry abo

[R] Testing the significance of gradients in trends

2009-12-01 Thread Steve Murray
Dear all, I want to determine if the slopes of the trends I have in my plot are significantly different from each other (I have 2 time-series trends). What statistical test is most suitable for this purpose and is it available in the R base package? Many thanks, Steve

[R] Testing for strength of fit using R

2009-11-26 Thread Steve Murray
Dear all, I am trying to validate a model by comparing simulated output values against observed values. I have produced a simple X-y scatter plot with a 1:1 line, so that the closer the points fall to this line, the better the 'fit' between the modelled data and the observation data. I am now

Re: [R] Subset returning unexpected result

2009-10-15 Thread Steve Murray
Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable MIME-Version: 1.0 Bill=2C It seems to be 'character' - odd...! > str(int1901$Latitude) =A0chr [1:61537] "5.75" "6.25" "6.75" "7.25" "7.75" "8.25" ... Thanks again=2C Steve > What does str(int1901)

Re: [R] Subset returning unexpected result

2009-10-15 Thread Steve Murray
<77eb52c6dd32ba4d87471dcd70c8d70001f11...@na-pa-vbe03.na.tibco.com> Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable MIME-Version: 1.0 Dear Bill and all=2C Yep you were right - for some strange reason (I'm not sure how...)=2C the l= atitude data were o

[R] Subset returning unexpected result

2009-10-15 Thread Steve Murray
Dear all, I am attempting to subset a data frame based on a range of latitude values. I want to extract the values of 'interception' where latitude ranges between 50 and 60. I am doing this using the following code, yet it doesn't return the results I expected: > test <- subset(int1901, Lati

Re: [R] Trendline for a subset of data

2009-10-09 Thread Steve Murray
Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable MIME-Version: 1.0 Thanks Mark=2C the reg.line trick seemed to work really well. David - hopefully the hex-text will have gone now - if not=2C please accept= my apologies as=2C this is=2C as far as I kn

Re: [R] Trendline for a subset of data

2009-10-09 Thread Steve Murray
; From: dwinsem...@comcast.net > To: smurray...@hotmail.com > Subject: Re: [R] Trendline for a subset of data > Date: Fri=2C 9 Oct 2009 09:27:43 -0400 > > > On Oct 9=2C 2009=2C at 5:50 AM=2C Steve Murray wrote: > >> >> Dear all=2C >> >> I am using abl

[R] Trendline for a subset of data

2009-10-09 Thread Steve Murray
Dear all, I am using abline(lm ...) to insert a linear trendline through a portion of my data (e.g. dataset[,36:45]). However, I am finding that whilst the trendline is correctly displayed and representative of the data portion I've chosen, the line continues to run beyond this data segment an

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread Steve Murray
month.abb should do the trick _ View your other email accounts from your Hotmail inbox. Add them now. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEA

[R] Using 'unlist' (incorrectly?!) to collate a series of objects

2009-08-20 Thread Steve Murray
Dear R Users, I am attempting to write a new netCDF file (using the ncdf) package, using 120 grids I've created and which are held in R's memory. I am reaching the point where I try to put the data into the newly created file, but receive the following error: > put.var.ncdf(evap_file, evap_di

[R] Creating a list of combinations

2009-08-20 Thread Steve Murray
Dear R Users, I have 120 objects stored in R's memory and I want to pass the names of these many objects to be held as just one single object. The naming convention is month, year in sequence for all months between January 1986 to December 1995 (e.g. Jan86, Feb86, Mar86... through to Dec95). I

[R] Replacing NA values in one column of a data.frame

2009-08-18 Thread Steve Murray
Dear all, I'm trying to replace NA values with - in one column of a data frame. I've tried using is.na and the testdata[testdata$onecolumn==NA] <- approach, but whilst neither generate errors, neither result in -s appearing - the NAs remain there! I'd be grateful for any advice o

[R] Reshape package: Casting data to form a grid

2009-08-18 Thread Steve Murray
Dear R Users, I'm trying to use the 'cast' function in the 'reshape' package to convert column-format data to gridded-format data. A sample of my dataset is as follows: head(finalframe) Latitude Longitude Temperature OrigLat p-value Blaney 1 -90-38.75 NA -87.75 17.10167

[R] Rounding to the nearest 5

2009-08-17 Thread Steve Murray
Dear all, A hopefully simple question: how do I round a series of values (held in an object) to the nearest 5? I've checked out trunc, round, floor and ceiling, but these appear to be more tailored towards rounding decimal places. Thanks, Steve _

Re: [R] Assigning values based on a separate reference (lookup) table

2009-08-14 Thread Steve Murray
Thanks to you both for your responses. 3 I think these approaches will *nearly* do the trick, however, the problem is that the reference/lookup table is based on 'bins' of latitude values, eg.>61, 60-56, 55-51, 50-46 etc. whereas the actual data (in my 720 x 360 data frame) are not binned, e.g.

[R] Assigning values based on a separate reference (lookup) table

2009-08-14 Thread Steve Murray
Dear R Users, I have a data frame of 360 rows by 720 columns (259200 values). For each value in this grid I am hoping to apply an equation to, to generate a new grid. One of the parts of the equation (called 'p') relies on reading from a separate reference table. This is Table 4 at: http://ww

[R] ncdf package problem - put.var.ncdf

2009-08-05 Thread Steve Murray
Dear all, I am attempting to convert 10 NetCDF files into a single NetCDF file, due to the data input requirements of a model I hope to use. I am using the ncdf package, version 1.6. The data are global-scale water values, on a monthly basis for 10 years (ie. 120 months of data in total; at p

[R] Writing a NetCDF file in R

2009-08-04 Thread Steve Murray
Dear all, I am attempting to convert 10 NetCDF files into a single NetCDF file, due to the data input requirements of a model I hope to use. I am using the ncdf package, version 1.6. The data are global-scale water values, on a monthly basis for 10 years (ie. 120 months of data in total; at pr

Re: [R] New line operator in mtext

2009-06-22 Thread Steve Murray
Thanks again for a very useful comment. That seems to have separated the text and put it onto separate lines. However, whilst this results in the text being centralised in relation to the axis, it means that the lower line is left-justified in relation to the upper line, rather than being cent

Re: [R] New line operator in mtext

2009-06-22 Thread Steve Murray
Thanks for the response, however, whilst this eliminates the 'new line' character from appearing, it doesn't actually cause a new line to be printed! Instead, I have the last few characters of the first character string overlapping with the first few characters of the next. How best should I c

[R] New line operator in mtext

2009-06-22 Thread Steve Murray
Dear R Users, I'm finding that when I execute the following bit of code, that the new line argument is actually displayed as text in the graphics device. How do I avoid this happening? mtext(side=2, line=5.5, expression(paste("Monthly Summed Runoff (mm/month)", "/n", "and Summed Monthly Preci

[R] Superscript in y-axis of plot

2009-06-16 Thread Steve Murray
Dear all, I've been trying to superscript the '2' in the following command (I don't want the '^' displayed), but as yet haven't had much luck. I've tried both the paste and expression commands, but neither have brought me any joy! mtext(side=2, line=5.5, "Monthly Precipitation (mm x 10^2/month

Re: [R] Aligning axis values when plotting more than one graph on same axes

2009-06-15 Thread Steve Murray
Thanks for the replies. This is a simple example which demonstrates exactly the problems I'm facing. As you can see, neither the x or y axes line up consistantly. > barplot(1:12, names.arg=substr(month.abb, 1,1)) > par(new=TRUE) > plot(1:12, 1:12, type="b", lwd=2, col="red", xaxt="n") I'd be

[R] Aligning axis values when plotting more than one graph on same axes

2009-06-15 Thread Steve Murray
Dear R Users, I am trying to plot a barchart with a line graph superimposed (using par(new=TRUE)). There are 12 bars and 12 corresponding points for the line graph. This is fine, except that I'm encountering two problems: 1) The position of the points (of the line graph) are not centred on the

[R] Memory problems when using ifelse

2009-05-25 Thread Steve Murray
Dear R Users, I have a data frame of 4 columns and ~58000 rows, the top of which looks like this: > head(max_out) Latitude Longitude Model Obs 1-0.25-49.25 4 4 2-0.25-50.25 4 5 3-0.25-50.75 4 4 4-0.25-51.25 311 5-0.25-5

[R] Changing a legend title to bold font

2009-05-22 Thread Steve Murray
Dear all, I have seen a response from Duncan Murdoch which comes close to solving this one, but I can't quite seem to tailor it to fit my needs! I am trying to make just the title in my legend as bold font, with the legend 'items' in normal typeface. I've tried setting par(font=2) external to

[R] Read.table problems

2009-05-18 Thread Steve Murray
Dear all, I have a file which I've converted from NetCDF (.nc) to text (.txt) using ncdump in Unix (as I had problems using the ncdf package to do this). The first few rows (as copied and pasted from the Unix console) of the file appear as follows: _, _, _, _, _, _, _, _, _, _, _, _, _, _, _

Re: [R] Simple plotting errors

2009-05-18 Thread Steve Murray
Many thanks once more for helping me to solve this. Gabor - I wasn't even aware of month.abb, so thanks for bringing this useful trick to my attention! Steve _ [[elided Hotmail spam]]

Re: [R] Simple plotting errors

2009-05-18 Thread Steve Murray
Thanks for all the useful information; use of 'c(...)' did the trick, although in future I'll try to hold the data in a more user-friendly setup. I've now got a plot, but have two issues that I can't seem to resolve: 1, The ylab is overlapping the y-axis tick mark values. I've tried using oma

[R] Simple plotting errors

2009-05-18 Thread Steve Murray
Dear R Users, I have 12 data frames, each of 12 rows and 2 columns. e.g. FeketeJAN MEANSUM_ AMAZON 144.4997874 68348.4 NILE 5.4701955 1394.9 CONGO71.3670036 21196.0 MISSISSIPPI 18.9273250 6511.0 AMUR 1.8426874 466.2 PARANA 58.38354

[R] Positioning a legend via X and Y coordinates

2009-05-06 Thread Steve Murray
Dear R Users, I'm able to display a legend using the following code: > legend("topright", c("Simulation", "Observation"), fill=2:3, bty="n") However, this causes the legend to be positioned too close to the bars in my barplot. I'd like to move the legend up slightly. I have been trying to det

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
Thanks Jim, that's great. Based on the information in the previous messages, is it possible to change the y-axis as I'd hoped? Thanks again, Steve _ [[elided Hotmail spam]] __ R-help@r-

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
> dput(total_sums) structure(c(17722202.6898231, 15276602.215475, 16875888.5155229, 14086271.625756, 18626581.9628846, 15387747.481166, 18428414.8535184, 15560882.404998, 17611181.5207881, 14905453.195546, 17290016.3934661, 14939493.120707, 16819288.8227961, 13979000.614402, 17657959.3656573,

Re: [R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
A very reasonable request! Sorry for not doing this initially, but please find below the data I am trying to plot: > total_sums sums86 sums87 sums88 sums89 sums90 sums91 sums92 Sim_1986 17722203 16875889 18626582 18428415 17611182 17290016 16819289 X198615276602 140862

[R] X and Y axis labels on a barplot

2009-05-06 Thread Steve Murray
Dear all, I have produced a barplot and wish to alter the axes a little. In place of the variable names appearing on the x-axis, I'd like to have the numbers 1986 to 1995. I have tried using the argument xlim=c(1986,1995) in the barplot command but receive: "Error in plot.window(xlim, ylim, lo

Re: [R] Plotting pairs of bars

2009-05-06 Thread Steve Murray
Jim and all, Thanks for the suggestion, however, I get the following error: > barplot(t(combine86[,1:2], beside = TRUE, las = 1)) Error in t(combine86[, 1:2], beside = TRUE, las = 1) : unused argument(s) (beside = TRUE, las = 1) I've looked up ?t and cannot see any extra arguments that I sho

Re: [R] Plotting pairs of bars

2009-05-06 Thread Steve Murray
Thanks for the reply - the 'beside' argument certainly looks useful, although I'm still not getting the output I'd hoped for. By doing: barplot(combine86[,1:2], beside = TRUE, las = 1, xlab=rownames(combine86)) ...I get all the bars for the 'Sim Mean' column plotted on the left side of the gra

[R] Plotting pairs of bars

2009-05-05 Thread Steve Murray
Dear all, I have a matrix called combine86 which looks as follows: > combine86 Sim Mean Obs Mean Sim Sum Obs Sum AMAZON 1172.0424 1394.44604 553204 659573 NILE 262.4440 164.23921 67973 41881 CONGO682.8007 722.63971 205523 214624 MISSISSIPPI 363.

[R] Producing a legend successfullly

2009-05-04 Thread Steve Murray
Dear all, I'm attempting to insert a legend into a line graph. I've sorted out the positioning, but I'm unable to display the sample line and associated colour to go within the legend box. Instead, under the variable names, the numbers 1, 2, 2, 3 are displayed in a column (with '2' repeated tw

[R] Reversing axis label order

2009-05-04 Thread Steve Murray
Dear R Users, I am executing the following command to produce a line graph: matplot(aggregate_1986[,1], aggregate_1986[,2:3], type="l", col=2:3) On the x-axis I have values of Latitude (in column 1) ranging from -60 to +80 (left to right on the x-axis). However, I wish to have these values sho

[R] Problems producing a simple plot

2009-05-02 Thread Steve Murray
Dear R Users, I have a data frame of the nature: > head(aggregate_1986) Latitude Mean Annual Simulated Runoff per 1° Latitudinal Band 1 -55574.09287 2 -54247.23078 3 -53

[R] Problems with 'valid columns' when using merge

2009-05-01 Thread Steve Murray
Dear all, I am trying to use 'merge' within a loop, however, I receive an error relating to the 'by' argument of the command, as follows: > merge_year <- 1986 > > for (i in 1:10) { # Number of file pairs + assign(paste("merged_arunfek_", merge_year, sep=""), merge(x=paste("arunoff_",start

Re: [R] Using 'aggregate' when dependent on row value increments

2009-04-30 Thread Steve Murray
Many thanks for the very useful responses in such a short time. I'm not a former SAS user - more a naive R user who didn't realise that a sort wasn't necessary! Jim, your solution worked really well - thanks. Thanks again for the great solutions. Steve ___

[R] Using 'aggregate' when dependent on row value increments

2009-04-30 Thread Steve Murray
Dear all, I have a data frame of three columns, which I have sorted by Latitude as follows: > test2[60:80,] Latitude Longitude Sim_1986 6194885.25-29.25 2.175345 6195785.25-28.75 8.750486 6196785.25-28.25 33.569305 6197785.25-27.75 23.702572 6198885.

Re: [R] Finding rows common to two datasets

2009-04-28 Thread Steve Murray
Dear Steve, > > Try > > ? intersect > > and see if that might help. > > Cheers, > Umesh > > On Tue, Apr 28, 2009 at 1:29 PM, Steve Murray> wrote: > > > > Dear all, > > > > I have 2 data frames, both with 14 columns of data and differing number

[R] Finding rows common to two datasets

2009-04-28 Thread Steve Murray
Dear all, I have 2 data frames, both with 14 columns of data and differing numbers of rows. The first two columns are 'Latitude' and 'Longitude'. I want to find the pairs of Latitude and Longitude coordinates which are common to both datasets, and output a new data frame which is composed of

Re: [R] Reshape - strange outputs

2009-04-10 Thread Steve Murray
Genius! Thanks very much Hadley - that was surprisingly easier to solve than I was anticipating! As a way of offering something back in return, I don't know if you plan to release a new version of the reshape package, but here's a suggestion to consider, just in case you do. On the basis of w

[R] Reshape - strange outputs

2009-04-08 Thread Steve Murray
Dear R Users, I am using the reshape package to reformat gridded data into column format using the code shown below. However, when I display the resulting object, a single column is fomed (instead of three) and all the latitude values (which should be in either column one or two) are collected

[R] Grouping of data frames to calculate means

2009-04-06 Thread Steve Murray
Dear R Users, I have 120 data frames of the format table_198601, table_198602... table_198612, table_198701, table_198702... table_198712 through to table_199512 (ie. the first 4 digits are years which vary from 1986 to 1995, and the final two digits are months which vary from 01 to 12 for eac

[R] Reshape: 'melt' numerous objects

2009-03-31 Thread Steve Murray
Dear R Users, I'm trying to use the reshape package to 'melt' my gridded data into column format. I've done this before on individual files, but this time I'm trying to do it on a directory of files (with variable file names) - therefore I have to also use the 'assign' command. I have come up

[R] Row/columns names within 'assign' command

2009-03-31 Thread Steve Murray
Dear all, I am attempting to add row and column names to a series of tables (120 in total) which have 2 variable parts to their name. These tables are created as follows: # Create table indexes index <- expand.grid(year = sprintf("%04d", seq(1986, 1995)), month = sprintf("%02d", 1:12)) #

Re: [R] Column name assignment problem

2009-03-30 Thread Steve Murray
Dear Peter, Jim and all, Thanks for the information regarding how to structure 'assign' commands. I've had a go at doing this, based on your advice, and although I feel I'm a lot closer now, I can't quite get it to work: rnames <- sprintf("%.2f", seq(from = -89.75, to = 89.75, length = 360)) c

Re: [R] Column name assignment problem

2009-03-30 Thread Steve Murray
Dear all, Apologies for yet another question (!). Hopefully it won't be too tricky to solve. I am attempting to add row and column names (these are in fact numbers) to each of the tables created by the code (120 in total). # Create index of file names files <- print(ls()[1:120], quote=FALSE)

Re: [R] Column name assignment problem

2009-03-30 Thread Steve Murray
Jim and all, Thanks - I managed to get it working based on your helpful advice. I'm now trying to do something very similar which simply involves changing the names of the variables in column 1 to make them more succinct. I'm trying to do this via the 'levels' command as I figured that I might

Re: [R] Column name assignment problem

2009-03-27 Thread Steve Murray
Dear all, I've been trying to implement the advice given to me, but without much success so far. I thought I'd provide the code in full in the hope that it might make more sense. Just to reiterate, I'm attempting to change the header of the 4th column of every table to "COUNT". year<- 1951:2

[R] Assignment to variables fails to loop

2009-03-27 Thread Steve Murray
Dear all, I think I'm nearly there in writing R code which will read in files with two variable parts to the file name and then assigning these file names to objects, which also have two variable parts. I have got the code running without encountering errors, however, I receive 50+ of the same

Re: [R] Reading in files with variable parts to names

2009-03-27 Thread Steve Murray
Thanks, that's great - just what I was looking for. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-conta

Re: [R] Reading in files with variable parts to names

2009-03-27 Thread Steve Murray
Dear all, Thanks for the help in the previous posts. I've considered each one and have nearly managed to get it working. The structure of the filelist being produced is correct, except for a single space which I can't seem to eradicate! This is my amended code, followed by the first twelve row

[R] Reading in files with variable parts to names

2009-03-26 Thread Steve Murray
Dear all, I'm trying to read in a whole directory of files which have two variable parts to the file name: year and month. E.g. comp198604.asc represents April of 1986 - 'comp' is fixed in each case. Years range between 1986 to 1995 and months are between 1 and 12. Just to be clear, there are

[R] Column name assignment problem

2009-03-26 Thread Steve Murray
Dear all, I'm trying to assign a name to the fourth column whilst using 'assign', but keep encountering errors. What have I done wrong?! > assign(colnames(c(paste("arunoff_",table_year, sep="")[4]), "COUNT")) Error in if (do.NULL) NULL else if (nc> 0) paste(prefix, seq_len(nc), : argument

Re: [R] Manual sort in a for loop

2009-03-26 Thread Steve Murray
Thanks all - I'm fairly new to R, so I was oblivious to the pros and cons of using a data frame as opposed to a list! The 'get' command also seemed to work successfully. Thanks again, Steve _ 25GB of FREE Online Storage – Find

[R] Manual sort in a for loop

2009-03-25 Thread Steve Murray
Dear all, I am trying to manually re-sort rows in a number of tables. The rows aren't sorted on any particular values but are simply ordered by user choice (as shown by the row numbers in the code). I have been able to carry out each re-arrangement without the use of the 'for' loop, but cannot

[R] Looping of read.table and assignment

2009-03-23 Thread Steve Murray
Dear all, I am trying to read in and assign data from 50 tables in an automated fashion. I have the following code, which I created with the help of textbooks and the internet, but it only seems to read in the final data file over and over again. For example, when I type:> table_1951 I get th

Re: [R] Working with duplicated rows

2009-01-12 Thread Steve Murray
Thanks - that's great! _ Choose the perfect PC or mobile phone for you __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide ht

[R] Working with duplicated rows

2009-01-12 Thread Steve Murray
Dear all, I have a dataframe of 3 columns, consisting of 'longitude', 'latitude' and a corresponding 'value'. Where identical 'longitude' and 'latitude' pairs occur more than once, I want their corresponding 'value' to be summed and the 'pair' to only appear once. For example: long lat

[R] Creating a Legend

2008-09-23 Thread Steve Murray
Dear all, I'm trying to create a legend for my graph. I hope to have the title as "Land Use Type" and the two elements being "Urban" and "Rural" with a red point and green point respectively. So far I have the following command, but obviously it isn't correct: > legend("topright", title="Land

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
say what the experiment died of. > ~ Sir Ronald Aylmer Fisher > > The plural of anecdote is not data. > ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answer can be extracted from a given body of &g

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
;> To call in the statistician after the experiment is done may be no more >> than asking him to perform a post-mortem examination: he may be able to >> say what the experiment died of. >> ~ Sir Ronald Aylmer Fisher >> >> The plural of anecdote is not data. >&g

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
rm a post-mortem examination: he may be able to > say what the experiment died of. > ~ Sir Ronald Aylmer Fisher > > The plural of anecdote is not data. > ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answ

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
dataset! Thanks, Steve > Date: Fri, 19 Sep 2008 21:57:03 +1000 > From: [EMAIL PROTECTED] > To: [EMAIL PROTECTED] > Subject: Re: [R] Use of colour in plots > > Steve Murray wrote: >> Greg, >> [[elided Hotmail spam]] >&g

Re: [R] Use of colour in plots

2008-09-18 Thread Steve Murray
gt; > Does this help, > > -- > Gregory (Greg) L. Snow Ph.D. > Statistical Data Center > Intermountain Healthcare > [EMAIL PROTECTED] > 801.408.8111 > > >> -Original Message- >> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] >> project.org] On Behalf

Re: [R] Use of colour in plots

2008-09-18 Thread Steve Murray
dote is not data. > ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answer can be extracted from a given body of > data. > ~ John Tukey > > -Oorspronkelijk bericht- > Van: [EMAIL PROTECTED] [ma

  1   2   >