similar to: How to produce rainfall maps

Hi, You might get more help from the R-sig-geo list, which is devoted to spatial topics. However. The *.asc file is an ArcGIS raster export format. You should use whatever the appropriate import commands are for your own gridded rainfall data. If you have a different format, you might or might not be able to import it directly with raster. ?raster will tell you more about the kinds of formats

Fwiw the engine behind geom_raster needs explicit observation-per-row form for input (with no structural normalization), so conversion to points is perfectly proper here, albeit confusing in context. (It's closer to what graphics devices actually use ultimately, but the expansion is laid out very early in ggplot2 because there's no standard for intermediate forms.) Cheers, Mike On Wed,

Thank you Sarah and Mike for your explanations. My final objective is to produce maps (png image or any kind of extension I can import in LaTeX) where rainfall data are interpolated, using the Inverse Distance method or Kriging. My input file (pointfile.csv in the reported example) reports the station code, lat and long of the meteorological station and the rainfall value (which might be the

Dear list users, let me ask you this trivial question. I worked on that for a long time, by now. Suppose to have a data frame with NAs and to sum some columns with rowSums: df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) df[1, ] <- NA rowSums(df[ , which(names(df) %in% c("A","B"))], na.rm=T) If all the elements of the selected columns are NA, rowSums

On 21/03/2018 11:44 AM, Stefano Sofia wrote: > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) > df[1, ] <- NA > rowSums(df[ , which(names(df) %in%

What do you mean by "should not"? NULL means "missing object" in R. The result of the sum function is always expected to be numeric... so NA_real or NA_integer could make sense as possible return values. But you cannot compute on NULL so no, that doesn't work. See the note under the "Value" section of ?sum as to why zero is returned when all inputs are removed.

Dear list users, I have to aggregate wind direction data (wd) using a function that requires also a second input variable, wind speed (ws). This is the function that I need to use: my_fun <- function(wd1, ws1){ u_component <- -ws1*sin(2*pi*wd1/360) v_component <- -ws1*cos(2*pi*wd1/360) mean_u <- mean(u_component, na.rm=T) mean_v <- mean(v_component, na.rm=T) mean_wd

No. The empty sum is zero. Adding it to another sum should not change it. Nothing audacious about that. This is consistent; other definitions just cause trouble. -pd > On 21 Mar 2018, at 18:05 , Boris Steipe <boris.steipe at utoronto.ca> wrote: > > Surely the result of summation of non-existent values is not defined, is it not? And since the NA values have been _removed_,

Should not the result be NULL if you have removed the NA with na.rm=TRUE ? B. > On Mar 21, 2018, at 11:44 AM, Stefano Sofia <stefano.sofia at regione.marche.it> wrote: > > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <-

Surely the result of summation of non-existent values is not defined, is it not? And since the NA values have been _removed_, there's nothing left to sum over. In fact, pretending the the result in that case is zero would appear audacious, no? Cheers, Boris > On Mar 21, 2018, at 12:58 PM, Jeff Newmiller <jdnewmil at dcn.davis.ca.us> wrote: > > What do you mean by

I see: consistency with additive identity. That makes sense. Thanks. B. > On Mar 21, 2018, at 1:22 PM, peter dalgaard <pdalgd at gmail.com> wrote: > > No. The empty sum is zero. Adding it to another sum should not change it. Nothing audacious about that. This is consistent; other definitions just cause trouble. > > -pd > >> On 21 Mar 2018, at 18:05 , Boris

Dear R-users, I am using R version 1.8.0-1 under Suse 8.2. I need to use the boxcox command because I want to apply a Cox-Box transformation to a vector of rainfall values. Within the libraries, the MASS library is present, but I don't know whether this means that is automatically installed or not. The command doesn't work. What do I have to do in order to make it working? Is just a

?s 15:51 de 13/05/2023, Stefano Sofia escreveu: > Dear list users, > > I have to aggregate wind direction data (wd) using a function that requires also a second input variable, wind speed (ws). > > This is the function that I need to use: > > > my_fun <- function(wd1, ws1){ > > u_component <- -ws1*sin(2*pi*wd1/360) > v_component <-

"I see: consistency with additive identity. " Ummm, well: > 1+NULL numeric(0) > sum(1,NULL) [1] 1 Of course, there could well be something here I don't get, but that doesn't look very consistent to me. However, as I said privately, so long as the corner case behavior is documented, which it is, I don't care. Cheers, Bert Bert Gunter "The trouble with

Hi all; I am running the following ggplot codes. Runs well. However, I need to reflect the numeric values of the log10_P to the point size in the graph. Your help highly appreciated, Regards, Greg p <- ggplot(mydata, aes(x = X, y = log10_P)) + theme_bw() +theme(panel.border=element_blank()) + theme(legend.position="top", axis.text=element_text(size = 8)) (p1 <- p +

I have code like this: data <- read.csv("test1.csv", stringsAsFactors=FALSE, header=TRUE) # Graph myplot=ggplot(data, aes(fill=condition, y=value, x=condition)) + geom_bar(position="dodge", stat="identity", width=0.5) + scale_fill_manual(values=c("#7b3294", "#c2a5cf", "#a6dba0", "#008837"))+

Dear R-users, suppose I have three dataframes like these df1: mydate min_temp 31032013 12 01042013 8 02042013 -999 df2: mydate min_temp 31032013 10 01042013 11 02042013 14 df3: mydate min_temp 31032013 4 01042013 3 02042013 5 where -999 means that the temperature data is not available (at the moment I cannot change it to NA because I am not the db administrator); suppose also that oggi is

You could add size = log10_P to the aes() inside geom_point(). Untested code below. See also http://ggplot2.tidyverse.org/reference/geom_point.html ggplot(mydata, aes(x = X, y = log10_P)) + theme_bw() + theme(panel.border=element_blank(), legend.position="top", axis.text=element_text(size = 8)) + geom_point(aes(color = Traits, size = log10_P)) Jean On Fri, Jun 16, 2017 at

Dear R-list users, I need to use strptime because I have to deal with date with hours and minutes. I read the manual for strptime and I also looked at many examples, but when I try to apply it to my code, I always encounter some problems. I try to change the default format, with no success. Why? How can I change the format? 1. init_day <- as.factor("2015-02-24-00-30")

Dear Jim and dear Enrico, thank you for your replies. Unfortunately your hints didn't solve my problem, and I am getting mad. Can I show you my whole process? I will be as quick as possible. I start from a data frame called Snow of the form year month day hh mm hs 2007 11 19 0 0 0.00 2007 11 19 0 30 0.00 2007 11 19 1 0 0.00 2007 11 19 1 30 0.00 2007 11 19 2 0 0.00 2007 11 19 2 30 0.00 2007 11