aggregate - Find value relative to mean for particular Day of Year in R -


i'm working meteorology data in r, , conceptually, i'm trying find out how day above/below average. this, want separate day of year, find average doy (e.g. average january 1 temperature?), , compare every date (e.g january 1, 2014 anomalously warm, how much?)

i can find 'mean' table every day of year using aggregate:

head(data)           x       date 1  5.072241 1970-01-01 2  6.517069 1970-01-02 3  4.413654 1970-01-03 4 11.129351 1970-01-04 5  9.331630 1970-01-05  library(lubridate) temp = aggregate(data$x, list(yday(data$date)), mean)

but i'm stuck how use aggregated table compare original data.frame, see how x @ 1970 jan 1 relates average jan 1 x.

we can remove 'year' part sub ('monthday'). use ave if mean variable needs created grouped 'monthday'. 

data$monthday <- sub('\\d+-', '', data$date) data$mean <- with(data, ave(x, monthday))

then, can compare 'x' variable, example

data$rel_temp <- with(data, x/mean)

-

Comments

Post a Comment

Popular posts from this blog

css - SVG using textPath a symbol not rendering in Firefox -

using textpath within symbol doesn't render in firefox. renders fine in latest chrome , ie, when try reference symbol, svg text doesn't render in firefox (37.0.1) - first box appears empty. code below (no external dependencies), there should 2 boxes word test flowing vertically in centre of each. edit: thought issue somehow involved flexbox layout, until paul pointed out issue exists without flexbox layout. the html is: <div> <svg id="not_working" viewbox="0 0 250 1200" preserveaspectratio="xmidymid meet"> <symbol id="test_symbol1" preserveaspectratio="xmidymid meet" viewbox="0 0 250 1200"> <path id="test_symbol_path" d="m 100 1200 l 100 0" /> <text font-size="100" fill="red"> <textpath text-anchor="middle" startoffset="50%" xlink:href="#test_symbol_path...
Read more

Java 8 + Maven Javadoc plugin: Error fetching URL -

i attempting generate javadoc links javadoc dependencies. have tried various means generate javadoc not produce qualified class names references classes dependencies. wanted links java doc simplified class names. however, java api classnames, no links , have qualified class names. working java 8. have following configuration: <plugin> <groupid>org.apache.maven.plugins</groupid> <artifactid>maven-javadoc-plugin</artifactid> <version>2.10.2</version> <configuration> <reportoutputdirectory>${project.basedir}/target</reportoutputdirectory> <destdir>javadoc</destdir> <windowtitle>epiphany</windowtitle> <doctitle>epiphany</doctitle> <show>private</show> <detectlinks>false</detectlinks> <detectoffline...
Read more

datatable - Matlab struct computations -

i download data yahoo! finance with: data = getyahoodailydata({'msft', 'axp'},'01/01/2000', '01/01/2015', 'dd/mm/yyyy'); and data stored 1x1 struct. want create tx2 matrix of daily adjusted closing prices msft , axp, column 7 in each table in struct. how can that? or better : there way make computations directly on information/prices in struct? you can access data in struct name. vecmsftadj = data.msft.adjclose; vecaxpadj = data.axp.adjclose; % if want n x 2 matrix madjclose = [vecmsftadj, vecaxpadj]; % personnaly prefer working table tadjclose = table; tadjclose.msftadj = vecmsftadj; tadjclose.axpadj = vecaxpadj;
Read more