r - Average count per group -


i real close missing last step, tell me need change below?

my dataset (mydata) looks (the first few entries missing):

date_daily station_date air_min air_avg air_max rain_mm temp_davg_c temp_dmax_c temp_dmin_c       date months years 1 2006/01/01         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-01     01  2006 2 2006/01/02         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-02     01  2006 3 2006/01/03         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-03     01  2006 4 2006/01/04         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-04     01  2006 5 2006/01/05         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-05     01  2006 6 2006/01/06         <na>    <na>    <na>    <na>      na          na          na          na 2006-01-06     01  2006

the daily data ranges 2006 2016, interested in average monthly rain average days per month rain > 2mm.

so tried following:

library(dplyr) stats <- mydata %>% group_by(months) %>%   summarise(days_w_rain = length(rain_mm[!is.na(rain_mm)]),             avg_rain = mean(rain_mm, na.rm=true),             days_2mm = sum(rain_mm > 2, na.rm=true),             p2mm = days_2mm / days_w_rain*100)

which gives me following:

# tibble: 12 x 5    months days_w_rain   avg_rain days_2mm      p2mm     <chr>       <int>      <dbl>    <int>     <dbl>  1     01         304 11.8078520      148 48.684211  2     02         310  7.4946903      125 40.322581  3     03         334  9.2680090      156 46.706587  4     04         327  3.7469113       92 28.134557  5     05         310  1.6073032       39 12.580645  6     06         300  1.1322100       39 13.000000  7     07         310  0.6212226       25  8.064516  8     08         333  1.6641862       47 14.114114  9     09         308  4.2796169       71 23.051948 10     10         310 10.7842742      155 50.000000 11     11         300 18.4062900      189 63.000000 12     12         287 20.6697038      155 54.006969

this correct number average rain per month gives me total number of days in januaries have rain (i.e. 304 days), instead of average per month of january, same amount of days have >2mm of rain (i.e. 148 days on months of january 2006 2016. want average per month , not total).

i sure easy fix don't know how it? suggestions appreciated cheers sandra

assuming have data years (2006 2016) give desired output:

 stats <- mydata %>% group_by(months) %>%           summarise(days_w_rain = length(rain_mm[!is.na(rain_mm)]/length(unique(years))),                     avg_rain = mean(rain_mm, na.rm=true),                     days_2mm = sum(rain_mm > 2, na.rm=true)/length(unique(years)),                     p2mm = days_2mm / days_w_rain*100)

this quick fix. dividing number of years , have average.


Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more